20.07.2013 Views

Software matematic

Software matematic

Software matematic

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

96 Algebră liniară în MATLAB<br />

>> [V,D]=eig([2, -1; 1,0])<br />

V =<br />

0.7071 0.7071<br />

0.7071 0.7071<br />

D =<br />

1 0<br />

0 1<br />

Vectorii proprii sunt scalat¸i astfel ca norma lor euclidiană să fie egală cu unu (lucru<br />

posibil, căci dacă x este un vector propriu, atunci orice multiplu al său este de<br />

asemenea vector propriu).<br />

Dacă A este hermitiană, atunci MATLAB returnează valorile proprii sortate<br />

crescător s¸i matricea vectorilor proprii unitară (în limita preciziei de lucru):<br />

>> [V,D]=eig([2,-1;-1,1])<br />

V =<br />

-0.5257 -0.8507<br />

-0.8507 0.5257<br />

D =<br />

0.3820 0<br />

0 2.6180<br />

>> norm(V’*V-eye(2))<br />

ans =<br />

2.2204e-016<br />

În exemplul următor vom calcula valorile proprii ale matricei lui Frank (nehermitiană):<br />

>> F = gallery(’frank’,5)<br />

F =<br />

5 4 3 2 1<br />

4 4 3 2 1<br />

0 3 3 2 1<br />

0 0 2 2 1<br />

0 0 0 1 1<br />

>> e = eig(F)’<br />

e =<br />

10.0629 3.5566 1.0000 0.0994 0.2812<br />

Dacăλeste valoare proprie a matricei F, atunci1/λ este de asemenea valoare proprie:<br />

>> 1./e<br />

ans =<br />

0.0994 0.2812 1.0000 10.0629 3.5566<br />

Motivul este acela că polinomul caracteristic este anti-palindromic, adică termenii<br />

egal depărtat¸i de extrem sunt numere opuse:

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

Saved successfully!

Ooh no, something went wrong!