18.01.2015 Views

Appunti di Calcolo Numerico - Esercizi e Dispense - Università degli ...

Appunti di Calcolo Numerico - Esercizi e Dispense - Università degli ...

Appunti di Calcolo Numerico - Esercizi e Dispense - Università degli ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

8. METODI ITERATIVI PER LA SOLUZIONE DI SISTEMI LINEARI<br />

Norma<br />

compatibile<br />

Norma<br />

naturale<br />

traccia <strong>di</strong> una<br />

matrice<br />

G ‖A‖ = 0 se e solo se A = 0<br />

G ‖αA‖ = |α|‖A‖<br />

G ‖A + B‖ ≤ ‖A‖ + ‖B‖<br />

G ‖AB‖ ≤ ‖A‖‖B‖<br />

Una proprietà importante che si richiede alle norme su matrici è che siano compatibili con norme vettoriali:<br />

la norma ‖A‖ <strong>di</strong> una matrice A si <strong>di</strong>ce compatibile con la norma ‖x‖ <strong>di</strong> un vettore x se vale la<br />

relazione<br />

‖Ax‖ ≤ ‖A‖‖x‖<br />

Alcune norme su matrici sono generate da norme su vettori: si parla allora <strong>di</strong> norma naturale o indotta<br />

dalla norma <strong>di</strong> vettori. In particolare, se ‖ · ‖ è una norma su vettori in R n , allora ‖A‖ = max ‖x‖=1 ‖Ax‖ è la<br />

norma naturale o indotta dalla norma ‖ · ‖ su vettori.<br />

Le norme <strong>di</strong> matrici indotte dalla norma 1 e dalla norma infinito su vettori sono:<br />

G Norma 1: ‖A‖ 1 = max j<br />

∑ n<br />

i=1 |a i j | (data dal massimo sulla somma delle colonne)<br />

G Norma infinito: ‖A‖ ∞ = max i<br />

∑ n<br />

j =1 |a i j | (data dal massimo sulla somma delle righe)<br />

La norma <strong>di</strong> matrice indotta dalla norma 2 è più complicata e vedremo in seguito come è definita.<br />

È facile vedere che le norme naturali sono norme compatibili con la norma <strong>di</strong> vettori da cui sono costruite.<br />

Una norma <strong>di</strong> matrici, che non è indotta, ma compatibile con la norma 2 è la cosiddetta norma euclidea<br />

(o <strong>di</strong> Frobenius). Tenendo presente che, data una matrice A, si chiama traccia della matrice o tr (A) la somma<br />

<strong>degli</strong> elementi della <strong>di</strong>agonale principale <strong>di</strong> A, la norma euclidea è data da<br />

G N (A) = √ tr (A T A) = √ √ ∑n<br />

tr (A A T ) =<br />

i=1<br />

|a i j | 2 .<br />

j =1<br />

8.5 Autovalori e autovettori<br />

Data una matrice quadrata A <strong>di</strong> or<strong>di</strong>ne n, se esiste un numero (reale o complesso) λ e un vettore x ≠ 0 tali<br />

che<br />

Ax = λx<br />

Autovalore e<br />

autovettore<br />

allora λ è un autovalore e x il corrispondente autovettore della matrice A.<br />

Scritta in maniera equivalente, la relazione definisce il sistema lineare<br />

(A − λI )x = 0<br />

Poichè x ≠ 0 e il termine noto del sistema è il vettore <strong>di</strong> tutti zeri, il determinante della matrice del sistema<br />

deve necessariamente essere uguale a zero, cioè det(A − λI ) = 0.<br />

Lo sviluppo del determinante porta a un polinomio <strong>di</strong> grado n nell’incognita λ:<br />

λ n − tr (A)λ n−1 + ... + (−1) n det(A) = 0<br />

Polinomio Questo polinomio si chiama polinomio caratteristico. Le sue n ra<strong>di</strong>ci, che chiamiamo λ 1 ,λ 2 ,...,λ n , sono gli<br />

caratteristico n autovalori della matrice A.<br />

Per le proprietà dei polinomi vale:<br />

n∑<br />

λ i = tr (A) = a 11 + a 22 + ... + a nn<br />

i=1<br />

e<br />

n∏<br />

λ i = det(A)<br />

i=1<br />

116<br />

Alcune proprietà sugli autovalori e autovettori sono le seguenti:<br />

G Se λ è autovalore della matrice A, allora λ k è autovalore della matrice potenza A k (cioè A · A ··· A k<br />

volte).<br />

G Se λ è autovalore della matrice A, e A è regolare, allora λ −1 è autovalore della matrice inversa A −1 .

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!