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.

426 CAPITOLUL 5. DESCOMPUNEREA VALORILOR SINGULARE<br />

DVS a matricei H, cu partiţiile din (5.139), avem U 2 Σ 2 = H 1 V 12 +H 2 V 22 . Acum,<br />

ţinând seama de con<strong>si</strong>stenţa normei spectrale, putem scrie următoarea secvenţă de<br />

inegalităţi<br />

σ n+1 = ‖Σ 2 ‖ = ‖U 2 ‖·‖Σ 2 ‖ ≥ ‖U 2 Σ 2 ‖ = ‖H 1 V 12 +H 2 V 22 ‖ =<br />

= max<br />

‖y‖=1 ‖(H 1V 12 +H 2 V 22 )y‖ ≥ ‖H 1 V 12 z‖ ≥ min<br />

‖w‖=1 ‖H 1w‖ = σ n (H 1 ),<br />

ceea ce contrazice ipoteza lemei. Deci V 22 este ne<strong>si</strong>ngulară.<br />

2 ◦ . Conform teoremei 5.11 (de separare a <strong>valorilor</strong> <strong>si</strong>ngulare) avem<br />

σ n = σ n (H) ≥ σ n (H [n+p−1] ) ≥ ... ≥ σ n (H [n] ) def<br />

= σ n (H 1 ), (5.141)<br />

de unde, în ipoteza lemei, rezultă σ n ≥ σ n (H 1 ) > σ n+1 , obţinându-se inegalitatea<br />

strictă din enunţ.<br />

✸<br />

Formulăm acum teorema de existenţă şi unicitate a soluţiei problemei CMMPT.<br />

Teorema 5.16 Utilizând notaţiile (5.138), (5.139), dacă σ n (H 1 ) > σ n+1 , atunci<br />

matricea G ∗ def<br />

= [ E ∗ R ∗ ] definită de<br />

G ∗ = [ E ∗ R ∗ ] = −C −1 U 2 Σ 2 [ V H<br />

12 V H<br />

22 ]D−1 (5.142)<br />

este o soluţie a problemei de minimizare CMMPT (5.137).<br />

În plus, dacă notăm<br />

D 1 = diag(d 1 ,d 2 ,...,d n ), D 2 = diag(d n+1 ,d n+2 ,...,d n+p ), (5.143)<br />

atunci matricea<br />

X ∗ = −D 1 V 12 V22 −1 D−1 2 (5.144)<br />

există şi este unica soluţie a <strong>si</strong>stemului<br />

(A+E ∗ )X = B +R ∗ , (5.145)<br />

i.e. este unica (pseudo)soluţie, în sens CMMPT, a <strong>si</strong>stemului liniar AX = B.<br />

Demonstraţie. Condiţia (B +R) ⊆ Im(A+E) este echivalentă cu existenţa unei<br />

matrice X ∈ IC n×p astfel încât (A + E)X = B + R. Cu notaţiile din (5.137) şi<br />

(5.138), ultima relaţie poate fi scrisă în următoarele forme echivalente<br />

([ ] ) [ ]<br />

[ ]<br />

X A B +G = 0 ⇔ (H +CGD)D −1 X<br />

= 0. (5.146)<br />

−I p −I p<br />

[ ]<br />

X<br />

Întrucât matricea este monică, din (5.146) rezultă rang(H +CGD) ≤ n.<br />

−I p<br />

NotândF def<br />

= CGD, încontinuareaplicămteorema5.15matricei−F = H−(H+F).<br />

Obţinem<br />

n+p<br />

min ‖F‖ 2 F = ∑<br />

σi 2 , (5.147)<br />

rang(H+F)≤n<br />

i=n+1

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

Saved successfully!

Ooh no, something went wrong!