Definizione
Data una matrice
[math]A \in \mathbb{R}^{n imes n}[/math]
, la potenza [math]k[/math]
-esima di [math]A[/math]
([math]k \in \mathbb{N}[/math]
), è definita ricorsivamente così come segue
[math]A^k = \egin{cases} I & \quad \text{se } k = 0 \\ A^{k-1} \cdot A & \quad \text{se } k > 0 \ \end{cases}[/math]
Potenza di una matrice diagonalizzabile
Se
[math]A \in \mathbb{R}^{n imes n}[/math]
è una matrice diagonalizzabile, il calcolo della potenza [math]k[/math]
-esima [math]A^k[/math]
può essere effettuato senza bisogno di ricorrere alla definizione. Detta [math]T[/math]
la matrice del cambio di coordinate, risulta
[math]A = T \Lambda T^{-1}[/math]
dove
[math]\Lambda[/math]
è una matrice in forma diagonale se tutti gli autovalori di [math]A[/math]
sono reali, mentre è in forma daigonale reali a blocchi se [math]A[/math]
ammette autovalori complessi. In entrambi i casi la potenza [math]k[/math]
-esima di [math]A[/math]
vale
[math]A^k = T \Lambda^k T^{-1}[/math]
Nel caso di matrici diagonalizzabili dunque è sufficiente calcolare, oltre alla matrice del calcolo di coordinate, la potenza della matrice diagonalizzata, e tale calcolo è relativamente semplice.
Nel calcolo di
[math]\Lambda^k[/math]
si distinguono i casi in cui [math]\Lambda[/math]
è una matrice in forma diagonale da quello in cui [math]\Lambda[/math]
è in forma diagonale reale a blocchi.
[math]A[/math] ha tutti gli autovalori reali
Se la matrice
[math]A[/math]
ha tutti gli autovalori reali, e questi sono pari a [math]\lambda_{1}, \lambda_{2}, \ldots, \lambda_n[/math]
, ed inoltre se [math]T[/math]
è la matrice del cambio di coordinate, allora
[math]A = T \Lambda T^{-1}[/math]
dove
[math]\Lambda = ((\lambda_{1}, \quad 0, \ldots, 0),(0, \lambda_{2}, \quad \ldots, 0),(vdots, \quad vdots, \quad ddots, \quad vdots),(0, \quad 0, \quad \ldots, \quad \lambda_n))[/math]
La potenza
[math]k[/math]
-esima di [math]\Lambda[/math]
è banalmente la matrice ottenuta calcolando la potenza [math]k[/math]
-esima degli elementi sulla diagonale principale
[math]\Lambda^k = ((\lambda_{1}^k, \quad 0, \ldots, 0),(0, \lambda_{2}^k, \quad \ldots, 0),(vdots, \quad vdots, \quad ddots, \quad vdots),(0, \quad 0, \quad \ldots, \quad \lambda_n^k))[/math]
Pertanto, in questo caso, la potenza
[math]k[/math]
-esima di [math]A[/math]
vale
[math]A^k = T \Lambda^k T^{-1}[/math]
[math]A[/math] ha almeno un autovalore complesso
Supponiamo ora che
[math]A[/math]
abbia autovalori complessi, e che i suoi autovalori siano [math]\lambda_{1}, \lambda_{2}, \ldots, \lambda_{r}, \lambda_{r+1}, \bar{\lambda_{r+1}}, \lambda_{r+2}, \bar{\lambda_{r+2}}, \ldots, \lambda_{c}, \bar{\lambda_{c}}[/math]
, dove i primi [math]r[/math]
sono reali, i restanti sono complessi. Scriviamo ogni autovalore complesso con parte immaginaria positiva come
[math]\lambda_i = \sigma_i + i \omega_i[/math]
Se
[math]T[/math]
è la matrice del cambio di coordinate, allora [math]A = T \Lambda T^{-1}[/math]
, dove [math]\Lambda[/math]
è una matrice in forma diagonale reale a blocchi, della forma
[math]\Lambda = [(\lambda_1, \quad 0, \quad \ldots, \quad \ldots, \quad \ldots, \quad \ldots, \quad \ldots, \quad \ldots, \quad 0),(0, \quad \lambda_2, \quad 0, \quad \ldots, \quad \ldots, \quad \ldots, \quad \ldots, \quad \ldots, \quad 0),(vdots, \quad 0, \quad ddots, \quad 0, \quad \ldots, \quad \ldots, \quad \ldots, \quad \ldots, \quad vdots),(vdots, \quad vdots, \quad 0, \quad \lambda_{r}, 0, \quad \ldots, \quad \ldots, \quad \ldots, \quad vdots),(vdots, \quad vdots, \quad vdots, \quad 0, \quad \sigma_{r+1}, \quad \omega_{r+1}, \quad 0, \quad \ldots, \quad vdots),(vdots, \quad vdots, \quad vdots, \quad ddots, \quad -\omega_{r+1}, \quad \sigma_{r+1}, \quad 0, \quad \ldots, \quad vdots),(vdots, \quad vdots, \quad vdots, \quad vdots, \quad ddots, \quad 0, \quad ddots, \quad 0, \quad 0),(vdots, \quad vdots, \quad vdots, \quad vdots, \quad vdots, \quad ddots, \quad vdots, \quad \sigma_c, \quad \omega_c),(0, \quad \ldots, \quad \ldots, \quad \ldots, \quad \ldots, \quad \ldots, \quad 0, \quad -\omega_c, \quad \sigma_c)][/math]
Per ogni autovalore complesso
[math]\lambda_i = \sigma_i + i \omega_i[/math]
(a cui corrisponde un blocco [math]2 imes 2[/math]
sulla diagonale di [math]\Lambda[/math]
), indichiamo con [math]\\rho_i[/math]
il suo modulo e con [math]\theta_i[/math]
la sua fase, ossia
[math]\\rho_i = \sqrt{\sigma_i^2 + \omega_i^2} qquad \theta_i = \text{arctg}{\frac{\omega_i}{\sigma_i}}[/math]
Fatto questo si può scrivere la potenza
[math]k[/math]
-esima di [math]\Lambda[/math]
, che è pari a
[math]\Lambda^k = [(\lambda_1^k, \quad 0, \quad \ldots, \quad \ldots, \quad \ldots, \quad \ldots, \quad \ldots, \quad \ldots, \quad 0),(0, \quad \lambda_2^k, \quad 0, \quad \ldots, \quad \ldots, \quad \ldots, \quad \ldots, \quad \ldots, \quad 0),(vdots, \quad 0, \quad ddots, \quad 0, \quad \ldots, \quad \ldots, \quad \ldots, \quad \ldots, \quad vdots),(vdots, \quad vdots, \quad 0, \quad \lambda_{r}^k, \quad 0, \quad \ldots, \quad \ldots, \quad \ldots, \quad vdots),(vdots, \quad vdots, \quad vdots, \quad 0, \quad \\rho_{r+1}^k \\cos(k \theta_{r+1}),\quad \\rho_{r+1}^k \\sin(k \theta_{r+1}), \quad 0, \quad \ldots, \quad vdots),(vdots, \quad vdots, \quad vdots, \quad ddots, \quad -\\rho_{r+1}^k \\sin(k \theta_{r+1}), \quad \\rho_{r+1}^k \\cos(k \theta_{r+1}), \quad 0, \quad \ldots, \quad vdots),(vdots, \quad vdots, \quad vdots, \quad vdots, \quad ddots, \quad 0, \quad ddots, \quad 0, \quad 0),(vdots, \quad vdots, \quad vdots, \quad vdots, \quad vdots, \quad ddots, \quad vdots, \quad \\rho_{c}^k \\cos(k \theta_{c}), \quad \\rho_{c}^k \\sin(k \theta_{c})),(0, \quad \ldots, \quad \ldots, \quad \ldots, \quad \ldots, \quad \ldots, \quad 0, \quad -\\rho_{c}^k \\sin(k \theta_{c}), \quad \\rho_{c}^k \\cos(k \theta_{c}))][/math]
Una volta calcolata
[math]\Lambda^k[/math]
è possibile calcolare anche la potenza [math]k[/math]
-esima di [math]A[/math]
, dato che
[math]A^k = T \Lambda^k T^{-1}[/math]
