Diagonalization, Jordan Form, and Linear Systems
mathdifferential-equationslinear-systemslecture
1Purpose of the Lecture
For the constant-coefficient system
\boldsymbol{x}'=A\boldsymbol{x},\qquad
\boldsymbol{x}(t_0)=\boldsymbol{x}_0,
Jordan form translates generalized eigenvector chains into solutions consisting of an exponential multiplied by a finite-degree polynomial. We work over \mathbb C. A real matrix may be complexified and placed in complex Jordan form; conjugate complex solutions then yield real solutions through their real and imaginary parts. Real Jordan form is outside the scope of this lecture.
data/lecture/math/linear-algebra/introduction-to-jordan-canonical-form.lecture.n.md
2Exponential of a Jordan Block
For a block of size k, write
J_k(\lambda)=\lambda I+N,\qquad
N=\begin{pmatrix}
0&1&0&\cdots&0\\
0&0&1&\cdots&0\\
\vdots&\vdots&\vdots&\ddots&\vdots\\
0&0&0&\cdots&1\\
0&0&0&\cdots&0
\end{pmatrix}.
The matrices \lambda I and N commute. Hence, with \tau=t-t_0,
e^{J_k(\lambda)\tau}
=e^{\lambda\tau}\sum_{r=0}^{k-1}\frac{\tau^r}{r!}N^r.
Here N^k=0 and, for k>1, N^{k-1}\ne0. Nilpotency terminates the series and produces polynomial factors through degree k-1. This factorization relies on commutativity and is not valid for arbitrary matrices X,Y.
3Return to the Original Coordinates
If
P^{-1}AP=J,
then
e^{A\tau}=Pe^{J\tau}P^{-1},\qquad
\boldsymbol{x}(t)=Pe^{J(t-t_0)}P^{-1}\boldsymbol{x}_0.
Each Jordan block may be exponentiated independently.
4Generalized Eigenvector Chains
Let
(A-\lambda I)v_1=0,\qquad
(A-\lambda I)v_j=v_{j-1}\quad(2\le j\le k)
be a Jordan chain. The solution with initial value v_j is
\boldsymbol{x}_j(t)=e^{\lambda\tau}
\sum_{\ell=0}^{j-1}\frac{\tau^\ell}{\ell!}v_{j-\ell}.
In particular,
\boldsymbol{x}_1=e^{\lambda\tau}v_1,\qquad
\boldsymbol{x}_2=e^{\lambda\tau}(v_2+\tau v_1),\qquad
\boldsymbol{x}_3=e^{\lambda\tau}\left(v_3+\tau v_2+\frac{\tau^2}{2}v_1\right).
Differentiation and the chain relations verify the system directly. The exponential derivative contributes the term with \lambda, while differentiating the polynomial lowers the chain index by one; (A-\lambda I)v_j=v_{j-1} then gives \boldsymbol{x}_j'=A\boldsymbol{x}_j.
5Example in Non-Jordan Coordinates
For
A=\begin{pmatrix}2&1\\-1&0\end{pmatrix},
\quad
v_1=\begin{pmatrix}1\\-1\end{pmatrix},
\quad
v_2=\begin{pmatrix}1\\0\end{pmatrix},
we have (A-I)v_1=0 and (A-I)v_2=v_1. Moreover,
\det(\mu I-A)=(\mu-1)^2,\qquad
P=\begin{pmatrix}1&1\\-1&0\end{pmatrix},\qquad
J=P^{-1}AP=\begin{pmatrix}1&1\\0&1\end{pmatrix}.
Consequently,
e^{A\tau}=e^\tau
\begin{pmatrix}
1+\tau&\tau\\
-\tau&1-\tau
\end{pmatrix}.
6Asymptotic Behavior
For \tau\ge0, a Jordan block satisfies
\|e^{J_k(\lambda)\tau}\|
\le C(1+\tau^{k-1})e^{\operatorname{Re}\lambda\,\tau}.
Thus, negative real part implies decay despite polynomial factors, while positive real part produces an exponentially growing direction. When the real part is zero, blocks of size one are bounded, whereas larger blocks admit polynomially unbounded solutions. Both eigenvalue real parts and Jordan block sizes are therefore required at the spectral boundary.
7Numerical Boundary
Jordan structure is sensitive to perturbations and is primarily a theoretical classification. For example,
\begin{pmatrix}\lambda&1\\\varepsilon&\lambda\end{pmatrix}
is one Jordan block when \varepsilon=0, whereas its eigenvalues split into \lambda\pm\sqrt{\varepsilon} when \varepsilon\ne0. Numerical computation of e^{At} should not generally assume that an exact Jordan block structure can be identified reliably in floating-point arithmetic.
8Previous Lecture
data/lecture/math/differential-equations/first-order-systems-and-matrix-exponentials.lecture.n.md
9Next Lecture
data/lecture/math/differential-equations/linearization-and-eigenvalue-stability.lecture.n.md
10Exercises
data/exercise/math/differential-equations/jordan-form-and-linear-systems.exercise.n.md