Matrices de Jordan
Qué hacer cuando una matriz no tiene bastantes vectores propios para diagonalizarse: se busca la matriz «casi diagonal» más sencilla que le es semejante, la forma de Jordan $J=P^{-1}AP$ (valores propios en la diagonal y unos justo encima). Aprenderás a leer cuántos bloques hay y de qué tamaño (con $\alpha$, $d$ y los rangos de $(A-\lambda I)^k$), a escribir $J$ y a construir la matriz de paso $P$ con cadenas de vectores. Al final, el método de las potencias (§4.4.3) para estimar el valor propio dominante, y un mapa que resume todo el tema 4.
0Antes de empezar U p.167-168
En el tema 4.3 vimos que una matriz $A$ de orden $n$ se diagonaliza cuando tiene $n$ vectores propios independientes, es decir, $n$ «ejes propios» que forman una base. Pero hay matrices a las que les faltan ejes: por ejemplo, un valor propio doble con un solo eje propio. En ese caso no se puede llegar a una matriz diagonal, pero sí a una casi diagonal: la que tiene los valores propios en la diagonal y, como mucho, unos «1» justo encima. Esa matriz se llama forma de Jordan de $A$. Cada «1» marca un eje que falta: en lugar de un vector propio, se usa un vector que $A$ no solo estira, sino que además empuja hacia un eje propio (como una baraja de cartas que se desliza al inclinarla).
Qué necesitas saber ya
- Del tema 4.2 (imprescindible): valores propios como raíces de $|A-\lambda I|$, multiplicidad algebraica $\alpha_i$, subespacio propio $L_i=\operatorname{Nuc}(A-\lambda_iI)$, multiplicidad geométrica $d_i=\dim L_i=n-\operatorname{rang}(A-\lambda_iI)$ y $1\le d_i\le\alpha_i$ (Teorema 4.5).
- Del tema 4.3: diagonalizable $\iff d_i=\alpha_i$ para todo $i$ (Teorema 4.7), la receta para $P$ y $\Lambda$ y la comprobación $AP=P\Lambda$ sin invertir $P$.
- Del tema 4.1: matrices semejantes ($B=P^{-1}AP$) tienen el mismo polinomio característico, traza, determinante y rango (Teorema 4.1, U p.146) y el mismo polinomio característico y traza (Teorema 4.4, U p.153).
- De los temas 2.3 y 3.2: base, independencia lineal, núcleo $\operatorname{Nuc}$ y $\dim\operatorname{Nuc}=n-\operatorname{rang}$. Del tema 1.2: rango y producto de matrices por columnas.
Mapa del tema
| Apartado | Qué contiene | Dónde |
|---|---|---|
| 1 La idea | Faltan ejes propios; ejemplo $2\times2$ y figura 3D | E 4.36 · P |
| 2 Teorema de Jordan | Bloques, $J_i=\lambda I+N$, nilpotencia, propiedades y avisos | U p.168-170 · E p.150 |
| 3 Qué $J$ son posibles | Tabla por $(\alpha,d)$; Ej. 4.16, 4.17, 4.22 | U p.169-172, 175-176, 179 |
| 4 Leer $J$ en los rangos | Regla $r_{k-1}-r_k$; subespacios propios generalizados (Def. 4.13); Ej. 4.20, 4.21, 4.23 | U p.173, 176-180 · P |
| 5 Matriz de paso | Cadenas, guion en 4 pasos, orden de columnas, Ej. 4.18-4.19, animación | U p.173-177 · E 4.38, 4.40 |
| 6 Método de potencias | Def. 4.14, algoritmo, Ej. 4.24-4.30 con tablas | U p.180-190 · L §5.8 |
| 7 Tabla Lay ↔ U | Lay apenas trata Jordan | L |
| 8 Ejercicios del libro | E 4.33–4.41 + 11 ejemplos propios + autoevaluación 9 y 10 | E pp.150-160 · U p.190 |
| 9 Ejemplos de examen | 4 ejercicios de exámenes anteriores + 3 propios + enlaces | X |
| 10 Chuleta | Fórmulas con hipótesis, lista de comprobación, errores | — |
| 11 Mapa del tema 4 | Qué hacer ante una matriz: de los valores propios a $\Lambda$ o $J$ y $P$ | P |
(a) Preguntas de teoría de verdadero/falso sobre $J$ (E 4.33-4.35, 4.37; autoevaluación U p.190, ítem 10). (b) Hallar $J$ de una matriz concreta, razonando con $\alpha$, $d$ y, si hace falta, los rangos de $(A-\lambda I)^k$ (E 4.36, 4.39; X Ene. 2020, Ene. 2022, Feb. 2025, Sept. 2024). (c) Hallar una matriz de paso $P$ con las cadenas (E 4.40; X Sept. 2024) o comprobar si una $P$ dada vale: $AP=PJ$ con $P$ regular (E 4.38). (d) Estimar el valor propio dominante con el método de potencias (E 4.41). Ningún examen escaneado (24 convocatorias) pregunta por el método de potencias ni por el valor propio dominante: es el apartado de menor prioridad. El examen exige desarrollar el razonamiento («$\alpha=\ldots$, $d=n-\operatorname{rang}(A-\lambda I)=\ldots$, luego hay $d$ bloques…»); Maxima no entra.
- U p.167-168: sugiere que la forma de Jordan resuelve también el caso «sin raíces reales» («la ecuación característica no tiene solución en $\mathbb R$… siempre se puede encontrar una matriz semejante casi diagonal»): es falso en $\mathbb R$ (giro de $90^\circ$: no hay $J$ real); solo vale sobre $\mathbb C$.
- U p.168, Teorema 4.8: enuncia la existencia de $J$ sin decir que todas las raíces de $|A-\lambda I|$ deben ser reales (si no, hay forma de Jordan solo sobre $\mathbb C$). U lo admite de pasada en la p.167 («no se puede diagonalizar»). Aquí se usa la forma segura.
- U p.170: los unos van en la «diagonal secundaria»; debe decir «línea inmediatamente superior a la diagonal principal» (E p.150 lo dice bien). La diagonal secundaria es la otra, la que va de la esquina superior derecha a la inferior izquierda.
- U p.173, Definición 4.13: los subespacios propios generalizados son los núcleos de $(A-\lambda I)^i$ con «$i\ge2$»; debería ser $i\ge1$ (la cadena empieza en $\operatorname{Nuc}(A-\lambda I)$, el subespacio propio).
- U p.175 (Ej. 4.18): «$\bar v_2$ pertenece a $(A-2I)^2=0$» debe leerse $\bar v_2\in\operatorname{Nuc}(A-2I)^2$ (la matriz $(A-2I)^2$ no es cero).
- U no define «base de Jordan» ni el orden de las columnas de $P$: se deduce de los Ej. 4.18, 4.19 y E 4.40 (vector propio primero, el resto de la cadena después). Es el convenio de aquí.
- U p.164, Teor. 4.7 («de rango $n$») debe leerse «de orden $n$» (ya avisado en 4.3). U p.169: escribe «nihilpotente» por «nilpotente».
- E 4.35: el libro da B «correcta» y D «falsa»; es al revés (B falsa, D correcta). E 4.36: «asociado a $\lambda=0$» debe ser $\lambda=2$. E 4.37-B: frase imprecisa. E 4.39: «$f(\bar e_4)=\bar e_3-2\bar e_4$» debe ser $+2\bar e_4$. E 4.40: escribe $P^{-1}=(\ldots)^{-1}$ con la matriz de $A$ (debe ser la de $P$). E 4.41: «$\bar x_4=(63/125,1)$» debe ser $\bar x_5$.
- U p.186 (Ej. 4.26): dice «$\bar x_k=(-1,1)$ para $k\ge1$» pero la tabla da $(1,-1)$ (el mismo vector con otro signo, sin consecuencias), y escribe «$\alpha_0=0$» para decir que $\bar x_0$ no tiene componente en $\bar v_1$ (se refiere al coeficiente de $\bar v_1$).
- Convenio: $J=P^{-1}AP\iff AP=PJ$ con $P$ regular; las columnas de $P$ forman la base de Jordan. Maxima está fuera de examen (E 4.36-D, 4.38-D); aquí no se usa.
1La idea: cuando faltan ejes propios U p.167 · E 4.36 · P
Una matriz diagonalizable solo estira cada eje propio. Una matriz con un valor propio doble y un solo eje propio hace algo más: sobre un eje estira, pero sobre otro vector estira y además lo empuja hacia el primer eje. Es una cizalla: imagina una baraja de cartas apilada y la empujas por arriba; cada carta se desplaza lateralmente una cantidad proporcional a su altura. La carta de abajo (el eje propio) no se mueve; las de arriba se desplazan hacia ella. Ese «empuje» es exactamente el «1» de la matriz de Jordan.
$A=\begin{pmatrix}3&1\cr-1&1\end{pmatrix}$. ¿Se puede diagonalizar? Si no, ¿cuánto se parece a una diagonal?
- Valores propios. $|A-\lambda I|=(3-\lambda)(1-\lambda)+1=\lambda^2-4\lambda+4=(\lambda-2)^2$: $\lambda=2$ con $\alpha=2$.
- Ejes propios. $A-2I=\begin{pmatrix}1&1\cr-1&-1\end{pmatrix}$ tiene rango $1$, luego $d=2-1=1$: un solo eje, $L=\langle p_1\rangle$ con $p_1=(1,-1)$. Como $d=1\lt2=\alpha$, no es diagonalizable (Teorema 4.7). Falta un eje.
- Completar la base con un vector cualquiera fuera del eje, por ejemplo $v=(1,0)$, y mirar qué hace $A$ con cada vector de la base $B'=\{p_1,v\}$: $$Ap_1=(2,-2)=2\,p_1,\qquad Av=(3,-1)=1\cdot(1,-1)+2\cdot(1,0)=1\cdot p_1+2\cdot v .$$ Con $v$ no hay solo estiramiento: $Av$ es «$2v$ más un empujón $p_1$».
- Matriz de $A$ en la base $B'$ (columnas = coordenadas de $Ap_1$ y de $Av$ en $B'$): $$J=\begin{pmatrix}2&1\cr0&2\end{pmatrix}.$$ Es lo más cerca de una diagonal que se puede estar. El «1» es el empujón.
- Con $P$. $P=(p_1\mid v)=\begin{pmatrix}1&1\cr-1&0\end{pmatrix}$ ($\det P=1\neq0$). Comprobación sin invertir: $AP=\begin{pmatrix}2&3\cr-2&-1\end{pmatrix}$ y $PJ=\begin{pmatrix}1\cdot2&1\cdot1+1\cdot2\cr-1\cdot2&-1\cdot1+0\end{pmatrix}=\begin{pmatrix}2&3\cr-2&-1\end{pmatrix}$ ✓, luego $J=P^{-1}AP$.
Aquí vale cualquier $v$ fuera de $L$ porque $(A-2I)^2=0$. En tamaños mayores habrá que elegir $v$ con cuidado (apartado 5).
E p.153 U p.167 propioJusto al revés: la forma de Jordan existe precisamente para las matrices que no se diagonalizan (siempre que sus raíces características sean reales). Y si $A$ sí es diagonalizable, su forma de Jordan es la matriz diagonal $\Lambda$ (todos los bloques $1\times1$, ningún «1»). Preguntar «halla la matriz de Jordan» en un examen puede tener como respuesta una diagonal (X Ene. 2020 y Ene. 2022, ejercicio 4).
La esfera unidad (gris) se transforma en su imagen (azul) por $A$. Las rectas y planos de colores son los subespacios propios. Prueba los casos: con «Doble con $d=2$» (diagonalizable) las rectas y el plano bastan para una base; con «Doble con $d=1$» o «Triple, $d=1$» faltan ejes y aparecen vectores «generalizados» (los que no son propios) de las cadenas (los comprenderás en el apartado 5). Con «Complejos» no hay ejes propios en un plano (giro): ni siquiera existe $J$ real. El deslizador $t$ deforma la esfera poco a poco. Si tu navegador no admite 3D, queda el texto con el veredicto.
P2Teorema de Jordan: bloques y nilpotencia U p.168-170 · E p.150
Un bloque de Jordan de tamaño $s$ y valor propio $\lambda$ es la matriz $s\times s$ $$J_s(\lambda)=\begin{pmatrix}\lambda&1&&\cr&\lambda&\ddots&\cr&&\ddots&1\cr&&&\lambda\end{pmatrix}$$ (triangular superior, $\lambda$ en toda la diagonal, $1$ en la línea inmediatamente superior y $0$ en el resto). Una matriz de Jordan es una matriz diagonal por bloques $J=\operatorname{diag}(J_1,\dots,J_s)$ formada por bloques de Jordan. Notación: $A\oplus B$ es la matriz diagonal por bloques $\operatorname{diag}(A,B)$ y $X\setminus Y$ son los elementos de $X$ que no están en $Y$. Ejemplo con $J_3(2)$ y $J_1(5)$: $$J=\begin{pmatrix}2&1&0&0\cr0&2&1&0\cr0&0&2&0\cr0&0&0&5\end{pmatrix}.$$
U p.168-170 E p.150Sea $A$ de orden $n$ con todas sus raíces características reales (hipótesis que U da por supuesta) y sea $s$ el número máximo de vectores propios linealmente independientes de $A$. Entonces $A$ es semejante a una matriz diagonal por bloques $$J=\operatorname{diag}(J_1,\dots,J_s),$$ donde cada $J_i$ es triangular superior con un único valor propio $\lambda_i$ en la diagonal y unos en la línea superior. $J$ es la forma canónica de Jordan de $A$. Es decir, existe $P$ regular con $J=P^{-1}AP$.
U p.168- $A$ y $J$ tienen los mismos valores propios (son semejantes: Teorema 4.4).
- $J$ es única salvo el orden de los bloques: si se cambian los bloques entre sí se obtiene otra matriz semejante a $A$ (corresponde a reordenar columnas de $P$ por grupos).
- Si todos los valores propios son distintos (o, en general, $d_i=\alpha_i$ para todo $i$), $J=\Lambda$ es diagonal. La forma de Jordan extiende la diagonalización.
- Cada bloque es $J_i=\lambda_iI+N$ con $N$ nilpotente: $N^m=0$ para algún $m$.
Para un bloque $3\times3$, $J_3(\lambda)=\lambda I+N$ con $N=\begin{pmatrix}0&1&0\cr0&0&1\cr0&0&0\end{pmatrix}$. Al multiplicar, los unos se alejan de la diagonal: $$N^2=\begin{pmatrix}0&0&1\cr0&0&0\cr0&0&0\end{pmatrix},\qquad N^3=0.$$ En general, un bloque de tamaño $s$ cumple $N^s=0$ y $\operatorname{rang}N^k=\max(s-k,0)$. Este dato es la clave del apartado 4: los rangos de las potencias «cuentan» los bloques. Ejemplo numérico: $4\begin{pmatrix}1&0\cr0&1\end{pmatrix}+\begin{pmatrix}0&1\cr0&0\end{pmatrix}=\begin{pmatrix}4&1\cr0&4\end{pmatrix}$ y $\begin{pmatrix}0&1\cr0&0\end{pmatrix}^2=0$.
U p.169 propioSi $|A-\lambda I|$ tiene raíces no reales, no hay forma de Jordan real. Ejemplo: $A=\begin{pmatrix}1&-2\cr1&-1\end{pmatrix}$ tiene $|A-\lambda I|=\lambda^2+1$, sin raíces reales. Si existiera $J$ real triangular con $J=P^{-1}AP$, tendría el mismo polinomio característico (Teorema 4.4) $(\lambda-j_{11})(\lambda-j_{22})$ con $j_{ii}$ reales: contradicción. Sobre $\mathbb C$ sí hay: $\lambda=\pm i$ distintos, luego $J=\operatorname{diag}(i,-i)$. Otro ejemplo: el giro de $90^\circ$. En los exámenes siempre se ha dado el caso de raíces reales; la nota oficial de Ene. 2020 («la matriz de Jordan siempre existe») solo es cierta con esa condición.
U p.167 propio- Los unos van encima de la diagonal (U, E, este tema). Con la base en orden contrario saldrían debajo: es la traspuesta, no la $J$ de U (error de la solución oficial de Sept. 2024, problema 8, 2.º apartado).
- «Diagonal secundaria» (U p.170) es una errata: se trata de la línea justo encima de la principal.
- Nº de bloques $\neq$ nº de valores propios distintos: se cuenta con vectores propios independientes, $\sum d_i$ (E 4.33, 4.37).
3¿Qué matrices $J$ son posibles? Número y tamaño de los bloques U p.169-172, 175-176, 179
Un valor propio $\lambda$ con multiplicidad algebraica $\alpha$ debe «ocupar» $\alpha$ posiciones de la diagonal de $J$. Las reparte en bloques: tantos como vectores propios independientes tenga, es decir, $d$. Piensa en $\alpha$ como el número de «personas» que hay que sentar y en $d$ como el número de «mesas» que hay: los tamaños de los bloques son cuántos se sientan en cada mesa. Cuantas más mesas ($d$ grande), menos unos hay. Con $d=\alpha$ cada uno en su mesa: diagonal.
- El nº de bloques de $\lambda_i$ es su multiplicidad geométrica $d_i=n-\operatorname{rang}(A-\lambda_iI)$; el nº total de bloques es $\sum d_i$ = nº máximo de vectores propios independientes.
- Los tamaños de los bloques de $\lambda_i$ suman $\alpha_i$. Los de un mismo $\lambda_i$ van seguidos y de mayor a menor.
- El nº de unos de $J$ es $n-\sum d_i$: tantos como vectores «no propios» hacen falta para completar una base (E 4.34).
Un valor propio con $\alpha=4$. Los bloques tienen tamaños que suman $4$ (en orden no creciente) y hay $d$ de ellos:
| $d$ | Tamaños de los bloques | Matriz $J$ (ese valor propio) |
|---|---|---|
| $1$ | $(4)$ | $J_4(\lambda)$ (tres unos) |
| $2$ | $(3,1)$ o $(2,2)$ | $J_3\oplus J_1$ o $J_2\oplus J_2$ (dos unos) |
| $3$ | $(2,1,1)$ | $J_2\oplus J_1\oplus J_1$ (un uno) |
| $4$ | $(1,1,1,1)$ | $\lambda I_4$ (diagonal) |
Con $\alpha=4$ y $d=2$ no basta con $d$: hay dos candidatas. Se decide con los rangos de las potencias (apartado 4).
U p.170-171| $\alpha$ | $d$ | Bloques | ¿Basta $d$ para decidir? |
|---|---|---|---|
| $1$ | $1$ | $(1)$ | sí |
| $2$ | $1$ / $2$ | $(2)$ / $(1,1)$ | sí |
| $3$ | $1$ / $2$ / $3$ | $(3)$ / $(2,1)$ / $(1,1,1)$ | sí (U p.172 y p.175-176) |
| $4$ | $1$ / $2$ / $3$ / $4$ | $(4)$ / $(3,1)$ ó $(2,2)$ / $(2,1,1)$ / $(1,1,1,1)$ | no, si $d=2$ |
- Tres valores propios distintos $\Rightarrow J=\Lambda$. Ejemplo 4.17: $A=\begin{pmatrix}1&2&3\cr0&4&5\cr0&0&6\end{pmatrix}$ es triangular, así que sus valores propios son la diagonal, $1,4,6$, distintos: $J=\operatorname{diag}(1,4,6)$ y $A$ es diagonalizable.
- Un valor propio doble $\lambda_2$ y uno simple: si $\operatorname{rang}(A-\lambda_2I)=1$ ($d=2$), $J$ es diagonal; si $\operatorname{rang}(A-\lambda_2I)=2$ ($d=1$), $J=\begin{pmatrix}\lambda_1&0&0\cr0&\lambda_2&1\cr0&0&\lambda_2\end{pmatrix}$ (un bloque $2\times2$). Se usa que $J-\lambda I=P^{-1}(A-\lambda I)P$ y las matrices semejantes tienen el mismo rango.
- Un valor propio triple $\lambda$: $\operatorname{rang}(A-\lambda I)=0,\,1,\,2$ da, respectivamente, $J=\lambda I$, $J'=\begin{pmatrix}\lambda&1&0\cr0&\lambda&0\cr0&0&\lambda\end{pmatrix}$ (bloques $2+1$) y $J''=\begin{pmatrix}\lambda&1&0\cr0&\lambda&1\cr0&0&\lambda\end{pmatrix}$ (un bloque $3\times3$).
Una matriz de orden $7$ con $|A-\lambda I|=(\lambda-4)^4(\lambda+3)^3$ (raíces $4$ con $\alpha=4$ y $-3$ con $\alpha=3$). El valor $4$ tiene entre $1$ y $4$ bloques y el $-3$ entre $1$ y $3$, así que el nº de bloques de $J$ está entre $1+1=2$ y $4+3=7$. Contando las formas: $4$ tiene $5$ posibilidades (las de arriba) y $-3$ tiene $3$ ($(3)$, $(2,1)$, $(1,1,1)$), luego hay $5\cdot3=15$ matrices $J$ distintas (salvo el orden de los bloques) compatibles con ese polinomio.
U p.179 propio4Leer $J$ en los rangos de $(A-\lambda I)^k$ U p.173, 176-180
Sea $N=A-\lambda I$. Dentro de un bloque de tamaño $s$, $N$ funciona como un ascensor que baja un piso: lleva la cabeza de la cadena $v$ al piso de abajo $Nv$, luego a $N^2v$, y así hasta el vector propio, que $N$ manda a $0$. Cuántas veces hay que aplicar $N$ para llegar a $0$ es el tamaño del bloque. Así, $\operatorname{Nuc}N$ (los vectores que caen a $0$ de golpe) son los del piso 1; $\operatorname{Nuc}N^2$ añade los del piso 2; etc.
Sea $N=A-\lambda I$ con $\lambda$ de multiplicidad algebraica $\alpha$ y $r_k=\operatorname{rang}N^k$ (con $r_0=n$). Se calculan las potencias hasta que $\dim\operatorname{Nuc}N^k=n-r_k=\alpha$. Entonces $$\#\{\text{bloques de tamaño}\ge k\}=r_{k-1}-r_k=\dim\operatorname{Nuc}N^k-\dim\operatorname{Nuc}N^{k-1}.$$ Los bloques de tamaño exactamente $k$ son (los $\ge k$) menos (los $\ge k+1$). El caso $k=1$ es $d=n-r_1$.
Por qué. En $J$, un bloque de tamaño $s$ aporta $\min(k,s)$ a $\dim\operatorname{Nuc}(J-\lambda I)^k$ (el rango de su potencia es $\max(s-k,0)$, apartado 2). Y $(A-\lambda I)^k=P(J-\lambda I)^kP^{-1}$ tiene el mismo rango que $(J-\lambda I)^k$ (matrices semejantes, U Teor. 4.1, p.146; ese es el argumento de U en la p.176).
propio U p.176Las dimensiones de $\operatorname{Nuc}N^k$ para $k=1,2,3,\dots$ son: para $(3,1)$: $2,\ 3,\ 4$ (el bloque de $3$ aporta $1,2,3$ y el de $1$ aporta $1,1,1$); para $(2,2)$: $2,\ 4,\ 4$. Con los dos casos $\operatorname{rang}N=2$ ($d=2$), pero $\operatorname{rang}N^2=1$ en el primero y $\operatorname{rang}N^2=0$ en el segundo. Es lo que ocurre con $A$ de E 4.39 (rangos $2,1,0$: bloques $(3,1)$) y con la matriz $B$ del ejemplo propio 10 (rangos $2,0$: bloques $(2,2)$; mismo polinomio, misma $d$ y no semejantes).
E 4.39 p.156 propioLos subespacios propios generalizados asociados a $\lambda$ son los núcleos $\operatorname{Nuc}(A-\lambda I)^i$, $i\ge1$ (U escribe $i\ge2$; con $i=1$ es el subespacio propio, el primer escalón). Forman una cadena creciente $$\operatorname{Nuc}(A-\lambda I)\subset\operatorname{Nuc}(A-\lambda I)^2\subset\cdots\subset\operatorname{Nuc}(A-\lambda I)^k,$$ que se estabiliza cuando la dimensión llega a $\alpha$ (no crece más). Sus vectores no nulos fuera del escalón anterior son los «vectores propios generalizados».
U p.173- Ej. 4.20. $A=\begin{pmatrix}0&1&2\cr0&0&1\cr0&0&0\end{pmatrix}$, triangular: $\lambda=0$ triple. $\operatorname{rang}(A-0I)=2$, $d=1$: un solo bloque, $J=\begin{pmatrix}0&1&0\cr0&0&1\cr0&0&0\end{pmatrix}$ (es $J''$; $A$ es nilpotente).
- Ej. 4.21. $A=\begin{pmatrix}0&-3&-1\cr2&5&1\cr-2&-3&1\end{pmatrix}$, $|A-\lambda I|=-(\lambda-2)^3$. $A-2I=\begin{pmatrix}-2&-3&-1\cr2&3&1\cr-2&-3&-1\end{pmatrix}$ tiene rango $1$ ($d=2$) y $(A-2I)^2=0$: bloques $(2,1)$, $J=\begin{pmatrix}2&1&0\cr0&2&0\cr0&0&2\end{pmatrix}$ (es $J'$).
- Ej. 4.23. Matriz de orden $5$ con $\lambda=0$ triple y $d=2$ (bloques $(2,1)$) y $\lambda=8$ doble y $d=1$ (bloque $2$): $J=\begin{pmatrix}0&1&0&0&0\cr0&0&0&0&0\cr0&0&0&0&0\cr0&0&0&8&1\cr0&0&0&0&8\end{pmatrix}$, es decir, $J_2(0)\oplus J_1(0)\oplus J_2(8)$.
5La matriz de paso: cadenas de vectores U p.173-177 · E 4.38, 4.40
Para cada bloque de tamaño $s$ hay que encontrar $s$ vectores tales que $A$ los vaya «empujando» hacia el vector propio. Se empieza por el más alto (la cabeza $v$), que debe estar «a $s$ pisos» del suelo, y se baja con $N$: $v$, $Nv$, $N^2v$, …, hasta el vector propio. Si eliges mal la cabeza (demasiado abajo) la cadena se queda corta y no llena el bloque.
Sea $N=A-\lambda I$ y sea $v$ con $N^3v=0$, $N^2v\neq0$ (bloque de tamaño $3$). Con las columnas $p_1=N^2v$, $p_2=Nv$, $p_3=v$: $$Ap_1=\lambda p_1,\qquad Ap_2=p_1+\lambda p_2,\qquad Ap_3=p_2+\lambda p_3 ,$$ porque $Ap_2=(N+\lambda I)Nv=N^2v+\lambda Nv$ y $Ap_3=Nv+\lambda v$. Las coordenadas de esos tres vectores en la base $\{p_1,p_2,p_3\}$ son justo las columnas $(\lambda,0,0)$, $(1,\lambda,0)$ y $(0,1,\lambda)$ de $J_3(\lambda)$. Por eso el vector propio va primero y la cabeza $v$ al final. En general, para un bloque de tamaño $k$: $P=(N^{k-1}v\mid\cdots\mid Nv\mid v)$.
propio E p.158-159- Tamaños de los bloques (apartado 4): $d$ y, si hace falta, los rangos de $N^k$, con $N=A-\lambda I$ para cada valor propio por separado.
- Cadena de núcleos: $\operatorname{Nuc}N\subset\operatorname{Nuc}N^2\subset\cdots$ hasta que la dimensión sea $\alpha$.
- Para cada bloque de tamaño $k$ (empezando por el mayor): elegir $v\in\operatorname{Nuc}N^k$ que sea independiente de $\operatorname{Nuc}N^{k-1}$ junto con las cabezas ya elegidas y los vectores de nivel $k$ de las cadenas más largas (no basta $v\notin\operatorname{Nuc}N^{k-1}$: con bloques $(2,2)$, tomar $h_2=h_1+u$ con $u$ propio da $P$ singular); calcular $Nv,\dots,N^{k-1}v$ (de derecha a izquierda). Los bloques de tamaño $1$ se rellenan con vectores propios independientes de los anteriores.
- Colocar las columnas bloque a bloque, en el mismo orden que en $J$, cada bloque como $(N^{k-1}v\mid\cdots\mid Nv\mid v)$. Comprobar $\det P\neq0$ y $AP=PJ$.
$A=\begin{pmatrix}1&0&3\cr-1&2&0\cr0&0&2\end{pmatrix}$.
- Polinomio (por la 3.ª fila): $|A-\lambda I|=(2-\lambda)\begin{vmatrix}1-\lambda&0\cr-1&2-\lambda\end{vmatrix}=(1-\lambda)(2-\lambda)^2$. $\lambda_1=1$ ($\alpha=1$) y $\lambda_2=2$ ($\alpha=2$).
- $\lambda_1=1$: $A-I=\begin{pmatrix}0&0&3\cr-1&1&0\cr0&0&1\end{pmatrix}\Rightarrow x_3=0,\ x_1=x_2$: vector propio $(1,1,0)$.
- $\lambda_2=2$: $N=A-2I=\begin{pmatrix}-1&0&3\cr-1&0&0\cr0&0&0\end{pmatrix}$, rango $2$, $d=1\lt2$: un bloque $2\times2$. $\operatorname{Nuc}N$: $x_1=0,\ x_3=0$, es $\langle(0,1,0)\rangle$ (dim $1$).
- Cadena de núcleos: $N^2=\begin{pmatrix}1&0&-3\cr1&0&-3\cr0&0&0\end{pmatrix}$, $\operatorname{Nuc}N^2$: $x_1=3x_3$, es $\langle(3,0,1),(0,1,0)\rangle$ (dim $2=\alpha$). Se para.
- Cabeza: $v_2=(3,0,1)\in\operatorname{Nuc}N^2$ y $\notin\operatorname{Nuc}N$ (U escribe «pertenece a $(A-2I)^2=0$»: errata de notación). Se baja: $Nv_2=(-3+3,\,-3,\,0)=(0,-3,0)$, vector propio ✓.
- Matriz de paso (propio de $\lambda=1$; luego el bloque: $Nv_2$, $v_2$): $$P=\begin{pmatrix}1&0&3\cr1&-3&0\cr0&0&1\end{pmatrix},\qquad J=\begin{pmatrix}1&0&0\cr0&2&1\cr0&0&2\end{pmatrix}.$$ $\det P=-3$ y $AP=PJ$ (columna 2: $A(0,-3,0)=(0,-6,0)=2(0,-3,0)$; columna 3: $A(3,0,1)=(6,-3,2)=(0,-3,0)+2(3,0,1)$ ✓).
$A=\begin{pmatrix}-1&3&0\cr0&-1&0\cr0&2&-1\end{pmatrix}$, triangular por bloques con $\lambda=-1$ triple.
- $N=A+I=\begin{pmatrix}0&3&0\cr0&0&0\cr0&2&0\end{pmatrix}$, rango $1\Rightarrow d=3-1=2$: bloques $(2,1)$ (orden 3, rango $1$: $J'$). $\operatorname{Nuc}N$: $x_2=0$, es $\langle(1,0,0),(0,0,1)\rangle$.
- $N^2=0$, luego $\operatorname{Nuc}N^2=\mathbb R^3$ (dim $3=\alpha$).
- Cabeza del bloque de $2$: $v_1=(0,1,0)\notin\operatorname{Nuc}N$. Se baja: $v_2=Nv_1=(3,0,2)$, vector propio.
- Vector del bloque $1\times1$: un propio independiente de $(3,0,2)$: $(1,0,0)$.
- Columnas $(Nv_1\mid v_1\mid(1,0,0))$: $$P=\begin{pmatrix}3&0&1\cr0&1&0\cr2&0&0\end{pmatrix},\qquad J=\begin{pmatrix}-1&1&0\cr0&-1&0\cr0&0&-1\end{pmatrix}.$$ $\det P=-2\neq0$ y $AP=PJ$.
- Cabeza en el núcleo equivocado. Si $v\in\operatorname{Nuc}N$ ($Nv=0$), la cadena no baja. Si $v\notin\operatorname{Nuc}N^k$, $N^{k-1}v$ no es propio. Ejemplo: $A=\begin{pmatrix}2&1&1\cr0&2&0\cr0&1&1\end{pmatrix}$ ($\lambda=2$ doble, $\lambda=1$): $v=(0,1,0)$ no está en $\operatorname{Nuc}N^2$ y da una $P$ con $P^{-1}AP$ no de Jordan; $v=(0,1,1)$ sí vale (ejemplo propio de examen 13).
- Orden de las columnas invertido: $(v\mid Nv)$ en lugar de $(Nv\mid v)$ da los unos debajo de la diagonal (error de la oficial de Sept. 2024).
- Dependencia: un vector propio extra proporcional a uno ya usado deja $P$ singular; comprueba $\det P\neq0$.
- Escalar mal: dentro de una cadena hay que escalar todos los vectores juntos ($p_2=Np_3$); no se puede cambiar solo uno (E 4.38).
- Mezclar valores propios: cada $\lambda$ se trabaja con su $A-\lambda I$.
Se pueden cambiar la cabeza $v$ por cualquier otro vector de $\operatorname{Nuc}N^k\setminus\operatorname{Nuc}N^{k-1}$ (y los de los bloques de tamaño $1$ por otros propios independientes), escalar y reordenar bloques (con su $J$). Todas las respuestas correctas valen en el examen. Para comprobar una $P$ dada: (1) $\det P\neq0$; (2) $AP=PJ$ columna a columna (E 4.38). En E 4.38, con $J=\operatorname{diag}(J_1(2),J_2(-2))$, todas las $P=(a\,u\mid b\,w\mid b\,v+c\,w)$ con $a,b\neq0$ y $c$ real son válidas (tres parámetros libres).
E p.155Escribe una matriz $A$ (de $2\times2$ a $4\times4$, enteros) o elige un caso del libro y recorre la receta: polinomio característico y valores propios, rangos de $(A-\lambda I)^k$ y tamaños de los bloques, matriz $J$, cadenas de vectores con la cabeza resaltada y la matriz de paso con la comprobación $AP=PJ$ en fracciones exactas. Prueba: «U 4.18» (bloques 1 y 2), «E 4.39» (bloques 3 y 1), «4×4: bloques 2 y 2» (mismo polinomio y $d$ que E 4.39, pero otra $J$), «Giro 90°» y «Raíces complejas» (no hay $J$ real) y «raíces irracionales» (diagonalizable con valores irracionales).
PIgual que en 4.3 ($A^k=P\Lambda^kP^{-1}$), con $J=\lambda I+N$ y $\lambda I$ conmutando con $N$ se puede usar el binomio de Newton, y como $N^s=0$ solo sobreviven los primeros términos: $$J_2(\lambda)^k=\begin{pmatrix}\lambda^k&k\lambda^{k-1}\cr0&\lambda^k\end{pmatrix},\qquad J_3(\lambda)^k=\begin{pmatrix}\lambda^k&k\lambda^{k-1}&\binom k2\lambda^{k-2}\cr0&\lambda^k&k\lambda^{k-1}\cr0&0&\lambda^k\end{pmatrix}.$$ Comprobación: $J_2(2)^3=\begin{pmatrix}8&12\cr0&8\end{pmatrix}$. Aplicación completa en el ejemplo propio 8. Ningún examen escaneado lo pide; es una aplicación útil (sistemas de ecuaciones recurrentes y diferenciales).
propio6Técnicas iterativas: el método de potencias U p.180-190 · L p.319-323
Hallar los valores propios exactos de una matriz grande es difícil (no hay fórmula general para el grado $\ge5$). Si solo interesa el mayor en valor absoluto, basta con una receta muy simple: coger un vector, multiplicarlo por $A$ una y otra vez. Como cada multiplicación estira cada eje propio por su valor propio, el eje con el valor propio más grande crece más rápido que los demás y acaba dominando la dirección del vector: es como una población en la que una especie se reproduce más deprisa que las otras. Se normaliza en cada paso (dividiendo por la coordenada mayor) para que los números no se disparen, y esa coordenada va dando el valor propio.
Si $A$ es de orden $n$ con valores propios $\lambda_1,\dots,\lambda_n$, $\lambda_1$ es estrictamente dominante si $|\lambda_1|\gt|\lambda_2|\ge\cdots\ge|\lambda_n|$: tiene mayor valor absoluto que cualquier otro.
U p.181 L §5.8 p.319Si $A$ es diagonalizable con base de vectores propios $v_1,\dots,v_n$ y $x_0=c_1v_1+\cdots+c_nv_n$ (aquí $c_i$ son coeficientes, no multiplicidades), entonces $A^kx_0=c_1\lambda_1^kv_1+\cdots+c_n\lambda_n^kv_n$. Dividiendo por $\lambda_1^k$: $$\frac1{\lambda_1^k}A^kx_0=c_1v_1+c_2\Big(\frac{\lambda_2}{\lambda_1}\Big)^kv_2+\cdots+c_n\Big(\frac{\lambda_n}{\lambda_1}\Big)^kv_n\ \longrightarrow\ c_1v_1,$$ porque $|\lambda_i/\lambda_1|\lt1$. Si $c_1\neq0$, la dirección de $A^kx_0$ tiende a la de $v_1$. La rapidez depende de $|\lambda_2/\lambda_1|$: si es casi $1$, converge muy despacio. (U lo supone diagonalizable «para simplificar»; no es necesario para aplicar el método.)
U p.181-185 L §5.8 p.319-322- Elegir un vector inicial $x_0\neq0$.
- Para $k=0,1,2,\dots$: calcular $Ax_k$; sea $\mu_k$ la coordenada de $Ax_k$ de mayor valor absoluto (con su signo); calcular $x_{k+1}=\dfrac1{\mu_k}Ax_k$ (así $x_{k+1}$ tiene una coordenada igual a $1$).
- Si $x_0$ es adecuado: $\mu_k\to\lambda_1$ (el valor propio estrictamente dominante) y $x_k\to v_1$ (un vector propio suyo). Si no converge, probar otro $x_0$.
Solo estima el valor propio dominante. Para el de menor módulo se aplica el método a $A^{-1}$ (Ej. 4.30); la variante general (desplazamientos $(A-\alpha I)^{-1}$) es el «método de potencias inverso» (U p.190; L p.322-323). No sale en ningún examen escaneado.
U p.185, 190 L §5.8 p.321-323Valores propios $7$ y $1$, con vectores $(5,1)$ y $(1,-1)$: el dominante es $7$. Sin escalar (Ej. 4.24) se ve cómo gira la dirección ($x_k=A^kx_0$); con escalado (Ej. 4.25) se lee el valor propio:
| $k$ | $x_k=A^kx_0$ (4.24) | pendiente $y/x$ | $Ax_k$ (con escalado, 4.25) | $\mu_k$ | $x_{k+1}$ |
|---|---|---|---|---|---|
| $0$ | $(0,1)$ | — | $(5,\,2)$ | $5$ | $(1,\,0.4)$ |
| $1$ | $(5,2)$ | $0.4$ | $(8,\,1.8)$ | $8$ | $(1,\,0.225)$ |
| $2$ | $(40,9)$ | $0.225$ | $(7.125,\,1.45)$ | $7.125$ | $(1,\,0.2035)$ |
| $3$ | $(285,58)$ | $0.2035$ | $(7.0175,\,1.4070)$ | $7.0175$ | $(1,\,0.2005)$ |
| $4$ | $(2000,401)$ | $0.2005$ | $(7.0025,\,1.4010)$ | $7.0025$ | $(1,\,0.20007)$ |
La pendiente tiende a $1/5=0.2$ (la de $(5,1)$) y $\mu_k\to7$. Convergencia rápida porque $|\lambda_2/\lambda_1|=1/7$. Comprobación: $A(1,0.2)=(7,1.4)=7(1,0.2)$.
U p.182-185 L §5.8 Ej.2 p.321- Ej. 4.26 (p.186): misma $A$ con $x_0=(2,-2)$, que es vector propio del valor $1$ (no dominante, no tiene componente en $v_1$). La tabla: $Ax_0=(2,-2)$, $\mu_0=2$, $x_1=(1,-1)$; después $Ax_k=(1,-1)$, $\mu_k=1$ para siempre. El método converge a un vector propio y a $1$, no a $7$: mal $x_0$.
- Ej. 4.27 (p.187): con $x_0=(1,-5)$ (sí tiene componente en $v_1$) vuelve a converger: $\mu_k=-19$, $8.3684$, $7.1635$, $7.0228$, $7.0033$ y $\to7$ y $x_k\to(1,0.2)$. Fíjate que $\mu_0=-19$ es negativo y aun así el método converge: lo importante es la componente en $v_1$.
Ej. 4.28 (pp.187-188) — $A=\begin{pmatrix}4&-1&1\cr-1&3&-2\cr1&-2&3\end{pmatrix}$, $x_0=(1,0,0)$ (valores propios $1,3,6$; el dominante es $6$, con $(1,-1,1)$):
| $k$ | $Ax_k$ | $\mu_k$ | $x_{k+1}$ |
|---|---|---|---|
| $0$ | $(4,-1,1)$ | $4$ | $(1,-0.25,0.25)$ |
| $1$ | $(4.5,-2.25,2.25)$ | $4.5$ | $(1,-0.5,0.5)$ |
| $2$ | $(5,-3.5,3.5)$ | $5$ | $(1,-0.7,0.7)$ |
| $3$ | $(5.4,-4.5,4.5)$ | $5.4$ | $(1,-0.8333,0.8333)$ |
| $4$ | $(5.667,-5.167,5.167)$ | $5.667$ | $(1,-0.9118,0.9118)$ |
| $5$ | $(5.824,-5.559,5.559)$ | $5.824$ | $(1,-0.9545,0.9545)$ (*) |
| $6$ | $(5.909,-5.773,5.773)$ | $5.909$ | $(1,-0.9769,0.9769)$ |
(*) U escribe $-0.9546$ y $0.9545$ en esta fila; el valor exacto es $0.95455\ldots$ en ambas coordenadas. Tras 7 iteraciones $\lambda\approx5.909$ y $v\approx(1,-0.9769,0.9769)$; la tabla sugiere $6$ y $(1,-1,1)$, que se confirma: $A(1,-1,1)^t=(6,-6,6)^t=6(1,-1,1)^t$. Aquí $|\lambda_2/\lambda_1|=3/6$: el error se reduce a la mitad en cada paso.
- Ej. 4.29 (p.188): $A=\begin{pmatrix}1&2\cr4&3\end{pmatrix}$, valores propios $-1$ y $5$ (el dominante es $5$, con vector propio $(1,2)$), $x_0=(-1,1)$, que es vector propio de $-1$: $Ax_0=(1,-1)$, y en $Ax_0$ las dos coordenadas tienen el mismo valor absoluto (empate); con $\mu_k=1$ los $x_k$ oscilan entre $(1,-1)$ y $(-1,1)$ y no convergen a ningún vector. Conclusión de U: $x_0$ no es adecuado, hay que elegir otro (p. ej. $(0,1)$).
- Ej. 4.30 (p.189): $A=\begin{pmatrix}4&3\cr1&2\end{pmatrix}$ (valores propios $1$ y $5$): el método con $A$ estima el $5$. Para estimar el $1$ se aplica a $A^{-1}=\begin{pmatrix}2/5&-3/5\cr-1/5&4/5\end{pmatrix}$ (existe porque $|A|=5\neq0$), cuyos valores propios son $1$ y $1/5$ (si $Ax=\lambda x$, $A^{-1}x=\lambda^{-1}x$; los vectores propios son los mismos). Con $x_0=(0,1)$: $A^{-1}x_0=(-0.6,0.8)$, $\mu_0=0.8$, $x_1=(-0.75,1)$; luego $\mu_k=0.95,\ 0.9895,\ 0.9979,\ 0.9996\to1$ y $x_k\to(-1,1)$. Comprobación: $A(-1,1)=(-1,1)$.
- No hay valor propio estrictamente dominante ($|\lambda_1|=|\lambda_2|$, o valores complejos conjugados): $\mu_k$ no converge.
- $x_0$ sin componente en $v_1$ (Ej. 4.26, 4.29): converge a otro vector propio o oscila. El redondeo a veces «arregla» el problema, pero no cuentes con ello.
- Empates en $\mu_k$ (dos coordenadas con el mismo valor absoluto): el signo escogido puede cambiar el resultado (Ej. 4.29).
- Lentitud: si $|\lambda_2/\lambda_1|$ es casi $1$, hacen falta muchas iteraciones.
- Signo: $\mu_k$ es la coordenada con su signo, no su valor absoluto (E 4.41): si no, $x_k$ alterna de signo y $\mu_k$ converge al valor opuesto.
- El método estima valores propios sin necesidad de que $A$ sea diagonalizable (U lo supone solo para explicar la idea).
7Tabla Lay ↔ U: vocabulario y alcance L §5.4 p.292-293 · §5.8 p.319-323
Lay casi no trata la forma de Jordan: la menciona en un ejemplo y una nota al pie, sin teorema ni método de construcción. Para el método de potencias, en cambio, Lay (§5.8) es una buena explicación alternativa; cuando algo difiere de U, manda U.
| Concepto | U (oficial) | Lay |
|---|---|---|
| Forma de Jordan | Teor. 4.8 (p.168); $J=P^{-1}AP$ con unos encima de la diagonal | §5.4, Ej. 4 (p.292): matriz con $(\lambda+2)^2$ y $d=1$; la $B$-matriz de $x\mapsto Ax$ es triangular, «forma de Jordan»; nota al pie: «toda matriz cuadrada es similar a una en forma de Jordan; la base consiste en vectores propios y vectores propios generalizados» (remite a Noble-Daniel, cap. 9) |
| Vectores propios generalizados | Def. 4.13 (p.173) y guion en cadenas | Solo el nombre, en la nota al pie (p.292) |
| Matriz de paso por cadenas | Guion p.173; Ej. 4.18-4.19 | No. Útil: nota numérica (p.293): $P^{-1}AP$ se obtiene reduciendo $[P\ \ AP]$ a $[I\ \ P^{-1}AP]$ sin calcular $P^{-1}$ |
| Semejanza y base | $J$ es la matriz de $f$ en la base de Jordan | «$B$-matriz» de $T$ cuando $B$ son las columnas de $P$: $[T]_B=P^{-1}AP$ (§5.4, Teor. 8 y p.292) |
| Método de potencias | §4.4.3: Def. 4.14, algoritmo p.185, Ej. 4.24-4.30 | §5.8 (p.319-323): mismo método con $x_0$ de entrada mayor $1$ y $\mu_k$ la entrada mayor; su Ejemplo 2 es el 4.24-4.25 de U (misma $A$ y $x_0$) |
| Valor propio menor | Ej. 4.30: método sobre $A^{-1}$ (p.189) | Método de potencias inverso con desplazamiento $(A-\alpha I)^{-1}$ (p.322-323); con $\alpha=0$ coincide con U; sugiere resolver $(A-\alpha I)y_k=x_k$ en lugar de invertir |
No cites matrices de Lay para $J$: no es el libro del curso y no desarrolla el tema. En el examen se usa la notación y el convenio de U (unos encima, vector propio primero en cada bloque).
8Ejercicios del libro, resueltos E p.150-160
Los Ejercicios 4.33 a 4.41 de E son los de esta sección: del 4.33 al 4.40 tratan la forma de Jordan y la matriz de paso; el 4.41 es del método de potencias. Los apartados con Maxima (4.36-D, 4.38-D) están fuera de examen; se razonan a mano. Orden de estudio recomendado: 4.33 → 4.34 → 4.35 (conceptos), 4.36 → 4.37 (hallar $J$), 4.38 (comprobar una $P$), 4.39 → 4.40 ($J$ y $P$ de una $4\times4$) y 4.41 (potencias). Pulsa «Mostrar soluciones» arriba para abrirlas todas.
Bloque I — Teoría de $J$ (E 4.33–4.35)
Ejercicio 4.33 — verdadero o falso sobre bloques y semejanza E p.151
Si $J$ es la matriz de Jordan de una matriz $A$ cualquiera, responder verdadero o falso: A) El número de bloques de $J$ es superior al de vectores propios linealmente independientes. B) Al cambiar de orden los bloques de $J$ se obtiene una matriz que no es semejante a $A$. C) Cada bloque de $J$ es la suma de un múltiplo de la matriz unidad y una matriz nihilpotente.
Qué se pide. Traducir cada frase a «$d_i$ = nº de bloques de $\lambda_i$» y a $J_i=\lambda I+N$.
- Qué me piden
- Verdadero o falso sobre bloques y semejanza de la matriz de Jordan.
- Cómo lo pienso
- Cada bloque corresponde a un vector propio independiente: número de bloques $=$ número de vectores propios independientes. Reordenar los bloques no cambia la clase de semejanza (es una matriz de paso con una permutación).
- Qué uso
- Forma de Jordan (Def. 4.13-4.14, U p.167): si todos los valores propios son reales, $A$ es semejante a $J$, diagonal por bloques $J_s(\lambda)$ (valor propio en la diagonal y unos justo encima). Para cada $\lambda$: el nº de bloques es $d(\lambda)=n-\operatorname{rang}(A-\lambda I)$ y la suma de sus tamaños es $\alpha(\lambda)$.
- Por dónde empiezo
- Contrasto cada frase con la definición: bloques, orden, estructura $\lambda I+N$.
A) Falsa. Por la construcción (U p.169) el nº de bloques de $\lambda_i$ es $d_i$ y el total $\sum d_i$ es el nº máximo de vectores propios independientes: es igual, no superior. Ejemplo:
$$A=\begin{pmatrix}2&1\cr0&2\end{pmatrix}$$tiene $d=2-\operatorname{rang}(A-2I)=2-1=1$: un bloque y un vector propio independiente.
- B) Falsa. U p.169: al permutar bloques la matriz sigue siendo de la misma clase de semejanza que $A$. Reordenar los bloques equivale a reordenar por grupos las columnas de $P$: si $J=P^{-1}AP$ y $\tilde P$ es $P$ con los grupos permutados, $\tilde J=\tilde P^{-1}A\tilde P$ es $J$ con los bloques permutados, semejante a $A$ por definición. Ejemplo: $\operatorname{diag}(J_2(3),5)$ y $\operatorname{diag}(5,J_2(3))$ son semejantes.
C) Verdadera. $J_i=\lambda_iI+N$ con $N$ con unos encima de la diagonal y ceros en el resto, nilpotente porque cada potencia desplaza los unos y, para un bloque de tamaño $s$, $N^s=0$. Ejemplo:
$$\begin{pmatrix}4&1\cr0&4\end{pmatrix}=4I+\begin{pmatrix}0&1\cr0&0\end{pmatrix}$$ $$\begin{pmatrix}0&1\cr0&0\end{pmatrix}^2=0.$$
Respuesta: A falsa, B falsa, C verdadera (coincide con E).
Error típico: pensar que el orden de los bloques cambia la clase de semejanza: solo cambia el orden de las columnas de $P$. Teoría: U Teor. 4.8 p.168; propiedades p.169; construcción pp.169-170. (E escribe «de orden» repetido en B: errata tipográfica.)
Ejercicio 4.34 — ¿a qué es igual el número de unos de $J$? E p.151
Si $J$ es la matriz de Jordan de una matriz $A$ cualquiera, responder verdadero o falso a si el número de unos que hay por encima de la diagonal principal de $J$ es: A) Igual al número de vectores propios linealmente independientes. B) Igual al número de bloques. C) Igual al número de vectores no propios que se necesitan para formar una base de vectores que determinan la matriz de paso. D) Igual al rango de $A$.
Qué se pide. Contar unos, bloques y vectores.
- Qué me piden
- A qué es igual el número de unos de $J$.
- Cómo lo pienso
- Un bloque de tamaño $s$ tiene $s-1$ unos. Sumando los tamaños da el orden $n$: unos $=n-$ nº de bloques.
- Qué uso
- Forma de Jordan (Def. 4.13-4.14, U p.167): si todos los valores propios son reales, $A$ es semejante a $J$, diagonal por bloques $J_s(\lambda)$ (valor propio en la diagonal y unos justo encima). Para cada $\lambda$: el nº de bloques es $d(\lambda)=n-\operatorname{rang}(A-\lambda I)$ y la suma de sus tamaños es $\alpha(\lambda)$. Un bloque de tamaño $s$ aporta $s-1$ unos.
- Por dónde empiezo
- Cuento unos por bloques y relaciono con cadenas de vectores.
- Cuenta clave. Un bloque de tamaño $s$ tiene $s-1$ unos. Con $b$ bloques de tamaños $s_1,\dots,s_b$ y $\sum s_i=n$, el número de unos es $\sum(s_i-1)=n-b$. Y $b=\sum d_i$ (U p.169).
- A) Falsa en general: unos $=n-b$ y vectores propios independientes $=b$; solo coinciden si $n=2b$. Contraejemplo: $A=I_2$, $0$ unos y $2$ vectores propios.
- B) Falsa: el nº de bloques es $b$, mismo contraejemplo ($b=2$, $0$ unos).
- C) Verdadera. En la base de Jordan cada bloque aporta un vector propio (el primero de su cadena) y $s_i-1$ vectores no propios; total de no propios $\sum(s_i-1)=n-b$, igual al nº de unos. Es lo que dice U p.170: «tantos 1 como vectores no propios se necesitan hasta completar una base con los vectores propios».
- D) Falsa. $A=I_n$ tiene rango $n$ y $J=I_n$ no tiene unos. (Para una matriz nilpotente sí coinciden, porque $\operatorname{rang}A=n-d=n-b$; el enunciado dice «cualquiera».)
Respuesta: A falsa, B falsa, C verdadera, D falsa (coincide con E).
Error típico: confundir «nº de unos» con «nº de bloques»: son complementarios hasta $n$. Teoría: U pp.169-170.
Ejercicio 4.35 — ¿cuándo es diagonal $J$? (errata en B y D) E p.152
Si $J$ es la matriz de Jordan de una matriz $A$ cualquiera, justificar si: A) Los elementos de $J$ que están encima de la diagonal principal son siempre 1. B) $J$ es diagonal cuando el número de valores propios distintos es igual al rango de $A$. C) Los elementos que están encima de la diagonal principal pueden ser valores distintos a 0 y distintos a 1 para alguna matriz $A$. D) Ninguna de las anteriores es correcta.
Qué se pide. Refutar con contraejemplos las tres primeras.
- Qué me piden
- Cuándo es diagonal $J$ y qué se puede decir de los elementos fuera de la diagonal.
- Cómo lo pienso
- $J$ es diagonal si todos los bloques son de tamaño 1, es decir, si $A$ es diagonalizable. Los elementos de encima de la diagonal son 0 o 1, no valores propios.
- Qué uso
- Forma de Jordan (Def. 4.13-4.14, U p.167): si todos los valores propios son reales, $A$ es semejante a $J$, diagonal por bloques $J_s(\lambda)$ (valor propio en la diagonal y unos justo encima). Para cada $\lambda$: el nº de bloques es $d(\lambda)=n-\operatorname{rang}(A-\lambda I)$ y la suma de sus tamaños es $\alpha(\lambda)$.
- Por dónde empiezo
- Contrasto cada opción con esa idea.
- A) Falsa. Si $J$ es diagonal hay ceros encima de la diagonal: por ejemplo $A=\operatorname{diag}(1,2)$, valores propios distintos, $J=\Lambda$. Además entre dos bloques distintos siempre hay un $0$.
B) Falsa (E dice «correcta»: errata). Contraejemplo
$$A=\begin{pmatrix}0&1\cr0&0\end{pmatrix}:$$$|A-\lambda I|=\lambda^2$, un único valor propio $0$, y $\operatorname{rang}A=1$: «nº de valores propios distintos = rango» se cumple ($1=1$). Pero $d=2-\operatorname{rang}A=1\lt2=\alpha$, no es diagonalizable y $J=A$ no es diagonal. Tampoco es necesario: $\operatorname{diag}(1,2,2)$ es diagonal, con $2$ valores distintos y rango $3$. La justificación de E no se apoya en ningún teorema; lo cierto (U p.169) es que $n$ valores propios distintos en una matriz de orden $n$ da diagonalizable (de ahí probablemente la confusión con «de rango $n$» de U Teor. 4.7).
- C) Falsa. Por la construcción (U pp.168-170) los elementos encima de la diagonal solo pueden ser $0$ o $1$.
- D) Verdadera, porque A, B y C son falsas.
Respuesta: D (y no B, como dice E).
Error típico: mezclar «rango» y «orden». Teoría: U Teor. 4.7 p.164; pp.168-170.
Bloque II — Hallar $J$ (E 4.36–4.38)
Ejercicio 4.36 — $2\times2$ con un valor propio doble (errata en el valor propio) E p.153
Si $A$ es la matriz de filas $(3,1)$, $(-1,1)$, justificar si: A) Es semejante a una matriz diagonal. B) Existe una base de $\mathbb R^2$ formada por vectores propios de $A$. C) Se puede calcular su matriz $J$ de Jordan. D) Para obtener $J$ con MAXIMA hay que evaluar load(diag); jordan(); dispJordan();.
Qué se pide. $A=\begin{pmatrix}3&1\cr-1&1\end{pmatrix}$: decidir con $\alpha$ y $d$, y hallar $J$ y $P$.
- Qué me piden
- Para una $2\times2$ con valor doble: semejanza a diagonal, base de vectores propios y existencia de $J$.
- Cómo lo pienso
- Calculo si el valor doble tiene dos vectores propios independientes; si no, hay un único bloque $2\times2$ y la matriz no es diagonalizable.
- Qué uso
- Forma de Jordan (Def. 4.13-4.14, U p.167): si todos los valores propios son reales, $A$ es semejante a $J$, diagonal por bloques $J_s(\lambda)$ (valor propio en la diagonal y unos justo encima). Para cada $\lambda$: el nº de bloques es $d(\lambda)=n-\operatorname{rang}(A-\lambda I)$ y la suma de sus tamaños es $\alpha(\lambda)$.
- Por dónde empiezo
- Polinomio y rango de $A-\lambda I$.
Polinomio:
$$|A-\lambda I|=(3-\lambda)(1-\lambda)+1=\lambda^2-4\lambda+4=(\lambda-2)^2:$$$\lambda=2$, $\alpha=2$.
$d$:
$$A-2I=\begin{pmatrix}1&1\cr-1&-1\end{pmatrix}$$de rango $1$, $d=1$. $L=\langle(1,-1)\rangle$, asociado a $\lambda=2$ (E escribe $\lambda=0$: errata).
- A) Falsa: $d=1\neq\alpha=2$ (Teor. 4.7). B) Falsa: por el Teor. 4.6 diagonalizable equivale a base de vectores propios; todos los vectores propios son múltiplos de $(1,-1)$.
C) Verdadera. Las raíces son reales, luego existe $J$ (Teor. 4.8). $d=1$: un solo bloque de tamaño $2$:
$$J=\begin{pmatrix}2&1\cr0&2\end{pmatrix}.$$Matriz de paso (guion de U p.173): $(A-2I)^2=0$, así que se elige $v=(1,0)\notin\operatorname{Nuc}(A-2I)$ y $Nv=(1,-1)$:
$$P=\begin{pmatrix}1&1\cr-1&0\end{pmatrix},$$$\det P=1$,
$$AP=\begin{pmatrix}2&3\cr-2&-1\end{pmatrix}=PJ$$(columna a columna: $A(1,-1)=2(1,-1)$ y $A(1,0)=(3,-1)=(1,-1)+2(1,0)$).
- D) Maxima, fuera de examen. E la da por correcta; U (p.169) escribe
jordan(A); dispJordan(%);: los paréntesis vacíos del enunciado no son sintaxis literal.
Respuesta: A falsa, B falsa, C verdadera, D verdadera (como E).
Error típico: concluir «no tiene forma de Jordan» porque no es diagonalizable: el teorema de Jordan existe para ese caso. Teoría: U Teor. 4.4 p.153, 4.6 p.162, 4.7 p.164, 4.8 p.168, guion p.173. (Verificado con sympy.)
Ejercicio 4.37 — $3\times3$ diagonalizable: «dos bloques» es falso E p.154
Si $A$ es la matriz de filas $(2,2,-5)$, $(3,7,-15)$, $(1,2,-4)$, justificar si: A) La matriz de Jordan de $A$ posee solo dos bloques. B) $S=\{x=\alpha,y=\beta,z=0\mid\alpha,\beta\in\mathbb R\}$ es un subespacio de vectores propios de $A$. C) $A$ es semejante a la matriz $B$ de filas $(3,0,0)$, $(0,1,0)$ y $(0,0,3)$. D) Ninguna de las anteriores.
Qué se pide. Hallar los valores propios y comprobar $d_i=\alpha_i$ antes de buscar bloques.
- Qué me piden
- Para una $3\times3$, cuántos bloques tiene $J$ y si es semejante a una diagonal dada.
- Cómo lo pienso
- El número de bloques por valor propio es $d(\lambda)$. Calculo valores propios y rangos de $A-\lambda I$.
- Qué uso
- Forma de Jordan (Def. 4.13-4.14, U p.167): si todos los valores propios son reales, $A$ es semejante a $J$, diagonal por bloques $J_s(\lambda)$ (valor propio en la diagonal y unos justo encima). Para cada $\lambda$: el nº de bloques es $d(\lambda)=n-\operatorname{rang}(A-\lambda I)$ y la suma de sus tamaños es $\alpha(\lambda)$. Tamaños de los bloques: el nº de bloques de tamaño $\ge k$ es $\operatorname{rang}(A-\lambda I)^{k-1}-\operatorname{rang}(A-\lambda I)^{k}$.
- Por dónde empiezo
- Polinomio y $d$ de cada valor.
- Polinomio (Sarrus): $p(\lambda)=-\lambda^3+5\lambda^2-7\lambda+3$. Comprobación rápida: coeficiente de $\lambda^2$ = traza $2+7-4=5$ y $p(0)=|A|=3$. Como $p(1)=0$, Ruffini: $p(\lambda)=-(\lambda-1)^2(\lambda-3)$. $\lambda_1=3$ ($\alpha=1$), $\lambda_2=1$ ($\alpha=2$).
$d$:
$$A-I=\begin{pmatrix}1&2&-5\cr3&6&-15\cr1&2&-5\end{pmatrix}$$con las tres filas proporcionales: rango $1$, $d(1)=3-1=2=\alpha$. $L_1:\ x+2y-5z=0$. Y $\lambda=3$: $d=1=\alpha$ ($1\le d\le\alpha=1$).
- A) Falsa. $d_i=\alpha_i$ para los dos valores propios $\Rightarrow$ diagonalizable (Teor. 4.7) $\Rightarrow J=\Lambda=\operatorname{diag}(3,1,1)$, con tres bloques $1\times1$ (nº de bloques $=d(3)+d(1)=1+2=3$).
- B) Falsa. «Subespacio de vectores propios» significa que todos sus vectores no nulos son propios de un mismo valor propio, es decir $S\subset L_1$ o $S\subset L_3$. Pero $(1,0,0)\in S$ y $A(1,0,0)=(2,3,1)$ no es múltiplo de $(1,0,0)$. Matiz (imprecisión de E; en E, $L_1$ es el subespacio de $\lambda_1=3$ y $L_2$ el de $\lambda_2=1$): E dice que los vectores $(\alpha,\beta,0)$ «no pertenecen a $L_1$ ni a $L_2$»; no es cierto para todos: $(-2,1,0)\in S$ cumple $A(-2,1,0)=(-2,1,0)$, es propio de $\lambda=1$. La conclusión es correcta; la frase debe ser «no todos».
- C) Falsa. Semejantes tienen el mismo determinante y la misma traza: $|A|=3\neq|B|=3\cdot1\cdot3=9$ (y $5\neq7$). $A$ es semejante a $\operatorname{diag}(3,1,1)$, no a $\operatorname{diag}(3,1,3)$.
- D) Verdadera.
Respuesta: D.
Error típico: responder «dos bloques» porque hay dos valores propios distintos: los bloques se cuentan por vectores propios independientes. Teoría: U p.153; Teor. 4.5 p.158; Teor. 4.7 p.164; p.169. (Verificado con sympy: jordan_form da $\operatorname{diag}(1,1,3)$.)
Ejercicio 4.38 — ¿es $P$ una matriz de paso? E p.155
Si $A$ es la matriz de filas $(0,3,1)$, $(2,-1,-1)$, $(-2,-1,-1)$ y $J$ la matriz de Jordan asociada, de filas $(2,0,0)$, $(0,-2,1)$, $(0,0,-2)$, justificar si la matriz de paso $P$ verifica: A) Tiene por filas $(1,1,0)$, $(1,-1,0)$, $(-1,1,1)$. B) Tiene por filas $(1,-1,1)$, $(1,1,-1)$, $(-1,-1,0)$. C) No es única. D) En MAXIMA hay instrucciones específicas para calcularla.
Qué se pide. Comprobar candidatas sin calcular $P$ de cero: $P$ vale si y solo si es regular y $AP=PJ$.
- Qué me piden
- Comprobar si dos matrices candidatas son matriz de paso $P$.
- Cómo lo pienso
- La matriz de paso cumple $AP=PJ$ columna a columna: cada columna es un vector propio o el siguiente de una cadena. Compruebo columna a columna.
- Qué uso
- Matriz de paso: por cada bloque, una cadena $v_1,\dots,v_s$ con $(A-\lambda I)v_1=0$ y $(A-\lambda I)v_{k}=v_{k-1}$; las cadenas por columnas forman $P$ y $AP=PJ$.
- Por dónde empiezo
- Calculo $AP$ y $PJ$ y comparo.
- Criterio. $J=P^{-1}AP\iff AP=PJ$ (multiplicando por $P$) con $P$ regular. Comprobar $AP=PJ$ evita invertir; hace falta $\det P\neq0$ porque con $P=0$ también se cumple.
A) Correcta.
$$P_A=\begin{pmatrix}1&1&0\cr1&-1&0\cr-1&1&1\end{pmatrix},$$$\det P_A=-2\neq0$:
$$AP_A=\begin{pmatrix}2&-2&1\cr2&2&-1\cr-2&-2&-1\end{pmatrix}=P_AJ.$$ Lectura por columnas: $u=(1,1,-1)$: $Au=2u$ (propio de $\lambda=2$); $w=(1,-1,1)$: $(A+2I)w=0$ (propio de $\lambda=-2$, primera columna del bloque); $v=(0,0,1)$: $(A+2I)v=(1,-1,1)=w$, es decir $Av=w-2v$ (tercera columna de $J$: $(0,1,-2)$).B) Correcta.
$$P_B=\begin{pmatrix}1&-1&1\cr1&1&-1\cr-1&-1&0\end{pmatrix},$$$\det P_B=-2\neq0$ y
$$AP_B=\begin{pmatrix}2&2&-3\cr2&-2&3\cr-2&2&-1\end{pmatrix}=P_BJ.$$C) Correcta — y por qué hay tantas. Con $P=(p_1\,p_2\,p_3)$, $AP=PJ$ dice: $(A-2I)p_1=0$ ($p_1=a\,u$, $a\neq0$); $(A+2I)p_2=0$ y $(A+2I)p_3=p_2$. Como
$$(A+2I)^2=\begin{pmatrix}8&8&0\cr8&8&0\cr-8&-8&0\end{pmatrix}:$$$\operatorname{Nuc}(A+2I)^2=\{x+y=0\}=\langle w,v\rangle$ y $\operatorname{Nuc}(A+2I)=\langle w\rangle$. Luego $p_3=b\,v+c\,w$ ($b\neq0$), $p_2=(A+2I)p_3=b\,w$ y
$$P=\big(a\,u\mid b\,w\mid b\,v+c\,w\big),\qquad a,b\neq0,\ c\in\mathbb R:$$ tres parámetros libres. $P_A$ es $a=b=1,c=0$; $P_B$ es $a=1,b=-1,c=1$. Ojo: $b$ es el mismo en $p_2$ y $p_3$.- D) Maxima (
ModeMatrix(A,l)), fuera de examen; E la da por correcta.
Respuesta: A, B, C y D correctas (como E).
Error típico: escalar solo la última columna de una cadena: dentro de un bloque $p_2=(A-\lambda I)p_3$ obliga a escalarlas juntas. Teoría: U Teor. 4.8 p.168; guion p.173; Ej. 4.18 pp.173-175. (Verificado con sympy; jordan_form da la misma $J$ con los bloques en otro orden, válido por U p.169.)
Bloque III — $J$ y $P$ de una $4\times4$ con cadenas (E 4.39–4.40)
Ejercicio 4.39 — forma de Jordan de una $4\times4$ con $\alpha=4$, $d=2$ E p.156-157
Si $(2,0,3,0)$, $(1,2,0,3)$, $(0,0,2,0)$, $(0,0,1,2)$ son las filas de $A$, se pide obtener la matriz de Jordan $J$ asociada.
Qué se pide. $A=\begin{pmatrix}2&0&3&0\cr1&2&0&3\cr0&0&2&0\cr0&0&1&2\end{pmatrix}$. Con $\alpha=4$ y $d=2$ hay dos candidatas, $(3,1)$ y $(2,2)$: hay que mirar los rangos de las potencias.
- Qué me piden
- Matriz de Jordan de una $4\times4$ con $\alpha=4$.
- Cómo lo pienso
- Un único valor propio. El número de bloques es $4-\operatorname{rang}(A-2I)$ y para saber sus tamaños miro también el rango de $(A-2I)^2$.
- Qué uso
- Forma de Jordan (Def. 4.13-4.14, U p.167): si todos los valores propios son reales, $A$ es semejante a $J$, diagonal por bloques $J_s(\lambda)$ (valor propio en la diagonal y unos justo encima). Para cada $\lambda$: el nº de bloques es $d(\lambda)=n-\operatorname{rang}(A-\lambda I)$ y la suma de sus tamaños es $\alpha(\lambda)$. Tamaños de los bloques: el nº de bloques de tamaño $\ge k$ es $\operatorname{rang}(A-\lambda I)^{k-1}-\operatorname{rang}(A-\lambda I)^{k}$.
- Por dónde empiezo
- Rangos de $A-2I$ y $(A-2I)^2$.
Polinomio (por la 3.ª fila, que solo tiene $2-\lambda$ en la columna 3; luego por la 1.ª fila):
$$|A-\lambda I|=(2-\lambda)\begin{vmatrix}2-\lambda&0&0\cr1&2-\lambda&3\cr0&0&2-\lambda\end{vmatrix}=(2-\lambda)^2\begin{vmatrix}2-\lambda&3\cr0&2-\lambda\end{vmatrix}=(2-\lambda)^4.$$Único valor propio $\lambda=2$, $\alpha=4$.
$d$:
$$N=A-2I=\begin{pmatrix}0&0&3&0\cr1&0&0&3\cr0&0&0&0\cr0&0&1&0\end{pmatrix}.$$La fila 1 es el triple de la fila 4 y la fila 2 es independiente: $\operatorname{rang}N=2$, $d=4-2=2\lt4$: no es diagonalizable, hay $2$ bloques. $\operatorname{Nuc}N$: $x_3=0,\ x_1+3x_4=0$.
- Candidatas: tamaños no crecientes que suman $4$ con $d=2$ bloques: $(3,1)$ o $(2,2)$ (U Ej. 4.16).
Cadena de núcleos:
$$N^2=\begin{pmatrix}0&0&0&0\cr0&0&6&0\cr0&0&0&0\cr0&0&0&0\end{pmatrix},$$$\operatorname{rang}N^2=1$; $N^3=0$. Dimensiones: $\dim\operatorname{Nuc}N=2$, $\dim\operatorname{Nuc}N^2=3$, $\dim\operatorname{Nuc}N^3=4=\alpha$ (se para).
- Decisión. Si fuera $(2,2)$, las dimensiones serían $2,4$; pero $\dim\operatorname{Nuc}N^2=3$ $\Rightarrow$ descartado. Con $(3,1)$: $1+1=2,\ 2+1=3,\ 3+1=4$ ✓. Equivalente: nº de bloques de tamaño $\ge k$ $=\dim\operatorname{Nuc}N^k-\dim\operatorname{Nuc}N^{k-1}$: $\ge1$: $2$; $\ge2$: $1$; $\ge3$: $1$ $\Rightarrow$ tamaños $3$ y $1$.
Resultado:
$$J=\begin{pmatrix}2&1&0&0\cr0&2&1&0\cr0&0&2&0\cr0&0&0&2\end{pmatrix}$$(sympy: la misma; rangos de $N^k$: $2,1,0$).
Errata (E p.157): en el caso hipotético $J_2$ escribe «$f(\bar e_4)=\bar e_3-2\bar e_4$»; debe ser $f(\bar e_4)=\bar e_3+2\bar e_4$.
Error típico: decidir los bloques solo con $d$ cuando $\alpha=4$ y $d=2$ (ejemplo propio 10). Teoría: U Teor. 4.1 p.146, Teor. 4.7 p.164, pp.169-170 (Ej. 4.16), guion p.173, p.176.
Ejercicio 4.40 — matriz de paso para el ejercicio 4.39 (errata en la comprobación) E p.158-159
Si $A$ y $J$ son las matrices del Ejercicio 4.39, se pide obtener una matriz de paso.
Qué se pide. Aplicar el guion de U p.173: $P=(e_1\,e_2\,e_3\,e_4)$ con $N=A-2I$ y $J=\operatorname{diag}(J_3(2),J_1(2))$.
- Qué me piden
- Matriz de paso para el $J$ del 4.39.
- Cómo lo pienso
- Para cada bloque necesito una cadena: un vector propio y los vectores que lo van generando con $(A-\lambda I)v=\text{anterior}$.
- Qué uso
- Matriz de paso: por cada bloque, una cadena $v_1,\dots,v_s$ con $(A-\lambda I)v_1=0$ y $(A-\lambda I)v_{k}=v_{k-1}$; las cadenas por columnas forman $P$ y $AP=PJ$.
- Por dónde empiezo
- Resuelvo por cadenas y compruebo $AP=PJ$.
- Qué base buscamos. $AP=PJ$ columna a columna: $Ae_1=2e_1$, $Ae_2=e_1+2e_2$, $Ae_3=e_2+2e_3$, $Ae_4=2e_4$; es decir $Ne_1=0$, $Ne_2=e_1$, $Ne_3=e_2$, $Ne_4=0$. Así $(e_1,e_2,e_3)=(N^2v,\,Nv,\,v)$ con $v=e_3$: una cadena que empieza en $v$ y termina en un vector propio. Se calcula de derecha a izquierda.
- Cadena de núcleos (de 4.39): $\operatorname{Nuc}N:\{x_3=0,x_1+3x_4=0\}$ (dim $2$) $\subset\operatorname{Nuc}N^2:\{x_3=0\}$ (dim $3$) $\subset\operatorname{Nuc}N^3=\mathbb R^4$ (dim $4=\alpha$, $k=3$).
- Elegir $v\in\operatorname{Nuc}N^3$ que no esté en $\operatorname{Nuc}N^2$, es decir con $x_3\neq0$: lo más sencillo, $v=(0,0,1,0)$. ¿Por qué fuera de $\operatorname{Nuc}N^2$? Para que $N^2v\neq0$: así la cadena tiene longitud $3$.
- Bajar: $e_2=Nv=$ tercera columna de $N=(3,0,0,1)$ (está en $\operatorname{Nuc}N^2$, no en $\operatorname{Nuc}N$: $3+3\neq0$). $e_1=N(3,0,0,1)=(0,6,0,0)$ (en $\operatorname{Nuc}N$ ✓: es propio). Hay $3$ vectores; faltan $\alpha-3=1$.
- Vector del bloque $1\times1$: propio ($Ne_4=0$) e independiente de $e_1$: $\operatorname{Nuc}N=\langle(0,1,0,0),(-3,0,0,1)\rangle$; vale $a(0,1,0,0)+b(-3,0,0,1)$ con $b\neq0$. E toma $e_4=(1,0,0,-\tfrac13)$; también vale $(-3,0,0,1)$, que evita fracciones.
- Matriz de paso (columnas $e_1,e_2,e_3,e_4$): $$P=\begin{pmatrix}0&3&0&1\cr6&0&0&0\cr0&0&1&0\cr0&1&0&-\tfrac13\end{pmatrix},\qquad\det P=12\neq0,\qquad AP=PJ.$$ Con $e_4=(-3,0,0,1)$: $\det P=-36$ y también $AP=PJ$.
- Si se ordenan mal las columnas $(v,Nv,N^2v)$: $P^{-1}AP$ tiene en ese bloque los unos debajo de la diagonal (la traspuesta), no la $J$ de U.
Errata (E p.159): tras comprobar $P^{-1}AP=J$ escribe «siendo $P^{-1}=\begin{pmatrix}2&0&3&0\cr1&2&0&3\cr0&0&2&0\cr0&0&1&2\end{pmatrix}^{-1}$», es decir, iguala $P^{-1}$ con la inversa de $A$; debe ser la inversa de $P$ (la matriz numérica $P^{-1}$ que da el libro sí es correcta).
Error típico: (a) tomar $v\in\operatorname{Nuc}N^2$ (p. ej. $(1,0,0,0)$): la cadena se corta antes; (b) $e_4$ proporcional a $e_1$: $P$ singular; (c) invertir el orden de la cadena. Teoría: U Def. 4.13 y guion p.173; Ej. 4.18-4.19 pp.173-177. (Verificado con sympy.)
Bloque IV — Método de potencias (E 4.41)
Ejercicio 4.41 — aproximar el vector propio estrictamente dominante E p.159-160
Se pide aproximar un vector propio estrictamente dominante de la matriz $\begin{pmatrix}-4&1\cr-6&1\end{pmatrix}$ mediante el método de potencias.
Qué se pide. Aplicar el algoritmo de U p.185 con $x_0=(1,1)$ y leer el resultado.
- Qué me piden
- Aproximar un vector propio dominante con el método de las potencias.
- Cómo lo pienso
- Multiplicar una y otra vez por $A$ hace que domine la dirección del valor propio de mayor módulo.
- Qué uso
- Método de potencias: la iteración $x_{k+1}=Ax_k$ normalizada converge hacia el vector propio del valor propio de mayor módulo.
- Por dónde empiezo
- Itero con un vector inicial y normalizo.
- Algoritmo: $\mu_k$ = coordenada de $Ax_k$ de mayor valor absoluto, con su signo; $x_{k+1}=\tfrac1{\mu_k}Ax_k$.
- Iteraciones a mano (fracciones exactas): $k=0$: $Ax_0=(-4+1,\,-6+1)=(-3,-5)$, $\mu_0=-5$, $x_1=(\tfrac35,1)$. $k=1$: $Ax_1=(-\tfrac75,-\tfrac{13}5)$, $\mu_1=-\tfrac{13}5$, $x_2=(\tfrac7{13},1)$. $k=2$: $Ax_2=(-\tfrac{15}{13},-\tfrac{29}{13})$, $\mu_2=-\tfrac{29}{13}$, $x_3=(\tfrac{15}{29},1)$. $k=3$: $\mu_3=-\tfrac{61}{29}$, $x_4=(\tfrac{31}{61},1)$. $k=4$: $\mu_4=-\tfrac{125}{61}$, $x_5=(\tfrac{63}{125},1)$. Atajo: con $x_k=(a,1)$ y $a\gt0$, $Ax_k=(1-4a,\,1-6a)$ y la 2.ª coordenada siempre gana en valor absoluto.
- Tabla:
(Coincide con la tabla de Maxima de E p.160, $k=7$: $\mu\approx-2.00589$, $x\approx(0.50049,1)$.)$k$ $\mu_k$ $\mu_k$ (dec.) $x_{k+1}$ 1.ª coord. $0$ $-5$ $-5$ $(\tfrac35,1)$ $0.6$ $1$ $-\tfrac{13}5$ $-2.6$ $(\tfrac7{13},1)$ $0.5385$ $2$ $-\tfrac{29}{13}$ $-2.2308$ $(\tfrac{15}{29},1)$ $0.5172$ $3$ $-\tfrac{61}{29}$ $-2.1034$ $(\tfrac{31}{61},1)$ $0.5082$ $4$ $-\tfrac{125}{61}$ $-2.0492$ $(\tfrac{63}{125},1)$ $0.504$ $5$ $-\tfrac{253}{125}$ $-2.024$ $(\tfrac{127}{253},1)$ $0.5020$ $6$ $-\tfrac{509}{253}$ $-2.0119$ $(\tfrac{255}{509},1)$ $0.5010$ $7$ $-\tfrac{1021}{509}$ $-2.0059$ $(\tfrac{511}{1021},1)$ $0.5005$ - Lectura: $\mu_k\to-2$ y $x_k\to(\tfrac12,1)\propto(1,2)$. Estimación: valor propio estrictamente dominante $\lambda_1\approx-2$, vector propio $\approx(1,2)$.
Comprobación exacta:
$$|A-\lambda I|=(-4-\lambda)(1-\lambda)+6=\lambda^2+3\lambda+2=(\lambda+1)(\lambda+2):$$valores propios $-2$ y $-1$; $|-2|\gt|-1|$, luego $-2$ es estrictamente dominante (U Def. 4.14).
$$A+2I=\begin{pmatrix}-2&1\cr-6&3\end{pmatrix}\Rightarrow(1,2);$$ $$A+I=\begin{pmatrix}-3&1\cr-6&2\end{pmatrix}\Rightarrow(1,3).$$- Por qué converge y a qué ritmo. $x_0=(1,1)=2(1,2)-(1,3)$, luego $A^kx_0=2(-2)^k(1,2)-(-1)^k(1,3)$ y $x_k=\Big(\dfrac{2^{k+1}-1}{2^{k+2}-3},1\Big)\to(\tfrac12,1)$, $\mu_k=-\dfrac{2^{k+3}-3}{2^{k+2}-3}\to-2$. El término «molesto» pesa como $(|\lambda_2|/|\lambda_1|)^k=(1/2)^k$: el error se divide por $2$ en cada paso (se ve en la tabla: $0.1,\ 0.038,\ 0.017,\ 0.008$…).
- Signo: aunque $\lambda_1\lt0$, $x_k$ no cambia de signo porque se divide por $\mu_k$ con su signo; por eso $\mu_k$ estima $-2$ y no $2$.
Errata (E p.160): en la última línea a mano, «$\mu_4=-\frac{125}{61},\ \bar x_4=(\frac{63}{125},1)$» debe decir $\bar x_5$ (la tabla de Maxima de esa página está bien).
Error típico: tomar $\mu_k$ en valor absoluto ($+5$ en vez de $-5$): entonces $x_k$ alterna de signo y $\mu_k\to2$, valor equivocado. Teoría: U §4.4.3, Def. 4.14 p.181, Ej. 4.24 p.182, algoritmo p.185. (Verificado con sympy en fracciones.)
Ejemplos propios (dificultad creciente)
Ejemplo propio 1 (fácil) — $J$ de una $2\times2$ no diagonalizable propio
Halla $J$ y una matriz de paso de $A=\begin{pmatrix}1&1\cr-1&3\end{pmatrix}$.
- Qué me piden
- Jordan y matriz de paso de una $2\times2$.
- Cómo lo pienso
- Valor doble con un único vector propio: un bloque $2\times2$. La cadena es $v_1\in L_\lambda$, $v_2$ con $(A-\lambda I)v_2=v_1$.
- Qué uso
- Forma de Jordan (Def. 4.13-4.14, U p.167): si todos los valores propios son reales, $A$ es semejante a $J$, diagonal por bloques $J_s(\lambda)$ (valor propio en la diagonal y unos justo encima). Para cada $\lambda$: el nº de bloques es $d(\lambda)=n-\operatorname{rang}(A-\lambda I)$ y la suma de sus tamaños es $\alpha(\lambda)$. Matriz de paso: por cada bloque, una cadena $v_1,\dots,v_s$ con $(A-\lambda I)v_1=0$ y $(A-\lambda I)v_{k}=v_{k-1}$; las cadenas por columnas forman $P$ y $AP=PJ$.
- Por dónde empiezo
- Polinomio, $v_1$ y $v_2$.
- $$|A-\lambda I|=(1-\lambda)(3-\lambda)+1=\lambda^2-4\lambda+4=(\lambda-2)^2:$$
$\lambda=2$, $\alpha=2$.
- $$A-2I=\begin{pmatrix}-1&1\cr-1&1\end{pmatrix},$$
rango $1$, $d=1\lt2$: no diagonalizable; un bloque:
$$J=\begin{pmatrix}2&1\cr0&2\end{pmatrix}.$$ - $(A-2I)^2=0$: sirve cualquier $v$ con $(A-2I)v\neq0$: $v=(1,0)$, $(A-2I)v=(-1,-1)$ (vector propio).
- $$P=\begin{pmatrix}-1&1\cr-1&0\end{pmatrix}$$
(vector propio primero), $\det P=1$. Comprobación: $A(-1,-1)=(-2,-2)$ ✓; $A(1,0)=(1,-1)=(-1,-1)+2(1,0)$ ✓. Es decir $AP=PJ$.
Ejemplo propio 2 (fácil-medio) — $3\times3$, $\lambda$ triple con $\operatorname{rang}(A-\lambda I)=1$ propio
Halla $J$ y $P$ de $A=\begin{pmatrix}3&-1&1\cr1&1&1\cr0&0&2\end{pmatrix}$.
- Qué me piden
- $J$ y $P$ de una $3\times3$ con valor triple y $\operatorname{rang}(A-\lambda I)=1$.
- Cómo lo pienso
- Rango 1 en dimensión 3: $d=2$, dos bloques que suman 3. Solo cabe un bloque de 2 y otro de 1.
- Qué uso
- Forma de Jordan (Def. 4.13-4.14, U p.167): si todos los valores propios son reales, $A$ es semejante a $J$, diagonal por bloques $J_s(\lambda)$ (valor propio en la diagonal y unos justo encima). Para cada $\lambda$: el nº de bloques es $d(\lambda)=n-\operatorname{rang}(A-\lambda I)$ y la suma de sus tamaños es $\alpha(\lambda)$. Tamaños de los bloques: el nº de bloques de tamaño $\ge k$ es $\operatorname{rang}(A-\lambda I)^{k-1}-\operatorname{rang}(A-\lambda I)^{k}$.
- Por dónde empiezo
- Cadena de longitud 2 y un vector propio independiente.
Por la 3.ª fila:
$$|A-\lambda I|=(2-\lambda)\big[(3-\lambda)(1-\lambda)+1\big]=(2-\lambda)(\lambda-2)^2=-(\lambda-2)^3.$$- $$N=A-2I=\begin{pmatrix}1&-1&1\cr1&-1&1\cr0&0&0\end{pmatrix},$$
rango $1$, $d=2$: bloques $(2,1)$,
$$J'=\begin{pmatrix}2&1&0\cr0&2&0\cr0&0&2\end{pmatrix}.$$ - $N^2=0$. $\operatorname{Nuc}N:\ x-y+z=0$, es $\langle(1,1,0),(0,1,1)\rangle$.
- Cabeza $v=(1,0,0)\notin\operatorname{Nuc}N$; $Nv=(1,1,0)$. Tercer vector: propio e independiente de $(1,1,0)$: $(0,1,1)$.
- $$P=\begin{pmatrix}1&1&0\cr1&0&1\cr0&0&1\end{pmatrix}$$
(columnas $Nv,v,(0,1,1)$), $\det P=-1$, $AP=PJ'$. Trampa: como 3.ª columna $(1,1,0)$ ($P$ singular) o $(0,0,1)$ (no es propio).
Ejemplo propio 3 (medio) — $\lambda$ triple con $\operatorname{rang}(A-\lambda I)=2$: un bloque $3\times3$ propio
Halla $J$ y $P$ de $A=\begin{pmatrix}0&1&0\cr0&0&1\cr1&-3&3\end{pmatrix}$.
- Qué me piden
- $J$ y $P$ con valor triple y rango 2.
- Cómo lo pienso
- Rango 2 ⇒ $d=1$ ⇒ un único bloque $3\times3$. Hace falta una cadena de tres vectores.
- Qué uso
- Forma de Jordan (Def. 4.13-4.14, U p.167): si todos los valores propios son reales, $A$ es semejante a $J$, diagonal por bloques $J_s(\lambda)$ (valor propio en la diagonal y unos justo encima). Para cada $\lambda$: el nº de bloques es $d(\lambda)=n-\operatorname{rang}(A-\lambda I)$ y la suma de sus tamaños es $\alpha(\lambda)$. Matriz de paso: por cada bloque, una cadena $v_1,\dots,v_s$ con $(A-\lambda I)v_1=0$ y $(A-\lambda I)v_{k}=v_{k-1}$; las cadenas por columnas forman $P$ y $AP=PJ$.
- Por dónde empiezo
- Resuelvo $(A-\lambda I)v_3=v_2$, $(A-\lambda I)v_2=v_1$.
- $$|A-\lambda I|=-\lambda\big[-\lambda(3-\lambda)+3\big]-1\cdot(0-1)=-\lambda^3+3\lambda^2-3\lambda+1=-(\lambda-1)^3.$$
- $$N=A-I=\begin{pmatrix}-1&1&0\cr0&-1&1\cr1&-3&2\end{pmatrix}:$$
fila 3 $=-$fila 1 $+2\cdot$fila 2, rango $2$, $d=1$: un bloque,
$$J''=\begin{pmatrix}1&1&0\cr0&1&1\cr0&0&1\end{pmatrix}.$$ - $$N^2=\begin{pmatrix}1&-2&1\cr1&-2&1\cr1&-2&1\end{pmatrix}$$
(rango $1$), $N^3=0$: dimensiones $1,2,3$.
- $v\notin\operatorname{Nuc}N^2=\{x-2y+z=0\}$: $v=(1,0,0)$; $Nv=(-1,0,1)$; $N^2v=(1,1,1)$ (propio: $N(1,1,1)=0$ ✓).
- $$P=(N^2v\mid Nv\mid v)=\begin{pmatrix}1&-1&1\cr1&0&0\cr1&1&0\end{pmatrix},$$
$\det P=1$, $AP=PJ''$.
Ejemplo propio 4 (medio) — $4\times4$ con dos valores propios dobles propio
Halla $J$ y $P$ de $A=\begin{pmatrix}1&1&1&0\cr0&1&0&1\cr0&0&2&0\cr0&0&0&2\end{pmatrix}$.
- Qué me piden
- $J$ y $P$ de una $4\times4$ con dos valores propios dobles.
- Cómo lo pienso
- Cada valor propio se trata por separado.
- Qué uso
- Forma de Jordan (Def. 4.13-4.14, U p.167): si todos los valores propios son reales, $A$ es semejante a $J$, diagonal por bloques $J_s(\lambda)$ (valor propio en la diagonal y unos justo encima). Para cada $\lambda$: el nº de bloques es $d(\lambda)=n-\operatorname{rang}(A-\lambda I)$ y la suma de sus tamaños es $\alpha(\lambda)$. Matriz de paso: por cada bloque, una cadena $v_1,\dots,v_s$ con $(A-\lambda I)v_1=0$ y $(A-\lambda I)v_{k}=v_{k-1}$; las cadenas por columnas forman $P$ y $AP=PJ$.
- Por dónde empiezo
- Rangos de $A-\lambda I$ y cadenas.
- Triangular: valores propios $1,1,2,2$; $|A-\lambda I|=(\lambda-1)^2(\lambda-2)^2$.
$\lambda=1$:
$$A-I=\begin{pmatrix}0&1&1&0\cr0&0&0&1\cr0&0&1&0\cr0&0&0&1\end{pmatrix},$$rango $3$, $d=1\lt2$: un bloque $2\times2$.
$\lambda=2$:
$$A-2I=\begin{pmatrix}-1&1&1&0\cr0&-1&0&1\cr0&0&0&0\cr0&0&0&0\end{pmatrix},$$rango $2$, $d=2=\alpha$: dos bloques $1\times1$.
- $$J=\begin{pmatrix}1&1&0&0\cr0&1&0&0\cr0&0&2&0\cr0&0&0&2\end{pmatrix}$$
(nº de bloques $3=1+2$; unos: $4-3=1$).
Cadena de $\lambda=1$:
$$(A-I)^2=\begin{pmatrix}0&0&1&1\cr0&0&0&1\cr0&0&1&0\cr0&0&0&1\end{pmatrix},$$$\operatorname{Nuc}=\langle e_1,e_2\rangle$ (dim $2=\alpha$) y $\operatorname{Nuc}(A-I)=\langle e_1\rangle$: $v=e_2$, $(A-I)e_2=e_1$.
- Propios de $\lambda=2$: $-x+y+z=0,\ -y+w=0$: $(1,0,1,0)$, $(1,1,0,1)$.
- $$P=\begin{pmatrix}1&0&1&1\cr0&1&0&1\cr0&0&1&0\cr0&0&0&1\end{pmatrix},$$
$\det P=1$, $AP=PJ$. Cada valor propio se trabaja por separado con su $A-\lambda I$.
Ejemplo propio 5 (medio) — todas las $J$ posibles para un polinomio dado (como U Ej. 4.16 y 4.22) propio
$A$ es de orden $5$ con $|A-\lambda I|=-(\lambda-2)^3(\lambda+1)^2$. ¿Qué formas de Jordan son posibles? ($J_s(\lambda)$ = bloque de tamaño $s$.)
- Qué me piden
- Enumerar las formas de Jordan posibles dado el polinomio característico.
- Cómo lo pienso
- Para cada valor propio, los tamaños de sus bloques forman una partición de su multiplicidad.
- Qué uso
- Forma de Jordan (Def. 4.13-4.14, U p.167): si todos los valores propios son reales, $A$ es semejante a $J$, diagonal por bloques $J_s(\lambda)$ (valor propio en la diagonal y unos justo encima). Para cada $\lambda$: el nº de bloques es $d(\lambda)=n-\operatorname{rang}(A-\lambda I)$ y la suma de sus tamaños es $\alpha(\lambda)$.
- Por dónde empiezo
- Listo particiones de 3 y de 2 y las combino.
- $\lambda=2$ ($\alpha=3$): $d=1$: $J_3(2)$; $d=2$: $J_2(2)\oplus J_1(2)$; $d=3$: $2I_3$.
- $\lambda=-1$ ($\alpha=2$): $d=1$: $J_2(-1)$; $d=2$: $-I_2$.
Combinaciones: $3\cdot2=6$ formas de Jordan (salvo orden). Nº de bloques $d(2)+d(-1)$, entre $2$ y $5$. Por ejemplo, con $d(2)=2$ y $d(-1)=1$:
$$J=\operatorname{diag}\big(J_2(2),\,J_1(2),\,J_2(-1)\big)=\begin{pmatrix}2&1&0&0&0\cr0&2&0&0&0\cr0&0&2&0&0\cr0&0&0&-1&1\cr0&0&0&0&-1\end{pmatrix}.$$- Solo una es diagonal ($d=\alpha$ en ambos). Con $\alpha\le3$, $d$ determina los bloques; con $\alpha=4$ y $d=2$ ya no (ejemplo propio 10).
Ejemplo propio 6 (medio) — deducir $J$ a partir de los rangos de $(A-\lambda I)^k$ propio
$A$ de orden $5$ con único valor propio $\lambda=3$, y $\operatorname{rang}(A-3I)=3$, $\operatorname{rang}(A-3I)^2=1$, $\operatorname{rang}(A-3I)^3=0$. Halla $J$.
- Qué me piden
- Deducir $J$ a partir de los rangos de $(A-3I)^k$.
- Cómo lo pienso
- Las diferencias entre rangos consecutivos cuentan cuántos bloques hay de cada tamaño.
- Qué uso
- Forma de Jordan (Def. 4.13-4.14, U p.167): si todos los valores propios son reales, $A$ es semejante a $J$, diagonal por bloques $J_s(\lambda)$ (valor propio en la diagonal y unos justo encima). Para cada $\lambda$: el nº de bloques es $d(\lambda)=n-\operatorname{rang}(A-\lambda I)$ y la suma de sus tamaños es $\alpha(\lambda)$. Tamaños de los bloques: el nº de bloques de tamaño $\ge k$ es $\operatorname{rang}(A-\lambda I)^{k-1}-\operatorname{rang}(A-\lambda I)^{k}$.
- Por dónde empiezo
- Resto rangos consecutivos.
- Dimensiones de los núcleos: $5-3=2$, $5-1=4$, $5-0=5=\alpha$ (se para).
- Nº de bloques de tamaño $\ge k$ $=\dim\operatorname{Nuc}_k-\dim\operatorname{Nuc}_{k-1}$: $\ge1$: $2$; $\ge2$: $4-2=2$; $\ge3$: $5-4=1$.
- Dos bloques, ambos de tamaño $\ge2$ y uno de ellos $\ge3$, con suma $5$: tamaños $(3,2)$. $J=\operatorname{diag}(J_3(3),J_2(3))$.
Justificación: $(J_s-\lambda I)^k$ tiene rango $\max(s-k,0)$ y
$$\operatorname{rang}(A-\lambda I)^k=\operatorname{rang}(J-\lambda I)^k$$por semejanza (U Teor. 4.1, p.146). sympy: rangos $3,1,0$ para esa $J$ ✓.
- Variante rápida: con $\alpha=4$, $\operatorname{rang}(A-\lambda I)=2$ y $(A-\lambda I)^2=0$ $\Rightarrow$ dimensiones $2,4$ $\Rightarrow(2,2)$; si $\operatorname{rang}(A-\lambda I)^2=1$ $\Rightarrow(3,1)$.
Ejemplo propio 7 (medio) — matriz nilpotente propio
Halla $J$ y $P$ de $A=\begin{pmatrix}1&1&0\cr-1&-1&1\cr0&0&0\end{pmatrix}$.
- Qué me piden
- $J$ y $P$ de una matriz nilpotente.
- Cómo lo pienso
- Único valor propio $0$; el tamaño máximo de bloque es el menor $k$ con $A^k=0$.
- Qué uso
- Forma de Jordan (Def. 4.13-4.14, U p.167): si todos los valores propios son reales, $A$ es semejante a $J$, diagonal por bloques $J_s(\lambda)$ (valor propio en la diagonal y unos justo encima). Para cada $\lambda$: el nº de bloques es $d(\lambda)=n-\operatorname{rang}(A-\lambda I)$ y la suma de sus tamaños es $\alpha(\lambda)$. Matriz de paso: por cada bloque, una cadena $v_1,\dots,v_s$ con $(A-\lambda I)v_1=0$ y $(A-\lambda I)v_{k}=v_{k-1}$; las cadenas por columnas forman $P$ y $AP=PJ$.
- Por dónde empiezo
- Potencias de $A$ y cadena en el núcleo.
- $$A^2=\begin{pmatrix}0&0&1\cr0&0&-1\cr0&0&0\end{pmatrix}$$
y $A^3=0$: nilpotente (U p.169).
- Único valor propio $0$: si $Av=\lambda v$ con $v\neq0$, entonces $0=A^3v=\lambda^3v$, luego $\lambda=0$. En efecto $|A-\lambda I|=-\lambda^3$.
$\operatorname{rang}A=2$, $d=1$: un bloque $3\times3$ con $\lambda=0$:
$$J=\begin{pmatrix}0&1&0\cr0&0&1\cr0&0&0\end{pmatrix}.$$- $v\notin\operatorname{Nuc}A^2$ ($A^2v=(z,-z,0)$, basta $z\neq0$): $v=e_3$; $Av=(0,1,0)$; $A^2v=(1,-1,0)$.
- $$P=(A^2v\mid Av\mid v)=\begin{pmatrix}1&0&0\cr-1&1&0\cr0&0&1\end{pmatrix},$$
$\det P=1$ y $AP=PJ$.
- Consecuencia: una nilpotente no nula nunca es diagonalizable (si lo fuera, $\Lambda=0$ y $A=P\,0\,P^{-1}=0$).
Ejemplo propio 8 (medio-difícil) — potencias de un bloque y $A^n=PJ^nP^{-1}$ propio
Con $A=\begin{pmatrix}1&1\cr-1&3\end{pmatrix}$ (ejemplo propio 1), halla $A^n$.
- Qué me piden
- Calcular $A^n$ con un bloque de Jordan.
- Cómo lo pienso
- Un bloque $\lambda I+N$ con $N^2=0$ se eleva con el binomio, que solo tiene dos términos.
- Qué uso
- Potencia de un bloque: $(\lambda I+N)^n=\lambda^nI+n\lambda^{n-1}N$ si $N^2=0$.
- Por dónde empiezo
- $J^n$ y $PJ^nP^{-1}$.
$J=2I+N$ con
$$N=\begin{pmatrix}0&1\cr0&0\end{pmatrix},$$$N^2=0$. Como $2I$ conmuta con $N$:
$$J^n=2^nI+n2^{n-1}N=\begin{pmatrix}2^n&n2^{n-1}\cr0&2^n\end{pmatrix}.$$(Calcular $J^n$ elevando cada entrada a $n$ solo vale para diagonales.)
De $A=PJP^{-1}$ con
$$P=\begin{pmatrix}-1&1\cr-1&0\end{pmatrix},$$ $$P^{-1}=\begin{pmatrix}0&-1\cr1&-1\end{pmatrix}:$$$A^n=PJ^nP^{-1}$ (los $P^{-1}P$ intermedios se cancelan).
$$PJ^n=\begin{pmatrix}-2^n&2^n-n2^{n-1}\cr-2^n&-n2^{n-1}\end{pmatrix}$$$$A^n=PJ^nP^{-1}=2^{n-1}\begin{pmatrix}2-n&n\cr-n&n+2\end{pmatrix}.$$Comprobaciones: $n=1$ da $A$; $n=2$ da
$$\begin{pmatrix}0&4\cr-4&8\end{pmatrix}=A^2$$(sympy: correcto para $n=1,2,3,5$).
Ejemplo propio 9 (medio-difícil) — raíces complejas: no hay $J$ real propio
Estudia si $A=\begin{pmatrix}1&-2\cr1&-1\end{pmatrix}$ tiene forma de Jordan real.
- Qué me piden
- Decidir si una matriz tiene forma de Jordan real.
- Cómo lo pienso
- Si el polinomio tiene raíces complejas, no hay $J$ real.
- Qué uso
- Forma de Jordan (Def. 4.13-4.14, U p.167): si todos los valores propios son reales, $A$ es semejante a $J$, diagonal por bloques $J_s(\lambda)$ (valor propio en la diagonal y unos justo encima). Para cada $\lambda$: el nº de bloques es $d(\lambda)=n-\operatorname{rang}(A-\lambda I)$ y la suma de sus tamaños es $\alpha(\lambda)$.
- Por dónde empiezo
- Raíces del polinomio.
- $|A-\lambda I|=(1-\lambda)(-1-\lambda)+2=\lambda^2+1$: sin raíces reales.
- Si existiera $J$ real triangular con $J=P^{-1}AP$, tendría el mismo polinomio característico (U Teor. 4.4), que sería $(\lambda-j_{11})(\lambda-j_{22})$ con $j_{ii}$ reales: contradicción. No hay forma de Jordan real; el Teor. 4.8 requiere tácitamente raíces reales.
- Sobre $\mathbb C$: $\lambda=\pm i$ distintos, diagonalizable, $J=\operatorname{diag}(i,-i)$, con vectores $(1+i,1)$ y $(1-i,1)$ (sympy).
Ejemplo propio 10 (difícil) — mismo polinomio y misma $d$, pero no semejantes propio
$A=\begin{pmatrix}2&0&3&0\cr1&2&0&3\cr0&0&2&0\cr0&0&1&2\end{pmatrix}$ (la de E 4.39) y $B=\begin{pmatrix}2&0&1&0\cr0&2&0&1\cr0&0&2&0\cr0&0&0&2\end{pmatrix}$. ¿Son semejantes?
- Qué me piden
- Decidir si dos matrices con igual polinomio y $d$ son semejantes.
- Cómo lo pienso
- Si los rangos de $(A-\lambda I)^2$ difieren, los tamaños de bloques difieren.
- Qué uso
- Forma de Jordan (Def. 4.13-4.14, U p.167): si todos los valores propios son reales, $A$ es semejante a $J$, diagonal por bloques $J_s(\lambda)$ (valor propio en la diagonal y unos justo encima). Para cada $\lambda$: el nº de bloques es $d(\lambda)=n-\operatorname{rang}(A-\lambda I)$ y la suma de sus tamaños es $\alpha(\lambda)$. Tamaños de los bloques: el nº de bloques de tamaño $\ge k$ es $\operatorname{rang}(A-\lambda I)^{k-1}-\operatorname{rang}(A-\lambda I)^{k}$.
- Por dónde empiezo
- Comparo rangos de potencias.
- Ambas tienen el único valor propio $2$ y $|A-\lambda I|=(\lambda-2)^4$ ($B$ es triangular; $A$ en 4.39).
$\operatorname{rang}(A-2I)=2$ y
$$B-2I=\begin{pmatrix}0&0&1&0\cr0&0&0&1\cr0&0&0&0\cr0&0&0&0\end{pmatrix}$$también tiene rango $2$: $d=2$ en las dos.
- Distinguirlas: $(A-2I)^2$ tiene rango $1$ (4.39), pero $(B-2I)^2=0$ (rango $0$).
- Si fueran semejantes, $B=Q^{-1}AQ$ y $(B-2I)^2=Q^{-1}(A-2I)^2Q$ tendría el mismo rango (U Teor. 4.1, p.146): $1\neq0$. No son semejantes.
- Sus formas de Jordan: $A\sim\operatorname{diag}(J_3(2),J_1(2))$ y $B\sim\operatorname{diag}(J_2(2),J_2(2))$ (sympy). Moraleja: polinomio característico, traza, determinante y $d$ son invariantes, pero no bastan cuando $\alpha\ge4$.
Ejemplo propio 11 (medio) — método de potencias con $A^{-1}$: el valor propio de menor módulo propio
$A=\begin{pmatrix}2&1\cr1&2\end{pmatrix}$ ($|A-\lambda I|=(\lambda-1)(\lambda-3)$; vectores $(1,1)$ para $3$ y $(1,-1)$ para $1$). Estima el valor propio menor con $x_0=(1,0)$.
- Qué me piden
- Estimar el valor propio de menor módulo con el método de potencias en $A^{-1}$.
- Cómo lo pienso
- Los valores propios de $A^{-1}$ son los inversos; el dominante de $A^{-1}$ da el menor de $A$.
- Qué uso
- Los valores propios de $A^{-1}$ son los inversos; el dominante de $A^{-1}$ da el de menor módulo de $A$.
- Por dónde empiezo
- Itero con $A^{-1}$.
- Idea (U Ej. 4.30): si $Av=\lambda v$ con $\lambda\neq0$, $A^{-1}v=\lambda^{-1}v$. Los valores de $A^{-1}$ son $\tfrac13$ y $1$; el dominante es $1$ (del $\lambda=1$ de $A$).
- $$A^{-1}=\tfrac13\begin{pmatrix}2&-1\cr-1&2\end{pmatrix}$$
($|A|=3$).
- Iteraciones:
$k$ $A^{-1}x_k$ $\mu_k$ $x_{k+1}$ $0$ $(\tfrac23,-\tfrac13)$ $0.6667$ $(1,-\tfrac12)$ $1$ $(\tfrac56,-\tfrac23)$ $0.8333$ $(1,-\tfrac45)$ $2$ $(\tfrac{14}{15},-\tfrac{13}{15})$ $0.9333$ $(1,-\tfrac{13}{14})$ $3$ $(\tfrac{41}{42},-\tfrac{20}{21})$ $0.9762$ $(1,-\tfrac{40}{41})$ $4$ $(\tfrac{122}{123},-\tfrac{121}{123})$ $0.9919$ $(1,-\tfrac{121}{122})$ Ejemplo de cálculo ($k=1$):
$$A^{-1}(1,-\tfrac12)=\tfrac13(2+\tfrac12,\,-1-1)=(\tfrac56,-\tfrac23).$$ - $\mu_k\to1\Rightarrow$ valor propio de $A$: $1/1=1$; $x_k\to(1,-1)$ ✓. Ritmo $\tfrac{1/3}{1}=\tfrac13$ por paso.
Truco sin fracciones: iterar con
$$3A^{-1}=\begin{pmatrix}2&-1\cr-1&2\end{pmatrix}$$(los mismos vectores; $\mu_k\to3$ y el valor propio de $A$ es $3/3=1$).
Autoevaluación de U (capítulo 4) — ítems de este tema
- Ítem 9: «Si $A$ es una matriz no diagonalizable de orden $4$ tal que posee dos valores propios $\lambda_1,\lambda_2$ distintos, puede suceder que $d_1=2$ y $\alpha_2=3$.» Falso. Si $\alpha_2=3$, como $\alpha_1+\alpha_2\le4$ y $\alpha_1\ge1$, necesariamente $\alpha_1=1$, y entonces $d_1\le\alpha_1=1$ (Teor. 4.5). Ya razonado en el tema 4.3.
- Ítem 10: «Si $A$ es una matriz de orden $8$ y solo posee un valor propio $\lambda\neq1$ con dimensión geométrica $7$, entonces la matriz de Jordan de $A$ tiene solo un uno.» Verdadero. Único valor propio y orden $8$: $\alpha=8$ (todas las raíces reales, como supone el enunciado). Con $d=7$ hay $7$ bloques cuyos tamaños suman $8$ y no crecen: $(2,1,1,1,1,1,1)$. Un bloque $2\times2$ tiene un único $1$ y los otros seis son $1\times1$. Nº de unos $=n-$ nº de bloques $=8-7=1$ (E 4.34). Equivalentemente $\operatorname{rang}(A-\lambda I)=8-7=1$, y U añade que basta buscar un solo vector del subespacio propio generalizado para hallar la matriz de paso. El valor concreto de $\lambda$ es irrelevante.
9Ejemplos de examen X
De los 24 exámenes escaneados, cuatro ejercicios (uno de ellos un apartado suelto) tienen la forma de Jordan como contenido principal; ninguno pregunta por el método de potencias. Se agrupan por método, de menos a más trabajo: A si la matriz es diagonalizable, su Jordan es la diagonal; B no diagonalizable: se leen los bloques en el rango; C Jordan con matriz de paso por cadenas. Las soluciones oficiales son breves: aquí se desarrolla el razonamiento completo y se avisa de sus errores. Todo verificado con sympy.
Método A — Diagonalizable $\Rightarrow$ su matriz de Jordan es la diagonal
Ene. 2022 (tipo A), ejercicio 4 — «una posible matriz de Jordan» de una matriz diagonalizable X Ene. 2022 (A) ej. 4
Si $A$ es la matriz de filas $(2,2,-5)$, $(3,7,-15)$, $(1,2,-4)$, obtener una posible matriz de Jordan suya.
Qué se pide. Es la matriz de E 4.37. Aquí lo que enseña es que si $A$ es diagonalizable, su $J$ es la diagonal (U p.169; Ej. 4.17, p.171-172). El trabajo de fondo ($d=\alpha$) es el de 4.3.
- Qué me piden
- Dar una forma de Jordan de una matriz diagonalizable.
- Cómo lo pienso
- Si es diagonalizable, $J$ es diagonal.
- Qué uso
- Forma de Jordan (Def. 4.13-4.14, U p.167): si todos los valores propios son reales, $A$ es semejante a $J$, diagonal por bloques $J_s(\lambda)$ (valor propio en la diagonal y unos justo encima). Para cada $\lambda$: el nº de bloques es $d(\lambda)=n-\operatorname{rang}(A-\lambda I)$ y la suma de sus tamaños es $\alpha(\lambda)$.
- Por dónde empiezo
- Valores propios y diagonal.
Escribir la matriz: cada terna es una fila,
$$A=\begin{pmatrix}2&2&-5\cr3&7&-15\cr1&2&-4\end{pmatrix}.$$Polinomio (Sarrus): términos $+$:
$$(2-\lambda)(7-\lambda)(-4-\lambda)+2\cdot(-15)\cdot1+(-5)\cdot3\cdot2=-\lambda^3+5\lambda^2+22\lambda-116;$$términos $-$:
$$(-5)(7-\lambda)\cdot1+(2-\lambda)(-15)\cdot2+2\cdot3\cdot(-4-\lambda)=-119+29\lambda.$$Restando: $p(\lambda)=-\lambda^3+5\lambda^2-7\lambda+3$.
- Raíces: $p(1)=0$; Ruffini con $1$ sobre $(-1,5,-7,3)$ da $(-1,4,-3\mid0)$: $p(\lambda)=-(\lambda-1)^2(\lambda-3)$. $\lambda_1=1$ ($\alpha_1=2$), $\lambda_2=3$ ($\alpha_2=1$).
$d$ del doble:
$$A-I=\begin{pmatrix}1&2&-5\cr3&6&-15\cr1&2&-5\end{pmatrix}:$$la 2.ª fila es el triple de la 1.ª y la 3.ª es igual a la 1.ª, rango $1$, $d_1=3-1=2=\alpha_1$. El simple: $d_2=1=\alpha_2$ (Teor. 4.5).
Conclusión: $d_i=\alpha_i$ y las raíces son reales $\Rightarrow$ diagonalizable (Teor. 4.7). Su forma de Jordan es la diagonal, con tres bloques $1\times1$:
$$J=\begin{pmatrix}3&0&0\cr0&1&0\cr0&0&1\end{pmatrix}$$(o $\operatorname{diag}(1,1,3)$: única salvo el orden de los bloques).
Matriz de paso (no se pide, sirve de comprobación): $L_1:\ x+2y-5z=0$, $L_1=\langle(-2,1,0),(5,0,1)\rangle$; $L_3=\langle(1,3,1)\rangle$. Con
$$P=\begin{pmatrix}1&-2&5\cr3&1&0\cr1&0&1\end{pmatrix}$$($\det P=2$): $P^{-1}AP=\operatorname{diag}(3,1,1)$.
Sobre la oficial. Da $p$ y los valores propios y concluye «es diagonalizable» sin justificar $d_1=2$ y sin escribir la matriz de Jordan, que es lo que se pide. Con el criterio del examen («lo que no se justifica no se puntúa») no llegaría a la nota completa.
Error típico: «tiene un valor propio doble $\Rightarrow$ no diagonalizable $\Rightarrow$ $J$ con un uno». Lo que decide es el rango de $A-\lambda I$: aquí es $1$ y no hay ningún uno.
Ene. 2020 (tipo A), ejercicio 4 — tres valores propios reales distintos (irracionales) X Ene. 2020 (A) ej. 4
Hallar una matriz de Jordan de la matriz $\begin{pmatrix}1&-1&3\cr3&1&2\cr4&0&1\end{pmatrix}$.
Qué se pide. Solo $J$. Con tres raíces reales distintas no hace falta ni un vector propio.
- Qué me piden
- Forma de Jordan con tres valores reales distintos.
- Cómo lo pienso
- Tres raíces reales distintas: $J$ diagonal.
- Qué uso
- Forma de Jordan (Def. 4.13-4.14, U p.167): si todos los valores propios son reales, $A$ es semejante a $J$, diagonal por bloques $J_s(\lambda)$ (valor propio en la diagonal y unos justo encima). Para cada $\lambda$: el nº de bloques es $d(\lambda)=n-\operatorname{rang}(A-\lambda I)$ y la suma de sus tamaños es $\alpha(\lambda)$.
- Por dónde empiezo
- Polinomio y raíces.
Polinomio (adjuntos por la 2.ª columna, que tiene un cero; $(-1)^{1+2}=-1$):
$$|A-\lambda I|=\begin{vmatrix}1-\lambda&-1&3\cr3&1-\lambda&2\cr4&0&1-\lambda\end{vmatrix}=\big(3(1-\lambda)-8\big)+(1-\lambda)\big((1-\lambda)^2-12\big)=(1-\lambda)^3+9\lambda-17=-\lambda^3+3\lambda^2+6\lambda-16.$$- Raíces: con $\lambda=2$: $-8+12+12-16=0$. Ruffini sobre $(-1,3,6,-16)$ da $(-1,1,8\mid0)$: $p(\lambda)=-(\lambda-2)(\lambda^2-\lambda-8)$, y $\lambda^2-\lambda-8=0\iff\lambda=\frac{1\pm\sqrt{33}}2$ (discriminante $1+32=33\gt0$: reales y distintas; ninguna vale $2$ porque $4-2-8=-6\neq0$).
- Tres valores propios reales distintos $\Rightarrow$ diagonalizable (corolario de U p.163), y por U p.169 («si todos los autovalores son diferentes, $\Lambda$ coincide con $J$») y el Ej. 4.17: $$J=\begin{pmatrix}2&0&0\cr0&\frac{1+\sqrt{33}}2&0\cr0&0&\frac{1-\sqrt{33}}2\end{pmatrix}.$$ Si se quisiera $P$: para $\lambda=2$ sale $(1,11,4)$ y $A(1,11,4)=(2,22,8)$.
Sobre la oficial. El resultado es correcto, pero la expresión intermedia «$-9(x-1)+(1-x)^3$» está mal (se desarrolla a $-x^3+3x^2-12x+10$); lo correcto es $(1-x)^3+9(x-1)-8=(1-x)^3+9x-17=-x^3+3x^2+6x-16$, como dice la propia oficial a continuación. Su «Nota: la matriz de Jordan siempre existe» es cierta solo si las raíces de $p$ son reales (aquí lo son).
Error típico: desanimarse por $\sqrt{33}$ y pensar «no es diagonalizable porque los valores propios no son enteros»: basta con que sean reales y distintos.
Método B — No diagonalizable: los bloques se leen en el rango
Feb. 2025 (1.ª sem.), pregunta 4 — matriz de Jordan de $\begin{pmatrix}0&1\cr-1&2\end{pmatrix}$ X Feb. 2025 (1.ª) P4
Calcule la matriz de Jordan correspondiente a $\begin{pmatrix}0&1\cr-1&2\end{pmatrix}$.
Qué se pide. $2\times2$ con un valor propio doble: hay que decidir entre $A=\lambda I$ y un bloque $2\times2$.
- Qué me piden
- Forma de Jordan de una $2\times2$.
- Cómo lo pienso
- Hay que ver si el valor doble tiene uno o dos vectores propios.
- Qué uso
- Forma de Jordan (Def. 4.13-4.14, U p.167): si todos los valores propios son reales, $A$ es semejante a $J$, diagonal por bloques $J_s(\lambda)$ (valor propio en la diagonal y unos justo encima). Para cada $\lambda$: el nº de bloques es $d(\lambda)=n-\operatorname{rang}(A-\lambda I)$ y la suma de sus tamaños es $\alpha(\lambda)$. Tamaños de los bloques: el nº de bloques de tamaño $\ge k$ es $\operatorname{rang}(A-\lambda I)^{k-1}-\operatorname{rang}(A-\lambda I)^{k}$.
- Por dónde empiezo
- Polinomio y rango.
Polinomio:
$$|A-\lambda I|=-\lambda(2-\lambda)+1=\lambda^2-2\lambda+1=(\lambda-1)^2:$$$\lambda=1$, $\alpha=2$. Comprobación: $\operatorname{tr}A=2=1+1$ y $|A|=1=1\cdot1$.
$d$:
$$A-I=\begin{pmatrix}-1&1\cr-1&1\end{pmatrix}$$tiene rango $1$ (filas iguales, no nulas), $d=2-1=1\lt2$: no es diagonalizable (Teor. 4.7). (Otra forma: si lo fuera, $\Lambda=I$ y $A=PIP^{-1}=I$, falso.)
Forma de Jordan: el nº de bloques de $\lambda=1$ es $d=1$: un único bloque, de tamaño $\alpha=2$:
$$J=\begin{pmatrix}1&1\cr0&1\end{pmatrix}.$$Un solo «1» porque falta un vector no propio.
Matriz de paso (no se pide): $(A-I)^2=0$; $v=(1,0)\notin\operatorname{Nuc}(A-I)=\langle(1,1)\rangle$, $(A-I)v=(-1,-1)$.
$$P=\begin{pmatrix}-1&1\cr-1&0\end{pmatrix}$$ $$AP=\begin{pmatrix}-1&0\cr-1&-1\end{pmatrix}=PJ.$$
Sobre la oficial. Dice «la matriz tiene un único autovalor de multiplicidad algebraica $2$. Por tanto $A$ no es diagonalizable». El «por tanto» no está justificado: $I$ también tiene un único autovalor doble y es diagonal. Falta comprobar $d=1$ (o el argumento $A\neq I$). El $J$ final es correcto.
Error típico: saltarse la comprobación de $d$ (como la oficial) o poner el $1$ debajo de la diagonal (U lo coloca encima). Práctica relacionada: es el mismo razonamiento que 4.3, Feb. 2023 (2.ª), pregunta 4 ($\begin{pmatrix}9&9\cr0&9\end{pmatrix}$, cuya $J$ es $\begin{pmatrix}9&1\cr0&9\end{pmatrix}$).
Método C — Jordan con matriz de paso: cadena de núcleos
Sept. 2024, problema 8 (2.º apartado) — $J$ y matriz de paso de $\begin{pmatrix}3&0&0\cr0&4&-1\cr0&1&2\end{pmatrix}$ X Sept. 2024 P8 (2.º)
Considere $A=\begin{pmatrix}3&0&0\cr0&4&-1\cr0&1&2\end{pmatrix}$. Demuestre que $A$ no es diagonalizable (1 punto). Calcule la matriz de Jordan asociada a $A$ (1 punto).
Qué se pide. El 1.er apartado está desarrollado en 4.3 (Sept. 2024, P8, 1.º). Aquí, el 2.º: hallar $J$; la oficial construye además $P$.
- Qué me piden
- $J$ y matriz de paso de una $3\times3$.
- Cómo lo pienso
- Un valor propio triple; un bloque $1\times1$ y uno $2\times2$.
- Qué uso
- Forma de Jordan (Def. 4.13-4.14, U p.167): si todos los valores propios son reales, $A$ es semejante a $J$, diagonal por bloques $J_s(\lambda)$ (valor propio en la diagonal y unos justo encima). Para cada $\lambda$: el nº de bloques es $d(\lambda)=n-\operatorname{rang}(A-\lambda I)$ y la suma de sus tamaños es $\alpha(\lambda)$. Matriz de paso: por cada bloque, una cadena $v_1,\dots,v_s$ con $(A-\lambda I)v_1=0$ y $(A-\lambda I)v_{k}=v_{k-1}$; las cadenas por columnas forman $P$ y $AP=PJ$.
- Por dónde empiezo
- Cadena de dos vectores y un propio.
Lo ya sabido:
$$p(\lambda)=(3-\lambda)\big[(4-\lambda)(2-\lambda)+1\big]=(3-\lambda)(\lambda^2-6\lambda+9)=(3-\lambda)^3:$$$\lambda=3$, $\alpha=3$.
$$N:=A-3I=\begin{pmatrix}0&0&0\cr0&1&-1\cr0&1&-1\end{pmatrix}$$tiene rango $1$: $d=3-1=2\lt3$.
- Forma de los bloques (U p.169-170 y p.175-176): $d=2$ bloques de $\lambda=3$ con tamaños no crecientes que suman $3$: un bloque $2\times2$ y uno $1\times1$. Equivalentemente $\operatorname{rang}(A-3I)=1\Rightarrow J'$: $$J=\begin{pmatrix}3&1&0\cr0&3&0\cr0&0&3\end{pmatrix}.$$ Con esto ya está lo que se pide; la matriz de paso sirve de comprobación.
- Cadena de núcleos (guion p.173): $\operatorname{Nuc}N=\{y-z=0\}=\langle(1,0,0),(0,1,1)\rangle$ (dim $2$). Se calcula $N^2$: $$N^2=\begin{pmatrix}0&0&0\cr0&1\cdot1+(-1)\cdot1&1\cdot(-1)+(-1)(-1)\cr0&0&0\end{pmatrix}=0,$$ luego $\operatorname{Nuc}N^2=\mathbb R^3$ (dim $3=\alpha$). Se para.
- Cabeza: $v\in\operatorname{Nuc}N^2$ (cualquiera) y $\notin\operatorname{Nuc}N$ ($y\neq z$): $v=(0,0,1)$. Baja: $v_1:=Nv=(0,-1,-1)$, vector propio ($Nv_1=N^2v=0$).
- Completar: un propio independiente de $v_1$ para el bloque $1\times1$: $v_3=(1,0,0)$ (o el de la oficial, $(1,2,2)$, que cumple $y=z$).
- Orden de columnas: para que el $1$ quede encima, primero $v_1=(A-3I)v$ y después $v$ (U Ej. 4.18-4.19: «los vectores de cada bloque se calculan de derecha a izquierda»): la 2.ª columna de $J$, $(1,3,0)^t$, significa $Av=1\cdot v_1+3v$, es decir $Nv=v_1$. $$P=(v_1\mid v\mid v_3)=\begin{pmatrix}0&0&1\cr-1&0&0\cr-1&1&0\end{pmatrix},\qquad\det P=-1\neq0.$$
Comprobación sin invertir:
$$AP=\begin{pmatrix}0&0&3\cr-3&-1&0\cr-3&2&0\end{pmatrix}=PJ;$$ $$P^{-1}=\begin{pmatrix}0&-1&0\cr0&-1&1\cr1&0&0\end{pmatrix},$$$P^{-1}AP=J$.
Errores de la solución oficial (verificados): (a) factorización (aparece en su «Fallo recurrente»): escribe $(3-\lambda)(4-\lambda)(2-\lambda)+(3-\lambda)=(3-\lambda)\big((2-\lambda)+(3-\lambda)+1\big)$; lo correcto es $(3-\lambda)\big[(4-\lambda)(2-\lambda)+1\big]=(3-\lambda)^3$. (b) Escribe $(A+3I)$ y $(A+3I)^2$ donde debe decir $(A-3I)$, y «un vector de $(A-3I)$» por «del núcleo de $(A-3I)$». (c) En la matriz de paso pone «$-11$» por $-1$. (d) Orden de columnas: usa $P=(v\mid(A-3I)v\mid v_3)$ con $v=(0,0,1)$, es decir $v$ antes que $(A-3I)v$; con ese orden $P^{-1}AP=\begin{pmatrix}3&0&0\cr1&3&0\cr0&0&3\end{pmatrix}$ (el $1$ debajo), no la $J$ que escribe; hay que intercambiar las dos primeras columnas. La $J$ final es correcta; la igualdad $J=P^{-1}AP$ que escribe, no.
Error típico: poner $v$ antes que $Nv$ (el uno queda debajo); elegir $v\in\operatorname{Nuc}N$ (entonces $Nv=0$, no sirve); tomar el tercer vector proporcional a $Nv$ ($P$ singular).
Ejemplos propios con estos esquemas
Ejemplo propio 12 (examen) — $2\times2$: $J$, $P$ y $A^n$ propio
Halla la forma de Jordan $J$, una matriz de paso $P$ y $A^n$ para $A=\begin{pmatrix}3&1\cr-1&1\end{pmatrix}$ (la de E 4.36).
- Qué me piden
- $J$, $P$ y $A^n$ de una $2\times2$.
- Cómo lo pienso
- Cadena de longitud 2 y binomio del bloque.
- Qué uso
- Forma de Jordan (Def. 4.13-4.14, U p.167): si todos los valores propios son reales, $A$ es semejante a $J$, diagonal por bloques $J_s(\lambda)$ (valor propio en la diagonal y unos justo encima). Para cada $\lambda$: el nº de bloques es $d(\lambda)=n-\operatorname{rang}(A-\lambda I)$ y la suma de sus tamaños es $\alpha(\lambda)$. Matriz de paso: por cada bloque, una cadena $v_1,\dots,v_s$ con $(A-\lambda I)v_1=0$ y $(A-\lambda I)v_{k}=v_{k-1}$; las cadenas por columnas forman $P$ y $AP=PJ$.
- Por dónde empiezo
- $P$, $J$, $J^n$.
$\operatorname{tr}A=4$ y $|A|=3+1=4$: $p(\lambda)=\lambda^2-4\lambda+4=(\lambda-2)^2$. $A-2I$ tiene rango $1$: $d=1$,
$$J=\begin{pmatrix}2&1\cr0&2\end{pmatrix}.$$$(A-2I)^2=0$. $v=(1,0)\notin\operatorname{Nuc}(A-2I)=\langle(1,-1)\rangle$ y $Nv=(1,-1)$.
$$P=\begin{pmatrix}1&1\cr-1&0\end{pmatrix},$$$\det P=1$,
$$P^{-1}=\begin{pmatrix}0&-1\cr1&1\end{pmatrix},$$ $$AP=\begin{pmatrix}2&3\cr-2&-1\end{pmatrix}=PJ.$$Potencias:
$$J^n=\begin{pmatrix}2^n&n2^{n-1}\cr0&2^n\end{pmatrix}$$ $$A^n=PJ^nP^{-1}=2^{n-1}\begin{pmatrix}n+2&n\cr-n&2-n\end{pmatrix}.$$Comprobaciones: $n=1$ da $A$; $n=2$ da
$$\begin{pmatrix}8&4\cr-4&0\end{pmatrix}=A^2.$$
Ejemplo propio 13 (examen) — valor simple y valor doble con $d=1$ (como U Ej. 4.18) propio
Halla $J$ y $P$ para $A=\begin{pmatrix}2&1&1\cr0&2&0\cr0&1&1\end{pmatrix}$.
- Qué me piden
- $J$ y $P$ de una matriz con valor simple y valor doble con $d=1$.
- Cómo lo pienso
- Un bloque $1\times1$ y uno $2\times2$.
- Qué uso
- Forma de Jordan (Def. 4.13-4.14, U p.167): si todos los valores propios son reales, $A$ es semejante a $J$, diagonal por bloques $J_s(\lambda)$ (valor propio en la diagonal y unos justo encima). Para cada $\lambda$: el nº de bloques es $d(\lambda)=n-\operatorname{rang}(A-\lambda I)$ y la suma de sus tamaños es $\alpha(\lambda)$. Matriz de paso: por cada bloque, una cadena $v_1,\dots,v_s$ con $(A-\lambda I)v_1=0$ y $(A-\lambda I)v_{k}=v_{k-1}$; las cadenas por columnas forman $P$ y $AP=PJ$.
- Por dónde empiezo
- Cadena de longitud 2 y vector propio del simple.
Por la 1.ª columna:
$$|A-\lambda I|=(2-\lambda)\begin{vmatrix}2-\lambda&0\cr1&1-\lambda\end{vmatrix}=(2-\lambda)^2(1-\lambda).$$$\lambda=1$ simple y $\lambda=2$ con $\alpha=2$.
- $$N=A-2I=\begin{pmatrix}0&1&1\cr0&0&0\cr0&1&-1\end{pmatrix}$$
tiene rango $2$, $d_2=1\lt2$:
$$J=\begin{pmatrix}1&0&0\cr0&2&1\cr0&0&2\end{pmatrix}.$$ Cadena para $\lambda=2$: $\operatorname{Nuc}N$: $y+z=0,\ y-z=0\Rightarrow y=z=0$, $\langle(1,0,0)\rangle$.
$$N^2=\begin{pmatrix}0&1&-1\cr0&0&0\cr0&-1&1\end{pmatrix}:$$$\operatorname{Nuc}N^2=\{y=z\}=\langle(1,0,0),(0,1,1)\rangle$ (dim $2=\alpha$).
Cabeza
$$v=(0,1,1)\in\operatorname{Nuc}N^2\setminus\operatorname{Nuc}N;$$$Nv=(2,0,0)$.
$\lambda=1$:
$$A-I=\begin{pmatrix}1&1&1\cr0&1&0\cr0&1&0\end{pmatrix},$$$y=0,\ x=-z$: $u=(-1,0,1)$.
- $$P=(u\mid Nv\mid v)=\begin{pmatrix}-1&2&0\cr0&0&1\cr1&0&1\end{pmatrix},$$
$\det P=2$; el orden sigue el de $J$ (primero el bloque de $1$, luego el de $2$, con $Nv$ antes que $v$).
$$AP=\begin{pmatrix}-1&4&2\cr0&0&2\cr1&0&2\end{pmatrix}=PJ.$$
Error típico (comprobado con sympy). Tomar $v=(0,1,0)$ «porque no es propio»: es cierto que $v\notin\operatorname{Nuc}N$, pero tampoco está en $\operatorname{Nuc}N^2$ ($N^2v=(1,0,-1)\neq0$). Entonces $Nv=(1,0,1)$ no es propio y con $P=(u\mid Nv\mid v)$ sale $P^{-1}AP=\begin{pmatrix}1&-1&0\cr0&2&1\cr0&0&2\end{pmatrix}$, que no es de Jordan (el $-1$ no es $0$ ni $1$).
Ejemplo propio 14 (examen) — valor triple con un solo bloque $3\times3$ y $J^n$ propio
a) Sin calcular $P$, ¿qué formas de Jordan son posibles para una $3\times3$ con $|A-\lambda I|=-(\lambda-1)^3$? b) Para $A=\begin{pmatrix}0&1&0\cr0&0&1\cr1&-3&3\end{pmatrix}$ halla $J$ y $P$.
- Qué me piden
- Posibles $J$ con polinomio $-(\lambda-1)^3$ y cálculo para una matriz concreta.
- Cómo lo pienso
- Particiones de 3: tres formas posibles.
- Qué uso
- Forma de Jordan (Def. 4.13-4.14, U p.167): si todos los valores propios son reales, $A$ es semejante a $J$, diagonal por bloques $J_s(\lambda)$ (valor propio en la diagonal y unos justo encima). Para cada $\lambda$: el nº de bloques es $d(\lambda)=n-\operatorname{rang}(A-\lambda I)$ y la suma de sus tamaños es $\alpha(\lambda)$. Matriz de paso: por cada bloque, una cadena $v_1,\dots,v_s$ con $(A-\lambda I)v_1=0$ y $(A-\lambda I)v_{k}=v_{k-1}$; las cadenas por columnas forman $P$ y $AP=PJ$.
- Por dónde empiezo
- Rango de $A-I$ decide.
a) Según U p.175-176: $J=I$,
$$J'=\begin{pmatrix}1&1&0\cr0&1&0\cr0&0&1\end{pmatrix}$$ $$J''=\begin{pmatrix}1&1&0\cr0&1&1\cr0&0&1\end{pmatrix},$$según $\operatorname{rang}(A-I)=0,1,2$.
b) Por Sarrus
$$|A-\lambda I|=\lambda^2(3-\lambda)+1-3\lambda=-(\lambda-1)^3.$$ $$B=A-I=\begin{pmatrix}-1&1&0\cr0&-1&1\cr1&-3&2\end{pmatrix}:$$fila 3 $=-$fila 1 $+2\,$fila 2, rango $2$, $d=1$: $J=J''$.
Cadena: $\operatorname{Nuc}B=\langle(1,1,1)\rangle$;
$$B^2=\begin{pmatrix}1&-2&1\cr1&-2&1\cr1&-2&1\end{pmatrix}$$(rango $1$, $\operatorname{Nuc}B^2=\{x-2y+z=0\}$, dim $2$); $B^3=0$ (dim $3=\alpha$).
- Cabeza $v\notin\operatorname{Nuc}B^2$: $v=(0,0,1)$ ($0-0+1\neq0$). $Bv=(0,1,2)$ y $B^2v=(1,1,1)$ (propio).
- $$P=(B^2v\mid Bv\mid v)=\begin{pmatrix}1&0&0\cr1&1&0\cr1&2&1\end{pmatrix},$$
$\det P=1$,
$$P^{-1}=\begin{pmatrix}1&0&0\cr-1&1&0\cr1&-2&1\end{pmatrix},$$ $$AP=\begin{pmatrix}1&1&0\cr1&2&1\cr1&3&3\end{pmatrix}=PJ.$$ (Ampliación.)
$$J^n=\begin{pmatrix}1&n&\binom n2\cr0&1&n\cr0&0&1\end{pmatrix}$$(porque $J=I+N$ con $N^3=0$); p. ej.
$$J^5=\begin{pmatrix}1&5&10\cr0&1&5\cr0&0&1\end{pmatrix},$$y $A^n=PJ^nP^{-1}$ (verificado con sympy para $n=0,\dots,6$).
Práctica: «calcula ahora su $J$» (ya resueltos en 4.3)
En 4.3 se demostró que no son diagonalizables; aquí se pide el paso siguiente. Resultados (verificados con sympy, con la regla de U p.172 y p.176):
| Ejercicio (4.3) | Matriz y dato clave | Forma de Jordan |
|---|---|---|
| Feb. 2023 (2.ª), P4 | $\begin{pmatrix}9&9\cr0&9\end{pmatrix}$ | $\begin{pmatrix}9&1\cr0&9\end{pmatrix}$ |
| Feb. 2019 (B), P a) (= E 4.18) | $\begin{pmatrix}2&2&-2\cr0&3&-1\cr0&1&1\end{pmatrix}$, $\operatorname{rang}(A-2I)=1$ | $\begin{pmatrix}2&1&0\cr0&2&0\cr0&0&2\end{pmatrix}$ |
| Feb. 2026 (1.ª), P8, 2.º | caso $a=b$: $\begin{pmatrix}a&0&0\cr1&a&2\cr0&0&a\end{pmatrix}$, $\operatorname{rang}(A-aI)=1$ | $\begin{pmatrix}a&1&0\cr0&a&0\cr0&0&a\end{pmatrix}$ |
| Feb. 2024 (2.ª), P8 | caso $a=-2$ ($b\neq0$): $\operatorname{rang}(A+2I)=2$ | $\begin{pmatrix}-2&1&0\cr0&-2&0\cr0&0&2\end{pmatrix}$ |
| Feb. 2025 (1.ª), P8 | $b=5$ (cualquier $a$): $\operatorname{rang}(A-5I)=2$; $b=-1$, $a\neq0$: $\operatorname{rang}(A+I)=2$ | $\operatorname{diag}\big(-1,\begin{pmatrix}5&1\cr0&5\end{pmatrix}\big)$; $\operatorname{diag}\big(\begin{pmatrix}-1&1\cr0&-1\end{pmatrix},5\big)$ |
| Sept. 2024, P8, 1.º | su 2.º apartado está desarrollado arriba (ex-sep24-p8b) | $\begin{pmatrix}3&1&0\cr0&3&0\cr0&0&3\end{pmatrix}$ |
En los temas 5 y 6 aparecen «formas canónicas» de formas cuadráticas (X Feb. 2018 A, Feb. 2022 A2 ej. 6, Feb. 2024 A p.6, Sept. 2025 P8). No son formas de Jordan: allí se busca congruencia $P^tAP$, no semejanza $P^{-1}AP$. Y recuerda: una matriz simétrica real (tema 5.3, p. ej. X Ene. 2020 A ej. 5) siempre es diagonalizable (incluso con $P$ ortogonal), así que su forma de Jordan es siempre diagonal.
10Chuleta, lista de comprobación y errores frecuentes
| Resultado | Fórmula / idea (hipótesis) |
|---|---|
| Teor. 4.8 (U p.168) | $\exists P$ regular con $J=P^{-1}AP$, $J=\operatorname{diag}(J_1,\dots,J_s)$ (todas las raíces de $|A-\lambda I|$ reales; si no, solo sobre $\mathbb C$) |
| Bloque $J_s(\lambda)$ | $\lambda$ en la diagonal, $1$ en la línea justo encima; $=\lambda I+N$ con $N^s=0$ |
| Unicidad | $J$ única salvo el orden de los bloques; $P$ no es única |
| Nº de bloques | de $\lambda_i$: $d_i=n-\operatorname{rang}(A-\lambda_iI)$; total $\sum d_i$ |
| Tamaños | suman $\alpha_i$; con $\alpha\le3$ los fija $d$; con $\alpha\ge4$ hacen falta los rangos de $N^k$ |
| Nº de unos | $n-\sum d_i$ = nº de vectores no propios de la base |
| Regla de los rangos | nº de bloques de tamaño $\ge k$ $=r_{k-1}-r_k$, $r_k=\operatorname{rang}(A-\lambda I)^k$, hasta $n-r_k=\alpha$ (por $\lambda$) |
| $A$ diagonalizable | $J=\Lambda$ (sin unos); $n$ valores distintos $\Rightarrow$ diagonal |
| Un único valor propio $\lambda$ | $J=\lambda I$ solo si $A=\lambda I$; si no, hay algún bloque de tamaño $\ge2$ |
| Subespacios propios generalizados (Def. 4.13) | $\operatorname{Nuc}(A-\lambda I)^i$, $i\ge1$; cadena creciente que se estabiliza en dim $=\alpha$ |
| Cadena de un bloque de tamaño $k$ | $v\in\operatorname{Nuc}N^k\setminus\operatorname{Nuc}N^{k-1}$; columnas $(N^{k-1}v\mid\cdots\mid Nv\mid v)$: propio primero |
| Comprobar | $\det P\neq0$ y $AP=PJ$ (sin invertir) |
| Potencias de un bloque | $J_2(\lambda)^n=\begin{pmatrix}\lambda^n&n\lambda^{n-1}\cr0&\lambda^n\end{pmatrix}$; $A^n=PJ^nP^{-1}$ (no está en U) |
| Método de potencias (U p.185) | $x_{k+1}=Ax_k/\mu_k$, $\mu_k$ = coordenada de mayor $|\cdot|$ con su signo; $\mu_k\to\lambda_1$, $x_k\to v_1$ (valor estrictamente dominante, $x_0$ con componente en $v_1$); con $A^{-1}$, el de menor módulo |
| Autoevaluación (U p.312-313) | 9 F, 10 V |
- Sé enunciar el Teorema de Jordan con su hipótesis (raíces reales) y explicar por qué cada bloque es $\lambda I+N$ con $N$ nilpotente (U Teor. 4.8, p.169; E 4.33).
- Sé calcular el nº de bloques de cada valor propio ($d_i$) y el nº de unos de $J$, y distinguir «nº de bloques» de «nº de valores propios» (E 4.33-4.35, 4.37; autoevaluación 10).
- Sé escribir todas las $J$ posibles para $\alpha\le4$ y reconocer cuándo $d$ no basta (U Ej. 4.16, 4.22; ejemplo propio 5).
- Sé usar los rangos de $(A-\lambda I)^k$ para decidir los tamaños de los bloques cuando $\alpha\ge4$ (E 4.39; ejemplos propios 6 y 10).
- Sé hallar $J$ de una $2\times2$ y de una $3\times3$ (E 4.36, U Ej. 4.17-4.21; X Feb. 2025, Ene. 2020, Ene. 2022).
- Sé construir la matriz de paso con cadenas: núcleos, cabeza $v\in\operatorname{Nuc}N^k\setminus\operatorname{Nuc}N^{k-1}$, propio primero (U Ej. 4.18-4.19; E 4.40; X Sept. 2024).
- Sé comprobar si una $P$ dada es matriz de paso ($\det P\neq0$ y $AP=PJ$) y por qué no es única (E 4.38).
- Sé reconocer cuándo no hay $J$ real (raíces complejas) y que una diagonalizable tiene $J=\Lambda$ (ejemplo propio 9; X Ene. 2020, Ene. 2022).
- Sé aplicar el método de potencias con $\mu_k$ con su signo y justificar la convergencia por $|\lambda_2/\lambda_1|$ (E 4.41; U Ej. 4.24-4.30). No cae en examen.
- «No diagonalizable ⇒ no hay $J$»: la forma de Jordan existe siempre (raíces reales). Recuerda también que «diagonalizable ⇒ su $J$ es diagonal»; y «hay unos en $J$ ⇒ $A$ no es diagonalizable» es verdadero.
- Contar bloques por valores propios distintos en lugar de por vectores propios independientes $d_i$ (E 4.37).
- Decidir $(3,1)$ o $(2,2)$ solo con $d$ cuando $\alpha=4$ y $d=2$: hacen falta los rangos de $(A-\lambda I)^2$ (E 4.39).
- Olvidar que las raíces deben ser reales (ejemplo propio 9; giro de $90^\circ$).
- Los unos debajo de la diagonal (por el orden de las columnas de $P$) en lugar de encima (X Sept. 2024 oficial).
- Cabeza mal elegida: $v\in\operatorname{Nuc}N^{k-1}$ (la cadena se corta) o $v\notin\operatorname{Nuc}N^k$ (no es propio el último). Ambas condiciones a la vez.
- $P$ singular por un propio extra dependiente; comprobar siempre $\det P\neq0$.
- Mezclar $A-\lambda_1I$ y $A-\lambda_2I$ en la misma cadena: cada valor propio con la suya.
- Dar por válida una solución oficial sin completarla: Ene. 2022 no llega a escribir $J$, Feb. 2025 no comprueba $d$ y Sept. 2024 tiene erratas de orden de columnas y de factorización.
- Método de potencias: usar $|\mu_k|$ (se pierde el signo), elegir un $x_0$ que sea propio de un valor no dominante, o confiar sin $|\lambda_1|\gt|\lambda_2|$.
11Mapa del tema 4: qué hacer ante una matriz propio
- Valores propios (tema 4.2): $|A-\lambda I|=0$, con sus $\alpha_i$. Si alguna raíz no es real: no es diagonalizable en $\mathbb R$ y no hay $J$ real.
- Para cada $\lambda_i$ repetido: $d_i=n-\operatorname{rang}(A-\lambda_iI)$ y comparar con $\alpha_i$ (los simples cumplen $d=\alpha=1$). Es el veredicto.
- $d_i=\alpha_i$ para todo $i$ $\Rightarrow$ diagonalizable (tema 4.3): $\Lambda$ con los $\lambda_i$ repetidos $\alpha_i$ veces; $P$ = bases de los $L_i$ en columnas, en el mismo orden. $J=\Lambda$.
- Algún $d_i\lt\alpha_i$ $\Rightarrow$ Jordan (este tema): tamaños de los bloques por $(\alpha,d)$ y, si $\alpha\ge4$, por los rangos de $(A-\lambda I)^k$. Escribe $J$.
- $P$ por cadenas: para cada bloque, $v\in\operatorname{Nuc}N^k\setminus\operatorname{Nuc}N^{k-1}$, y columnas $(N^{k-1}v\mid\cdots\mid v)$; completa con propios independientes.
- Comprueba $\det P\neq0$ y $AP=P\Lambda$ (o $AP=PJ$). Si te piden potencias: $A^k=P\Lambda^kP^{-1}$ o $PJ^kP^{-1}$.
| Situación | Qué calcular | Resultado |
|---|---|---|
| $A$ y $B$ semejantes (o no) | invariantes: polinomio característico, traza, det, rango de $A-\lambda I$ y de sus potencias (4.1, 4.4) | un invariante distinto prueba que no lo son; polinomio, traza, det y $d$ iguales no bastan si $\alpha\ge4$; si además coinciden $\operatorname{rang}(A-\lambda I)^k$ para todo $\lambda$ y $k$ (raíces reales), sí son semejantes (misma $J$) |
| $n$ valores propios reales distintos | nada más (4.3) | $A\sim\Lambda$ y $J=\Lambda$ |
| Hay repetidos, $d_i=\alpha_i$ siempre | bases de los $L_i$ (4.3) | $\Lambda$, $P$ de vectores propios |
| Hay repetidos, algún $d_i\lt\alpha_i$ | tamaños de bloques, $J$, cadenas (4.4) | $J$ con unos, $P$ de cadenas |
| Raíces no reales | nada en $\mathbb R$ | sin $\Lambda$ ni $J$ real |
| $A$ simétrica real | (tema 5.3) | siempre diagonalizable, con $P$ ortogonal |
| Valor propio dominante, sin resolver $|A-\lambda I|$ | método de potencias (4.4.3) | $\mu_k\to\lambda_1$, $x_k\to v_1$ |