In matematica, e in particolare in probabilità, con matrice aleatoria si intende una matrice in cui una parte o tutti i suoi ingressi sono variabili aleatorie estratte in modo casuale da una distribuzione di probabilità. In altre parole, è una variabile aleatoria a valori matriciali. Talvolta è utilizzata anche l’espressione matrice casuale come sinonimo, tuttavia tale dizione è non preferibile in contesti matematici.
La teoria delle matrici aleatorie studia le proprietà delle matrici aleatorie, focalizzandosi sul calcolo e la stima di quantità, come traccia, determinante, norma e spettro, in particolare per il limite di matrici a dimensionalità elevata. La teoria si avvale di strumenti come la teoria di campo medio, i metodi diagrammatici e il metodo delle cavità. Getta le sue basi nella seconda metà del XX secolo, grazie ai fondamentali contributi di Wigner, Dyson e Lal Mehta. L'importanza dell'argomento è giustificata dal fatto che numerosi fenomeni fisici e non, come lo spettro dei nuclei di atomi pesanti, la conducibilità termica di un reticolo o l’emergere del caos quantistico, i mercati azionari e lo studio della forma delle proteine possono essere modellati matematicamente come problemi relativi a matrici grandi e aleatorie.
Enrico Fermi e collaboratori furono i primi che, tramite esperimenti, mostrarono che i nucleoni all’interno di un nucleo atomico non possono essere modellati come particelle che si muovono indipendentemente. Dunque questa osservazione indusse Niels Bohr alla formulazione del concetto di nucleo composto, secondo il quale le interazioni tra nucleoni, prima di decadere, danno luogo a stati intermedi. Fu questa teoria, nell’ambito della fisica nucleare, la prima ad utilizzare gli strumenti delle matrici aleatorie, al di fuori della letteratura matematica.
In assenza di una conoscenza dettagliata delle interazioni dirette tra nucleoni, Eugene Wigner propose di modellare l’Hamiltoniana nucleare come una matrice aleatoria. Per atomi con un numero elevato di nucleoni, la distribuzione degli autovalori dell’Hamiltoniana poteva essere utilizzata per stimare le sezioni d’urto dei processi di scattering, facendo ricorso alla distribuzione di Wishart.

Negli anni successivi, Freeman Dyson formalizzò ulteriormente il quadro teorico introducendo il concetto di ensemble di matrici aleatorie, in particolare classificando i cosiddetti ensemble classici (gaussiano ortogonale, unitario e simplettico), noti con l’acronimo inglese GOE, GUE e GSE. Tali ensemble descrivono differenti simmetrie fisiche associate all’Hamiltoniana, e si distinguono per il comportamento statistico degli autovalori.
Oltre alla fisica nucleare, la teoria delle matrici aleatorie ha trovato applicazioni in numerosi altri contesti, tra cui la fisica delle particelle, i sistemi disordinati (come la teoria della localizzazione di Anderson), la teoria dei numeri, la finanza quantitativa, e più recentemente, l’apprendimento automatico.
In particolare, sviluppi recenti nel campo del deep learning hanno mostrato che modelli di reti neurali profonde possono essere analizzati utilizzando strumenti della teoria delle matrici aleatorie. Le matrici dei pesi nelle reti neurali, soprattutto dopo l’addestramento, presentano spettro e proprietà statistiche che possono essere descritte in termini di distribuzioni di Wigner o Marčenko–Pastur. Tali osservazioni hanno portato alla nascita di una nuova disciplina, talvolta chiamata “Random Matrix Theory for Deep Learning”, che fornisce una cornice teorica per analizzare generalizzazione, dinamiche di apprendimento e stabilità nei modelli neurali complessi.
Una matrice aleatoria è una funzione misurabile definita su uno spazio di probabilità, i cui valori sono matrici a coefficienti reali o complessi. Formalmente, data una tripla che rappresenta uno spazio di probabilità, una matrice aleatoria è una funzione
dove oppure , e ogni elemento della matrice è una variabile aleatoria. La distribuzione congiunta degli elementi (oppure degli autovalori) della matrice è specificata da un modello probabilistico, detto ensemble.
Tra i casi più studiati si trovano:
In molti casi, l'interesse principale non risiede negli elementi della matrice, ma nel comportamento statistico degli autovalori, ad esempio mediante la misura spettrale empirica
dove sono gli autovalori della matrice . Tale approccio è alla base della moderna teoria delle matrici aleatorie e delle sue applicazioni in fisica teorica, statistica multivariata, teoria dei numeri e apprendimento automatico.
Nel contesto della teoria delle matrici aleatorie, un ensemble è una famiglia di matrici aleatorie dotata di una distribuzione di probabilità specifica, spesso determinata da proprietà di simmetria. Gli ensemble classici più studiati sono tre, introdotti e classificati da Freeman J. Dyson nel 1962 sulla base delle simmetrie fondamentali dei sistemi fisici e del comportamento sotto trasformazioni ortogonali, unitarie o simplettiche.
I principali ensemble gaussiani sono:
Questi ensemble condividono una densità di probabilità del tipo:
dove (GOE), (GUE), (GSE), e indica la traccia della matrice. è una matrice aleatoria e è la traccia del quadrato della matrice, ovvero la somma dei quadrati dei suoi autovalori. Il parametro codifica invece il tipo di simmetria. Dunque, è una densità: descrive quanto è probabile che la matrice aleatoria assuma certi valori spettrali. Il parametro ha anche interpretazione fisica in termini di temperatura nella meccanica statistica.
| Sigla | Denominazione | Tipologia | Sistema | |
|---|---|---|---|---|
| GOE | Ortogonale | Reale simmetrica | 1 | Sistemi con simmetria temporale |
| GUE | Unitaria | Complessa hermitiana | 2 | Sistemi senza simmetria temporale |
| GSE | Simplettica | Quaternionica hermitiana | 4 | Sistemi con spin e simmetria temporale |
Un'importante classe di matrici aleatorie non gaussiane è costituita dalle matrici di Wishart. Introdotte da John Wishart nel 1928 in ambito statistico, esse emergono naturalmente nel contesto dell’analisi multivariata, della teoria dell’informazione e della fisica statistica.
Sia una matrice rettangolare di dimensione i cui elementi sono variabili aleatorie gaussiane reali, indipendenti e identicamente distribuite, con media zero e varianza unitaria. La matrice
è detta matrice di Wishart. Si tratta di una matrice simmetrica semi-definita positiva di dimensione , che segue una distribuzione di Wishart. Più in generale, se le righe di sono distribuite secondo una legge normale multivariata , allora segue la distribuzione di Wishart con parametro di scala e gradi di libertà.
Le matrici di Wishart svolgono un ruolo cruciale in numerosi contesti. Per esempio:
Un ensemble importante associato a queste matrici è detto ensemble di Wishart–Laguerre, in quanto la densità spettrale degli autovalori segue, nel limite asintotico, la legge di Marchenko–Pastur:
dove e nel limite con costante.

Al di fuori del regime classico degli ensemble gaussiani densi, esiste un ampio interesse per le matrici aleatorie sparse e non gaussiane, che modellano sistemi fisici e reti con connettività limitata o distribuzioni di probabilità non standard.
Una matrice aleatoria sparsa è una matrice in cui la maggior parte degli elementi è nulla, e solo una frazione di essi è diversa da zero. Tipicamente, gli indici non nulli sono scelti in modo casuale, e i valori associati seguono una certa distribuzione (gaussiana, di Bernoulli, uniforme, ecc.).
Queste matrici emergono in:
Diversamente dagli ensemble gaussiani, gli autovalori delle matrici sparse possono mostrare localizzazione spettrale, comportamento non universale nei bordi dello spettro e assenza di legge semicircolare nei limiti di densità spettrale. Le distribuzioni degli elementi non devono essere gaussiane: studi recenti hanno esaminato matrici con elementi di tipo sub-gaussiano, di Lévy, o appartenenti a leggi a coda pesante (heavy-tailed).
Questi ensemble richiedono metodi differenti per l’analisi teorica, come il metodo delle cavità, il formalismo di replica e la teoria dei grafi aleatori.
Una delle principali questioni nello studio delle matrici aleatorie riguarda il comportamento dello spettro, ovvero l'insieme degli autovalori, e in particolare la loro distribuzione in limiti asintotici. Per matrici di grandi dimensioni, si osservano regolarità notevoli e universali nel modo in cui gli autovalori si distribuiscono.
Sia una matrice casuale con autovalori reali . La densità empirica degli autovalori è definita come
Nel limite , la densità converge (in senso debole) a una densità di probabilità deterministica , che dipende dalla classe di ensemble considerato.
Le proprietà dello spettro includono:
Questi risultati sono fondamentali sia in matematica pura sia in applicazioni a fisica, statistica, finanza e machine learning.
Nel contesto degli ensemble gaussiani simmetrici (come il GOE), la distribuzione limite degli autovalori è descritta dalla legge semicircolare di Wigner.
Sia una matrice simmetrica di dimensione i cui elementi sono variabili aleatorie indipendenti (fatta eccezione per la simmetria), con media nulla e varianza . Allora, nel limite , la densità degli autovalori converge a
Questa distribuzione è supportata sull'intervallo ed è detta legge semicircolare per via della forma del suo grafico. Il risultato, scoperto da Eugene Wigner negli anni '50, è uno dei primi esempi di comportamento universale nella teoria delle matrici aleatorie: la forma semicircolare emerge indipendentemente dalla distribuzione precisa degli elementi (purché soddisfino ipotesi di momenti).

È possibile derivare la distribuzione congiunta degli autovalori reali di una matrice casuale (GOE, GUE, GSE) mediante un cambiamento di base. La densità di probabilità congiunta ha la forma:
dove sono gli autovalori, è una costante di normalizzazione, un parametro di scala, e corrisponde rispettivamente agli ensemble GOE, GUE e GSE. Questa distribuzione implica una repulsione tra gli autovalori: la probabilità che due autovalori siano molto vicini () tende a zero, in netto contrasto con il comportamento di un processo di Poisson. Questa proprietà è nota come repulsione dei livelli.
La legge di Marčenko–Pastur descrive il comportamento asintotico dello spettro delle matrici di covarianza campionarie, ed è la controparte della legge semicircolare per matrici rettangolari.
Sia una matrice con elementi reali indipendenti e identicamente distribuiti, a media zero e varianza . Si considera la matrice di covarianza non normalizzata
che è una matrice di Wishart. Se con rapporto fisso, allora la densità empirica degli autovalori di converge (quasi sicuramente) alla distribuzione di Marčenko–Pastur:
dove i bordi dello spettro sono
Per , la distribuzione presenta un atomo di Dirac in zero con peso .
La legge di Marčenko–Pastur trova applicazioni in statistica multivariata (analisi delle componenti principali), fisica teorica (spettro delle matrici di tipo Wishart) ed apprendimento automatico e random feature models.
Oltre alla distribuzione globale degli autovalori, grande interesse teorico e applicativo riguarda il comportamento degli autovalori estremi, in particolare il massimo (o minimo) autovalore di matrici aleatorie di grande dimensione.

Per molte classi di matrici (come l’ensemble gaussiano ortogonale – GOE, unitario – GUE, o simmetrico reale), si dimostra che l’autovalore più grande non cresce indefinitamente con , ma fluttua intorno a un valore deterministico, con fluttuazioni descritte da una legge non gaussiana, detta legge di Tracy–Widom.
Sia una matrice del GUE di dimensione , allora vale:
dove è la distribuzione di Tracy–Widom di tipo 2. Analogamente, per il GOE si ottiene . Queste distribuzioni sono legate alle soluzioni dell’equazione di Painlevé II e non hanno espressioni chiuse elementari. Tuttavia, esse descrivono in modo universale le fluttuazioni dell’autovalore massimo per un’ampia classe di matrici, anche al di fuori dei modelli gaussiani.
La legge di Tracy–Widom compare in numerosi ambiti:
Nello studio delle matrici aleatorie non solo gli autovalori sono oggetto di analisi, ma anche le proprietà strutturali degli autovettori associati. Una questione centrale è se un autovettore associato a un certo autovalore sia localizzato (concentrato su poche componenti) oppure delocalizzato (distribuito uniformemente su tutte le componenti).
Formalmente, dati gli autovettori di una matrice aleatoria simmetrica , si definisce un indice di partecipazione come:
Un valore piccolo di indica che il vettore è distribuito su molte componenti (delocalizzato), mentre un valore grande segnala localizzazione. Per le matrici gaussiane dense (GOE, GUE), è noto che gli autovettori sono fortemente delocalizzati: ciascuna componente ha modulo quadrato dell’ordine di con alta probabilità. Viceversa, in modelli sparsi o con disordine forte (es. matrici a banda o matrici con struttura casuale), può emergere localizzazione. In tali casi, solo un sottoinsieme ristretto delle componenti dell'autovettore è significativamente diverso da zero.
La transizione tra comportamento localizzato e delocalizzato è legata a fenomeni di transizione di Anderson, con implicazioni nella:

I fondamenti della teoria delle matrici aleatorie sono stati gettati nella seconda metà del XX secolo attraverso una serie di lavori pionieristici che ne hanno definito il quadro concettuale e metodologico. Tra i contributi più influenti si annoverano:
| Controllo di autorità | LCCN (EN) sh86001920 · J9U (EN, HE) 987007543966505171 |
|---|