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:VVT: 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:VVT: V \to V and some λF\lambda \in F. Then:

T has JNF [Jn1(λ1)00Jnk(λk)]      TλIdimV has JNF [Jn1(λ1λ)00Jnk(λ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:VVT: V \to V as a JNF matrix, it suffices to write T:=TλIdimVT' := T - \lambda I_{\dim V} as a JNF matrix for any λ\lambda of our choice. We'll pick λ\lambda so that T=TλIdimVT' = T - \lambda I_{\dim V} has zero as an eigenvalue.

§ Proving JNF: Factoring Out a Nilpotent Matrix

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

VTVT2V.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:=imTkU := \im T^k be the subspace to which we stabilize, and also consider W:=kerTkW := \ker T^k.

Theorem. (Vector Space Decomposition) Given U=imTkU = \im T^k and W=kerTkW = \ker T^k as defined above, we must have UWVU \oplus W \cong V. Equivalently, dimU+dimW=dimV\dim U + \dim W = \dim V and UW={0}U \cap W = \{0\}.

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

dimU+dimW=dim(imTk)+dim(kerTk)=dimV.\dim U + \dim W = \dim (\im T^k) + \dim (\ker T^k) = \dim V.

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

Now consider some vUWv \in U \cap W. Since vW=kerTkv \in W = \ker T^k, it must be that Tkv=0T^kv = 0. But since vUv \in U as well, and TU:UUT_{|U}: U \to U is injective—hence so is TUkT^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 TU:UUT_{|U}: U \to U is a bijective map and TWk:WWT^k_{|W}: W \to W annihilates WW, the matrix MTM_T from this choice of basis must look like:

MT=[MTU00MTW] where MTU is an invertible matrix and MTW 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:VVT: V \to V rewrites MTM_T as a JNF matrix by choosing an appropriate basis for VV.

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

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

  1. Shift the not-necessarily-nilpotent T:VVT: V \to V by λ\lambda so that T=TλIdimVT' = T - \lambda I_{\dim V} has a zero eigenvalue.

  2. Decompose MTM_{T'} as the direct sum of submatrices MTUM_{T'_{|U}} and MTWM_{T'_{|W}}.

  3. Apply the hammer to MTWM_{T'_{|W}} so it becomes a JNF matrix.

  4. Zoom into MTUM_{T'_{|U}}. Repeat this algorithm, starting from Step #1, for TU:UUT'_{|U}: U \to U, so that MTUM_{T'_{|U}} is written as a JNF matrix.

  5. Now both MTUM_{T'_{|U}} and MTWM_{T'_{|W}} are JNF matrices. So T=T+λIdimVT = 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 TTT \mapsto T' comes in! Because TT' has a zero eigenvalue, WW contains the nonzero subspace kerT\ker T', so the linear operator TU:UUT'_{|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:VVT: 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:VVT: 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:

e3Te2Te1T0    and    e5Te4T0.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:VVT: 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,k1TTe1,2Te1,1T0e2,k2TTe2,2Te2,1T0   e,kTTe,2Te,1T0\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:VVT: 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 dimW<dimV\dim W < \dim V, the hammer exists for any nilpotent linear operator T:WWT: W \to W.

For the inductive step, take W:=imTW := \im T. Since TT is nilpotent, dimW<dimV\dim W < \dim V by the Dimension Theorem. Then by the inductive hypothesis on the linear operator TW:WWT_{|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 TWT_{|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,k1TWTWe1,2TWe1,1TW0   e,kTWTWe,2TWe,1TW0 extending basis    {e1,k1+1Te1,k1TTe1,2Te1,1T0   e,k+1Te,kTTe,2Te,1T0e+1,1T0   edim(kerT),1T0\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(kerT)\dim(\ker T) chains with lengths kik_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(kerT)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 MCn×nM \in \mathbb{C}^{n \times n}. Consider a degree nn polynomial p:CCp: \mathbb{C} \to \mathbb{C} defined by:

p(λ):=det(λInM)=λn+c1λn1++cn1λ+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+c1Mn1++cn1M+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(MInM)=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:=σSni=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(λInM)\det(\lambda I_n - M) is replaced with perm(λInM)\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.