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.

INDICAŢII, RĂSPUNSURI, SOLUŢII 519<br />

cazuri particulare de matrice normale, demonstraţia poate fi mai directă. De exemplu,<br />

dacă matricea unitară Q este, e.g. superior triunghiulară, atunci inversa ei Q −1 = Q H<br />

este <strong>si</strong>multan superior [ şi inferior ] triunghiulară, i.e. diagonală.<br />

α β<br />

P4.21 Fie A = . Din A T A = AA T rezultă β 2 = γ 2 . Dacă γ = β matricea<br />

γ δ<br />

este <strong>si</strong>metrică, iar dacă γ = −β rezultă δ = α.<br />

P4.22 Fie A normală şi S = Q H AQ o formă Schur reală a lui A, în care, fără a reduce<br />

generalitatea, putem presupune că valorile <strong>proprii</strong> reale (în număr de q) sunt <strong>si</strong>tuate în<br />

primele q poziţii diagonale. Deci S este normală şi are structura<br />

⎡<br />

⎤<br />

R A 1,q+1 ··· A 1p<br />

A q+1,q+1 ··· A q+1,p<br />

S = ⎢<br />

⎣ . ..<br />

. ⎥ .. ⎦ ,<br />

A pp<br />

cu R ∈ IR q×q superior triunghiulară. Din (S T S = SS T ) 11 rezultă R T R = RR T +<br />

+ ∑ p<br />

Bj, unde j=q+1 Bj = A1jAT 1j, j = q + 1 : p, sunt matrice <strong>si</strong>metrice, pozitiv semidefinite.<br />

Cum însă tr(R T R) = tr(RR T ), rezultă ∑ p<br />

trBj = 0. În continuare, din<br />

j=q+1<br />

faptul că λ i(B j) ≥ 0 pentru toţi i, rezultă trB j = ∑ λi(Bj) ≥ 0. Deci, trBj = 0 pentru<br />

i<br />

toţi j şi, prin urmare, λ i(B j) = 0 pentru toţi i şi j. Cum însă o matrice <strong>si</strong>metrică având<br />

toate valorile <strong>proprii</strong> nule este nulă (demonstraţi!) B j = 0 şi, de aici, A 1j = 0 pentru toţi<br />

j. Acum R T R = RR T , i.e. R este normală, şi cum este triunghiulară, este diagonală (v.<br />

problema 4.20). În continuare se procedează <strong>si</strong>milar. Din (S T S = SS T ) q+1,q+1 rezultă<br />

A q+1,j = 0, j = q+2 : p şi că blocul 2 × 2 A q+1,q+1 este normal. Având valori <strong>proprii</strong><br />

complexe, conform problemei 4.21, are structura din teoremă etc. Reciproca este imediată.<br />

P4.23 c) Fie B = 1 (A + 2 AH ) şi C = 1 (A − 2 AH ). Atunci A = B + C şi, conform<br />

punctelor a), b), B este hermitică iar C este antihermitică. Presupunem că avem şi A =<br />

= ˜B+ ˜C cu ˜B hermitică şi ˜C antihermitică. Atunci2B = A+A H = ˜B+ ˜C+ ˜B H + ˜C H = 2˜B.<br />

Deci B = ˜B. Analog, 2C = A−A H = 2˜C, i.e. C = ˜C. Deci descompunerea este unică.<br />

d) Se utilizează c) cu S = B şi T = −iC.<br />

P4.24 Se utilizează relaţiile din definiţii.<br />

P4.25 Se utilizează relaţiile din definiţii şi expre<strong>si</strong>ile părţilor hermitică şi antihermitică<br />

(v. soluţia problemei 4.23).<br />

P4.26 Se con<strong>si</strong>deră un set de n vectori liniar independenţi, e.g. ortogonali.<br />

P4.27 Fie P o matrice de permutare. Întrucât P este unitară, A este normală (hermitică,<br />

antihermitică, <strong>si</strong>metrică, anti<strong>si</strong>metrică) dacă şi numai dacă la fel este şi matricea<br />

C = P T AP. Putem alege P astfel încât B = C(1 : k,1 : k). a) Evident. b) Nu. De<br />

[ ] 6 3 1<br />

exemplu, matricea A = 1 3 3 este normală, dar B = A(1:2,1:2) nu este.<br />

3 −1 3<br />

P4.28 Conform teoremei 4.3, λ min(A) ≤ µ ≤ λ max(A).<br />

P4.29 Dacă B = Q H AQ ∈ IC p×p , unde Q H Q = I p, atunci aplicând matricei B teorema<br />

Courant-Fisher avem µ k = max dimV=k min x∈VS x H Q H AQx, unde V S este mulţimea<br />

<strong>vectorilor</strong> de normă euclidiană unitară din subspaţiul V ⊂ IC p . Acum, este uşor de constatat<br />

că Ṽ = {y ∈ IC n | y = Qx, x ∈ V} este un subspaţiu liniar al lui IC n , de aceeaşi<br />

dimen<strong>si</strong>une cu dimen<strong>si</strong>unea lui V (i.e. k) şi că mulţimea tuturor subspaţiilor Ṽ este<br />

numai o parte a mulţimii tuturor subspaţiilor de dimen<strong>si</strong>une k din IC n . Prin urmare,<br />

λ k = max dim Ṽ=k min x∈ṼS xH Ax ≥ µ k . Pentru cel de al doilea set de inegalităţi se utilizează<br />

cealaltă caracterizare minimax a <strong>valorilor</strong> <strong>proprii</strong> din teorema Courant-Fisher.

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

Saved successfully!

Ooh no, something went wrong!