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.

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

astfel încât<br />

⎡<br />

D = X −1 TX = ⎢<br />

⎣<br />

⎤<br />

T 11 0 ··· 0<br />

0 T 22 ··· 0<br />

.<br />

. . ..<br />

⎥<br />

. ⎦ . (4.233)<br />

0 0 ··· T qq<br />

Demonstraţie. Dovadase obţine imediat prinaplicarearepetată alemei 4.5pentru<br />

a proba existenţa şi pentru a calcula submatricele X ij care definesc matricea de<br />

transformare X. Procedura are q −1 paşi.<br />

Pasul 1 ◦ . Fie partiţia<br />

T =<br />

⎡<br />

[ ]<br />

T11 ˜T12<br />

, unde ˜T12 = [ ] ⎢<br />

T 12 ··· T 1q , ˜T22 = ⎣<br />

0 ˜T22<br />

⎤<br />

T 22 ··· T 2q<br />

.<br />

. .. . 0 ··· T qq<br />

Din (4.231) rezultă λ(T 11 )∩λ(˜T 22 ) = ∅. Prin urmare, conform lemei 4.5, transformarea<br />

definită de T ← T (1) = X1 −1 TX 1 cu<br />

[ ]<br />

In1 ˜X12<br />

X 1 = ,<br />

0 I n−n1<br />

unde ˜X 12 este soluţia ecuaţiei Sylvester<br />

T 11 ˜X12 − ˜X 12˜T22 + ˜T 12 = 0<br />

a<strong>si</strong>gură anularea blocurilor extradiagonale de pe prima bloc-linie a matricei T.<br />

Pasul k ◦ . Presupunem că la primii k − 1 paşi am realizat bloc-diagonalizarea<br />

din primele bloc linii, i.e.<br />

T ← T (k−1) = X −1<br />

k−1···X−1 2 X−1 1 TX 1X 2···X k−1 =<br />

⎡<br />

⎣<br />

⎤<br />

˜T 11 0 0<br />

0 T kk<br />

˜Tk,k+1 ⎦,<br />

0 0 ˜Tk+1,k+1<br />

⎥<br />

⎦.<br />

unde<br />

˜T 11 =<br />

⎡<br />

⎢<br />

⎣<br />

⎤<br />

T 11 ··· 0<br />

.<br />

. .. .<br />

⎥<br />

⎦, ˜Tk,k+1 = [ ]<br />

T k,k+1 ··· T kq ,<br />

0 ··· T k−1,k−1<br />

˜T k+1,k+1 =<br />

⎡<br />

⎢<br />

⎣<br />

T k+1,k+1<br />

⎤<br />

··· T k+1,q<br />

.<br />

. .. .<br />

⎥<br />

⎦.<br />

0 ··· T qq<br />

Din aceleaşi motive ca la pasul 1 ◦ , dacă ˜X k,k+1 este soluţia ecuaţiei Sylvester<br />

T kk ˜Xk,k+1 − ˜X k,k+1 ˜Tk+1,k+1 + ˜T k,k+1 = 0, (4.234)

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

Saved successfully!

Ooh no, something went wrong!