Funzione di matrice

In algebra lineare si può estendere il concetto di funzione alle matrici quadrate di qualsiasi ordine n attraverso l'associazione di una serie di Maclaurin ad ogni funzione, riducendola a una somma infinita di potenze di matrici:

Da cui risulta già chiaro che una funzione di matrice quadrata è una matrice dello stesso ordine i cui elementi sono costituiti da una combinazione lineare della funzione degli elementi della matrice di partenza, mentre in generale non risultano semplicemente le funzioni dell'elemento corrispondente della matrice di partenza. Le funzioni di matrice sono impiegate in particolare per risolvere i sistemi differenziali, di cui i più semplici sono i sistemi differenziali del prim'ordine, nella cui soluzione compare in particolare la matrice esponenziale e la matrice potenza.

Limitazione della serieModifica

Grazie al teorema di Hamilton-Cayley nella forma:

 

si può ridurre la procedura dal calcolo delle potenze di matrice da quello infinito fornito dalla definizione a quello di   potenze (l'identità e la matrice stessa banalmente non si calcolano), pur complicandone i coefficienti: moltiplicando più volte a destra e a sinistra del segno di uguale per la matrice  , è facile verificare che ogni potenza   può essere espressa come combinazione lineare delle sole prime   matrici.

 

Individuazione dei coefficienti della serie finitaModifica

Si può osservare che ciascun autovalore   della matrice di partenza   annulla, per definizione, il polinomio caratteristico, quindi analogamente a quanto avviene per la matrice:

 

Quindi scrivendo questa relazione per ciascun autovalore si ottiene un sistema lineare la cui matrice è quadrata di tipo Vandermonde di   righe e   colonne, che però non è invertibile quando esistono   autovalori distinti perché alcuni hanno molteplicità   con  , in quanto le righe corrispondenti risultano ripetute:

 

Si ricorre quindi al metodo dell'interpolazione polinomiale per vincolare sufficientemente il sistema, ricorrendo alle derivate successive fino alla  -esima:

 

Questa è una matrice invertibile di   righe e   colonne che definisce univocamente i coefficienti e ne permette il calcolo. Per maggiore semplicità di comprensione se ne dà la forma espansa:

 

ProceduraModifica

In base alle considerazioni fatte, il calcolo di una funzione di matrice si compone dei seguenti sei passaggi elementari:

  • individuazione della serie di MacLaurin associata alla funzione;
  • calcolo degli autovalori della matrice originaria;
  • calcolo della funzione di questi autovalori, e nel caso di autovalori con molteplicità nj dei valori assunti anche dalle derivate fino alla  -esima compresa;
  • calcolo dei coefficienti delle potenze come soluzione del sistema lineare di Vandermonde non omogeneo di cui sopra;
  • calcolo delle potenze della matrice originaria come prodotto o molto più velocemente sfruttando il teorema di Hamilton-Cayley;
  • calcolo della combinazione lineare di queste potenze con i coefficienti già calcolati.

Esempio applicativoModifica

Si voglia calcolare il seno di matrice:

 

Gli autovalori risultano essere   con molteplicità   e   con molteplicità   perciò i coefficienti sono:

 

perciò risulta che:

 

Potenza di matriceModifica

Anche la potenza di matrice necessaria per il calcolo di qualsiasi altra funzione diventa ottenibile molto più velocemente in base alle considerazioni precedenti, tanto più quanto minore è l'ordine della matrice   rispetto all'esponente  , e in generale risulta conveniente rispetto allo svolgimento di tutte le semplici moltiplicazioni necessarie quando  .

In realtà per esponente intero risulta ancora di gran lunga più veloce sfruttare direttamente il teorema di Hamilton-Cayley; tuttavia questo sistema permette di generalizzare la definizione di potenza fino ad ammettere un esponente complesso:

 

dove in questo caso i coefficienti risultano:

 

Per esempio si voglia calcolare:

 

poiché l'autovalore   ha molteplicità  , la matrice di Vandermonde risulta coincidere con la matrice originaria:

 

e quindi:

 

Per cui per esempio possiamo ammettere senza difficoltà che:

 

Voci correlateModifica


Controllo di autoritàThesaurus BNCF 54100
  Portale Matematica: accedi alle voci di Wikipedia che trattano di matematica