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.

5.6. APLICAŢIILE DVS 425<br />

[ (A(i,:))<br />

T<br />

normei ‖ · ‖ D ) de la punctele z =<br />

b(i)<br />

apropiate ale lui P x ∗ (vezi (5.136)) este minimă.<br />

]<br />

∈ IR n+1 la punctele cele mai<br />

✸<br />

A treia treaptă de generalizare o introducem con<strong>si</strong>derând un membru drept<br />

multiplu, i.e. formulând problema CMMPT pentru <strong>si</strong>stemul liniar matriceal AX =<br />

= B cu B ∈ IC m×p , respectiv problema de minimizare<br />

‖CG ∗ D‖ = min<br />

(B +R) ⊆ Im(A+E) ‖CGD‖ F, E ∈ IC m×n , R ∈ IC m×p , G = [E R],<br />

(5.137)<br />

unde matricele A ∈ IC m×n , B ∈ IC m×p , cu m ≥ n+p, precum şi matricele diagonale<br />

ne<strong>si</strong>ngulare C ∈ IC m×m şi D ∈ IC (n+p)×(n+p) sunt date. La fel ca şi până acum,<br />

dacă (E ∗ ,R ∗ ) este o soluţie a problemei de minimizare (5.137), atunci orice matrice<br />

X ∗ ∈ IC n×p care satisface <strong>si</strong>stemul liniar matriceal (A+E ∗ )X = B+R ∗ va fi numită<br />

(pseudo)soluţia, în sens CMMPT, a <strong>si</strong>stemului AX = B.<br />

Pentru a formula mai concis rezultatul referitor la existenţa şi unicitatea soluţiei<br />

problemelor de minimizare ce definesc CMMPT, vom introduce unele notaţii şi vom<br />

stabili un rezultat preliminar. Fie<br />

H def<br />

= C[ A B ]D =<br />

n p<br />

{}}{ {}}{<br />

[ ]<br />

H1 H 2<br />

(5.138)<br />

şi H = UΣV H DVS a matricei H, cu următoarele partiţii ale matricelor U, V şi Σ<br />

impuse de structura lui H<br />

U =<br />

n<br />

p<br />

{}}{ m−n−p<br />

{}}{ {}}{<br />

[ ]<br />

U1 U 2 U 3<br />

, V =<br />

n p<br />

{}}{ {}}{<br />

[ ]<br />

V11 V 12 }n<br />

V 21 V 22 }p<br />

Σ =<br />

n<br />

p<br />

{}}{ {}}{<br />

⎡<br />

⎣ Σ ⎤<br />

1 0<br />

0 Σ 2<br />

⎦ }n<br />

}p ,<br />

0 0 }m−n−p<br />

Introducem următoarea lemă.<br />

Σ 1 = diag(σ 1 ,σ 2 ,...,σ n )<br />

Σ 2 = diag(σ n+1 ,σ n+2 ,...,σ n+p ). (5.139)<br />

Lema 5.1 Dacă σ n (H 1 ) > σ n+1 , atunci<br />

1 ◦ . Matricea V 22 din (5.139) este ne<strong>si</strong>ngulară.<br />

2 ◦ . Cu notaţiile din (5.139), avem inegalitatea strictă<br />

σ n > σ n+1 . (5.140)<br />

Demonstraţie. 1 ◦ . Presupunem că matricea V 22 este <strong>si</strong>ngulară. Atunci există un<br />

vector z ∈ IC p nenul, pe care îl putem con<strong>si</strong>dera de normă euclidiană [ ] unitară, astfel<br />

def V12<br />

încât V 22 z = 0. Mai departe, din faptul că matricea V 2 = are coloanele<br />

V 22<br />

ortogonale, i.e. V2 H V 2 = I p , obţinem ‖V 2 z‖ = ‖V 12 z‖ = 1. Pe de altă parte din

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

Saved successfully!

Ooh no, something went wrong!