12.07.2015 Views

confronto fattorizzazioni - Università degli studi di Cagliari.

confronto fattorizzazioni - Università degli studi di Cagliari.

confronto fattorizzazioni - Università degli studi di Cagliari.

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

1.3.1 FATTORIZZAZIONE DI HOUSEHOLDERLa fattorizzazione <strong>di</strong> Householder consiste nella decomposizione della matrice A indue matrici Q ed R a partire dalle matrici elementari <strong>di</strong> Householder.Una matrice elementare <strong>di</strong> Householder reale ha la seguente forma:H=I−2 w w T ,w∈R n ,∥w∥=1,H è una matrice ortogonale e simmetrica.Questo metodo presenta alcuni vantaggi pratici, in particolare esso non richiede ilcalcolo dell’inversa <strong>di</strong> una matrice triangolare superiore.Poiché ogni trasformazione <strong>di</strong> Householder è ortogonale, ci basta far vedere cheesiste un prodotto <strong>di</strong> matrici del tipo H(x i ; e j ) che trasforma A in una matricetriangolare superiore con entrate sulla <strong>di</strong>agonale principale tutte positive. A questopunto, l’unicità della fattorizzazione garantisce che quanto abbiamo ottenuto siaeffettivamente la decomposizione QR cercata, che è detta matrice <strong>di</strong> Householder(alternativamente riflessione <strong>di</strong> Householder, trasformazione <strong>di</strong> Householder)associata a w.H 1 a 1(1)=k 1 e 1.Dobbiamo creare una successione <strong>di</strong> matrici A (i) , con i=1,...n in modo che A (n) siatriangolare superiore. Moltiplichiamo la matrice <strong>di</strong> Householder H 1 a sinistra dellamatrice A (1) e otteniamo:A (2) =H 1 A (1) =[a 1(2)a 2(2)… a n(2)]in cuia1(2)=k1e1,mentre i restanti elementiaj(2)=H 1 aj(1) con j=2,...,n.La matrice A(2) ha la seguente struttura:A (2) =(k(2) (2)1 a 12 ⋯ a 1n(2) (2)0 a 22⋯ a 2n⋮ ⋮ ⋮(2)0 a n2 ⋯ a nn(2))= (k 10Tv 1 Â , 13

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!