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→V has a full set of eigenvalues to begin with. Thus, from here onward, we will assume the characteristic polynomial of T splits into linear factors.
This can be achieved, say, by demanding that the field of V is algebraically closed. So Jordan Normal Form does not necessarily exist for vector spaces over R, but it certainly exists for vector spaces over C.
Of course, we will also restrict our analysis to finite-dimensional V.
§ Proving JNF: Assuming Zero is an Eigenvalue
Theorem. (Shifting JNFs) Consider some linear operator T:V→V and some λ∈F. Then:
T has JNF [Jn1(λ1)⋮0⋯⋱⋯0⋮Jnk(λk)]⟺T−λIdimV has JNF [Jn1(λ1−λ)⋮0⋯⋱⋯0⋮Jnk(λk−λ)]
Proof: No matter what basis we choose, λI must look like a diagonal matrix with λ's along the diagonal. So no matter what basis we choose, the matrix T−λI will just be T with λ's subtracted along the diagonal. ■
This tells us that to write T:V→V as a JNF matrix, it suffices to write T′:=T−λIdimV as a JNF matrix for any λ of our choice. We'll pick λ so that T′=T−λIdimV has zero as an eigenvalue.
§ Proving JNF: Factoring Out a Nilpotent Matrix
Consider some linear operator T:V→V, and assume without loss of generality that T has zero as an eigenvalue. Now consider the infinite chain:
V⊇TV⊇T2V⊇⋯.
Since V is finite-dimensional, this chain must eventually stabilize; that is, there must exist some k for which:
TkV=Tk+1V=Tk+2V=⋯.
Let U:=imTk be the subspace to which we stabilize, and also consider W:=kerTk.
Theorem. (Vector Space Decomposition) Given U=imTk and W=kerTk as defined above, we must have U⊕W≅V. Equivalently, dimU+dimW=dimV and U∩W={0}.
Proof: The former condition is easy; the Dimension Theorem gives it to us for free!
dimU+dimW=dim(imTk)+dim(kerTk)=dimV.
It remains to show U∩W={0}. Since U=TU, the restriction map T∣U:U→U is surjective. But T∣U is a linear operator, so by the Dimension Theorem, T∣U must then be injective.
Now consider some v∈U∩W. Since v∈W=kerTk, it must be that Tkv=0. But since v∈U as well, and T∣U:U→U is injective—hence so is T∣Uk—the equation Tkv=Tk(0)=0 implies v=0, as desired. ■
Thus, we may choose a basis of V to be the union of a basis of U and a basis of W. Now, recalling that T∣U:U→U is a bijective map and T∣Wk:W→W annihilates W, the matrix MT from this choice of basis must look like:
MT=[MT∣U00MT∣W] where MT∣U is an invertible matrix and MT∣W is nilpotent.
Now suppose, hypothetically, we had the following tool:
Tool (Hammer). Using the hammer on a nilpotent linear operator T:V→V rewrites MT as a JNF matrix by choosing an appropriate basis for V.
Then given any linear operator T:V→V, we propose the following algorithm to write T as a JNF matrix.
Algorithm for writing T:V→V as a JNF Matrix.
Shift the not-necessarily-nilpotent T:V→V by λ so that T′=T−λIdimV has a zero eigenvalue.
Decompose MT′ as the direct sum of submatrices MT∣U′ and MT∣W′.
Apply the hammer to MT∣W′ so it becomes a JNF matrix.
Zoom into MT∣U′. Repeat this algorithm, starting from Step #1, for T∣U′:U→U, so that MT∣U′ is written as a JNF matrix.
Now both MT∣U′ and MT∣W′ are JNF matrices. So T=T′+λIdimV 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′ comes in! Because T′ has a zero eigenvalue, W contains the nonzero subspace kerT′, so the linear operator T∣U′:U→U acts on a space of strictly smaller dimension than V. But since V is finite-dimensional, our “zooming in” must eventually yield 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→V such that MT looks like the direct sum of Jordan Blocks Jni(0). What can we say about the hammer's choice of basis?
Example. Suppose our nilpotent linear operator T:V→V, with our hammer's choice of basis {e1,e2,e3,e4,e5}, appears in JNF as:
Then this choice of basis {e1,e2,e3,e4,e5} has the following property:
e3↦Te2↦Te1↦T0 and e5↦Te4↦T0.
That's interesting—the hammer chooses a basis that forms chains descending down to 0 when operated on by T. So we can rewrite the “docstring” for our hammer as follows:
Tool (Hammer, Revised). Using the hammer on a nilpotent linear operator T:V→V chooses a basis {{ei,j}j=1ki}i=1ℓ of V that forms ℓ chains, the i-th of length ki, as shown below:
To build our hammer, it suffices to be able to choose a basis of any nilpotent T:V→V satisfying this chain-like structure. We'll prove this with a proof by induction on the dimension of V.
Inductive Hypothesis. For all vector spaces W with dimW<dimV, the hammer exists for any nilpotent linear operator T:W→W.
For the inductive step, take W:=imT. Since T is nilpotent, dimW<dimV by the Dimension Theorem. Then by the inductive hypothesis on the linear operator T∣W:W→W, the hammer gives us a basis {ei,j} of W with a chain-like structure with respect to T∣W.
We need to extend {ei,j}⊆W to a basis of V that preserves our chain-like structure. Here are two extensions:
Extension #1. Extend each individual chain by one: more precisely, for each chain i=1,…,ℓ, extend it by a vector called ei,ki+1∈V that satisfies ei,ki+1↦Tei,ki. Such an ei,ki+1∈V exists because W=imT.
Thus, we can extend {ei,j}⊆W by adding {ei,ki+1}i=1ℓ.
Extension #2. Recall the following fact:
Fact. For any vector space B, a subspace A⊆B, and any basis of A, we may extend the basis of A by adding dimB−dimA vectors to form a basis of B.
Using the above fact for B=kerT and A=W∩(kerT)=span({e1,1,…,eℓ,1}), there exist dim(kerT)−ℓ basis vectors {eℓ+1,1,…,edim(kerT),1} that extends the basis of W∩(kerT) to a basis of kerT.
Thus, we can extend {ei,j}⊆W by adding {ei,1}i=ℓ+1dim(kerT).
Extending our basis by both Extension #1 and Extension #2 indeed preserves the chain-like structure:
There are now dim(kerT) chains with lengths ki′, where ki′=ki+1 for all i=1,…,ℓ, and ki′=1 for all i=ℓ+1,…,dim(kerT). There's just one more thing left to check now: is our extended basis really a valid basis?
Linear Independence. Suppose there were some nontrivial linear dependence on our extended basis. Then the idea is to apply T to this linear dependence:
When we apply T, the linear dependence on our extended basis becomes a linear dependence on the original basis. But the original basis was properly chosen by our hammer to be a linearly independent basis of W, so every ci,j with j≥2 must vanish. All that survives is ∑ici,1ei,1=0, and {ei,1}i=1dim(kerT) is a basis of kerT by Extension #2, so those coefficients vanish too—there cannot exist such a linear dependence!
Spanning. We know our extended basis is linearly independent. Furthermore, the original basis contained dimW vectors, and our new basis extended the original basis by dim(kerT) vectors.
So our new basis contains dimW+dim(kerT)=dim(imT)+dim(kerT)=dimV vectors, and it's linearly independent, so it must be spanning!
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×n. Consider a degree n polynomial p:C→C defined by:
p(λ):=det(λIn−M)=λn+c1λn−1+⋯+cn−1λ+cn.
Then p(M):=Mn+c1Mn−1+⋯+cn−1M+cnIn must equal 0.
Importantly, here is an incorrect proof of this theorem.
Incorrect Proof. By definition, p(M)=det(MIn−M)=det(0)=0. Done. ■
This doesn't work, of course.
Remark. The permanent is a variant of the determinant defined by:
permA:=σ∈Sn∑i=1∏nAi,σ(i).
In other words, the permanent is a determinant without the signσ factor. Then the Cayley-Hamilton theorem is not true if det(λIn−M) is replaced with perm(λIn−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.