Matrix Exponentials: A, Eigenvalues, Jordan Form, Fundamental Matrix T

[Somefantastik]

A = \[ \left( \begin{array}{ccc}<br /> 0 &amp; 0 &amp; 1 \\<br /> 0 &amp; 1 &amp; 0 \\<br /> 1 &amp; 0 &amp; 0 \end{array} \right)\]

eigenvalues are \lambda_{1} = -1, \ \lambda_{2} = \lambda_{3} = 1

(A-\lambda_{1}I)u^{(1)} = 0 \ =&gt; \ u^{(1)} = \[ \left( \begin{array}{c}<br /> -1 \\<br /> 0 \\<br /> 1 \end{array} \right)\]

(A-\lambda_{2}I)u^{(2)} = 0 \ =&gt; u^{(1)} = \[ \left( \begin{array}{c}<br /> 1 \\<br /> 0 \\<br /> 1 \end{array} \right)\]

(A -\lambda_{3}I)u^{(3)} = u^{(2)} \ =&gt; \ <br /> \left(<br /> \begin{array}{ccc|c}<br /> 1&amp;0&amp;-1&amp;-1\\<br /> 0&amp;0 &amp;0&amp;1\\<br /> 0&amp;0&amp;0&amp;0<br /> \end{array}<br /> \right)<br />

since we cannot have 0 = 1, we can say that there is only one eigenvector for

\lambda = 1

which means that the Jordan form will be

\[ \left( \begin{array}{ccc}<br /> 1 &amp; 1 &amp; 0 \\<br /> 0 &amp; 1 &amp; 0 \\<br /> 0 &amp; 0 &amp; -1 \end{array} \right)\]

Am I correct?

Now I need to find exp(Jt) and I'm not sure how.

If I only have 2 eigenvectors, how can I find the fund. matrix T such that

e^{At} = Te^{Jt}T^{-1}?

 

[HallsofIvy]

In general, to find such a function of other than a "number", use the Taylor's series

e^A= I+ A+ (1/2)A²+ (1/3!)A³+ ...

If A happens to be diagonal, that's easy: Aⁿ is just the diagonal matrix with exponentials of the diagonal elements of A on its diagonal.

For Jordan form, its a little more complicated. Try calculating J², J³, ... for this particular J yourself and see if you spot a pattern.

 

[Somefantastik]

I noticed that A^2k = I and A^(2k+1)'s = A, but I don't know how to get e^A from that.

 

[Hurkyl]

I noticed that A^2k = I and A^(2k+1)'s = A, but I don't know how to get e^A from that.

How exactly are you having trouble summing the series?

 

[Somefantastik]

I'm just an idiot. I didn't realize that the Taylor expansions' limits came out to be cosh and sinh. I think I finished this problem. Thanks for the help so far.

 

[Hurkyl]

For the record, there's another method that can be useful. Note that you can separate that Jordan block into

<br /> \left[ \begin{array}{cc} 1 &amp; 1 \\ 0 &amp; 1 \end{array} \right]<br /> =<br /> \left[ \begin{array}{cc} 1 &amp; 0 \\ 0 &amp; 1 \end{array} \right]<br /> +<br /> \left[ \begin{array}{cc} 0 &amp; 1 \\ 0 &amp; 0 \end{array} \right]<br />

...

 

[Somefantastik]

Ah, I see, and the second matrix is nilpotent.

But still, if I got that, how would I compute

<br /> e^{At} = Te^{Jt}T^{-1} <br />

Not sure what T would be in this case.

 

[Hurkyl]

Right. And you'll find the results of doing arithmetic look a lot like derivatives -- a nilpotent matrix behaves a little bit like an infinitessimal.

(Again, look to the series to get started)