You can find the canonical form of Jordan for and a non singular matrix such that . Then, .
Dear MHF members,
my problem is explained below.
Problem. Compute the matrix exponential of the matrix
I can compute the matrix exponential by using the series expansion
However, I could not figure out how to compute it by using the eigenvalues.
I would be glad if you show me the solution in this way.
Yes, and so the matrix may not be diagonalizable. But the crucial part is not the number of eigenvalues but the number of independent eigenvectors. Here, 0 is a single eigenvalue and 1 is triple eigenvalue. An eigenvector corresponding to eigenvalue 0 (any multiple of <1, 0, 0, -1> is such an eigenvalue) is independent of any other eigenvalues but the eigenvalue 1 may have 3, 2, or only 1 independent eigenvector. In fact, looking for eigenvectors corresponding to eigenvalue 1 means looking for solutions to
So we must have x- z+ w= x, y+ z= y, z= z, and z= w y+ z= y gives immediately z= 0 so then w= 0. x- z+ w= x becomes x= x and y+ z= 7 becomes y= y both of which are satisfied by all x and y. That is, any eigenvector corresponding to eigenvalue 1 is of the form <x, y, 0, 0>= <x, 0, 0, 0>+ <0, y, 0, 0>= x<1, 0, 0, 0>+ y<0, 1, 0, 0>. That is, <1, 0, 0, 0> and <0, 1, 0, 0> are independent eigenvectors corresponding to eigenvalue 1. Because there are only 3 independent eigenvectors and not 4, we cannot "diagoalize" this matrix, but we can write it in "Jordan Normal form' as Fernando Revilla suggested. To do that, we need to find a "generalized eigenvector"- a vector, v, such that but . That, of course, is the same as saying which, in turn, means that is an eigenvector. That is, we must have
must be equal to either or .
The first of those is impossible but the second is satisfied by .
Constructing a matrix, P, having those vectors as columns, say,
then will be the Jordan Normal form.
(That's the idea. The calculations were done quickly and I don't guarentee them.)