Numero di Condizionamento

Indice dei contenuti

    Il numero di condizionamento di una matrice misura quanto la soluzione di un sistema lineare è sensibile alle perturbazioni dei dati: un numero di condizionamento elevato indica che piccoli errori nei dati (o negli arrotondamenti) possono produrre grandi errori nella soluzione.

    Vedi anche: Virgola Mobile, Norma, Decomposizione SVD.

    Norme Vettoriali e Matriciali

    Una norma matriciale indotta dalla norma vettoriale \|\cdot\| è:

    A=supx0Axx\|A\| = \sup_{\vec{x} \neq 0} \frac{\|A\vec{x}\|}{\|\vec{x}\|}

    Norme matriciali comuni:

    NormaDefinizioneCalcolo
    A1\|A\|_1massima somma di colonna$\max_j \sum_i
    A\|A\|_\inftymassima somma di riga$\max_i \sum_j
    A2\|A\|_2 (spettrale)massimo valore singolareσmax(A)\sigma_{\max}(A)
    AF\|A\|_F (Frobenius)radice somma quadratiijaij2\sqrt{\sum_{ij} a_{ij}^2}

    Definizione del Numero di Condizionamento

    Per una matrice invertibile AA:

    κ(A)=AA1\kappa(A) = \|A\| \cdot \|A^{-1}\|

    Con la norma spettrale: κ2(A)=σmax(A)/σmin(A)\kappa_2(A) = \sigma_{\max}(A) / \sigma_{\min}(A) (rapporto tra il massimo e il minimo valore singolare). Vedi: Decomposizione SVD.

    Per una matrice ortogonale: κ2(Q)=1\kappa_2(Q) = 1 (perfettamente condizionata). Per una matrice singolare: κ=\kappa = \infty.

    Sensibilità dei Sistemi Lineari

    Dato il sistema Ax=bA\vec{x} = \vec{b}, se i dati sono perturbati di δA\delta A e δb\delta \vec{b}, l’errore relativo sulla soluzione soddisfa:

    δxxκ(A)(δAA+δbb)\frac{\|\delta\vec{x}\|}{\|\vec{x}\|} \leq \kappa(A) \left(\frac{\|\delta A\|}{\|A\|} + \frac{\|\delta\vec{b}\|}{\|\vec{b}\|}\right)

    Un sistema con κ(A)=10k\kappa(A) = 10^k può perdere fino a kk cifre decimali significative per effetto degli errori di arrotondamento.

    Interpretazione Geometrica

    κ2(A)\kappa_2(A) è il rapporto tra il semiasse maggiore e il semiasse minore dell’ellissoide immagine della sfera unitaria. Una matrice mal condizionata «schiaccia» lo spazio in alcune direzioni più di altre.

    Miglioramento del Condizionamento: Precondizionamento

    Per sistemi mal condizionati si usa un precondizionatore PA1P \approx A^{-1} per risolvere P1Ax=P1bP^{-1}A\vec{x} = P^{-1}\vec{b}, che ha numero di condizionamento κ(P1A)κ(A)\kappa(P^{-1}A) \ll \kappa(A). Vedi: Gradiente Coniugato.

    Applicazioni ingegneristiche

    • FEM: le matrici di rigidezza di strutture con elementi molto diversi per dimensione o rigidezza hanno κ1\kappa \gg 1; il precondizionamento è essenziale per la convergenza dei solvitori iterativi.
    • Regolarizzazione: nei problemi inversi (tomografia, ricostruzione di immagini) le matrici di osservazione sono mal condizionate; la regolarizzazione di Tikhonov aggiunge λI\lambda I per limitare κ\kappa.
    • Reti neurali: il condizionamento della matrice hessiana della funzione di perdita influenza la velocità di convergenza della discesa del gradiente; metodi del secondo ordine (Adam, L-BFGS) stimano e compensano il mal condizionamento.

    Ultimo aggiornamento: