Matrice esponenziale — ingegnerismo.it

La matrice esponenziale estende l’esponenziale reale al caso di una matrice quadrata. Se $A\in\mathbb{R}^{n\times n}$ oppure $A\in\mathbb{C}^{n\times n}$ , si definisce tramite la serie di potenze

e^{At}=\sum_{k=0}^{\infty}\dfrac{A^k t^k}{k!}.

La serie converge per ogni matrice quadrata e per ogni valore reale o complesso di $t$ . La definizione non richiede che $A$ sia diagonalizzabile: anche matrici difettive, triangolari o non simmetriche hanno una matrice esponenziale ben definita.

Proprietà principali

Per ogni matrice quadrata $A$ :

e^{0}=I, \qquad (e^A)^{-1}=e^{-A}, \qquad \det(e^A)=e^{\operatorname{tr}A}.

Se $A$ e $B$ commutano, cioè $AB=BA$ , allora

e^{A+B}=e^Ae^B.

La condizione di commutazione è essenziale: in generale l’identità non vale. Questo è uno degli errori più frequenti quando si manipolano sistemi accoppiati, operatori non commutativi o modelli a coefficienti variabili.

Significato dinamico

La proprietà fondamentale è

\dfrac{d}{dt}e^{At}=Ae^{At}, \qquad e^{A0}=I.

Per il sistema lineare omogeneo $x'=Ax$ , la soluzione è

x(t)=e^{A(t-t_0)}x_0.

Per questo la matrice esponenziale è centrale in analisi dei sistemi di EDO, automazione e teoria dei sistemi lineari. La matrice $e^{At}$ è la matrice di transizione di stato: trasporta lo stato iniziale $x_0$ fino allo stato al tempo $t$ .

Nel caso non omogeneo

x'(t)=Ax(t)+b(t)

la formula di variazione delle costanti diventa

x(t)=e^{A(t-t_0)}x_0+ \int_{t_0}^{t}e^{A(t-\tau)}b(\tau)\,d\tau.

Calcolo e casi particolari

Se $A$ è diagonalizzabile, con $A=PDP^{-1}$ , allora

e^{At}=Pe^{Dt}P^{-1},

dove $e^{Dt}$ si ottiene esponenziando gli elementi diagonali di $D$ . Se $A$ non è diagonalizzabile, si ricorre alla forma di Jordan, alla decomposizione di Schur o a metodi numerici come scaling and squaring con approssimanti di Padé.

Per un blocco di forma canonica di Jordan

J_k(\lambda)=\lambda I+N,

dove $N$ è nilpotente, vale

e^{J_k(\lambda)t}=e^{\lambda t} \left(I+Nt+\dfrac{N^2t^2}{2!}+\cdots+\dfrac{N^{k-1}t^{k-1}}{(k-1)!}\right).

Il teorema di Cayley-Hamilton permette inoltre, per matrici di ordine $n$ , di esprimere $e^{At}$ come combinazione di

I,\ A,\ A^2,\ \dots,\ A^{n-1}.

Questa strada è utile nei casi piccoli e nei calcoli simbolici, mentre per applicazioni numeriche robuste si preferiscono metodi specializzati.

Applicazioni ingegneristiche

Controlli automatici: $e^{At}$ descrive l’evoluzione libera di un sistema LTI; la stabilità dipende dagli autovalori di $A$ .
Circuiti elettrici: transitori di reti RC, RL e RLC possono essere scritti in forma di stato e risolti con la matrice esponenziale.
Meccanica strutturale: sistemi smorzati a più gradi di libertà possono essere studiati tramite formulazione di stato e riduzione modale.
Robotica: rotazioni finite e moti rigidi si rappresentano con esponenziali di matrici antisimmetriche o twist.

Tavola di sintesi

Lettura	Formula chiave	Interpretazione operativa
Definizione per serie	$\displaystyle e^{At}=\sum_{k=0}^{\infty}\dfrac{A^k t^k}{k!}$	Definisce sempre la matrice esponenziale, anche quando $A$ non è diagonalizzabile. È la base teorica del concetto.
Calcolo spettrale o di Jordan	$\displaystyle A=PDP^{-1}\Rightarrow e^{At}=Pe^{Dt}P^{-1}$ ; per $\displaystyle J=\lambda I+N$ , $\displaystyle e^{Jt}=e^{\lambda t}e^{Nt}$	Trasforma il calcolo in un problema sugli autovalori e sui blocchi nilpotenti. È utile per calcoli simbolici, stabilità e interpretazione modale.
Propagatore di stato	$\displaystyle x(t)=e^{A(t-t_0)}x_0$	Descrive l’evoluzione libera di un sistema lineare. In automazione è la matrice che trasporta lo stato iniziale fino al tempo $t$ .

La stessa formula ha quindi tre ruoli complementari: definizione analitica, strumento di calcolo e operatore dinamico. Nelle applicazioni ingegneristiche il terzo punto è spesso il più importante, perché collega direttamente gli autovalori di $A$ alla risposta temporale del sistema.