Matrix Exponentials and Repeated Modes

Solve a constant system even when the matrix lacks an eigenvector basis.

Builds on Phase Portraits and Stability

The bigger question: How do coupled variables move together?

On this page

The idea

The scalar exponential generalizes to eAt=I+At+A2t2/2!+⋯e^{At}=I+At+A^2t^2/2!+\cdots. For a constant matrix, differentiation gives (eAt)′=AeAt(e^{At})'=Ae^{At} and eA0=Ie^{A0}=I. Thus x=eAtx0\mathbf x=e^{At}\mathbf x_0 solves any homogeneous constant system.

Visual guide

VISUAL GUIDEA stable system can have temporary component growth
For the defective matrix [[−1, 1], [0, −1]] and initial state (0, 1), x₁ = te⁻ᵗ first rises, then decays; x₂ = e⁻ᵗ decreases throughout. The polynomial factor comes from the repeated mode, not an unstable eigenvalue.001.250.2752.50.553.750.82551.1tstate
  • x₁ = te⁻ᵗ
  • x₂ = e⁻ᵗ
For the defective matrix [[−1, 1], [0, −1]] and initial state (0, 1), x₁ = te⁻ᵗ first rises, then decays; x₂ = e⁻ᵗ decreases throughout. The polynomial factor comes from the repeated mode, not an unstable eigenvalue.

Method and assumptions

If A=PDP−1A=PDP^{-1}, then eAt=PeDtP−1e^{At}=Pe^{Dt}P^{-1}. For a Jordan block A=λI+NA=\lambda I+N with nilpotent NN, factor eAt=eλteNte^{At}=e^{\lambda t}e^{Nt}; the nilpotent series terminates. In general, eA+B=eAeBe^{A+B}=e^Ae^B requires commuting matrices.

Worked example: a defective matrix

For A=(−110−1)A=\begin{pmatrix}-1&1\\0&-1\end{pmatrix}, write A=−I+NA=-I+N with N2=0N^2=0. Then eAt=e−t(1t01)e^{At}=e^{-t}\begin{pmatrix}1&t\\0&1\end{pmatrix}. Initial state (0,1)(0,1) yields (te−t,e−t)(te^{-t},e^{-t}), which cannot be produced by a single eigenvector alone.

PAUSE & THINKA quick check, not a grade

Try it yourself.

Choose an answer and explain your reasoning to yourself. Use a hint if you get stuck.

If N²=0, what is eᴺᵗ?

Hint 1 · Find a starting point

Expand the matrix exponential as a power series.

Hint 2 · Take the next step

Every power N² and higher is zero.

Show the reasoning

Answer: I+tN

Only I and tN survive. The identity is necessary to give eᴺ⁰=I.

Before moving on: what would make one of the other answers wrong? Saying why is part of understanding.

Worked example: rotation

For J=(0−110)J=\begin{pmatrix}0&-1\\1&0\end{pmatrix}, J2=−IJ^2=-I. Splitting the series into even and odd powers gives eJt=Icos⁡t+Jsin⁡te^{Jt}=I\cos t+J\sin t, the rotation matrix. State norm is preserved.

Interpreting the result

A negative repeated eigenvalue can produce a transient polynomial factor before eventual decay. A plot over a short interval may show growth in one component even though all states eventually approach zero.

Practice

  1. What is eA0e^{A0}?
  2. Compute eNte^{Nt} when N2=0N^2=0.
  3. Does te−tte^{-t} eventually decay?
Show worked solutions
  1. The identity matrix.
  2. I+tNI+tN, since every higher power vanishes.
  3. Yes. Its maximum for t≥0t\ge0 is at t=1t=1, after which it tends to zero.

Further study

MIT OpenCourseWare: Differential Equations provides a full university course with additional lectures and exercises.

MAKE IT YOURS

Pause before the next idea.

Can you explain how the visualization connects to this lesson’s goal? If a step still feels uncertain, put this lesson on your review list and try the check again another day.

Optional marks, not a grade. Saved in this browser only. Open notebook →