Scribd
Cargar
994 vistas

Ejemplo Krylov

Utilizando el método de Krylov, se obtiene el polinomio característico λ3 − 2λ + 3 = 0 de la matriz dada A = [[1, -1, 0], [2, 0, 1], [-2, 3, -1]]. El método construye la matriz C y el vector D que definen el sistema de ecuaciones lineales CW=D, cuya solución proporciona los coeficientes del polinomio característico.

Cargado por

Alfredo Flores
Derechos de autor
© © All Rights Reserved
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
994 vistas

Ejemplo Krylov

Utilizando el método de Krylov, se obtiene el polinomio característico λ3 − 2λ + 3 = 0 de la matriz dada A = [[1, -1, 0], [2, 0, 1], [-2, 3, -1]]. El método construye la matriz C y el vector D que definen el sistema de ecuaciones lineales CW=D, cuya solución proporciona los coeficientes del polinomio característico.

Cargado por

Alfredo Flores
Derechos de autor
© © All Rights Reserved
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
Está en la página 1/ 4

Análisis Numérico Apuntes

M. I. Jaime Alfonso Reyes Cortés

Ejemplo de Krylov

𝟏 −𝟏 𝟎
Dada la matriz 𝑨 = 𝟐 𝟎 𝟏 aplicar el método de Krylov para obtener su polinomio
−𝟐 𝟑 −𝟏
característico.

R.

Utilizando el algoritmo mejorado de Krylov de los apuntes del tema 3.

Paso 1.

𝟏 −𝟏 𝟎
𝟎 𝟏 𝟐 𝟏 𝟐 𝟎
𝒅𝒆𝒕 𝑨 = 𝟐 𝟎 𝟏 =𝟏 −𝟏 +𝟎
𝟑 −𝟏 −𝟐 −𝟏 −𝟐 𝟑
−𝟐 𝟑 −𝟏
= 𝟏 𝟎 − 𝟑 − 𝟏 −𝟐 − −𝟐 + 𝟎 𝟔 − 𝟎 = −𝟑 ∴ 𝒄𝒐𝒏𝒗𝒆𝒓𝒈𝒆 𝒚 𝒄𝒐𝒏𝒕𝒊𝒏ú𝒂

Paso 2.

𝟏 𝟎 𝟎 𝟎 𝟎
𝒚 = 𝟎 ,𝑪 = 𝟎 𝟎 𝟎 , 𝑫 = 𝟎 , 𝒏 = 𝟑
𝟎 𝟎 𝟎 𝟎 𝟎

Paso 3.

𝒊 = 𝟎, ¿ 𝒊 ≤ 𝒏? 𝑽

¿ 𝒊 = 𝟎? 𝑽

𝟏
𝑪 𝒄𝒋𝒏A𝒊 𝑪 𝒄𝒋𝟑
= =𝒔=𝒚= 𝟎
∀𝒋 ∀𝒋
𝟎

La última columna de C toma el valor del vector 𝒚, que es el mismo que el vector 𝒔, con lo
que quedan como:

𝟎 𝟎 𝟏 𝟏
𝑪= 𝟎 𝟎 𝟎 ,𝒔 = 𝟎
𝟎 𝟎 𝟎 𝟎

𝒊 = 𝟏, ¿ 𝒊 ≤ 𝒏? 𝑽

¿ 𝒊 = 𝟎? 𝑭

¿ 𝒊 ≥ 𝟏 𝑨𝑵𝑫 𝒊 ≤ 𝒏 − 𝟏 ? 𝑽
Análisis Numérico Apuntes

M. I. Jaime Alfonso Reyes Cortés

Ejemplo de Krylov

Ahora se hace

𝟏 −𝟏 𝟎 𝟏 𝟏
𝒔 = 𝑨𝒔 = 𝟐 𝟎 𝟏 𝟎 = 𝟐
−𝟐 𝟑 −𝟏 𝟎 −𝟐

𝟏
Este producto, con 𝒊 = 𝟏 equivale a 𝒔 = 𝑨𝒚 = 𝟐 con lo que
−𝟐

𝟏
𝑪 𝒄𝒋𝒏A𝒊 𝑪 𝒄𝒋𝟐
= =𝒔= 𝟐
∀𝒋 ∀𝒋
−𝟐

La segunda columna de C toma el valor de 𝑨𝒚 a través del vector 𝒔, con lo que quedan como

𝟎 𝟏 𝟏 𝟏
𝑪= 𝟎 𝟐 𝟎 ,𝒔 = 𝟐
𝟎 −𝟐 𝟎 −𝟐

𝒊 = 𝟐, ¿ 𝒊 ≤ 𝒏? 𝑽

¿ 𝒊 = 𝟎? 𝑭

¿ 𝒊 ≥ 𝟏 𝑨𝑵𝑫 𝒊 ≤ 𝒏 − 𝟏 ? 𝑽

Nuevamente, se hace

𝟏 −𝟏 𝟎 𝟏 −𝟏
𝒔 = 𝑨𝒔 = 𝑨𝑨𝒚 = 𝟐 𝟎 𝟏 𝟐 = 𝟎
−𝟐 𝟑 −𝟏 −𝟐 𝟔

−𝟏
Lo que, con 𝒊 = 𝟐 equivale a haber hecho 𝒔 = 𝑨𝟐 𝒚 = 𝟎
𝟔

con lo que, la primera columna de C toma el valor del producto 𝑨𝟐 𝒚 a través del vector 𝒔

−𝟏
𝑪 𝒄𝒋𝒏A𝒊 𝑪 𝒄𝒋𝟏
= =𝒔= 𝟎
∀𝒋 ∀𝒋
𝟔

quedando como:

−𝟏 𝟏 𝟏 −𝟏
𝑪= 𝟎 𝟐 𝟎 ,𝒔 = 𝟎
𝟔 −𝟐 𝟎 𝟔
Análisis Numérico Apuntes

M. I. Jaime Alfonso Reyes Cortés

Ejemplo de Krylov

Al llegar a la penúltima iteración, la matriz C queda llena, es decir, ya se tiene la matriz de


coeficientes del nuevo sistema de ecuaciones lineales.

𝒊 = 𝟑, ¿ 𝒊 ≤ 𝒏? 𝑽

¿ 𝒊 = 𝟎? 𝑭

¿ 𝒊 ≥ 𝟏 𝑨𝑵𝑫 𝒊 ≤ 𝒏 − 𝟏 ? 𝑭

¿ 𝒊 = 𝒏? 𝑽

Se vuelve a hacer

𝟏 −𝟏 𝟎 −𝟏 −𝟏
𝒔 = 𝑨𝒔 = 𝑨𝑨𝟐 𝒚 = 𝟐 𝟎 𝟏 𝟎 = 𝟒
−𝟐 𝟑 −𝟏 𝟔 −𝟒

−𝟏
que con 𝒊 = 𝒏 equivale a haber hecho 𝒔 = 𝑨𝟑 𝒚 = 𝟒
−𝟒

Luego, para formar el vector de términos independientes

−𝟏 𝟏
𝑫 = −𝒔 = − 𝟒 = −𝟒
−𝟒 𝟒

𝒊 = 𝟒, ¿ 𝒊 ≤ 𝒏? 𝑭

salimos del ciclo

Paso 4.

El sistema de ecuaciones lineales en forma matricial 𝑪𝑾 = 𝑫 queda como:

−𝟏 𝟏 𝟏 𝒃𝟏 𝟏
𝟎 𝟐 𝟎 𝒃𝟐 = −𝟒
𝟔 −𝟐 𝟎 𝒃𝟑 𝟒

que equivale al mismo sistema que ya habíamos resuelto donde

𝒃𝟏 = 𝟎, 𝒃𝟐 = −𝟐, 𝒃𝟑 = 𝟑

Por lo que la ecuación característica es

𝝀𝟑 + 𝒃𝟏 𝝀𝟐 + 𝒃𝟐 𝝀 + 𝒃𝟑 = 𝟎 ⟹ 𝝀𝟑 − 𝟐𝝀 + 𝟑 = 𝟎
Análisis Numérico Apuntes

M. I. Jaime Alfonso Reyes Cortés

Ejemplo de Krylov

Con lo que termina el ejemplo.

Queda pendiente hallar las raíces de la ecuación característica, para ello, se tendría que
aplicar uno de los métodos de obtención de raíces que hemos visto durante el curso.

También podría gustarte