Processing math: 100%

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