L'intento di questo articolo è di affrontare, si spera sempre nel modo più corretto e chiaro possibile, la trattazione di un argomento importante nell'ambito dell'algebra lineare quale quello della diagonalizzazione.
Cercherò di affrontare l'argomento in maniera graduale percorrendo le varie tappe necessarie alla comprensione dell'argomento nel suo insieme.
Data la lunghezza dell'articolo, per orientarsi nella sua consultazione, si consiglia di cliccare sul pulsante Mostra indice posto in alto.
Buona lettura.
Indice |
Polinomi di matrici
Consideriamo, su un campo K, il polinomio seguente:

Data una matrice A quadrata, essa può essere definita come:

dove I è la nota matrice identità.
In particolare, diremo che A è una radice, o uno zero, del polinomio f(t) se
.
Esempio
Abbiamo la matrice:
; si avrà quindi che
Date inoltre:
,
avremo quanto segue.
e
Si deduce quindi che A è uno zero di g(t).
Enunciamo ora il seguente teorema.
Siano f e g dei polinomi su K. Per ogni matrice A e per ogni scalare k,
-
;
-
;
-
;
-
(due polinomi in A commutano).
Le matrici e gli operatori lineari
Consideriamo ora un operatore lineare
su uno spazio vettoriale V. Le potenze di T sono definite mediante l'operazione di composizione:

Per ogni polinomio
definiamo f(T) in maniera uguale a quanto fatto per le matrici, cioè:

dove I è l'applicazione identità. Si dirà quindi che T èuna radice o zero di f(t) se
.
Va precisato che quanto esposto al teorema prima enunciato vale tanto per gli operatori lineari quanto per le matrici.
Polinomio caratteristico e teorema di Cayley - Hamilton
Prendiamo ora in considerazione una matrice quadrata A di ordine n.
La matrice
, dove In è la matrice quadrata identità di dimensione n e t è una incognita, può essere ottenuta se si sottrae t dalla diagonale A.
L'opposta della matrice M è la matrice
ed il suo determinante:

è un polinomio di grado n in t d è più precisamente detto polinomio caratteristico di A.
Passiamo ora ad enunciare il seguente importante teorema dovuto a Cayley ed Hamilton.
Il teorema di Cayley - Hamilton afferma che ogni matrice A è uno zero del suo polinomio caratteristico.
Esempio
Facciamo un piccolo esempio a riguardo.
Prendiamo in considerazione la matrice A (quadrata):

Scriviamone il polinomio caratteristico:

Possiamo verificare quanto dichiarato dal teorema di Cayley - Hamilton e cioè che A è uno zero di Δ(t):

Si dimostra inoltre che matrici simili hanno lo stesso polinomio caratteristico.
Alcune formule per il calcolo di polinomi caratteristici di matrici di ordine 2 e 3
Ordine 2
Data la matrice
, allora :

dove tr(A) è la traccia di A, cioè la somma degli elementi diagonali.
Ordine 3
Data la matrice
, allora :

dove
sono i cofattori degli elementi diagonali.
Possiamo enunciare quindi il seguente teorema.
Data una matrice quadrata di ordine n, il suo polinomio caratteristico è :

dove Sk è la somma dei minori principali di ordine k.
Esempio
Scriviamo il polinomio caratteristico della matrice :

Si ha
. I cofattori degli elementi sulla diagonale sono:

Quindi
. Inoltre
Conseguentemente:
Polinomio caratteristico per un operatore lineare
Consideriamo ancora un operatore lineare
su uno spazio vettoriale V di dimensione finita.
Il polinomio caratteristico Δ(t) di T è definito come il polinomio caratteristico di una qualunque matrice che rappresenta T.
Se A e B sono matrici di rappresentazione di T, si ha
, dove P è la matrice del cambiamento di base. A e B sono quindi simili e giacché due matrici simili hanno lo stesso polinomio caratteristico (lo abbiamo detto precedentemente), esse hanno lo stesso polinomio caratteristico.
Il polinomio caratteristico di T è quindi indipendente dalla base rispetto a cui è calcolata la sua matrice di rappresentazione.
se e solo se
, dove f(t) è un qualunque polinomio ed A è una matrice qualunque che rappresenta T: fatte salve tale considerazioni si può enunciare il seguente teorema per un operatore lineare.
Un'operatore lineare T è uno zero del suo polinomio caratteristico.
La diagonalizzazione: autovalori ed autovettori
Consideriamo la nostra matrice quadrata A di ordine n.
Diremo che è possibile rappresentare A con una matrice diagonale:

se e solo se esiste una base S costituita da vettori (colonna)
, tali che:

In questo caso A si dice diagonalizzabile.
Inoltre: sia A una matrice qualunque e sia λ uno scalare appartenente a K. Ebbene λ si dirà autovalore di A se esiste un vettore (colonna) non nullo
tale che :

Si dirà invece autovettore di A relativo all'autovalore λ ogni vettore che soddisfa tale relazione.
L'insieme Eλ di tutti gli autovettori relativi a λ è un sottospazio di V e si chiama autospazio di λ.
Enunciamo ora il seguente teorema.
Una matrice quadrata A di dimensione n è simile ad una matrice diagonale D se e solo se A ha n autovettori linearmente indipendenti. In tal caso gli elementi diagonali di D soon i corrispondenti autovalori e
, dove P è la matrice le cui colonne sono gli autovettori.
Se una matrice A può essere diagonalizzata si dice che A ammette una fattorizzazione diagonale

Con tale fattorizzazione l'algebra di A viene quindi ad essere l'algebra della matrice diagonale D, calcolabile facilmente.
Le proprietà degli autovalori e degli autovettori
Elenchiamo nel teorema seguente le proprietà che occorreranno per trovare una rappresentazione diagonale di una matrice quadrata.
Sia A una matrice quadrata. Sono equivalenti le seguenti affermazioni:
- lo scalare λ è un autovalore di A ;
- la matrice M = A − 2λI è singolare ;
- lo scalare λ è una radice del polinomio caratteristico Δ(t) di A.
L'autospazio Eλ è quindi lo spazio delle soluzioni del sistema omogeneo MX = 0, dove
, ottenuta cioè sottraendo λ dagli elementi sulla diagonale di A.
Alcune matrici possono non avere autovalori (e quindi autovettori). Ricorrendo però al Teorema fondamentale dell'algebra (ogni polinomio su
ha almeno una radice) e al teorema prima esposto si può ottenere il risultato che segue.
Sia A una matrice quadrata sul campo complesso
, allora A ha almeno un autovalore.
Consideriamo adesso la matrice A e supponiamo che λ sia un suo autovalore.
Per molteplicità algebrica di λ si intende la molteplicità di λ come polinomio caratteristico di A.
La molteplicità geometrica di λ è la dimensione del suo autospazio, in simboli : dimEλ.
Alla luce delle considerazioni finora fatte vale il seguente teorema.
Sia λ un autovalore di una matrice A: allora la molteplicità geometrica di λ non eccede la sua molteplicità algebrica.
Diagonalizzazione di operatori lineari
Cominciamo col dire che un operatore lineare T su uno spazio vettoriale V è detto diagonalizzabile se può essere rappresentato da una matrice diagonale D. Dunque T sarà diagonalizzabile se e solo se esiste una base di V

per la quale si ha:

In tal caso T è rappresentato dalla matrice diagonale:

relativa alla base S prima citata.
In maniera analoga a quanto fatto per le matrici procediamo col dare le seguenti definizioni.
Uno scalare
si dice autovalore di T se esiste un vettore
non nullo per il quale

Ogni vettore che soddisfa questa relazione è detto autovettore di T relativo all'autovalore λ.
L'insieme Eλ di tutti questi vettori è un sottospazio di V ed è detto autospazio di λ.
Senza ripetere quanto già detto per le matrici, anche per gli operatori lineari valgono le medesime definizioni di molteplicità algebrica e geometrica.
Per completare enunciamo infine i seguenti teoremi che si applicano ad un operatore lineare T su uno spazio vettoriale di dimensione finita V.
T si può rappresentare con una matrice diagonale D se e solo se esiste una base S di V costituita da autovettori di T. In tal caso gli elementi diagonali di D soon i corrispondenti autovalori.
Sia T un operatore lineare. Sono equivalenti le seguenti affermazioni:
- lo scalare λ è un autovalore di T ;
- l'operatore lineare λI − T è singolare ;
- lo scalare λ è una radice del polinomio caratteristico Δ(t) di T .
Dato uno spazio vettoriale complesso, allora T ha almeno un autovalore.
Supponiamo che
siano gli autovettori diversi da zero di T e relativi a distinti autovalori
, allora
sono linearmente indipendenti.
Supponiamo che il polinomio caratteristico di T sia prodotto di n fattori distinti

Ebbene T può essere rappresentato dalla matrice diagonale

La molteplicità geometrica di un autovalore λ di T non supera la sua molteplicità algebrica.
Infine rilevante è anche il teorema che segue che riduce la ricerca della diagonalizzazione di un operatore lineare T a quella di una matrice A.
Se A è la rappresentazione matriciale di T, T è diagonalizzabile se e solo se A è diagonalizzabile.
Calcolo di autovalori ed autovettori e diagonalizzazione delle matrici
L'algoritmo
Diamo ora un algoritmo di diagonalizzazione per una matrice quadrata di ordine n.
Passo 1. Trovare il polinomio caratteristico di A.
Passo 2. Trovare le radici del polinomio caratteristico di A, ossia trovarne gli autovalori.
Passo 3. Ripetere 1. e 2. per ogni autovalore λ di A.
- Determinare M = A − λI sottraendo λ dalla diagonale di A.
- Trovare una base per lo spazio delle soluzioni del sistema omogeneo MX = 0.
Passo 4. Consideriamo la collezione
di tutti gli autovettori che abbiamo ottenuto nel Passo 3.
- Se
A non è diagonalizzabile.
- Se m = n, A è diagonalizzabile. Detta P la matrice le cui colonne sono gli autovettori
, allora

dove λi è l'autovalore corrispondente all'autovettore vi.
Esempio
Applichiamo l'algoritmo di diagonalizzazione espresso ad una matrice A quadrata di ordine n:

Calcoliamo il polinomio caratteristico:

quindi

Le radici λ1 e λ2 sono rispettivamente 5 e − 2 e sono quindi gli autovalori di A.
Troviamo un autovettore v1 di A relativo a λ1 = 5, sottraendo quest'ultimo dalla diagonale di A ottenendo la matrice

Gli autovettori relativi a λ1 = 5 formano la soluzione del sistema omogeneo MX = 0:
oppure 
Il sistema ha una sola variabile libera. Una soluzione diversa da zero, ad esempio v1 = (2,1), è un autovettore che genera l'autospazio di λ1 = 5.
Troviamo un autovettore v2 di A relativo a λ2 = − 2, sottraendo quest'ultimo dalla diagonale di A ottenendo la matrice


Il sistema ha una sola variabile libera. Una soluzione diversa da zero, ad esempio v2 = ( − 1,3), è un autovettore che genera l'autospazio di λ2 = − 2.
Sia P la matrice le cui colonne sono gli autovettori trovati. Allora:
e quindi 
In definitiva D = P − 1AP è la matrice diagonale i cui elementi sono i rispettivi autovalori:

Diagonalizzazione di matrici reali simmetriche
Una generica matrice reale A può non essere diagonalizzabile. Se però A è una matrice simmetrica reale, valgono i seguenti teoremi.
Data una matrice A simmetrica reale, ogni radice λ del suo polinomio caratteristico è reale.
Data una matrice A simmetrica reale, detti u e v degli autovettori non di A, relativi ad autovalori distinti λ1 e λ2, allora u e v soon ortogonali.
A partire da questi due teoremi si arriva al teorema che segue.
Data una matrice A simmetrica reale, esiste una matrice ortogonale P, tale che D = P − 1AP è diagonale.
La matrice ortogonale P si ottiene normalizzando una base di autovettori ortogonali di A. A viene così detta ortogonalmente diagonalizzabile.
L'algoritmo
Esponiamo ora l'algoritmo di diagonalizzazione ortonormale:
Sia A una matrice reale simmetrica.
Passo 1. Trovare il polinomio caratteristico di A.
Passo 2. Trovare gli autovalori di A come radici del polinomio caratteristico.
Passo 3. Per ogni autovalore λ di A, trovare una base ortogonale del sottospazio relativo.
Passo 4. Normalizzare tutti gli autovalori ottenuti al Passo 3, formando quindi una base ortonormale di
.
Passo 5. Sia P la matrice le cui colonne sono gli autovettori nor,alizzati del Passo 4.
Esempio
Data la matrice simmetrica
,troviamo una matrice ortogonale P tale che P − 1AP è diagonale.
Procediamo come segue.
Troviamo il polinomio caratteristico di A:

Quindi λ1 = 6 e λ2 = 1 sono gli autovalori di A.
Sottraendo λ1 = 6 dalla diagonale di A si ha la matrice
ed il sistema omogeneo
oppure
Una soluzione diversa da zero è u1 = (1, − 2).
Sottraendo λ2 = 1 dalla diagonale di A si ha la matrice
ed il sistema omogeneo
La seconda equazione è eliminata in quanto multiplo della prima. Una soluzione diversa da zero è u2 = (2,1).
u1 e u2 sono ortogonali. Li normalizziamo e otteniamo i vettori ortonormali:
e 
Sia infine la matrice P di rispettive colonne u1 e u2, allora:
e 
Gli elementi diagonali di P − 1AP sono gli autovalori corrispondenti alle colonne di P.
Bibliografia
Algebra lineare - Lipschutz, Lipson.

Elettrotecnica e non solo (admin)
Un gatto tra gli elettroni (IsidoroKZ)
Esperienza e simulazioni (g.schgor)
Moleskine di un idraulico (RenzoDF)
Il Blog di ElectroYou (webmaster)
Idee microcontrollate (TardoFreak)
PICcoli grandi PICMicro (Paolino)
Il blog elettrico di carloc (carloc)
DirtEYblooog (dirtydeeds)
Di tutto... un po' (jordan20)
AK47 (lillo)
Esperienze elettroniche (marco438)
Telecomunicazioni musicali (clavicordo)
Automazione ed Elettronica (gustavo)
Direttive per la sicurezza (ErnestoCappelletti)
EYnfo dall'Alaska (mir)
Apriamo il quadro! (attilio)
H7-25 (asdf)
Passione Elettrica (massimob)
Elettroni a spasso (guidob)
Bloguerra (guerra)