MIT 18.701 — Lecture 8

Proving Jordan Normal Form, Cayley-Hamilton

The goal for today will be to prove the validity of Jordan Normal Form.

Remark. It only makes sense to discuss Jordan Normal Form if our linear operator T:V→VT: V \to V has a full set of eigenvalues to begin with. Thus, from here onward, we will assume the characteristic polynomial of TT splits into linear factors.

This can be achieved, say, by demanding that the field of VV is algebraically closed. So Jordan Normal Form does not necessarily exist for vector spaces over R\mathbb{R}, but it certainly exists for vector spaces over C\mathbb{C}.

Of course, we will also restrict our analysis to finite-dimensional VV.

§ Proving JNF: Assuming Zero is an Eigenvalue

Theorem. (Shifting JNFs) Consider some linear operator T:V→VT: V \to V and some λ∈F\lambda \in F. Then:

T has JNF [Jn1(λ1)⋯0⋮⋱⋮0⋯Jnk(λk)]   ⟺   T−λIdim⁡V has JNF [Jn1(λ1−λ)⋯0⋮⋱⋮0⋯Jnk(λk−λ)]T \text{ has JNF } \left[\begin{smallmatrix} J_{n_1}(\lambda_1) & \cdots & 0 \\ \vdots & \ddots & \vdots \\ 0 & \cdots & J_{n_k}(\lambda_k) \end{smallmatrix}\right] ~ \iff ~ T - \lambda I_{\dim V} \text{ has JNF } \left[\begin{smallmatrix} J_{n_1}(\lambda_1 - \lambda) & \cdots & 0 \\ \vdots & \ddots & \vdots \\ 0 & \cdots & J_{n_k}(\lambda_k - \lambda) \end{smallmatrix}\right]

Proof: No matter what basis we choose, λI\lambda I must look like a diagonal matrix with λ\lambda's along the diagonal. So no matter what basis we choose, the matrix T−λIT - \lambda I will just be TT with λ\lambda's subtracted along the diagonal.   ■\blacksquare

This tells us that to write T:V→VT: V \to V as a JNF matrix, it suffices to write T′:=T−λIdim⁡VT' := T - \lambda I_{\dim V} as a JNF matrix for any λ\lambda of our choice. We'll pick λ\lambda so that T′=T−λIdim⁡VT' = T - \lambda I_{\dim V} has zero as an eigenvalue.

§ Proving JNF: Factoring Out a Nilpotent Matrix

Consider some linear operator T:V→VT: V \to V, and assume without loss of generality that TT has zero as an eigenvalue. Now consider the infinite chain:

V⊇TV⊇T2V⊇⋯ .V \supseteq TV \supseteq T^2V \supseteq \cdots.

Since VV is finite-dimensional, this chain must eventually stabilize; that is, there must exist some kk for which:

TkV=Tk+1V=Tk+2V=⋯ .T^kV = T^{k + 1}V = T^{k + 2}V = \cdots.

Let U:=im⁡TkU := \im T^k be the subspace to which we stabilize, and also consider W:=ker⁡TkW := \ker T^k.

Theorem. (Vector Space Decomposition) Given U=im⁡TkU = \im T^k and W=ker⁡TkW = \ker T^k as defined above, we must have U⊕W≅VU \oplus W \cong V. Equivalently, dim⁡U+dim⁡W=dim⁡V\dim U + \dim W = \dim V and U∩W={0}U \cap W = \{0\}.

Proof: The former condition is easy; the Dimension Theorem gives it to us for free!

dim⁡U+dim⁡W=dim⁡(im⁡Tk)+dim⁡(ker⁡Tk)=dim⁡V.\dim U + \dim W = \dim (\im T^k) + \dim (\ker T^k) = \dim V.

It remains to show U∩W={0}U \cap W = \{0\}. Since U=TUU = TU, the restriction map T∣U:U→UT_{|U} : U \to U is surjective. But T∣UT_{|U} is a linear operator, so by the Dimension Theorem, T∣UT_{|U} must then be injective.

Now consider some v∈U∩Wv \in U \cap W. Since v∈W=ker⁡Tkv \in W = \ker T^k, it must be that Tkv=0T^kv = 0. But since v∈Uv \in U as well, and T∣U:U→UT_{|U}: U \to U is injective—hence so is T∣UkT^k_{|U}—the equation Tkv=Tk(0)=0T^kv = T^k(0) = 0 implies v=0v = 0, as desired.   ■\blacksquare

Thus, we may choose a basis of VV to be the union of a basis of UU and a basis of WW. Now, recalling that T∣U:U→UT_{|U}: U \to U is a bijective map and T∣Wk:W→WT^k_{|W}: W \to W annihilates WW, the matrix MTM_T from this choice of basis must look like:

MT=[MT∣U00MT∣W] where MT∣U is an invertible matrix and MT∣W is nilpotent.M_T = \begin{bmatrix} M_{T_{|U}} & 0 \\ 0 & M_{T_{|W}} \end{bmatrix} \text{ where } M_{T_{|U}} \text{ is an invertible matrix and } M_{T_{|W}} \text{ is nilpotent.}

Now suppose, hypothetically, we had the following tool:

Tool (Hammer). Using the hammer on a nilpotent linear operator T:V→VT: V \to V rewrites MTM_T as a JNF matrix by choosing an appropriate basis for VV.

Then given any linear operator T:V→VT: V \to V, we propose the following algorithm to write TT as a JNF matrix.

Algorithm for writing T:V→VT: V \to V as a JNF Matrix.

  1. Shift the not-necessarily-nilpotent T:V→VT: V \to V by λ\lambda so that T′=T−λIdim⁡VT' = T - \lambda I_{\dim V} has a zero eigenvalue.

  2. Decompose MT′M_{T'} as the direct sum of submatrices MT∣U′M_{T'_{|U}} and MT∣W′M_{T'_{|W}}.

  3. Apply the hammer to MT∣W′M_{T'_{|W}} so it becomes a JNF matrix.

  4. Zoom into MT∣U′M_{T'_{|U}}. Repeat this algorithm, starting from Step #1, for T∣U′:U→UT'_{|U}: U \to U, so that MT∣U′M_{T'_{|U}} is written as a JNF matrix.

  5. Now both MT∣U′M_{T'_{|U}} and MT∣W′M_{T'_{|W}} are JNF matrices. So T=T′+λIdim⁡VT = T' + \lambda I_{\dim V} is a JNF matrix, as desired.

The only way this algorithm could fail is if we “zoom in” indefinitely. We claim this is impossible; this is where our decision to shift T↦T′T \mapsto T' comes in! Because T′T' has a zero eigenvalue, WW contains the nonzero subspace ker⁡T′\ker T', so the linear operator T∣U′:U→UT'_{|U}: U \to U acts on a space of strictly smaller dimension than VV. But since VV is finite-dimensional, our “zooming in” must eventually yield U={0}U = \{0\}, at which point we're done!

Remark. Of course, the above algorithm can be phrased more formally using strong induction.

So for the remainder of this proof, it suffices to build a hammer.

§ Proving JNF: The Nilpotent Inductive Step

Recall that the hammer chooses a basis for a nilpotent linear operator T:V→VT: V \to V such that MTM_T looks like the direct sum of Jordan Blocks Jni(0)J_{n_i}(0). What can we say about the hammer's choice of basis?

Example. Suppose our nilpotent linear operator T:V→VT: V \to V, with our hammer's choice of basis {e1,e2,e3,e4,e5}\{e_1, e_2, e_3, e_4, e_5\}, appears in JNF as:

MT=[J3(0)00J2(0)]=[0100000100000000000100000]M_T = \begin{bmatrix} J_3(0) & 0 \\ 0 & J_2(0) \end{bmatrix} = \left[\begin{smallmatrix} 0 & 1 & 0 & 0 & 0 \\ 0 & 0 & 1 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 1 \\ 0 & 0 & 0 & 0 & 0 \end{smallmatrix}\right]

Then this choice of basis {e1,e2,e3,e4,e5}\{e_1, e_2, e_3, e_4, e_5\} has the following property:

e3↦Te2↦Te1↦T0    and    e5↦Te4↦T0.e_3 \stackrel{T}{\mapsto} e_2 \stackrel{T}{\mapsto} e_1 \stackrel{T}{\mapsto} 0 ~~~ \text{ and } ~~~ e_5 \stackrel{T}{\mapsto} e_4 \stackrel{T}{\mapsto} 0.

That's interesting—the hammer chooses a basis that forms chains descending down to 00 when operated on by TT. So we can rewrite the “docstring” for our hammer as follows:

Tool (Hammer, Revised). Using the hammer on a nilpotent linear operator T:V→VT: V \to V chooses a basis {{ei,j}j=1ki}i=1ℓ\{\{e_{i, j}\}_{j = 1}^{k_i}\}_{i = 1}^{\ell} of VV that forms ℓ\ell chains, the ii-th of length kik_i, as shown below:

{e1,k1↦T…↦Te1,2↦Te1,1↦T0e2,k2↦T…↦Te2,2↦Te2,1↦T0   ⋮eℓ,kℓ↦T…↦Teℓ,2↦Teℓ,1↦T0\begin{cases} e_{1, k_1} \stackrel{T}{\mapsto} \dots \stackrel{T}{\mapsto} e_{1, 2} \stackrel{T}{\mapsto} e_{1, 1} \stackrel{T}{\mapsto} 0 \\ e_{2, k_2} \stackrel{T}{\mapsto} \dots \stackrel{T}{\mapsto} e_{2, 2} \stackrel{T}{\mapsto} e_{2, 1} \stackrel{T}{\mapsto} 0 \\ ~~~\vdots \\ e_{\ell, k_{\ell}} \stackrel{T}{\mapsto} \dots \stackrel{T}{\mapsto} e_{\ell, 2} \stackrel{T}{\mapsto} e_{\ell, 1} \stackrel{T}{\mapsto} 0 \end{cases}

To build our hammer, it suffices to be able to choose a basis of any nilpotent T:V→VT: V \to V satisfying this chain-like structure. We'll prove this with a proof by induction on the dimension of VV.

Inductive Hypothesis. For all vector spaces WW with dim⁡W<dim⁡V\dim W < \dim V, the hammer exists for any nilpotent linear operator T:W→WT: W \to W.

For the inductive step, take W:=im⁡TW := \im T. Since TT is nilpotent, dim⁡W<dim⁡V\dim W < \dim V by the Dimension Theorem. Then by the inductive hypothesis on the linear operator T∣W:W→WT_{|W}: W \to W, the hammer gives us a basis {ei,j}\{e_{i, j}\} of WW with a chain-like structure with respect to T∣WT_{|W}.

We need to extend {ei,j}⊆W\{e_{i, j}\} \subseteq W to a basis of VV that preserves our chain-like structure. Here are two extensions:

Extending our basis by both Extension #1 and Extension #2 indeed preserves the chain-like structure:

{e1,k1↦T∣W…↦T∣We1,2↦T∣We1,1↦T∣W0   ⋮eℓ,kℓ↦T∣W…↦T∣Weℓ,2↦T∣Weℓ,1↦T∣W0 →extending basis    {e1,k1+1↦Te1,k1↦T…↦Te1,2↦Te1,1↦T0   ⋮eℓ,kℓ+1↦Teℓ,kℓ↦T…↦Teℓ,2↦Teℓ,1↦T0eℓ+1,1↦T0   ⋮edim⁡(ker⁡T),1↦T0\begin{cases} e_{1, k_1} \stackrel{T_{|W}}{\mapsto} \dots \stackrel{T_{|W}}{\mapsto} e_{1, 2} \stackrel{T_{|W}}{\mapsto} e_{1, 1} \stackrel{T_{|W}}{\mapsto} 0 \\ ~~~\vdots \\ e_{\ell, k_{\ell}} \stackrel{T_{|W}}{\mapsto} \dots \stackrel{T_{|W}}{\mapsto} e_{\ell, 2} \stackrel{T_{|W}}{\mapsto} e_{\ell, 1} \stackrel{T_{|W}}{\mapsto} 0 \end{cases} ~ \xrightarrow{\text{extending basis}} ~~~~ \begin{cases} e_{1, k_1 + 1} \stackrel{T}{\mapsto} e_{1, k_1} \stackrel{T}{\mapsto} \dots \stackrel{T}{\mapsto} e_{1, 2} \stackrel{T}{\mapsto} e_{1, 1} \stackrel{T}{\mapsto} 0 \\ ~~~\vdots \\ e_{\ell, k_{\ell} + 1} \stackrel{T}{\mapsto} e_{\ell, k_{\ell}} \stackrel{T}{\mapsto} \dots \stackrel{T}{\mapsto} e_{\ell, 2} \stackrel{T}{\mapsto} e_{\ell, 1} \stackrel{T}{\mapsto} 0 \\ e_{\ell + 1, 1} \stackrel{T}{\mapsto} 0 \\ ~~~ \vdots \\ e_{\dim(\ker T), 1} \stackrel{T}{\mapsto} 0 \end{cases}

There are now dim⁡(ker⁡T)\dim(\ker T) chains with lengths ki′k_i', where ki′=ki+1k_i' = k_i + 1 for all i=1,…,ℓi = 1, \dots, \ell, and ki′=1k_i' = 1 for all i=ℓ+1,…,dim⁡(ker⁡T)i = \ell + 1, \dots, \dim(\ker T). There's just one more thing left to check now: is our extended basis really a valid basis?

Thus, our extended basis really is a valid basis! So—at long last—we're done.

§ The Cayley-Hamilton Theorem

Theorem. (Cayley-Hamilton) Consider some M∈Cn×nM \in \mathbb{C}^{n \times n}. Consider a degree nn polynomial p:C→Cp: \mathbb{C} \to \mathbb{C} defined by:

p(λ):=det⁡(λIn−M)=λn+c1λn−1+⋯+cn−1λ+cn.p(\lambda) := \det(\lambda I_n - M) = \lambda^n + c_1 \lambda^{n - 1} + \dots + c_{n - 1} \lambda + c_n.

Then p(M):=Mn+c1Mn−1+⋯+cn−1M+cnInp(M) := M^n + c_1 M^{n - 1} + \dots + c_{n - 1} M + c_n I_n must equal 00.

Importantly, here is an incorrect proof of this theorem.

Incorrect Proof. By definition, p(M)=det⁡(MIn−M)=det⁡(0)=0p(M) = \det (M I_n - M) = \det (0) = 0. Done.   ■\blacksquare

This doesn't work, of course.

Remark. The permanent is a variant of the determinant defined by:

perm A:=∑σ∈Sn∏i=1nAi,σ(i).\mathrm{perm} \ A := \sum_{\sigma \in S_n} \prod_{i = 1}^n A_{i, \sigma(i)}.

In other words, the permanent is a determinant without the sign σ\mathrm{sign} \ \sigma factor. Then the Cayley-Hamilton theorem is not true if det⁡(λIn−M)\det(\lambda I_n - M) is replaced with perm(λIn−M)\mathrm{perm}(\lambda I_n - M).

A more correct proof is to just hit the theorem with our JNF sledgehammer: since the statement of our theorem is independent of basis, it suffices to prove the theorem for JNF matrices, at which point it's obvious.