Siano
λ
1
,
…
λ
n
{\displaystyle \lambda _{1},\dots \lambda _{n}}
autovalori (reali o complessi) di una matrice
A
∈
C
n
×
n
{\displaystyle A\in \mathbb {C} ^{n\times n}}
. Allora il suo raggio spettrale
ρ
(
A
)
{\displaystyle \rho (A)}
è definito come:
ρ
(
A
)
≡
max
i
(
|
λ
i
|
)
{\displaystyle \rho (A)\equiv \max _{i}(|\lambda _{i}|)}
Un limite superiore per il raggio spettrale è dato dal seguente lemma. Sia
A
∈
C
n
×
n
{\displaystyle A\in \mathbb {C} ^{n\times n}}
una matrice complessa,
ρ
(
A
)
{\displaystyle \rho (A)}
il suo raggio spettrale e
‖
⋅
‖
{\displaystyle \|\cdot \|}
una norma matriciale consistente . Allora per ogni
k
∈
N
{\displaystyle k\in \mathbb {N} }
si ha:
ρ
(
A
)
≤
‖
A
k
‖
1
/
k
{\displaystyle \rho (A)\leq \|A^{k}\|^{1/k}}
Infatti, sia
(
v
,
λ
)
{\displaystyle (\mathbf {v} ,\lambda )}
una coppia autovettore-autovalore relativi ad
A
{\displaystyle A}
. Per la proprietà sub-moltiplicativa della norma matriciale:
|
λ
|
k
‖
v
‖
=
‖
λ
k
v
‖
=
‖
A
k
v
‖
≤
‖
A
k
‖
⋅
‖
v
‖
{\displaystyle |\lambda |^{k}\|\mathbf {v} \|=\|\lambda ^{k}\mathbf {v} \|=\|A^{k}\mathbf {v} \|\leq \|A^{k}\|\cdot \|\mathbf {v} \|}
e dato che
v
≠
0
{\displaystyle \mathbf {v} \neq 0}
per ogni
λ
{\displaystyle \lambda }
si ha:
|
λ
|
k
≤
‖
A
k
‖
{\displaystyle |\lambda |^{k}\leq \|A^{k}\|}
e dunque:
ρ
(
A
)
≤
‖
A
k
‖
1
/
k
{\displaystyle \rho (A)\leq \|A^{k}\|^{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
∈
C
n
×
n
{\displaystyle A\in \mathbb {C} ^{n\times n}}
una matrice complessa e
ρ
(
A
)
{\displaystyle \rho (A)}
il suo raggio spettrale. Allora
lim
k
→
∞
A
k
=
0
{\displaystyle \lim _{k\to \infty }A^{k}=0}
se e solo se
ρ
(
A
)
<
1
{\displaystyle \rho (A)<1}
. Inoltre, se
ρ
(
A
)
>
1
{\displaystyle \rho (A)>1}
allora
‖
A
k
‖
{\displaystyle \|A^{k}\|}
non è limitato per valori di
k
{\displaystyle k}
crescenti.
Per mostrare che
lim
k
→
∞
A
k
=
0
{\displaystyle \lim _{k\to \infty }A^{k}=0}
implica
ρ
(
A
)
<
1
{\displaystyle \rho (A)<1}
, sia
(
v
,
λ
)
{\displaystyle (\mathbf {v} ,\lambda )}
una coppia autovettore-autovalore relativi ad
A
{\displaystyle A}
. Dato che:
A
k
v
=
λ
k
v
{\displaystyle A^{k}\mathbf {v} =\lambda ^{k}\mathbf {v} }
si ha:
0
=
(
lim
k
→
∞
A
k
)
v
=
lim
k
→
∞
A
k
v
=
lim
k
→
∞
λ
k
v
=
v
lim
k
→
∞
λ
k
{\displaystyle {\begin{aligned}0&=\left(\lim _{k\to \infty }A^{k}\right)\mathbf {v} \\&=\lim _{k\to \infty }A^{k}\mathbf {v} \\&=\lim _{k\to \infty }\lambda ^{k}\mathbf {v} \\&=\mathbf {v} \lim _{k\to \infty }\lambda ^{k}\end{aligned}}}
e dato che per ipotesi
v
≠
0
{\displaystyle \mathbf {v} \neq 0}
si verifica:
lim
k
→
∞
λ
k
=
0
{\displaystyle \lim _{k\to \infty }\lambda ^{k}=0}
che implica
|
λ
|
<
1
{\displaystyle |\lambda |<1}
. Poiché questo deve essere vero per ogni autovalore, succede che
ρ
(
A
)
<
1
{\displaystyle \rho (A)<1}
.
Per mostrare che
ρ
(
A
)
<
1
{\displaystyle \rho (A)<1}
implica
lim
k
→
∞
A
k
=
0
{\displaystyle \lim _{k\to \infty }A^{k}=0}
, dal teorema di Jordan segue che per ogni matrice a valori nel campo complesso
A
∈
C
n
×
n
{\displaystyle A\in \mathbb {C} ^{n\times n}}
esistono una matrice non singolare
V
∈
C
n
×
n
{\displaystyle V\in \mathbb {C} ^{n\times n}}
e una matrice diagonale a blocchi
J
∈
C
n
×
n
{\displaystyle J\in \mathbb {C} ^{n\times n}}
tali che:
A
=
V
J
V
−
1
{\displaystyle A=VJV^{-1}}
con:
J
=
[
J
m
1
(
λ
1
)
0
0
⋯
0
0
J
m
2
(
λ
2
)
0
⋯
0
⋮
⋯
⋱
⋯
⋮
0
⋯
0
J
m
s
−
1
(
λ
s
−
1
)
0
0
⋯
⋯
0
J
m
s
(
λ
s
)
]
{\displaystyle J={\begin{bmatrix}J_{m_{1}}(\lambda _{1})&0&0&\cdots &0\\0&J_{m_{2}}(\lambda _{2})&0&\cdots &0\\\vdots &\cdots &\ddots &\cdots &\vdots \\0&\cdots &0&J_{m_{s-1}}(\lambda _{s-1})&0\\0&\cdots &\cdots &0&J_{m_{s}}(\lambda _{s})\end{bmatrix}}}
dove:
J
m
i
(
λ
i
)
=
[
λ
i
1
0
⋯
0
0
λ
i
1
⋯
0
⋮
⋮
⋱
⋱
⋮
0
0
⋯
λ
i
1
0
0
⋯
0
λ
i
]
∈
C
m
i
,
m
i
1
≤
i
≤
s
{\displaystyle J_{m_{i}}(\lambda _{i})={\begin{bmatrix}\lambda _{i}&1&0&\cdots &0\\0&\lambda _{i}&1&\cdots &0\\\vdots &\vdots &\ddots &\ddots &\vdots \\0&0&\cdots &\lambda _{i}&1\\0&0&\cdots &0&\lambda _{i}\end{bmatrix}}\in \mathbb {C} ^{m_{i},m_{i}}\qquad 1\leq i\leq s}
Si vede facilmente che:
A
k
=
V
J
k
V
−
1
{\displaystyle A^{k}=VJ^{k}V^{-1}}
e dato che
J
{\displaystyle J}
è diagonale a blocchi:
J
k
=
[
J
m
1
k
(
λ
1
)
0
0
⋯
0
0
J
m
2
k
(
λ
2
)
0
⋯
0
⋮
⋯
⋱
⋯
⋮
0
⋯
0
J
m
s
−
1
k
(
λ
s
−
1
)
0
0
⋯
⋯
0
J
m
s
k
(
λ
s
)
]
{\displaystyle J^{k}={\begin{bmatrix}J_{m_{1}}^{k}(\lambda _{1})&0&0&\cdots &0\\0&J_{m_{2}}^{k}(\lambda _{2})&0&\cdots &0\\\vdots &\cdots &\ddots &\cdots &\vdots \\0&\cdots &0&J_{m_{s-1}}^{k}(\lambda _{s-1})&0\\0&\cdots &\cdots &0&J_{m_{s}}^{k}(\lambda _{s})\end{bmatrix}}}
Un noto risultato riguardante la k -esima potenza di un blocco di Jordan
m
i
×
m
i
{\displaystyle m_{i}\times m_{i}}
stabilisce che per
k
≥
m
i
−
1
{\displaystyle k\geq m_{i}-1}
si ha:
J
m
i
k
(
λ
i
)
=
[
λ
i
k
(
k
1
)
λ
i
k
−
1
(
k
2
)
λ
i
k
−
2
⋯
(
k
m
i
−
1
)
λ
i
k
−
m
i
+
1
0
λ
i
k
(
k
1
)
λ
i
k
−
1
⋯
(
k
m
i
−
2
)
λ
i
k
−
m
i
+
2
⋮
⋮
⋱
⋱
⋮
0
0
⋯
λ
i
k
(
k
1
)
λ
i
k
−
1
0
0
⋯
0
λ
i
k
]
{\displaystyle J_{m_{i}}^{k}(\lambda _{i})={\begin{bmatrix}\lambda _{i}^{k}&{k \choose 1}\lambda _{i}^{k-1}&{k \choose 2}\lambda _{i}^{k-2}&\cdots &{k \choose m_{i}-1}\lambda _{i}^{k-m_{i}+1}\\0&\lambda _{i}^{k}&{k \choose 1}\lambda _{i}^{k-1}&\cdots &{k \choose m_{i}-2}\lambda _{i}^{k-m_{i}+2}\\\vdots &\vdots &\ddots &\ddots &\vdots \\0&0&\cdots &\lambda _{i}^{k}&{k \choose 1}\lambda _{i}^{k-1}\\0&0&\cdots &0&\lambda _{i}^{k}\end{bmatrix}}}
In questo modo, se
ρ
(
A
)
<
1
{\displaystyle \rho (A)<1}
allora
|
λ
i
|
<
1
{\displaystyle |\lambda _{i}|<1}
per ogni
i
{\displaystyle i}
, sicché:
lim
k
→
∞
J
m
i
k
=
0
∀
i
{\displaystyle \lim _{k\to \infty }J_{m_{i}}^{k}=0\ \forall i}
e questo implica:
lim
k
→
∞
J
k
=
0
{\displaystyle \lim _{k\to \infty }J^{k}=0}
Quindi:
lim
k
→
∞
A
k
=
lim
k
→
∞
V
J
k
V
−
1
=
V
(
lim
k
→
∞
J
k
)
V
−
1
=
0
{\displaystyle \lim _{k\to \infty }A^{k}=\lim _{k\to \infty }VJ^{k}V^{-1}=V(\lim _{k\to \infty }J^{k})V^{-1}=0}
D'altra parte, se
ρ
(
A
)
>
1
{\displaystyle \rho (A)>1}
allora vi è almeno un elemento in
J
{\displaystyle J}
che non rimane limitato per
k
{\displaystyle k}
crescente, concludendo la dimostrazione.
Formula di Gel'fand
modifica
La formula di Gel'fand (1941) stabilisce che per ogni norma matriciale
‖
⋅
‖
{\displaystyle \|\cdot \|}
si ha:
ρ
(
A
)
=
lim
k
→
∞
‖
A
k
‖
1
/
k
{\displaystyle \rho (A)=\lim _{k\to \infty }\|A^{k}\|^{1/k}}
In altri termini, mostra come il raggio spettrale di
A
{\displaystyle A}
fornisca l'entità della crescita asintotica della norma di
A
k
{\displaystyle A^{k}}
, cioè:
‖
A
k
‖
∼
ρ
(
A
)
k
{\displaystyle \|A^{k}\|\sim \rho (A)^{k}}
per
k
→
∞
{\displaystyle k\to \infty }
.
Per la dimostrazione, si consideri la matrice:
A
~
=
(
ρ
(
A
)
+
ϵ
)
−
1
A
ϵ
>
0
{\displaystyle {\tilde {A}}=(\rho (A)+\epsilon )^{-1}A\qquad \epsilon >0}
Allora:
ρ
(
A
~
)
=
ρ
(
A
)
ρ
(
A
)
+
ϵ
<
1
{\displaystyle \rho ({\tilde {A}})={\frac {\rho (A)}{\rho (A)+\epsilon }}<1}
e per il teorema precedente:
lim
k
→
∞
A
~
k
=
0
{\displaystyle \lim _{k\to \infty }{\tilde {A}}^{k}=0}
Per la definizione di limite di una successione , esiste un numero naturale
N
1
∈
N
{\displaystyle N_{1}\in \mathbb {N} }
tale per cui:
‖
A
~
k
‖
<
1
∀
k
≥
N
1
{\displaystyle \|{\tilde {A}}^{k}\|<1\qquad \forall k\geq N_{1}}
che implica:
‖
A
k
‖
<
(
ρ
(
A
)
+
ϵ
)
k
∀
k
≥
N
1
{\displaystyle \|A^{k}\|<(\rho (A)+\epsilon )^{k}\qquad \forall k\geq N_{1}}
o equivalentemente:
‖
A
k
‖
1
/
k
<
(
ρ
(
A
)
+
ϵ
)
∀
k
≥
N
1
{\displaystyle \|A^{k}\|^{1/k}<(\rho (A)+\epsilon )\qquad \forall k\geq N_{1}}
Considerando ora la matrice:
A
ˇ
=
(
ρ
(
A
)
−
ϵ
)
−
1
A
{\displaystyle {\check {A}}=(\rho (A)-\epsilon )^{-1}A}
in modo analogo si ha:
ρ
(
A
ˇ
)
=
ρ
(
A
)
ρ
(
A
)
−
ϵ
>
1
{\displaystyle \rho ({\check {A}})={\frac {\rho (A)}{\rho (A)-\epsilon }}>1}
e per il teorema precedente
‖
A
ˇ
k
‖
{\displaystyle \|{\check {A}}^{k}\|}
non è limitata. Esiste quindi
N
2
∈
N
{\displaystyle N_{2}\in \mathbb {N} }
tale per cui:
‖
A
ˇ
k
‖
>
1
∀
k
≥
N
2
{\displaystyle \|{\check {A}}^{k}\|>1\qquad \forall k\geq N_{2}}
che implica:
‖
A
k
‖
>
(
ρ
(
A
)
−
ϵ
)
k
∀
k
≥
N
2
{\displaystyle \|A^{k}\|>(\rho (A)-\epsilon )^{k}\qquad \forall k\geq N_{2}}
o:
‖
A
k
‖
1
/
k
>
(
ρ
(
A
)
−
ϵ
)
∀
k
≥
N
2
{\displaystyle \|A^{k}\|^{1/k}>(\rho (A)-\epsilon )\qquad \forall k\geq N_{2}}
Considerando:
N
:=
max
(
N
1
,
N
2
)
{\displaystyle N:=\max(N_{1},N_{2})}
allora per ogni
ϵ
>
0
{\displaystyle \epsilon >0}
esiste
N
∈
N
{\displaystyle N\in \mathbb {N} }
tale che per ogni
k
≥
N
{\displaystyle k\geq N}
:
ρ
(
A
)
−
ϵ
<
‖
A
k
‖
1
/
k
<
ρ
(
A
)
+
ϵ
{\displaystyle \rho (A)-\epsilon <\|A^{k}\|^{1/k}<\rho (A)+\epsilon }
dunque:
lim
k
→
∞
‖
A
k
‖
1
/
k
=
ρ
(
A
)
{\displaystyle \lim _{k\to \infty }\|A^{k}\|^{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:
ρ
(
A
1
A
2
…
A
n
)
≤
ρ
(
A
1
)
ρ
(
A
2
)
…
ρ
(
A
n
)
{\displaystyle \rho (A_{1}A_{2}\ldots A_{n})\leq \rho (A_{1})\rho (A_{2})\ldots \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
{\displaystyle \epsilon >0}
esiste
N
∈
N
{\displaystyle N\in \mathbb {N} }
tale per cui:
ρ
(
A
)
≤
‖
A
k
‖
1
/
k
<
ρ
(
A
)
+
ϵ
∀
k
≥
N
{\displaystyle \rho (A)\leq \|A^{k}\|^{1/k}<\rho (A)+\epsilon \qquad \forall k\geq N}
e dunque:
lim
k
→
∞
‖
A
k
‖
1
/
k
=
ρ
(
A
)
+
{\displaystyle \lim _{k\to \infty }\|A^{k}\|^{1/k}=\rho (A)^{+}}
Operatori lineari limitati
modifica
Per un operatore lineare limitato
A
{\displaystyle A}
e una norma operatoriale
‖
⋅
‖
{\displaystyle \|\cdot \|}
, il raggio spettrale
ρ
(
A
)
{\displaystyle \rho (A)}
di
A
{\displaystyle A}
è dato dalla formula di Gel'fand.
(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.
Collegamenti esterni
modifica