Forma canonica di Jordan

Video brevi n. 51–53 della playlist ("Forma canonica di Jordan", parti 1–3). Il prof la presenta come complemento: all'esame si chiede di riconoscere se una matrice è diagonalizzabile, non di calcolarne la forma di Jordan. Prerequisiti: DiagonalizzazioneUna matrice è diagonalizzabile se è simile a una diagonale, cioè se esiste una base di autovettori: allora A = S D S⁻¹ con gli autovettori nelle colonne di S e gli autovalori in D. Criterio: tutti gli autovalori nel campo e molteplicità geometrica uguale a quella algebrica. Le matrici simmetriche reali hanno autovalori reali.Diagonalizzazione →, Autovalori e autovettoriUn autovettore è un vettore non nullo che una funzione lineare manda in un suo multiplo; si trovano gli autovalori come radici del polinomio caratteristico det(A − λI) e gli autovettori come nucleo di A − λI. Matrici simili hanno gli stessi autovalori.Autovalori e autovettori →. Applicazione: Esponenziale di una matrice e sistemi differenziali lineariIl sistema di equazioni differenziali y' = Ay ha soluzione y(x) = e^{Ax} ȳ. L'esponenziale di una matrice si definisce con la serie di e^x; per una matrice diagonalizzabile A = S D S⁻¹ vale e^{Ax} = S e^{Dx} S⁻¹, e e^{Dx} ha sulla diagonale gli esponenziali e^{λx} degli autovalori.Esponenziale di una matrice e sistemi differenziali lineari →.

Il problema

Una matrice è diagonalizzabile se e solo se ha tutti gli autovalori nel campo e, per ognuno, molteplicità geometrica uguale a quella algebrica. Se la seconda condizione fallisce, per esempio

$$A = \begin{pmatrix} 3 & 1 \\ -1 & 1 \end{pmatrix}, \qquad \det(A - \lambda I) = (3 - \lambda)(1 - \lambda) + 1 = (\lambda - 2)^2,$$

con $A - 2I = \begin{pmatrix} 1 & 1 \\ -1 & -1 \end{pmatrix}$ di rango $1$ (quindi un solo autovettore indipendente, $m_g = 1 < 2 = m_a$), non esiste una base di autovettori. Domanda: qual è la matrice "più semplice possibile" simile ad $A$?

Blocchi e matrici di Jordan

Definizione. Il blocco di Jordan di ordine $k$ relativo a $\lambda$ è la matrice $k \times k$ $$J_k(\lambda) = \begin{pmatrix} \lambda & 1 & & \\ & \lambda & \ddots & \\ & & \ddots & 1 \\ & & & \lambda \end{pmatrix}$$ con $\lambda$ sulla diagonale, $1$ subito sopra la diagonale e $0$ altrove. Una matrice di Jordan è una matrice diagonale a blocchi i cui blocchi sono blocchi di Jordan.

Esempi: $J_1(\lambda) = (\lambda)$; $J_2(2) = \begin{pmatrix} 2 & 1 \\ 0 & 2 \end{pmatrix}$; $J_3(0) = \begin{pmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ 0 & 0 & 0 \end{pmatrix}$. Una matrice diagonale è una matrice di Jordan con tutti i blocchi di ordine $1$.

Che cosa fa un blocco sulla base. Se in una base $\{v_1, \dots, v_k\}$ la matrice di $f$ è $J_k(\lambda)$, leggendo le colonne:

$$f(v_1) = \lambda v_1, \qquad f(v_2) = v_1 + \lambda v_2, \qquad \dots, \qquad f(v_k) = v_{k-1} + \lambda v_k ,$$

cioè, con $N = f - \lambda\,\mathrm{id}$:

$$N v_1 = 0, \qquad N v_2 = v_1, \qquad \dots, \qquad N v_k = v_{k-1} .$$

$v_1$ è un autovettore; gli altri non lo sono, ma $N$ li fa "scendere" uno sull'altro fino ad arrivare all'autovettore. Una sequenza così si chiama catena di Jordan, e i $v_i$ autovettori generalizzati (vale $N^i v_i = 0$).

Il teorema

Teorema (forma canonica di Jordan). Sia $A \in M_n(K)$ con tutti gli autovalori in $K$ (sempre vero per $K = \mathbb{C}$). Allora $A$ è simile a una matrice di Jordan $J$, unica a meno dell'ordine dei blocchi: $$P^{-1} A P = J .$$

Come si leggono i blocchi. Per ogni autovalore $\lambda$:

• numero di blocchi relativi a $\lambda$ = $m_g(\lambda)$ = $\dim \operatorname{Ker}(A - \lambda I)$ (ogni blocco contiene esattamente un autovettore indipendente, il primo vettore della catena);
• somma degli ordini dei blocchi relativi a $\lambda$ = $m_a(\lambda)$ (sulla diagonale di $J$, $\lambda$ compare $m_a$ volte);
• ordini dei singoli blocchi: si ricavano dalle dimensioni di $\operatorname{Ker}(A - \lambda I)^j$ per $j = 1, 2, \dots$. Il numero di blocchi di ordine $\ge j$ è $\dim \operatorname{Ker}(A - \lambda I)^j - \dim \operatorname{Ker}(A - \lambda I)^{j-1}$.

Fino a ordine $3$ bastano le prime due regole:

Da ordine $4$ in poi possono servire le potenze: con $m_a = 4$ e $m_g = 2$ i blocchi possono essere $J_3, J_1$ oppure $J_2, J_2$.

$A$ è diagonalizzabile esattamente quando tutti i blocchi hanno ordine $1$, cioè $m_g = m_a$ per ogni autovalore: è il criterio di diagonalizzabilità visto da un'altra angolazione.

Come si trova $P$

Per un blocco $J_k(\lambda)$ serve una catena $v_1, \dots, v_k$ con $N v_1 = 0$, $N v_{i} = v_{i-1}$ ($N = A - \lambda I$). Due modi:

• dal basso: si sceglie $v_1$ autovettore e si risolvono successivamente $N v_2 = v_1$, $N v_3 = v_2$, … (attenzione: non ogni autovettore va bene come $v_1$ quando $m_g \ge 2$, perché il sistema potrebbe non avere soluzione);
• dall'alto (più sicuro): si sceglie $v_k \in \operatorname{Ker} N^k$ ma $v_k \notin \operatorname{Ker} N^{k-1}$, e si pone $v_{k-1} = N v_k$, $v_{k-2} = N v_{k-1}$, …

Le colonne di $P$ sono le catene, una dopo l'altra, ciascuna in ordine $v_1, v_2, \dots, v_k$ (prima l'autovettore).

Esempio 1: la matrice $2 \times 2$ di sopra

$\lambda = 2$, $m_a = 2$, $m_g = 1$: un blocco $J_2(2)$.

• Autovettore: $(A - 2I)x = 0$ dà $x_1 + x_2 = 0$, quindi $v_1 = (1, -1)$.
• Secondo vettore: $(A - 2I) v_2 = v_1$, cioè $x_1 + x_2 = 1$ (la seconda equazione, $-x_1 - x_2 = -1$, è la stessa): per esempio $v_2 = (1, 0)$.

$$P = \begin{pmatrix} 1 & 1 \\ -1 & 0 \end{pmatrix}, \qquad P^{-1} A P = \begin{pmatrix} 2 & 1 \\ 0 & 2 \end{pmatrix}.$$

Verifica senza invertire $P$: $A v_1 = (3 - 1, \ -1 - 1) = (2, -2) = 2v_1$ ✓ e $A v_2 = (3, -1) = v_1 + 2v_2 = (1 + 2, \ -1 + 0)$ ✓.

Esempio 2: la matrice dell'Esercizio 62 · matrice 3×3 non diagonalizzabile

$$A = \begin{pmatrix} 2 & 0 & 1 \\ -2 & -4 & 0 \\ 0 & 0 & 2 \end{pmatrix}, \qquad \lambda = 2 \ (m_a = 2, \ m_g = 1), \quad \lambda = -4 \ (m_a = m_g = 1).$$

Forma di Jordan: un blocco $J_2(2)$ e un blocco $J_1(-4)$.

Per $\lambda = 2$, $N = A - 2I = \begin{pmatrix} 0 & 0 & 1 \\ -2 & -6 & 0 \\ 0 & 0 & 0 \end{pmatrix}$, $\operatorname{Ker} N = \langle (3, -1, 0) \rangle$. Si procede dall'alto: serve $v_2$ con $N^2 v_2 = 0$ ma $N v_2 \ne 0$. $N^2 = \begin{pmatrix} 0 & 0 & 0 \\ 12 & 36 & -2 \\ 0 & 0 & 0 \end{pmatrix}$, quindi $\operatorname{Ker} N^2$: $12x_1 + 36x_2 - 2x_3 = 0$. Si prende $v_2 = (1, 0, 6)$ (soddisfa $12 - 12 = 0$ e non sta in $\operatorname{Ker} N$). Allora $v_1 = N v_2 = (6, -2, 0)$, che è un autovettore (multiplo di $(3, -1, 0)$) ✓.

Per $\lambda = -4$: $A + 4I = \begin{pmatrix} 6 & 0 & 1 \\ -2 & 0 & 0 \\ 0 & 0 & 6 \end{pmatrix}$ dà $x_1 = x_3 = 0$: $u = (0, 1, 0)$.

$$P = \begin{pmatrix} 6 & 1 & 0 \\ -2 & 0 & 1 \\ 0 & 6 & 0 \end{pmatrix}, \qquad P^{-1} A P = \begin{pmatrix} 2 & 1 & 0 \\ 0 & 2 & 0 \\ 0 & 0 & -4 \end{pmatrix}.$$

Esempio 3: un solo blocco $3 \times 3$

La matrice dell'Esercizio 69 · diagonalizzabilità al variare di h e matrice P (lezione 31) per $h = -1$, $A = \begin{pmatrix} 0 & 1 & 0 \\ -1 & 0 & 1 \\ 0 & 1 & 0 \end{pmatrix}$, ha l'unico autovalore $0$ con $m_a = 3$ e $m_g = 1$: la forma di Jordan è $J_3(0)$. Coerentemente $A^2 = \begin{pmatrix} -1 & 0 & 1 \\ 0 & 0 & 0 \\ -1 & 0 & 1 \end{pmatrix} \ne 0$ ma $A^3 = 0$, come $J_3(0)^3 = 0$ (matrici simili hanno potenze simili).

A che cosa serve: potenze ed esponenziale

Si scrive $J_k(\lambda) = \lambda I + N$, con $N$ la matrice con gli $1$ sopra la diagonale. $\lambda I$ e $N$ commutano e $N^k = 0$ ($N$ sposta ogni vettore di base "di un posto", dopo $k$ passi non resta nulla). Allora la formula del binomio si tronca e, per l'esponenziale,

$$e^{J_k(\lambda) x} = e^{\lambda x}\left(I + N x + \frac{N^2 x^2}{2!} + \dots + \frac{N^{k-1} x^{k-1}}{(k-1)!}\right).$$

Per $k = 2$: $e^{J_2(\lambda) x} = e^{\lambda x}\begin{pmatrix} 1 & x \\ 0 & 1 \end{pmatrix}$. Ecco da dove vengono, nelle equazioni differenziali con radici multiple dell'equazione caratteristica, le soluzioni del tipo $x e^{\lambda x}$: è il caso non diagonalizzabile del metodo dell'esponenziale di una matriceIl sistema di equazioni differenziali y' = Ay ha soluzione y(x) = e^{Ax} ȳ. L'esponenziale di una matrice si definisce con la serie di e^x; per una matrice diagonalizzabile A = S D S⁻¹ vale e^{Ax} = S e^{Dx} S⁻¹, e e^{Dx} ha sulla diagonale gli esponenziali e^{λx} degli autovalori.Esponenziale di una matrice e sistemi differenziali lineari →, con $e^{Ax} = P\, e^{Jx} P^{-1}$.

Errori comuni

• Mettere gli $1$ sotto la diagonale invece che sopra: con la convenzione $Nv_i = v_{i-1}$ e le catene in colonna nell'ordine $v_1, \dots, v_k$, gli $1$ stanno sopra.
• Ordinare la catena al contrario in $P$ (prima $v_k$, poi l'autovettore): si ottengono gli $1$ sotto la diagonale.
• Prendere come $v_2$ un altro autovettore: $v_2$ deve risolvere $N v_2 = v_1$, quindi non è un autovettore.
• Credere che $m_a$ e $m_g$ bastino sempre: da ordine $4$ servono anche le dimensioni dei nuclei delle potenze.