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.

220 CAPITOLUL 4. VALORI ŞI VECTORI PROPRII<br />

Demonstraţie. Fie V un subspaţiu arbitrar de dimen<strong>si</strong>une k şi v j , j = 1 : k, o bază<br />

a lui V. Fie, de asemenea, w j , j = 1 : n−k, o bază a lui W. Notăm cu V ∈ IC k ,<br />

respectiv W ∈ IC n−k , matricele <strong>vectorilor</strong> care formează bazele celor două subspaţii<br />

complementare. Conform teoremei precedente<br />

λ n ≤ x H Ax ≤ λ 1 (4.34)<br />

pentru toţi x din S, i.e. funcţia µ este mărginită pe compactul V S şi, în consecinţă,<br />

îşi atinge marginile pe această mulţime. La fel ca în demonstraţia teoremei precedente,<br />

fie y = Q H x, unde Q este o matrice unitară de vectori <strong>proprii</strong>, ordonaţi<br />

conform (4.26). Avem, evident, ‖y‖ = ‖x‖ şi x = Qy ∈ V dacă şi numai dacă este<br />

ortogonal pe W, i.e.<br />

W H x = W H Qy = 0. (4.35)<br />

[ ]<br />

Întrucât W este monică, factorizarea QR a matricei ˜W = Q H W = ˜Q R<br />

are<br />

0<br />

matriceasuperior triunghiularăR ∈ IC (n−k)×(n−k) ne<strong>si</strong>ngulară. În consecinţă, (4.35)<br />

devine [<br />

R H 0 ] ˜QH y = 0. (4.36)<br />

Notând z def<br />

= ˜Q H y relaţia (4.36) impune z(1 : n−k) = 0. Notând, încă o dată,<br />

u def<br />

= z(n−k +1 : n) ∈ IC k şi ţinând seama de faptul că transformările unitare<br />

conservă norma euclidiană, din (4.35), (4.36) rezultă că x = Qy = Q˜Qz = ˆQu, unde<br />

ˆQ = Q˜Q(:,n−k+1 : n), aparţine mulţimii V S dacă şi numai dacă ‖u‖ = 1, fără<br />

nici o altă restricţie asupra lui u. Acum, putem alege u astfel încât y(1 : k−1) = 0.<br />

Într-adevăr, y = ˜Q(:,n−k+1: n)u şi orice soluţie normată(i.e. de normă euclidiană<br />

unitară)a<strong>si</strong>stemuluisubdeterminat ˆQ(1 : k−1,,n−k+1 : n)u = 0a<strong>si</strong>gurăsatisfacerea<br />

acestei condiţii. Cu această alegere a lui u, pentru vectorul corespunzător x din V S ,<br />

avem<br />

n∑<br />

µ(x) = x H Ax = y H Λy = λ k − (λ k −λ j )|y (j) | 2 ≤ λ k , (4.37)<br />

j=k+1<br />

unde am ţinut seama de faptul că ∑ n<br />

j=k |y(j) | 2 = ‖y‖ 2 = 1 şi de ordonarea descrescătoare<br />

a <strong>valorilor</strong> <strong>proprii</strong>. Natural, din (4.37) rezultă<br />

min x H Ax ≤ λ k (4.38)<br />

x ∈ V S<br />

şi, cum subspaţiul V, de dimen<strong>si</strong>une k, era arbitrar, inegalitatea (4.38) are loc în<br />

toate subspaţiile de aceeaşi dimen<strong>si</strong>une sau, altfel spus,<br />

max<br />

dimV = k<br />

min x H Ax ≤ λ k . (4.39)<br />

x ∈ V S<br />

Rămâne să arătăm că această margine este atinsă efectiv. Aceasta se întâmplă în<br />

subspaţiul A-invariant asociat primelor k valori <strong>proprii</strong> din secvenţa (4.26). Întradevăr,<br />

fie V = ImQ ′ k şi x = Q′ k z cu ‖z‖ = 1. Rezultă ‖x‖ = 1, i.e. x ∈ V S şi<br />

k−1<br />

∑<br />

µ(x) = x H Ax = (λ j −λ k )|z (j) | 2 +λ k ≥ λ k , (4.40)<br />

j=1

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

Saved successfully!

Ooh no, something went wrong!