12.01.2015 Views

Calculul valorilor si vectorilor proprii

Calculul valorilor si vectorilor proprii

Calculul valorilor si vectorilor proprii

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

452 CAPITOLUL 6. VALORI ŞI VECTORI PROPRII GENERALIZAŢI<br />

Deci, dimV ≤ k, adică S este un subspaţiu de deflaţie k-dimen<strong>si</strong>onal al perechii<br />

(A,B). Cu alte cuvinte, primele k coloane ale matricei de transformare Z, i.e.<br />

primii k vectori Schur la dreapta ai perechii (A,B), formează o bază ortogonală a<br />

subspaţiului de deflaţie k-dimen<strong>si</strong>onal asociat <strong>valorilor</strong> <strong>proprii</strong> generalizate definite<br />

de perechile (s ii ,t ii ), i = 1 : k. În acest mod, prin ordonarea corespunzătoare a<br />

elementelor diagonale ale matricelor S şi T, se pot construi baze ortogonale pentru<br />

subspaţiide deflaţieasociateunorgrupuriimpuse devalori<strong>proprii</strong>generalizate(vezi<br />

secţiunea 6.4).<br />

6.2.2 Forma Schur reală generalizată<br />

În cazul în care matricele A, B sunt reale se obţine un spor important de eficienţă<br />

dacă se utilizează în exclu<strong>si</strong>vitate o aritmetică reală. Corespondentul generalizat al<br />

formei Schur reale este introdus prin următoarea teoremă pe care o prezentăm fără<br />

demonstraţie.<br />

Teorema 6.2 (Forma Schur reală generalizată) Oricare ar fi perechea (A,B) ∈<br />

∈ IR n×n ×IR n×n există matricele ortogonale Q,Z ∈ IR n×n astfel încât<br />

Q T AZ = S, Q T BZ = T, (6.16)<br />

unde matricea S este în formă Schur reală iar matricea T este superior triunghiulară.<br />

Perechile de blocuri diagonale (S ii ,T ii ), i = 1:p, de dimen<strong>si</strong>uni 1×1 sau 2×2<br />

ale matricelor S şi T determină valorile <strong>proprii</strong> generalizate ale perechii (A,B), mai<br />

precis dacă blocul diagonal i al lui S are ordinul n i , atunci ∑ p<br />

i=1 n i = n şi<br />

λ(A,B) = ∪ p i=1 λ(S ii,T ii ). (6.17)<br />

Perechile de blocuri diagonale pot fi dispuse în orice ordine predeterminată.<br />

Perechea (S,T) se numeşte forma Schur reală generalizată (FSRG) a perechii<br />

(A,B), iar coloanele q i , respectiv z i , ale matricelor ortogonale de transformare Q<br />

şi Z se numesc vectori Schur generalizaţi ai perechii (A,B) la stânga, respectiv la<br />

dreapta, asociaţi FSRG.<br />

Conform(6.17), dacădispunemdeFSRGauneiperechi(A,B), calculul<strong>valorilor</strong><br />

<strong>proprii</strong> generalizate se reduce la rezolvarea ecuaţiilor algebrice<br />

det(S ii −λT ii ) = 0, i = 1:p, (6.18)<br />

de grad cel mult doi.<br />

Toate con<strong>si</strong>deraţiile făcute în legătură cu FSG au un corespondent transparent<br />

pentru FSRG. De exemplu, dacă dimen<strong>si</strong>unea cumulată a primelor l blocuri diagonale<br />

ale matricei S este k, atunci primele k coloane ale matricei ortogonalede transformare<br />

Z formează o bază ortogonală a unui subspaţiu de deflaţie k-dimen<strong>si</strong>onal<br />

(din IR n ) al perechii (A,B) asociat ”primelor” k valori <strong>proprii</strong> generalizate.<br />

Din cele de mai sus rezultă că problema de calcul a <strong>valorilor</strong> <strong>proprii</strong> generalizate<br />

se reduce, în esenţă, la obţinerea formei Schur (reale) generalizate. Cum acest lucru<br />

nu este po<strong>si</strong>bil, în cazul general, printr-o secvenţă finită de operaţii elementare,<br />

calculul va fi bazat, în mod necesar, pe trunchierea unui proces infinit, <strong>si</strong>milar<br />

algoritmuluiQR.Înformasaceamaiperformantăaceastăprocedurăestecunoscută<br />

sub numele de algoritm QZ şi este prezentată în secţiunea ce urmează.

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

Saved successfully!

Ooh no, something went wrong!