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 219<br />

Demonstraţie. Conform teoremei 4.2, matricea Q este unitară, A = QΛQ H unde<br />

Λ = diag(λ 1 ,λ 2 ,...,λ n ) şi, prin urmare,<br />

µ(x) = x H Ax = y H Λy =<br />

n∑<br />

λ k |y (k) | 2 , y = Q H x = [y (1) y (2) ··· y (n) ] T .<br />

k=1<br />

(4.30)<br />

Cum vectorii x şi y din (4.30) se află într-o relaţie biunivocă, iar transformările<br />

unitare conservă norma euclidiană, rezultă că extremele funcţiei µ coincid cu extremele<br />

funcţiei ν : S → IR, ν(y) = y H Λy. Din faptul că vectorii y sunt de normă<br />

unitară, i.e. ∑ n<br />

j=1 |y(j) | 2 = 1, rezultă<br />

ν(y) = λ 1 −<br />

n∑<br />

n−1<br />

∑<br />

(λ 1 −λ j )|y (j) | 2 = (λ j −λ n )|y (j) | 2 +λ n . (4.31)<br />

j=2<br />

Întrucât sumele din relaţia (4.31) sunt, datorită (4.26), nenegative, iar valoarea<br />

nulă a acestor sume se poate realiza, e.g. pentru y (j) = 0, j = 2 : n în primul<br />

caz şi j = 1 : n − 1 în cel de al doilea, avem egalităţile (4.28). Dacă valorile<br />

<strong>proprii</strong> maximă, respectiv minimă, sunt <strong>si</strong>mple, atunci valorile extreme ale funcţiei<br />

ν se ating pentru vectorii y de forma y 1 = [y (1) 0 ··· 0] T = e iθ1 e 1 , respectiv<br />

y n = [0 ··· 0 y (n) ] T = e iθn e n , cu θ 1 , θ n ∈ IR. Prin urmare, cele două extreme ale<br />

funcţiei µ se ating pentru vectorii x de forma x 1 = e iθ1 q 1 şi, respectiv x n = e iθn q n .<br />

Dacă λ 1 are multiplicitatea s, iar λ n multiplicitatea t, atunci maximul se atinge<br />

pentru orice vector x de normă unitară din V s = ImQ ′ s, i.e. subspaţiul A-invariant<br />

asociat <strong>valorilor</strong> <strong>proprii</strong> λ j , j = 1 : s, iar minimul se atinge pentru orice vector de<br />

normă unitară din W n−t .<br />

Pentru cea de a doua parte a teoremei, dacă x ∈ W k = V ⊥ k atunci xH Q ′ k = 0 şi<br />

y = Q H x = [0 ··· 0 y (k+1) ··· y (n) ] T . Prin urmare,<br />

µ(x) = ν(y) = λ k+1 −<br />

n∑<br />

j=k+2<br />

j=1<br />

de unde, cu aceleaşi argumente ca mai sus, se obţine (4.29).<br />

(λ k+1 −λ j )|y (j) | 2 , (4.32)<br />

Rezultatul următor prezintă o interesantă caracterizareminimax a <strong>valorilor</strong><strong>proprii</strong><br />

ale unei matrice hermitice (în cazul real, <strong>si</strong>metrice) şi este util prin consecinţele<br />

sale. Notăm, generic, cu V subspaţiile liniare ale spaţiului IC n şi cu W = V ⊥ complementele<br />

lor ortogonale în IC n . De asemenea, vom nota cu V S = V ∩S şi, respectiv,<br />

W S = W ∩S, mulţimile <strong>vectorilor</strong> de normă euclidiană unitară din V şi W.<br />

Teorema 4.4 (Courant – Fisher) Dacă matricea hermitică A ∈ IC n×n are valorile<br />

<strong>proprii</strong> ordonate ca în (4.26) atunci pentru toţi k ∈ 1 : n avem<br />

λ k = max<br />

dimV = k<br />

min x H Ax = min<br />

x ∈ V S dimV = k<br />

✸<br />

max x H Ax 8 . (4.33)<br />

x ∈ W S<br />

8 Întrucât oricărui subspaţiu n − k dimen<strong>si</strong>onal din IC n îi corespunde un complement ortogonal<br />

k dimen<strong>si</strong>onal, ultimul termen al egalităţilor (4.33) poate fi scris şi în forma λ k =<br />

= min dimV = n−k max x ∈ VS x H Ax.

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

Saved successfully!

Ooh no, something went wrong!