Trouble with Matrix Exponentials

  • Thread starter Thread starter Blanchdog95
  • Start date Start date
  • Tags Tags
    Matrix
Click For Summary

Homework Help Overview

The discussion revolves around the computation of matrix exponentials, specifically focusing on the matrix A and its decomposition into diagonal and nilpotent components. Participants are exploring the conditions under which the exponential of a sum of matrices can be expressed as the product of their exponentials.

Discussion Character

  • Conceptual clarification, Mathematical reasoning, Assumption checking

Approaches and Questions Raised

  • The original poster attempts to apply the matrix exponential formula by separating matrix A into a diagonal matrix D and a nilpotent matrix N. They question the validity of this method after finding a discrepancy in the result compared to the expected answer.
  • Some participants point out the importance of the commutation relationship between matrices D and N, suggesting that the failure of the original method is due to the non-commutativity in this case.
  • Others suggest verifying the work through an alternative approach involving the matrix diagonalization and the use of the matrix P.

Discussion Status

The discussion is active, with participants providing insights into the conditions necessary for the matrix exponential properties to hold. There is acknowledgment of the original poster's confusion and a recognition of the importance of understanding commutation in this context. Guidance has been offered regarding alternative methods to check the work.

Contextual Notes

Participants note that the textbook provided examples where the method worked, raising questions about the specific conditions under which the method fails for matrix A. There is also mention of an upcoming test, indicating a time constraint for the original poster.

Blanchdog95
Messages
2
Reaction score
0
Homework Statement
Given the matrix A = {{1, 0}, {1, 0}}, find the matrix exponential e^(At)
Relevant Equations
e^(Mt) = I + Mt + (Mt)^2/2! + ... + (Mt)^k/k!
e^t = 1 + t + t^2/2! + t^3/3! + ... + t^k/k!
I've attempted to solve this by separating A into a diagonal matrix D and nilpotent matrix N:
D = {{1, 0}, {0, 0}}
N = {{0, 0}, {1, 0}}

e^(At) = e^((D + N)t) = e^(Dt) * e^(Nt)

When N is raised to the second power, it becomes the zero matrix. Therefore,
e^(Nt) = I + Nt = {{1, 0}, {t, 1}}

Note that D^2 = D
e^(Dt) = I + Dt + Dt^2/2! + ... + Dt^k/k! = {{1, 0}, {0, 1}} + {{1, 0}, {0, 0}}t + {{1, 0}, {0, 0}}t^2/2! + ... + {{1, 0}, {0, 0}}t^k/k!
e^(Dt) = {{1 + t + t^2/2! + ...+ t^k/k!, 0}, {0, 1}} = {{e^t, 0}, {0, 1}}

e^(Dt) * e^(Nt) = {{e^t, 0}, {0, 1}} * {{1, 0}, {t, 1}} = {{e^t, 0}, {t, 1}}

e^(At) = {{e^t, 0}, {t, 1}}

This answer is incorrect. The correct answer is

e^(At) = {{e^t, 0}, {e^t-1, 1}}, which can be obtained without separating A into D and N by noting that A^2 = A and substituting the series representation of e^t. I can see how this is done, but what I can't see is why the method I used failed, since the textbook itself used this method to find the matrix exponential of B = {{2, 1}, {0, 2}} by separating it into a diagonal and nilpotent matrix. If someone could teach me why this method does not work for A, or if it does, where I made a mistake, I would appreciate it a lot. Thanks in advance!
 
Physics news on Phys.org
Blanchdog95 said:
e^((D + N)t) = e^(Dt) * e^(Nt)
This is your problem. The statement ##e^{X+Y}= e^{X}e^{Y}## holds if XY=YX. Otherwise, it is not necessarily true. As you have seen, it is not true in your case. Notice that in the case where
$$
B=\begin{pmatrix}
2 & 1\\
0 & 2
\end{pmatrix}=N+D=
\begin{pmatrix}
2 & 0\\
0 & 2
\end{pmatrix}=
\begin{pmatrix}
0 & 1\\
0 & 0
\end{pmatrix}$$
N and D commute (as N=2I, and I commutes with anything), and thus you can use ##e^{X+Y}= e^{X}e^{Y}##. In your current case, however, N and D do not commute, and thus you cannot use ##e^{X+Y}= e^{X}e^{Y}##.
 
  • Like
Likes   Reactions: Blanchdog95
Isaac0427 said:
This is your problem. The statement ##e^{X+Y}= e^{X}e^{Y}## holds if XY=YX. Otherwise, it is not necessarily true. As you have seen, it is not true in your case. Notice that in the case where
$$
B=\begin{pmatrix}
2 & 1\\
0 & 2
\end{pmatrix}=N+D=
\begin{pmatrix}
0 & 1\\
0 & 0
\end{pmatrix}+
\begin{pmatrix}
2 & 0\\
0 & 2
\end{pmatrix}$$
N and D commute (as D=2I, and I commutes with anything), and thus you can use ##e^{X+Y}= e^{X}e^{Y}##. In your current case, however, N and D do not commute, and thus you cannot use ##e^{X+Y}= e^{X}e^{Y}##.
It won't let me edit the post as I made a few mistakes, so see the above quote for my revised post.
 
  • Like
Likes   Reactions: Blanchdog95
Isaac0427 said:
Much needed help.
THANK YOU SO MUCH! Now that I think about it the textbook did mention something about commutation; this probably saved me half a dozen points on my test next week.
 
Well, good luck on your test!

Just a quick thought-- you can easily check your work by noting that ##A=PDP^{-1}## where
$$D=
\begin{pmatrix}
1 & 0\\
0 & 1
\end{pmatrix}$$
and
$$P=
\begin{pmatrix}
1 & 0\\
1 & 1
\end{pmatrix}.$$
It is easy to compute ##e^{Dt}##, and then you can just use ##e^{At}=Pe^{Dt}P^{-1}##.
 

Similar threads

Replies
12
Views
2K
Replies
6
Views
3K
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 3 ·
Replies
3
Views
1K
  • · Replies 2 ·
Replies
2
Views
1K
  • · Replies 6 ·
Replies
6
Views
1K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
Replies
2
Views
2K