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.

414 CAPITOLUL 5. DESCOMPUNEREA VALORILOR SINGULARE<br />

De asemenea, cum era de aşteptat, condiţionarea întregului ansamblu de valori<br />

<strong>si</strong>ngulare, privită ca o normă a vectorului condiţionărilor <strong>valorilor</strong> <strong>si</strong>ngulare, este<br />

cea mai bună po<strong>si</strong>bilă. În sprijinul acestei afirmaţii vine şi teorema 5.13 care arată<br />

că norma euclidiană a vectorului variaţiilor absolute ale <strong>valorilor</strong> <strong>si</strong>ngulare este<br />

inferioară normei Frobenius a matricei variaţiilor elementelor matricei iniţiale.<br />

[ ]<br />

1.60 0.36 0.48<br />

Exemplul 5.4 Fie matricea A =<br />

∈ IR 2×3 ale cărei valori<br />

−1.20 0.48 0.64<br />

<strong>si</strong>ngulare exacte sunt σ 1 [ = 2 şi σ 2 = 1. Valorile ] <strong>si</strong>ngulare ale matricei [ perturbate ]<br />

1.60 0.36 0.481<br />

0 0 1<br />

F = A+E = A+ǫG =<br />

, unde ǫ = 10<br />

−1.20 0.48 0.64<br />

−3 , G =<br />

0 0 0<br />

(cu ‖G‖ = ‖G‖ F = 1), sunt ˆσ 1 ≈ 2.0000002, ˆσ 2 ≈ 1.0004800. Se observă că<br />

nici variaţiile absolute ale <strong>valorilor</strong> <strong>si</strong>ngulare individuale, nici norma euclidiană a<br />

vectorului lor √ (ˆσ 1 −σ 1 ) 2 +(ˆσ 2 −σ 2 ) 2 ≈ 0.48·10 −4 nu depăşesc valoarea lui ǫ. ✸<br />

Condiţionarea excelentă a <strong>valorilor</strong> <strong>si</strong>ngulare este unul din argumentele fundamentale<br />

ale utilizării lor pentru rezolvarea numerică a unei multitudini de probleme<br />

aplicative de algebră liniară.<br />

5.4.3 Condiţionarea <strong>vectorilor</strong> <strong>si</strong>ngulari<br />

Similar cazului matricelor hermitice (v. §4.10), perfecta condiţionare a <strong>valorilor</strong><br />

<strong>si</strong>ngulare nu implică în mod necesar buna condiţionare a <strong>vectorilor</strong> <strong>si</strong>ngulari şi a<br />

subspaţiilor generate de aceştia.<br />

Fie A ∈ IC m×n şi A = U H ΣV descompunerea <strong>valorilor</strong> sale <strong>si</strong>ngulare. Aprecierea<br />

condiţionării <strong>vectorilor</strong> <strong>si</strong>ngulari, i.e a coloanelor matricelor unitare de transformare<br />

U şi V, se face prin evaluarea diferenţei unghiulare dintre vectorul exact şi cel<br />

perturbat, definită prin<br />

θ(u i ,û i ) = arccos|u H i û i |, i = 1 : m θ(v j ,ˆv j ) = arccos|v H j ˆv j |, j = 1 : n,<br />

(5.110)<br />

raportată la norma variaţiei matricei date.<br />

Această abordare poate fi extinsă la exprimarea condiţionării subspaţiilor generate<br />

de vectori <strong>si</strong>ngulari care este apreciată prin variaţia unghiulară (v. §4.10)<br />

a subspaţiilor respective. Concret, fiind dată o matrice A şi σ I ⊂ σ(A) un set de<br />

valori <strong>si</strong>ngulare ale acesteia, prin condiţionarea subspaţiului U, generat de vectorii<br />

<strong>si</strong>ngulari asociaţi setului λ I , vom înţelege variaţia unghiulară (sau o margine superioară<br />

a acesteia) a subspaţiului U raportată la nivelul perturbaţiilor în elementele<br />

matricei A.<br />

Condiţionarea subspaţiilor generate de vectori <strong>si</strong>ngulari este determinată în<br />

mod deci<strong>si</strong>v de localizarea <strong>valorilor</strong> <strong>si</strong>ngulare asociate. Este po<strong>si</strong>bil, şi aici, ca un<br />

subspaţiu generat de vectori <strong>si</strong>ngulari rău condiţionaţi să aibă o condiţionare foarte<br />

bună dacă grupul corespunzător de valori <strong>si</strong>ngulare este bine separat de celelalte.<br />

În contextul DVS, vom nota separarea unei valori <strong>si</strong>ngulare individuale σ i şi,<br />

respectiv a setului de valori <strong>si</strong>ngulare σ I , de celelalte valori <strong>si</strong>ngulare ale aceleiaşi

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

Saved successfully!

Ooh no, something went wrong!