Decomposizione QR
Da Wikipedia, l'enciclopedia libera.
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.
[modifica] Calcolo
Si può calcolare esplicitamente la fattorizzazione QR di una matrice data in tempo O(n3) operazioni aritmetiche attraverso l'uso delle trasformazioni di Householder o di quelle di Givens. Qui è riportato un algoritmo che utilizza le prime.
[modifica] 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(Qtb)
Il calcolo di c = Qtb richiede O(n2) operazioni, mentre il calcolo di x = R − 1b si può effettuare attraverso un algoritmo di sostituzione all'indietro sempre con O(n2) 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.
