Laboratorio #4 - Valores Propios
Laboratorio #4 - Valores Propios
Laboratorio #4 - Valores Propios
Fin-Desde
3. L ← 𝐿0 , 𝑉 = 𝑣0 ;
Programa#2
Programar (Potencia inversa: Devuelve el menor valor propio)
Datos de ingreso: Matriz A y 𝑢0: Vector inicial, n: # de iteraciones
Salida: V y L (vector y valor propio respectivamente)
𝑢0
4. 𝑣0 = ‖𝑢0‖ ; 𝐿0 = 1/((𝑣0)𝑡 𝐴(𝑣0));
Fin-Desde
6. L ← 𝐿0 , 𝑉 = 𝑣0 ;
𝟓 𝟔 −𝟒 𝟏 𝟎 𝟏
(𝟐 𝟔 𝟐 ) ( 𝟎 𝟏 𝟏)
−𝟐 𝟔 𝟏 𝟏 𝟏 𝟎
Programa#3
Programar (Método de Jacobi: Devuelve todos los valores propios de una matriz simétrica)
Paso #1
Hacer una función [f,c]=eligeMax(A) que devuelva la posición (f, c) del elemento mayor en valor absoluto
sobre la diagonal de A, A es simétrica.
Paso #2
Programe el algoritmo siguiente:
Datos de ingreso: Matriz simétrica A y N=# de iteraciones
Salida: Matriz D (valores propios en su diagonal) y V (matriz cuyas columnas son
vectores propios).
7. A1← A ; Dim ← dimension (A) ; 𝑃 ← 𝐼𝐷𝑖𝑚
8. Desde i=1 hasta N hacer
2.1 [f,c]=eligeMax(A1)
2.2. Si (A1 (f, f)=A1(c, c)) entonces t ← 𝜋/4
Sino 𝑡 = 0.5 ∗ atan ((2 ∗ 𝐴1(𝑓, 𝑐))/(𝐴1(𝑓, 𝑓) − 𝐴1(𝑐, 𝑐)))
2.3. 𝑃1 ← 𝐼𝐷𝑖𝑚 % matriz identidad del orden de la matriz A
2.4. 𝑃1(𝑓, 𝑓) ← cos (𝑡), 𝑃1(𝑐, 𝑐) ← cos (𝑡) , 𝑃1(𝑓, 𝑐) ← −sen(t) , 𝑃1(𝑐, 𝑓) ← sen (𝑡)
2.8. 𝐴1 ← 𝑃1𝑡 ∗ 𝐴1 ∗ 𝑃1; 𝑃 ← 𝑃 ∗ 𝑃1
Fin-Desde
9. D ← A1 , V ← P
10. Fin.
4 1 −1 1
( 1 4 −1 −1 )
−1 −1 5 1
1 −1 1 3