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.

Matrix Powers and Exponentials<br />

Y =<br />

145.0000 -207.0000 81.0000<br />

-207.0000 298.0000 -117.0000<br />

81.0000 -117.0000 46.0000<br />

Fractional powers, like A^(2/3), are also permitted; the results depend upon<br />

the distribution <strong>of</strong> the eigenvalues <strong>of</strong> the matrix.<br />

Element-by-Element Powers<br />

The .^ operator produces element-by-element powers. For example,<br />

X = A.^2<br />

A =<br />

1 1 1<br />

1 4 9<br />

1 9 36<br />

Exponentials<br />

The function<br />

sqrtm(A)<br />

computes A^(1/2) by a more accurate algorithm. The m in sqrtm distinguishes<br />

this function from sqrt(A) which, like A.^(1/2), does its job<br />

element-by-element.<br />

A system <strong>of</strong> linear, constant coefficient, ordinary differential equations can be<br />

written<br />

dx ⁄ dt = Ax<br />

where x = x(t) is a vector <strong>of</strong> functions <strong>of</strong> t and A is a matrix independent <strong>of</strong> t.<br />

The solution can be expressed in terms <strong>of</strong> the matrix exponential,<br />

xt () = e tA x( 0)<br />

The function<br />

expm(A)<br />

1-35

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

Saved successfully!

Ooh no, something went wrong!