Autor Tema: Cálculo de la matriz de Jordan.

0 Usuarios y 1 Visitante están viendo este tema.

28 Mayo, 2024, 11:01 am
Respuesta #10

josepapaiii

  • $$\Large \color{#5b61b3}\pi\,\pi\,\pi\,\pi\,\pi$$
  • Mensajes: 1,061
  • Karma: +0/-0
  • Sexo: Masculino
Hola

¿Sería esta parte correcta?

Si.

Citar
Esto funciona siempre; aunque a veces puedes acabarse más rápido. Por ejemplo si \( m.a=4 \) y \( d_1=3 \) la única posibilidad de tres cajas de Jordan de tamaño total \( 4 \) es \( 2+1+1+1 \). Esto correspondería a este esquema:

\( \begin{array}{cccc}
{4}&{1}\\
{\downarrow}&{\downarrow}\\
{*}&{*}\\
{*}&{}\\
{*}&{}\\
{*}&{}\\
\end{array} \)

Esto no lo he comprendido bien. Para el caso que pones, si \( m_a=4 \), entonces estamos diciendo que el tamaño total de los bloques es 4 y si \( d_1=3 \) estamos diciendo que habrá tres cajas de Jordan, ¿no sería entonces la distribución 2+1+1?

Perdona, fue una errata. Ya lo he corregido. Sería;

\( \begin{array}{cccc} {\color{red}3\color{black}}&{1}\\ {\downarrow}&{\downarrow}\\ {*}&{*}\\ {*}&{}\\ {*}&{}\\ \end{array} \)

Saludos.

Hola, Luis.

Al final me he dado cuenta que todo debe de encajar (como no puede ser de otra manera). Ciertamente (y me costó verlo) al final el número total de de las cajas debe ser igual a \( m_a \) como de igual manera pasa cuando una matriz es diagonalizable, vas a tener tantas veces el valor propio en la diagonal como su multiplicidad algebraica. Luego otro tema es saber cuántas cajas son y su tamaño. En una matriz diagonal, \( d_1 \) siempre va a ser igual la multiplicidad algebraica y geométrica, de ahí que las cajas sean de como mucho tamaño 1 (matriz diagonal). Los \( f_i \) nos darán el tamaño de cada una de las cajas.

Voy ahora con el tema de la matriz de paso. Miraré con detenimiento la teoría que me has puesto.

Muchas gracias por la ayuda.

Saludos.


28 Mayo, 2024, 11:53 am
Respuesta #11

feriva

  • $$\Large \color{#a53f54}\pi\,\pi\,\pi\,\pi\,\pi\,\pi\,\pi$$
  • Mensajes: 12,032
  • País: es
  • Karma: +1/-0
  • Sexo: Masculino

Hola, Josepapaiii, te lo hago a mi manera para la matriz B.

\( B=\left(\begin{array}{cccc}
0 & 1 & -2 & 1\\
-2 & 1 & -6 & 3\\
2 & -3 & 0 & 1\\
2 & -3 & -2 & 3
\end{array}\right) \)

Una vez hallado el polinomio característico se ve que tiene el autovalor 2 triple y el autovalor -2 simple. Y la multiplicidad geométrica del autvalor 2 es dos; con esto falta un autovector para que la matriz B sea diagonalizable y entonces es necesario completar la matriz “J” con algún “1” añadido en la caja de los doses positivos.

En principio, la disposición para la matriz de Jordan no es única, depende de cómo se coloquen los vectores de la base.

Se puede formar primero la matriz “J” y después hallar la base (junto con el orden apropiado de los vectores).

Si colocamos los autovalores en una matriz diagonal, por ejemplo, así, arbitrariamente con el -2 en la caja de arriba

\( \left(\begin{array}{cccc}
-2\\
 & 2\\
 &  & 2\\
 &  &  & 2
\end{array}\right) \)

el “1” que se añade lo podemos poner, por ejemplo, aquí

\( \left(\begin{array}{cccc}
-2 & 0 & 0 & 0\\
0 & 2 & 0 & 0\\
0 & 1 & 2 & 0\\
0 & 0 & 0 & 2
\end{array}\right) \)

o aquí

\( \left(\begin{array}{cccc}
-2 & 0 & 0 & 0\\
0 & 2 & 0 & 0\\
0 & 0 & 2 & 1\\
0 & 0 & 0 & 2
\end{array}\right) \)

...

*(En el libro de Álgebra lineal y Geometría de Alba Valverde Colmeiro, que es el que yo usé, los unos se colocaban debajo de los autovalores; pero veo que la tendencia más habitual es ponerlos por encima, así que me quedo con la última forma).

Es decir, una vez hallado el polinomio característico (y los autovalores y su multiplicadad) la matriz J ya está hallada, sólo hay que escribirla con una colocación arbitraria; respetando que los autovalores vayan colocados por cajas, sin mezclarse, eso sí, pero sin importar si tal caja está más arriba o más abajo.

Para hallar la base adaptada a una determinada colocación de la forma de Jordan, multiplicamos los vectores de la base canónica por “J”: \( J\vec{e_{1}} \), \( J\vec{e_{2}} \) etc. El subíndice se va a corresponder con los vectores de la base cuando los hallemos.

Por ejemplo, para el primer caso (no hace falta ni escribirlo, pero para que se vea) tenemos

\( \left(\begin{array}{cccc}
-2 & 0 & 0 & 0\\
0 & 2 & 0 & 0\\
0 & 0 & 2 & 1\\
0 & 0 & 0 & 2
\end{array}\right)\left(\begin{array}{c}
1\\
0\\
0\\
0
\end{array}\right)=\left(\begin{array}{c}
-2\\
0\\
0\\
0
\end{array}\right) \)

lo que se traduce en esta ecuación \( J\vec{e_{1}}=-2\vec{e_{1}} \).

Y así con el resto. Entonces quedan estas ecuaciones (quito las flechitas por comodidad, pero son vectores)

\( Je_{1}=-2e_{1} \)

\( Je_{2}=2e_{2} \)

\( Je_{3}=2e_{3} \)

\( Je_{4}=e_{3}+2e_{4} \).

Y ahora, despejando y sacando factor común a los vectores en la primera y siguientes, tienes

\( (J+2I)e_{1}=\overrightarrow{0} \)

\( (J-2I)e_{2}=\overrightarrow{0} \)

\( (J-2I)e_{3}=\overrightarrow{0} \)

\( (J-2I)e_{4}=e_{3} \)

Sustituyendo \( (J-2I)e_{4}=e_{3} \) en la tercera

\( (J-2I)^{2}e_{4}=\overrightarrow{0} \)

Esas ecuaciones se van a cumplir igual para la matriz B (en vez de J) respecto de la base de vectores que buscas ( en vez de respecto los vectores de la base canónica). Es decir:

Si \( (J+2I)e_{1}=\overrightarrow{0} \), entonces \( (B+2I)v_{1}=\overrightarrow{0} \), donde ese “v” es el primer vector que aparecerá (en columna) a la izquierda en la matriz base. Y así con los demás.

Y ya con eso supongo que ves cómo se termina.

De esta forma, hallas la base, como te decía, adaptada a la colocación arbitraria de la matriz J que hayas elegido.

Si no me equivoco, una base de vectores puede ser

\( \left(\begin{array}{cccc}
0 & 0 & 2 & 0\\
1 & 0 & 2 & 1\\
1 & 1 & -2 & 0\\
1 & 2 & -2 & 1
\end{array}\right) \)

De tal forma que tendría que funcionar

\( B=\left(\begin{array}{cccc}
0 & 1 & -2 & 1\\
-2 & 1 & -6 & 3\\
2 & -3 & 0 & 1\\
2 & -3 & -2 & 3
\end{array}\right)=\left(\begin{array}{cccc}
0 & 0 & 2 & 0\\
1 & 0 & 2 & 1\\
1 & 1 & -2 & 0\\
1 & 2 & -2 & 1
\end{array}\right)\left(\begin{array}{cccc}
-2 & 0 & 0 & 0\\
0 & 2 & 0 & 0\\
0 & 0 & 2 & 1\\
0 & 0 & 0 & 2
\end{array}\right)\left(\begin{array}{cccc}
0 & 0 & 2 & 0\\
1 & 0 & 2 & 1\\
1 & 1 & -2 & 0\\
1 & 2 & -2 & 1
\end{array}\right)^{-1} \)

Saludos.

28 Mayo, 2024, 12:36 pm
Respuesta #12

josepapaiii

  • $$\Large \color{#5b61b3}\pi\,\pi\,\pi\,\pi\,\pi$$
  • Mensajes: 1,061
  • Karma: +0/-0
  • Sexo: Masculino
Hola, Josepapaiii

Hola, feriva. Estoy muy agradecido de tu respuesta, me sirve muchísimo lo que has puesto.

Muchas gracias.

28 Mayo, 2024, 04:43 pm
Respuesta #13

josepapaiii

  • $$\Large \color{#5b61b3}\pi\,\pi\,\pi\,\pi\,\pi$$
  • Mensajes: 1,061
  • Karma: +0/-0
  • Sexo: Masculino
Ahora si además quieres saber la base en la cual se obtiene esa forma de Jordan reproduce el esquema anterior de la siguiente forma:

\( \begin{array}{cccc}
{ker(A-\lambda Id)}&{ker((A-\lambda Id)^2)}&{ker((A-\lambda Id)^3)}&\ldots&{ker((A-\lambda Id)^k)}\\
{\downarrow}&{\downarrow}&{\downarrow}&{\ldots}&{\downarrow}\\
{u_1}&{u_2}&{u_{3}}&\ldots&{u_k}\\
{v_1}&{v_2}&{v_3}&&{}\\
{w_1}&{w_2}&{}&&{}\\
\end{array} \)

Por cada asterisco vas a tener un vector. Comienzas hallando el último de cada fila, empezando por la primera fila, escogiendo \( u_k \) verificando:

\( (A-\lambda Id)^ku_k=0 \) (el exponente coincide con la columna donde está el vector)

y después directamente el resto de la fila es automático hacia atrás multiplicando por \( (A-\lambda Id) \) de cada vez:

\( u_{k-1}=(A-\lambda Id)U_k,\quad  u_{k-2}=(A-\lambda I)u_{k-1d} \)

Pasas a la segunda columna; calculas el último vector.

\( (A-\lambda Id)^3v_3=0 \) (el exponente coincide con la columna donde está el vector)

y OJO lo escoges que sea linealmente independiente de los que tiene encima de él en la misma columna.

Después el resto de la fila automático hacia atrás multiplicando por \( A-\lambda Id \)

Y así sucesivamente. La base asociada a la forma de Jordan es la formada por todos esos vectores de izquierda a derecha y de arriba a abajo:

\( u_1,u_2,\ldots,u_k,v_1,v_2,v_3,w_1,\ldots \)


Hola, Luis.

Por aquí voy:

En mi caso tenemos que calcular:

\( \begin{array}{cc}
{ker(A-\lambda Id)}&{ker((A-\lambda Id)^2)}\\
{\downarrow}&{\downarrow}\\
{u_1}&{u_2}\\
{v_1}&{v_2}\\
\end{array} \)

Donde \( ((A-\lambda Id)^2)u_2=0 \) y \( u_1=(A-\lambda Id)u_2 \)
De la misma manera: \( ((A-\lambda Id)^2)v_2=0 \) y \( v_1=(A-\lambda Id)v_2 \).

Calculamos \( u_2 \):

Como \( (A-\lambda Id)^2)=0 \), entonces un posible vector puede ser \( \{\alpha, \beta, \gamma, \delta \} \), para cualquier valor de \( \alpha, \beta, \gamma, \delta \in \mathbb{R} \). Tomamos por ejemplo \( u_2=\{1,0,0,0\} \)

Calculamos \( u_1 \):

\( u_1=\left(\begin{matrix}-1 & -1 & -1 & -1\\1 & 1 & 1 & 1 \\ 0 & 0 & -1 & -1 \\ 0 & 0 & 1 & 1\end{matrix}\right)\left(\begin{matrix} 1 \\ 0 \\ 0 \\ 0\end{matrix}\right)=\{-1,1,0,0\} \)

Calculamos \( v_2 \):

Podemos tomar por ejemplo el vector \( \{0,0,1,0\} \). Ya que tiene que ser linealmente independiente de los anteriores, y una base sería como la base de \( u_2 \) con los parámetros libres, basta con tomar  \( \gamma=1 \) y el resto a 0.

Calculamos \( v_1 \):

\( v_1=\left(\begin{matrix}-1 & -1 & -1 & -1\\1 & 1 & 1 & 1 \\ 0 & 0 & -1 & -1 \\ 0 & 0 & 1 & 1\end{matrix}\right)\left(\begin{matrix} 0 \\ 0 \\ 1 \\ 0\end{matrix}\right)=\{-1,1,-1,1\} \)

Siendo entonces una base \( \mathcal{B}=\{(-1,1,0,0)^t,(1,0,0,0)^t,(-1,1,-1,1)^t,(0,0,1,0)^t\}] \) y la matriz de paso esos vectores colocados por columnas:

\( P=\left(\begin{matrix}-1 & 1 & -1 & 0\\1 & 0 & 1 & 0 \\ 0 & 0 & -1 & 1 \\ 0 & 0 & 1 & 0\end{matrix}\right) \)

Cumpliéndose que: \( J=P^{-1}AP \).

Los cálculos me salen correctos, pero no sé si es por casualidad o porque está bien.

Saludos.




28 Mayo, 2024, 05:26 pm
Respuesta #14

Masacroso

  • “Lo que oigo, lo olvido; lo que veo, lo recuerdo; lo que hago, lo aprendo” (antiguo proverbio chino)
  • Moderador Global
  • Mensajes: 5,485
  • País: es
  • Karma: +0/-0
*(En el libro de Álgebra lineal y Geometría de Alba Valverde Colmeiro, que es el que yo usé, los unos se colocaban debajo de los autovalores; pero veo que la tendencia más habitual es ponerlos por encima, así que me quedo con la última forma).

Sí, así es. Y está muy bien que lo hayas mencionado por si a josepapaiii se lo han enseñado con los unos "por abajo".

28 Mayo, 2024, 05:53 pm
Respuesta #15

Luis Fuentes

  • el_manco
  • Administrador
  • Mensajes: 58,907
  • País: es
  • Karma: +0/-0
Hola

*(En el libro de Álgebra lineal y Geometría de Alba Valverde Colmeiro, que es el que yo usé, los unos se colocaban debajo de los autovalores; pero veo que la tendencia más habitual es ponerlos por encima, así que me quedo con la última forma).

Sí, así es. Y está muy bien que lo hayas mencionado por si a josepapaiii se lo han enseñado con los unos "por abajo".

Que se pongan arriba  o abajo tiene que ver con el convenio de si se trabajan con los vectores como filas o como columnas.

Lo más habitual es como columnas.

Si los unos están por encima de la diagonal en la forma de Jordan al formar la matriz \( P \) los vectores hallados deben de colocarse como las columnas de las mismas.

Si los unos están por debajo de la diagonal en la forma de Jordan al formar la matriz \( P \) los vectores hallados deben de colocarse como las filas de las mismas.

Saludos.

29 Mayo, 2024, 11:01 am
Respuesta #16

Luis Fuentes

  • el_manco
  • Administrador
  • Mensajes: 58,907
  • País: es
  • Karma: +0/-0
Hola

Los cálculos me salen correctos, pero no sé si es por casualidad o porque está bien.

Está bien.

Del procedimiento donde describía el método general, he corregido esta frase:

Citar
y OJO lo escoges que los vectores que se van obteniendo en la primera columna sean  linealmente independientes (y entre otras cosas no nulos) de los que tiene encima de él en la misma columna.

Es decir cada vez que escoges un vector en la última columna de cada fila (pongamos la \( k \)-ésima) debes de hacerlo de manera que al pasarlo a la primera multiplicando por \( (A-\lambda Id)^{k-1} \) de un vector independiente de los anteriormente hallados en esa primera columna (y en particular no nulo).

Saludos.

29 Mayo, 2024, 11:37 am
Respuesta #17

josepapaiii

  • $$\Large \color{#5b61b3}\pi\,\pi\,\pi\,\pi\,\pi$$
  • Mensajes: 1,061
  • Karma: +0/-0
  • Sexo: Masculino
Hola

Del procedimiento donde describía el método general, he corregido esta frase:

Citar
y OJO lo escoges que los vectores que se van obteniendo en la primera columna sean  linealmente independientes (y entre otras cosas no nulos) de los que tiene encima de él en la misma columna.

Es decir cada vez que escoges un vector en la última columna de cada fila (pongamos la \( k \)-ésima) debes de hacerlo de manera que al pasarlo a la primera multiplicando por \( (A-\lambda Id)^{k-1} \) de un vector independiente de los anteriormente hallados en esa primera columna (y en particular no nulo).

Saludos.

Hola, Luis. Gracias por la aclaración, la tendré en cuenta también.

Saludos.

29 Mayo, 2024, 10:50 pm
Respuesta #18

josepapaiii

  • $$\Large \color{#5b61b3}\pi\,\pi\,\pi\,\pi\,\pi$$
  • Mensajes: 1,061
  • Karma: +0/-0
  • Sexo: Masculino
Hola,

Para la matriz B, intentando replicar el ejercicio del forero feriva por el método de Luis, no me cuadran las cuentas:

Sea \( B=\left(\begin{matrix} 0 & 1 & -2 & 1 \\ -2 & 1 & -6 & 3 \\ 2 & -3 & 0 & 1 \\ 2 & -3 & -2 & 3 \end{matrix}\right). \)

El polinomio característico es: \( P(\lambda)=(\lambda-2)^3(\lambda+2) \), con autovalores \( \lambda_1=2 \) y \( \lambda_2=-2 \)

Para el autovalor \( \lambda_2=-2 \), se obtiene una única caja de Jordan y de tamaño 1. Y un vector \( w_1 \) que nos servirá para luego la matriz de paso, me coincide con el vector que obtiene feriva,  sacándolo del subespacio generado por el autovalor, es decir \( Ker (B+2I)w_1=0 \) (comprobado), siendo \( w_1=\{0,1,1,1\} \)

Para el autovalor \( \lambda_1=2 \) hago el método de las \( f_i \), obteniendo:

\( \begin{array}{cc}
{f_1}&{f_2}\\
{\downarrow}&{\downarrow}\\
{*}&{*}\\
{*}
\end{array} \)

Y la matriz de Jordan sería:

\( \left(\begin{matrix}-2 & 0 & 0 & 0\\0 & 2 & 0& 0 \\ 0 & 0 & 2 & 1 \\ 0 & 0 & 0 & 2\end{matrix}\right) \)

Ahora, necesitaría 3 vectores para formar la matriz de paso (que junto al vector del autovalor \( \lambda_2 \) ya la tendríamos) y utilizo el esquema:

\( \begin{array}{cc}
{ker(A-\lambda Id)}&{ker((A-\lambda Id)^2)}\\
{\downarrow}&{\downarrow}\\
{u_1}&{u_2}\\
{v_1}\\
\end{array} \)

Donde \( ((A-\lambda Id)^2)u_2=0 \) y \( u_1=(A-\lambda Id)u_2 \)
Y para \( v_1 \), en este caso uso: \( (A-\lambda Id)v_1)=0 \).

Calculamos \( u_2: \)

\( (B-2I)^2u_2=0 \Rightarrow \left(\begin{matrix}0 & 1 & 2 & -1\\0 & 0 & 0& 0 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0\end{matrix}\right)u_2=0  \), pudiendo obtener \( u_2=\{0,0,1,2\} \)

Calculamos \( u_1: \)

\( (B-2I)u_2=u_1 \Rightarrow \left(\begin{matrix}1 & 0 & 2 & -1\\0 & 1 & 2& -1 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0\end{matrix}\right)u_2=0  \). Es decir, \( u_1=\{0,0,0,0\} \) y eso no podría, ser.

¿sabríais dónde estoy cometiendo el error?.

Saludos.

29 Mayo, 2024, 11:15 pm
Respuesta #19

Luis Fuentes

  • el_manco
  • Administrador
  • Mensajes: 58,907
  • País: es
  • Karma: +0/-0
Hola

  El error está en que no has respetado la observación en rojo que te indiqué en mi último mensaje. Tienes que elegir otro vector \( \vec u_2 \).

Saludos.