Raggio spettrale

In matematica, il raggio spettrale di una matrice o di un operatore lineare limitato è l'estremo superiore della norma del modulo degli elementi del suo spettro. Spesso è denotato con $\rho (\cdot )$.

In analisi numerica il raggio spettrale viene utilizzato per determinare se un metodo iterativo è convergente verso la soluzione di un problema. È dimostrato infatti che un metodo iterativo per la risoluzione di un sistema lineare (come il metodo di Jacobi o quello di Gauss-Seidel) converge alla soluzione del sistema se e solo se il raggio spettrale della matrice di iterazione è strettamente minore di 1.

Matrici

Siano ${\lambda }_1,\dots {\lambda }_n$ autovalori (reali o complessi) di una matrice $A\in {ℂ}^{n\times n}$. Allora il suo raggio spettrale $\rho (A)$ è definito come:

$\rho (A)\equiv \underset{i}{\max }(|{\lambda }_i|)$

Un limite superiore per il raggio spettrale è dato dal seguente lemma. Sia $A\in {ℂ}^{n\times n}$ una matrice complessa, $\rho (A)$ il suo raggio spettrale e $\Vert \cdot \Vert$ una norma matriciale consistente. Allora per ogni $k\in \mathbb{N}$ si ha:

$\rho (A)\le \Vert A^k{\Vert }^{1/k}$

Infatti, sia $(𝐯,\lambda )$ una coppia autovettore-autovalore relativi ad $A$. Per la proprietà sub-moltiplicativa della norma matriciale:

$|\lambda {|}^k\Vert 𝐯\Vert =\Vert {\lambda }^k𝐯\Vert =\Vert A^k𝐯\Vert \le \Vert A^k\Vert \cdot \Vert 𝐯\Vert$

e dato che $𝐯\ne 0$ per ogni $\lambda$ si ha:

$|\lambda {|}^k\le \Vert A^k\Vert$

e dunque:

$\rho (A)\le \Vert A^k{\Vert }^{1/k}$

come si voleva mostrare.

Il raggio spettrale è strettamente legato al comportamento della convergenza della successione delle potenze di una matrice. In pratica, vale il seguente teorema. Sia $A\in {ℂ}^{n\times n}$ una matrice complessa e $\rho (A)$ il suo raggio spettrale. Allora $\underset{k\to \mathrm{\infty }}{\lim }{A}^k=0$ se e solo se $\rho (A)<1$. Inoltre, se $\rho (A)>1$ allora $\Vert A^k\Vert$ non è limitato per valori di $k$ crescenti.

Per mostrare che $\underset{k\to \mathrm{\infty }}{\lim }{A}^k=0$ implica $\rho (A)<1$, sia $(𝐯,\lambda )$ una coppia autovettore-autovalore relativi ad $A$. Dato che:

$A^k𝐯={\lambda }^k𝐯$

si ha:

$\begin{array}{cc}0& =\left(\underset{k\to \mathrm{\infty }}{\lim }{A}^k\right)𝐯\\ & =\underset{k\to \mathrm{\infty }}{\lim }{A}^k𝐯\\ & =\underset{k\to \mathrm{\infty }}{\lim }{\lambda }^k𝐯\\ & =𝐯\underset{k\to \mathrm{\infty }}{\lim }{\lambda }^k\end{array}$

e dato che per ipotesi $𝐯\ne 0$ si verifica:

$\underset{k\to \mathrm{\infty }}{\lim }{\lambda }^k=0$

che implica $|\lambda |<1$. Poiché questo deve essere vero per ogni autovalore, succede che $\rho (A)<1$.

Per mostrare che $\rho (A)<1$ implica $\underset{k\to \mathrm{\infty }}{\lim }{A}^k=0$, dal teorema di Jordan segue che per ogni matrice a valori nel campo complesso $A\in {ℂ}^{n\times n}$ esistono una matrice non singolare $V\in {ℂ}^{n\times n}$ e una matrice diagonale a blocchi $J\in {ℂ}^{n\times n}$ tali che:

$A=VJ{V}^{-1}$

con:

$J=\left[\begin{array}{ccccc}{J}_{m_1}({\lambda }_1)& 0& 0& \cdots & 0\\ 0& J_{m_2}({\lambda }_2)& 0& \cdots & 0\\ ⋮& \cdots & \ddots & \cdots & ⋮\\ 0& \cdots & 0& J_{m_{s-1}}({\lambda }_{s-1})& 0\\ 0& \cdots & \cdots & 0& J_{m_s}({\lambda }_s)\end{array}\right]$

dove:

$J_{m_i}({\lambda }_i)=\left[\begin{array}{ccccc}{\lambda }_i& 1& 0& \cdots & 0\\ 0& {\lambda }_i& 1& \cdots & 0\\ ⋮& ⋮& \ddots & \ddots & ⋮\\ 0& 0& \cdots & {\lambda }_i& 1\\ 0& 0& \cdots & 0& {\lambda }_i\end{array}\right]\in {ℂ}^{m_i,m_i}1\le i\le s$

Si vede facilmente che:

$A^k=V{J}^k{V}^{-1}$

e dato che $J$ è diagonale a blocchi:

$J^k=\left[\begin{array}{ccccc}{J}_{m_1}^k({\lambda }_1)& 0& 0& \cdots & 0\\ 0& J_{m_2}^k({\lambda }_2)& 0& \cdots & 0\\ ⋮& \cdots & \ddots & \cdots & ⋮\\ 0& \cdots & 0& J_{m_{s-1}}^k({\lambda }_{s-1})& 0\\ 0& \cdots & \cdots & 0& J_{m_s}^k({\lambda }_s)\end{array}\right]$

Un noto risultato riguardante la k-esima potenza di un blocco di Jordan $m_i\times m_i$ stabilisce che per $k\ge m_i-1$ si ha:

$J_{m_i}^k({\lambda }_i)=\left[\begin{array}{ccccc}{\lambda }_i^k& (\genfrac{}{}{0ex}{}{k}{1}){\lambda }_i^{k-1}& (\genfrac{}{}{0ex}{}{k}{2}){\lambda }_i^{k-2}& \cdots & (\genfrac{}{}{0ex}{}{k}{m_i-1}){\lambda }_i^{k-m_i+1}\\ 0& {\lambda }_i^k& (\genfrac{}{}{0ex}{}{k}{1}){\lambda }_i^{k-1}& \cdots & (\genfrac{}{}{0ex}{}{k}{m_i-2}){\lambda }_i^{k-m_i+2}\\ ⋮& ⋮& \ddots & \ddots & ⋮\\ 0& 0& \cdots & {\lambda }_i^k& (\genfrac{}{}{0ex}{}{k}{1}){\lambda }_i^{k-1}\\ 0& 0& \cdots & 0& {\lambda }_i^k\end{array}\right]$

In questo modo, se $\rho (A)<1$ allora $|{\lambda }_i|<1$ per ogni $i$, sicché:

$\underset{k\to \mathrm{\infty }}{\lim }{J}_{m_i}^k=0\mathrm{\forall }i$

e questo implica:

$\underset{k\to \mathrm{\infty }}{\lim }{J}^k=0$

Quindi:

$\underset{k\to \mathrm{\infty }}{\lim }{A}^k=\underset{k\to \mathrm{\infty }}{\lim }V{J}^k{V}^{-1}=V(\underset{k\to \mathrm{\infty }}{\lim }{J}^k)V^{-1}=0$

D'altra parte, se $\rho (A)>1$ allora vi è almeno un elemento in $J$ che non rimane limitato per $k$ crescente, concludendo la dimostrazione.

Formula di Gel'fand

La formula di Gel'fand (1941) stabilisce che per ogni norma matriciale $\Vert \cdot \Vert$ si ha:

$\rho (A)=\underset{k\to \mathrm{\infty }}{\lim }\Vert A^k{\Vert }^{1/k}$

In altri termini, mostra come il raggio spettrale di $A$ fornisca l'entità della crescita asintotica della norma di $A^k$, cioè:

$\Vert A^k\Vert \sim \rho (A{)}^k$

per $k\to \mathrm{\infty }$.

Per la dimostrazione, si consideri la matrice:

$\tilde{A}=(\rho (A)+ϵ{)}^{-1}Aϵ>0$

Allora:

$\rho (\tilde{A})=\frac{\rho (A)}{\rho (A)+ϵ}<1$

e per il teorema precedente:

$\underset{k\to \mathrm{\infty }}{\lim }{\tilde{A}}^k=0$

Per la definizione di limite di una successione, esiste un numero naturale $N_1\in \mathbb{N}$ tale per cui:

$\Vert {\tilde{A}}^k\Vert <1\mathrm{\forall }k\ge N_1$

che implica:

$\Vert A^k\Vert <(\rho (A)+ϵ{)}^k\mathrm{\forall }k\ge N_1$

o equivalentemente:

$\Vert A^k{\Vert }^{1/k}<(\rho (A)+ϵ)\mathrm{\forall }k\ge N_1$

Considerando ora la matrice:

$\check{A}=(\rho (A)-ϵ{)}^{-1}A$

in modo analogo si ha:

$\rho (\check{A})=\frac{\rho (A)}{\rho (A)-ϵ}>1$

e per il teorema precedente $\Vert {\check{A}}^k\Vert$ non è limitata. Esiste quindi $N_2\in \mathbb{N}$ tale per cui:

$\Vert {\check{A}}^k\Vert >1\mathrm{\forall }k\ge N_2$

che implica:

$\Vert A^k\Vert >(\rho (A)-ϵ{)}^k\mathrm{\forall }k\ge N_2$

o:

$\Vert A^k{\Vert }^{1/k}>(\rho (A)-ϵ)\mathrm{\forall }k\ge N_2$

Considerando:

$N:=\max (N_1,N_2)$

allora per ogni $ϵ>0$ esiste $N\in \mathbb{N}$ tale che per ogni $k\ge N$:

$\rho (A)-ϵ<\Vert A^k{\Vert }^{1/k}<\rho (A)+ϵ$

dunque:

$\underset{k\to \mathrm{\infty }}{\lim }\Vert A^k{\Vert }^{1/k}=\rho (A)$

come si voleva mostrare.

La formula di Gel'fand conduce direttamente ad un limite per il raggio spettrale del prodotto di infinite matrici. Nello specifico, assumendo che esse commutano reciprocamente:

$\rho (A_1{A}_2\dots A_n)\le \rho (A_1)\rho (A_2)\dots \rho (A_n)$

Inoltre, nel caso la norma matriciale sia consistente, grazie al lemma enunciato in precedenza si può rimpiazzare, nella definizione del limite, il limite inferiore sinistro con il raggio spettrale stesso. Quindi per ogni $ϵ>0$ esiste $N\in \mathbb{N}$ tale per cui:

$\rho (A)\le \Vert A^k{\Vert }^{1/k}<\rho (A)+ϵ\mathrm{\forall }k\ge N$

e dunque:

$\underset{k\to \mathrm{\infty }}{\lim }\Vert A^k{\Vert }^{1/k}=\rho (A{)}^{+}$

Operatori lineari limitati

Per un operatore lineare limitato $A$ e una norma operatoriale $\Vert \cdot \Vert$, il raggio spettrale $\rho (A)$ di $A$ è dato dalla formula di Gel'fand.

Esempio

Si consideri la matrice:

$A=\left[\begin{array}{ccc}9& -1& 2\\ -2& 8& 4\\ 1& 1& 8\end{array}\right]$

i cui autovalori sono 5, 10, 10. Per definizione, il suo raggio spettrale è $\rho (A)=10$. Nella tabella che segue sono riportati i valori di $\Vert A^k{\Vert }^{1/k}$ per le quattro norme più utilizzate, ordinati per $k$ crescente. Si nota che a causa della particolare forma della matrice $\Vert .{\Vert }_1=\Vert .{\Vert }_{\mathrm{\infty }}$.

k	$\Vert .{\Vert }_1=\Vert .{\Vert }_{\mathrm{\infty }}$	$\Vert .{\Vert }_F$	$\Vert .{\Vert }_2$
1	14	15.362291496	10.681145748
2	12.649110641	12.328294348	10.595665162
3	11.934831919	11.532450664	10.500980846
4	11.501633169	11.151002986	10.418165779
5	11.216043151	10.921242235	10.351918183
$⋮$	$⋮$	$⋮$	$⋮$
10	10.604944422	10.455910430	10.183690042
11	10.548677680	10.413702213	10.166990229
12	10.501921835	10.378620930	10.153031596
$⋮$	$⋮$	$⋮$	$⋮$
20	10.298254399	10.225504447	10.091577411
30	10.197860892	10.149776921	10.060958900
40	10.148031640	10.112123681	10.045684426
50	10.118251035	10.089598820	10.036530875
$⋮$	$⋮$	$⋮$	$⋮$
100	10.058951752	10.044699508	10.018248786
200	10.029432562	10.022324834	10.009120234
300	10.019612095	10.014877690	10.006079232
400	10.014705469	10.011156194	10.004559078
$⋮$	$⋮$	$⋮$	$⋮$
1000	10.005879594	10.004460985	10.001823382
2000	10.002939365	10.002230244	10.000911649
3000	10.001959481	10.001486774	10.000607757
$⋮$	$⋮$	$⋮$	$⋮$
10000	10.000587804	10.000446009	10.000182323
20000	10.000293898	10.000223002	10.000091161
30000	10.000195931	10.000148667	10.000060774
$⋮$	$⋮$	$⋮$	$⋮$
100000	10.000058779	10.000044600	10.000018232

Bibliografia

• Template:En Gradshteyn, I. S. and Ryzhik, I. M. Tables of Integrals, Series, and Products, 6th ed. San Diego, CA: Academic Press, pp. 1115-1116, 2000.

Voci correlate

• Autovettore e autovalore
• Norma (matematica)
• Operatore limitato
• Operatore lineare continuo
• Spettro (matematica)
• Spettro essenziale
• Teorema spettrale
• Teoria spettrale
• Trasformazione lineare

Collegamenti esterni

• Template:Collegamenti esterni