12.01.2015 Views

Calculul valorilor si vectorilor proprii

Calculul valorilor si vectorilor proprii

Calculul valorilor si vectorilor proprii

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

4.1. FORMULAREA PROBLEMEI 215<br />

În cazul general, structura ”fină” a unei matrice, care poate fi dezvăluită prin<br />

transformărideasemănarecorespunzătoare,estedatădeaşanumitaformă canonică<br />

Jordan. Deşi forma canonică Jordan joacă un rol esenţial în analiza matriceală,<br />

conţinând maximum de informaţie structurală privitor la o matrice dată, totuşi<br />

rolul ei în calculul numeric este mult diminuat de sen<strong>si</strong>bilitatea structurii Jordan la<br />

perturbaţii numerice în elementele matricei iniţiale, perturbaţii inerente în calcule<br />

efectuate pe un calculator datorită reprezentării informaţiei numerice în virgulă<br />

mobilă. Acesta este motivul pentru care în toate dezvoltările numerice se preferă<br />

o structură mult mai robustă şi anume forma Schur reală sau complexă prezentată<br />

într-una din secţiunile următoare 5 .<br />

4.1.4 Valorile <strong>proprii</strong> ale matricelor <strong>si</strong>metrice şi hermitice<br />

Prezentăm în continuare câteva rezultate referitoare la valorile şi vectorii <strong>proprii</strong><br />

pentru matricele hermitice (<strong>si</strong>metrice). Matricele hermitice (<strong>si</strong>metrice) se întâlnesc<br />

în numeroase aplicaţii şi prezintă particularităţi remarcabile.<br />

Definiţia 4.4 Fie A ∈ IC n×n . Matricea A se numeşte normală dacă<br />

A H A = AA H . (4.19)<br />

În cazul real, matricea A ∈ IR n×n este normală dacă<br />

A T A = AA T . (4.20)<br />

În acest context reamintim că matricea A se numeşte hermitică dacă A H = A<br />

şi <strong>si</strong>metrică dacă A T = A. De asemenea, o matrice A ∈ IC n×n se numeşte unitară<br />

dacă A H A = I n şi ortogonală dacă A T A = I n .<br />

Se constată imediat că matricele hermitice şi cele unitare sunt matrice normale.<br />

Matricele hermitice au proprietatea că elementele <strong>si</strong>metrice faţă de diagonala principală<br />

sunt complex conjugate, i.e. a ij = ā ji , i,j ∈ 1 : n, deci elementele diagonale<br />

ale matricelor hermitice sunt reale. O matrice hermitică reală este <strong>si</strong>metrică. O<br />

matrice unitară reală este ortogonală. Prin urmare matricele reale <strong>si</strong>metrice sau<br />

ortogonale sunt normale. Există [ matrice ] normale care nu sunt nici <strong>si</strong>metrice nici<br />

1 −1<br />

ortogonale, de exemplu A = .<br />

1 1<br />

Prezentăm în continuare câteva rezultate fundamentale, urmând ca aspectele<br />

specifice legate de calculul efectiv al <strong>valorilor</strong> şi <strong>vectorilor</strong> <strong>proprii</strong> pentru matrice<br />

hermitice (<strong>si</strong>metrice) să fie date în două secţiuni distincte (§4.8 şi §4.9), iar cele<br />

legate de condiţionare şi stabilitate în §4.10 şi §4.11.<br />

Teorema 4.1 O matrice n×n complexă A este normală dacă şi numai dacă admite<br />

un set complet de n vectori <strong>proprii</strong> ortogonali, adică există o matrice unitară Q ∈<br />

∈ IC n×n ale cărei coloane sunt vectori <strong>proprii</strong> ai matricei A astfel încât<br />

Q H AQ = Λ = diag(λ 1 ,λ 2 ,...,λ n ) ∈ IC n×n . (4.21)<br />

5 Algoritmii de reducere la forma canonică Jordan, prezentaţi în unele lucrări de matematică<br />

(vezi, e.g. [XVI]) nu prezintă interes practic decât în contextul unor medii de calcul exact. Pentru<br />

detalii privitoare la aspectele numerice şi algoritmice ale calculului formei canonice Jordan, vezi<br />

secţiunea 4.7.

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

Saved successfully!

Ooh no, something went wrong!