domingo, 18 de diciembre de 2022

Diagonalización de matrices con GNU Octave. Valores y vectores propios

Para obtener los valores propios y los vectores propios de una matriz $A$ existe también una instrucción en GNU Octave. A modo de ejemplo, voy encontrar los valores y vectores propios de la matriz de orden $3$ $$A=\begin{pmatrix}1&2&3\\ 4&5&6 \\ 7&1&9 \end{pmatrix}$$ utilizando GNU Octave (de manera automática). La instrucción nos da la matriz del cambio de base $C$, siendo por tanto los coeficiente de la misma, por columnas, las coordenadas de los vectores propios; además, $D$ representa la matriz diagonal: en cuya diagonal principal tendremos los valores propios correspondientes.

  >> A=[1,2,3;4,5,6;7,1,9]
A =

   1   2   3
   4   5   6
   7   1   9

>> [C,D]=eig(A)
C =

  -0.285862  -0.818536  -0.227349
  -0.667939  -0.027924  -0.877584
  -0.687125   0.573776   0.422088

D =

Diagonal Matrix

   12.8843         0         0
         0   -1.0347         0
         0         0    3.1505


>> eig(A)
ans =

   12.8843
   -1.0347
    3.1505

-oOo-

Referencias

[1] John W. Eaton; David Bateman; Søren Hauberg; Rik Wehbring, Free Your Numbers (Manual de GNU Octave), https://www.gnu.org/software/octave/octave.pdf

$\diamond$

No hay comentarios:

Publicar un comentario