domingo, 21 de junio de 2015

Diagonalización de matrices de 2x2, 3x3 y 4x4. Y aplicación que se le puede dar a la diagonalización.

  Antes de comenzar repasaremos algunas conceptos para la mejor comprensión del tema que nos ocupa, que en este caso es, la diagonalización.


Matriz identidad: Matriz cuadrada que tiene en cada una de las entradas de su diagonal principal el valor 1 y en el resto 0.




Ejemplos de matriz Identidad de orden 2 (2x2), orden 3 (3x3) y orden 4 (4x4).

Matriz Identidad de Orden 2 (2x2), Orden 3 (3x3) y Orden 4 (4x4)

Polinomio característico: Se obtendrá de aplicar el determinante l A-λI l=0 e igualando a cero para obtener las soluciones o raíces.

Matriz de Paso: Matriz que se forma de los autovectores o vectores propios puestos en columna.

Autovalores o valores propios: Son las soluciones que se obtienen del polinomio característico y pueden tener multiplicidad mayor de 1.


Autovectores o vectores propios: Son los vectores no nulos, que obtenemos a partir de los valores propios obtenidos en el polinomio característico. La forma de obtener los vectores es usando la forma:

Formula para obtener vectores propios o autovectores

Ejemplo de como podemos obtener un autovector a partir de un autovalor

Supongamos que tenemos el autovalor  λ = -1 de la matriz :

 Matriz simétrica de orden 2 (2x2)

Aplicando la formula Formula para obtener vectores propios o autovectorescitada anteriormente tendríamos:

Resolución de una autovector

  Como sabemos que λ = -1 sustituimos en la matriz y resolvemos, quedando un sistema de ecuaciones, donde se aprecia que son la misma ecuación:

2 V1 + 2 V2 = 0
2 V1 + 2 V2 = 0

  Eliminaremos una ecuación y resolvemos quedándonos que V1 = -V2, con lo que podemos determinar el vector propio asociado al valor propio λ = -1, siendo este:

                       (-V2 , V2 ) = V2 ( -1 , 1 ) , El vector propio es ( -1 , 1 )



Diagonalización de Matrices de orden 2 (2x2).

Vamos a comenzar diagonalizando una matriz no simétrica de orden 2, como por ejemplo:

Matriz de 2x2

Pasos a seguir:

  Obtendremos el polinomio característico de la matriz mediante l A-λI l=0

Cálculo de polinomio característico de matriz de 2x2



  Solucionamos el determinante obteniendo: 

(1-λ) (-2-λ) = 0  →  λ ² + λ - 2 = 0  

  El polinomio característico sería   λ ² + λ - 2 = 0  y nos da como resultado λ = 1 y  λ = -2, siendo estos los autovalores o valores propios asociados a la matriz A.

  Puesto que los autovalores son distintos, ya podríamos predecir que la matriz será diagonalizable. Ademas ya podríamos construir la matriz diagonal a partir de sus valores propios.
Matriz diagonal de 2x2

  Continuamos ahora con el cálculo de autovectores o vectores propios, para ello usaremos la forma:
Cálculo de autovectores

  Para el valor propio:  λ = 1.

Cálculo del autovector para el valor 1

  Como resultado obtenemos la ecuación: 

 x -3y = 0 →  x = 3y

sustituimos en ( x , y ) y obtenemos

(x,y) → (3y , y) = y (3 , 1)

  Siendo el autovector: (3 , 1)


  Para el valor propio:  λ = -2

   Como resultado obtenemos la ecuación: 

 3x = 0  x = 0

sustituimos en ( x , y ) y obtenemos
                           
(x,y) → (0 , y) = y (0 , 1)
    Siendo el autovector: (0 , 1)

  La matriz de paso y su inversa serán :   




  Teniendo la matriz diagonal, la de paso y su inversa, vemos si se cumple:




-----------------------------

  Como segunda matriz a diagonalizar será :  
Matriz 2x2


 Obtenemos el polinomio característico de la matriz mediante l A-λI l=0

Resolución de polinomio característico de 2x2

Solucionamos el determinante obteniendo: 

(1-λ) (1-λ) - 0 = 0  →  (1-λ) (1-λ) = 0, las soluciones en este caso se nos dan directas sin necesidad de hacer ningún calculo, siendo λ = 1 con multiplicidad = 2, al tener un autovalor que se repite en este caso no podemos predecir que la matriz será diagonalizable, sino que debemos calcular los vectores propios para saber, si es o no, diagonalizable.

   Construimos su matriz diagonal: 


Matriz diagonal 2x2

   Calculamos autovectores o vectores propios, usando la forma:


Cálculo de autovectores


   Para el valor propio:  λ = 1


Cálculo de autovectores para valor 1

  Como resultado, obtenemos la ecuación: 
 -x = 0 → x = 0

(x,y) → (0, y) = y (0 , 1)

  Puesto que solo hemos obtenido un vector propio, la matriz de paso no se puede construir y por lo tanto, la matriz A propuesta, no es diagonalizable.



  Diagonalización de Matrices de orden 3 (3x3).

  Primero elegimos una matriz cuadrada de 3x3, que llamaremos A.

  Los pasos a seguir serán, primero calcular su matriz diagonal a partir de sus autovalores o valores propios y encontrar la matriz de paso que se genera a partir de los autovectores o vectores propios.

  1. Cálculo de los valores propios mediante el polinomio característico,                                   usaremos  l A-λI l=0


Cálculo de polinomio característico de matriz 3x3

  Solucionando el determinante obtenemos que (1-λ) (2-λ) (3-λ) = 0  de aquí se desprenden de forma directa, los autovalores, que determinan la matriz diagonal y nos auguran que la matriz será diagonalizable, ya que todos sus autovalores poseen multiplicidad 1, o sea, que no se repiten.

Matriz diagonal 3x3 de A

  2. Cálculo de los vectores propios o autovectores.  De la misma forma que hicimos en la matriz de 2x2.





  Comenzamos a calcular los vectores propios con el autovector  λ=1


Cálculo del primer autovector

  Resolviendo obtenemos dos ecuaciones.

x + y + 3z = 0 → x = -y
2z = 0  → z = 0

sustituimos en ( x , y , z ) y obtenemos

  ( x , y , z ) → (-y , y , 0 ) → y (-1 , 1 , 0 )


  Ya hemos conseguido el primer autovector : ( -1 , 1 , 0 )


  Continuamos con el autovector o vector propio  λ = 2


Resolviendo obtenemos dos ecuaciones.

             -x    = 0     → x = 0
   x + 3z = 0    
                  z = 0     →     z = 0 

sustituimos en ( x , y , z ) y obtenemos

  ( x , y , z ) → (0 , y , 0 ) → y (0 , 1 , 0 )


el segundo vector propio : ( 0 , 1 , 0 )


  El último autovector o vector propio  λ = 3


Resolviendo obtenemos dos ecuaciones.

             -2x    = 0     → x = 0
   x -y + 3z = 0    → y = 3z
                 

sustituimos en ( x , y , z ) y obtenemos

  ( x , y , z ) → (0 , 3z , z ) → z (0 , 3 , 1 )



  El tercer vector propio sería: (0 , 3 , 1 )


  Al obtener los tres autovectores, podemos construir la matriz de paso, colocando los vectores propios en columna.


 Diagonalización de Matrices de orden 4 (4x4).



  Vamos a escoger una matriz 4x4, por ejemplo esta:

  
  Para el cálculo de autovectores y vectores propios, procederemos igual que en los casos de orden 2 y orden 3. Así que a partir de  l A-λI l=0  obtendrenos los valores propios y podremos construir la matriz diagonal. 





  Aplicando la propiedad de los determinantes, ya que la primera columna esta compuesta de ceros excepto un valor, de esta forma queda un determinante de 3x3 mas sencillo de resolver. Aunque podríamos reducirlo mas, lo solucionaremos con esta forma.



  Al solucionarlo nos queda (1-λ) (1-λ) (2-λ) (-1-λ) = 0, teniendo sus cuatro autovalores, ya podríamos construir la matriz diagonal. Como tenemos que λ = 1 tiene multiplicidad 2, este sera el primer valor propio en verificar, pues necesitamos obtener dos vectores de él.






 Para λ = 1



  Resolviendo obtenemos  tres ecuaciones.

2y - z + t = 0
y = 0
z = 0
z -2t = 0  →   t = 0                     

sustituimos en ( x , y , z , t ) y obtenemos



  ( x , y , z , t ) → (x , 0 , 0 , 0 ) → x (1 , 0 , 0 , 0 )

  Puesto que solo hemos obtenido un vector y necesitábamos dos la matriz no será diagonalizable. 

----------------------------------------------------




Vamos a diagonalizar otra matriz, por ejemplo.

  Calculamos sus autovalores a partir de  l A-λI l=0



 Aplicando la propiedad de los determinantes, nos queda:



 Al solucionar el determiante obtenemos (1-λ) (-1-λ) (2-λ) (-1-λ) = 0, teniendo sus cuatro autovalores, ya podríamos construir la matriz diagonal. Como tenemos que λ = -1 tiene multiplicidad 2, este sera el primer valor propio en verificar.



Para  λ = -1



Resolviendo obtenemos dos ecuaciones.

2x + 2y - z + t = 0   →  t = -2x -2y

                  3z = 0    →  z = 0


sustituimos en ( x , y , z , t ) y obtenemos



  ( x , y , z , t ) → (x , y , 0 , -2x -2y) →

→(x , 0 , 0 , -2x ) + ( 0 , y , 0 , -2y) = x (1 , 0 , 0 , -2 ) + y( 0 , 1 , 0 , -2)

  Hemos obtenido dos autovectores V2 (1 , 0 , 0 , -2 ) y V4 ( 0 , 1 , 0 , -2), el hecho de que sean V2 y V4 lo explico mas adelante.

  

  Para  λ = 2



Resolviendo obtenemos tres ecuaciones.

-x + 2y - z + t = 0   →  x = -z       

         -3y = 0  →  y = 0

         -3t = 0  →  t = 0


sustituimos en ( x , y , z , t ) y obtenemos



  ( x , y , z , t ) → (-z , 0 , z , 0) → z (-1 , 0 , 1 , 0)



  Hemos obtenido otro autovector  V3 (-1 , 0 , 1 , 0 )



  Para  λ = 1



Resolviendo obtenemos tres ecuaciones.

              2y - z + t = 0   

         -2y = 0  →  y = 0

           z = 0  →  z = 0

         -2t = 0  →  t = 0


sustituimos en (  x , y , z , t ) y obtenemos



  ( x , y , z , t ) → (x , 0 , 0 , 0) → x (1 , 0 , 0 , 0)



  Hemos obtenido otro autovector  V1 (1 , 0 , 0 , 0 )


  La matriz de paso será:

  Los vectores propios de la matriz de paso se colocan en la columna exacta donde están los valores propios que los generan, por este motivo que los dos primeros vectores que haya obtenido sean V2 y V4 no es al azar, es porque el autovalor que los genera esta en la columna dos y cuatro de la matriz diagonal. Os pongo la matriz diagonal para que veáis donde he situado los autovalores.

  Para comprobar si el ejercicio esta bien y poder asegurar que la matriz es diagonalizable necesitamos la matriz de paso inversa, que sería esta:

  Y por último la comprobación:


  Si hacéis las operaciones comprobareis que da como resultado la matriz M, puesto que esta matriz es diagonalizable.















10 comentarios:

  1. y si la matriz 4x4 tiene en la columna en vez de 0 tiene un dos y dos 0?

    ResponderEliminar
  2. Hola Borja, haces referencia al enunciado de la matriz 4x4, pero no se exactamente a que columna te refieres, si me pudieras especificar la columna que quieres modificar, podría ayudarte en tu consulta.

    ResponderEliminar
  3. Buena explicacion pero en la matriz de 3x3 como terminaria el ejercicio ...... muy bien explicado

    ResponderEliminar
  4. Hola Erick, gracias por tu comentario. Te comento, el último paso sería comprobar que A = P*D*P(inversa), donde A es la matriz que estamos diagonalizando, P es la matriz de paso que se genera a partir de los vectores propios, D es la matriz diagonal de A que se ha obtenido a partir de los autovalores y P(inversa) se obtiene haciendo la inversa a la matriz P. Si se cumple la igualdad, aseguramos que la matriz en cuestión es digonalizable, pero como en el caso 3x3 que he desarrollado, hemos obtenido todos sus autovalores distintos, ya podemos afirmar que la matriz es diagonalizable por este motivo no se ha realizado el último paso, pero si me lo pides y a ti te sirve, no tendré ningún problema en acabarlo para mejor comprensión del mismo. Saludos

    ResponderEliminar
  5. En la primer matriz 2x2 el resultado de A=PDP(INVERSA) el resultado me da en la segunda fila (0,-2) estoy mal en algo? Gracias

    ResponderEliminar
  6. Hola, el resultado debe de darte igual que la matriz que se diagonaliza, o sea, la matriz A, esto es solamente una comprobación para ver que la matriz diagonal y de paso son correctas.

    ResponderEliminar
  7. Propongo como ejemplo, esta matriz 4x4. Por favor, cuando puedas, dime como podria despejarla para hacer una columna o fila todo 'ceros' salvo un valor para hacer el ejercicio como tu lo has hecho:
    1 1 1 1
    1 -1 0 0
    1 0 -1 0
    1 0 0 -1

    ResponderEliminar
    Respuestas
    1. Este comentario ha sido eliminado por el autor.

      Eliminar
    2. Hola, para hacer ceros lo mas sencillo es hacer
      Fila2 = Fila1 - Fila2
      Fila3 = Fila1 - Fila3
      Fila4 = Fila1 - Fila4

      El resultado sería:
      1 1 1 1
      0 2 1 1
      0 1 2 1
      0 1 1 2

      Espero haberte ayudado. Un saludo.

      Eliminar
  8. Este comentario ha sido eliminado por el autor.

    ResponderEliminar