Decomposizione QR

La decomposizione QR o fattorizzazione QR di una matrice quadrata a coefficienti reali o complessi $M$ è una scomposizione del tipo

$M=QR,$

dove $Q$ è una matrice ortogonale, e $R$ è una matrice triangolare superiore. Si può dimostrare che tutte le matrici quadrate ammettono una decomposizione QR, anche se essa non è unica. Nel caso in cui la matrice $M$ sia a coefficienti complessi, allora $Q$ è una matrice unitaria.

Calcolo

Si può calcolare esplicitamente la fattorizzazione QR di una matrice data in tempo $O(n^3)$ operazioni aritmetiche attraverso l'uso delle trasformazioni di Householder o di quelle di Givens.

Applicazioni

L'applicazione principale della fattorizzazione QR è la soluzione di sistemi lineari: una volta fattorizzata la matrice $A=QR$ di un sistema lineare $Ax=b$, la soluzione del sistema è data da

$x=R^{-1} (Q^t b)$

Il calcolo di $c=Q^t b$ richiede $O(n^3)$ operazioni, mentre il calcolo di $x=R^{-1}c$ si può effettuare attraverso un algoritmo di sostituzione all'indietro sempre con $O(n^2)$ operazioni. Il costo dominante quindi è proprio quello della fattorizzazione.

La complessità computazionale di questo metodo risolutivo è quindi la stessa della soluzione mediante la fattorizzazione LU (o algoritmo di Gauss), ma questo algoritmo risulta avere una migliore stabilità numerica.

Inoltre, la fattorizzazione QR può essere utilizzata per il calcolo di basi ortonormali e per la risoluzione di un sistema ai minimi quadrati.

Portale Matematica: accedi alle voci di Wikipedia che trattano di matematica