15.11.2014 Views

MATLAB Mathematics - SERC - Index of

MATLAB Mathematics - SERC - Index of

MATLAB Mathematics - SERC - Index of

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.

Eigenvalues<br />

The real part <strong>of</strong> each <strong>of</strong> the eigenvalues is negative, so e λt approaches zero as<br />

t increases. The nonzero imaginary part <strong>of</strong> two <strong>of</strong> the eigenvalues, ± ω ,<br />

contributes the oscillatory component, sin( ωt)<br />

, to the solution <strong>of</strong> the<br />

differential equation.<br />

With two output arguments, eig computes the eigenvectors and stores the<br />

eigenvalues in a diagonal matrix:<br />

[V,D] = eig(A)<br />

V =<br />

-0.8326 0.2003 - 0.1394i 0.2003 + 0.1394i<br />

-0.3553 -0.2110 - 0.6447i -0.2110 + 0.6447i<br />

-0.4248 -0.6930 -0.6930<br />

D =<br />

-3.0710 0 0<br />

0 -2.4645+17.6008i 0<br />

0 0 -2.4645-17.6008i<br />

The first eigenvector is real and the other two vectors are complex conjugates<br />

<strong>of</strong> each other. All three vectors are normalized to have Euclidean length,<br />

norm(v,2), equal to one.<br />

The matrix V*D*inv(V), which can be written more succinctly as V*D/V, is<br />

within round<strong>of</strong>f error <strong>of</strong> A. And, inv(V)*A*V, or V\A*V, is within round<strong>of</strong>f error<br />

<strong>of</strong> D.<br />

Defective Matrices<br />

Some matrices do not have an eigenvector decomposition. These matrices are<br />

defective, or not diagonalizable. For example,<br />

A = [ 6 12 19<br />

-9 -20 -33<br />

4 9 15 ]<br />

For this matrix<br />

[V,D] = eig(A)<br />

1-39

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

Saved successfully!

Ooh no, something went wrong!