Mehmet Emre Tiryaki
Back to Notes

Jordan Canonical Form

Jordan Blocks

Jk(λ)=[λ10λ⋱⋱10λ]∈MkJ_k(\lambda) = \begin{bmatrix} \lambda & 1 & & 0 \\ & \lambda & \ddots & \\ & & \ddots & 1 \\ 0 & & & \lambda \end{bmatrix} \in M_k

has the single eigenvalue λ\lambda, and its eigenvectors are in Span{e1}\mathrm{Span}\{e_1\}.

Jk(0)J_k(0) shifts entries up by one, and its powers move the ones further up,

(Jk(0)L)ij={1,j−i=L0,else\big(J_k(0)^L\big)_{ij} = \begin{cases} 1, & j - i = L \\ 0, & \text{else} \end{cases}

so rank Jk(0)L=max⁡(k−L,0)\mathrm{rank}\, J_k(0)^L = \max(k - L, 0).

Jordan Matrices

A Jordan matrix is the direct sum of Jordan blocks,

J=[Jn1(λ1)0Jn2(λ2)0⋱]=⨁iJni(λi)J = \begin{bmatrix} J_{n_1}(\lambda_1) & & 0 \\ & J_{n_2}(\lambda_2) & \\ 0 & & \ddots \end{bmatrix} = \bigoplus_i J_{n_i}(\lambda_i)

Permuting the blocks of the direct sum gives a similar matrix, via PJPTPJP^T.

Jordan's Theorem For every A∈MnA \in M_n there is a Jordan matrix J∼AJ \sim A, unique up to rearrangement of blocks.

  • AA is diagonalizable   ⟺  \iff ni=1n_i = 1 for all blocks.
  • If A=SJS−1A = SJS^{-1}, rank(A)=rank(J)=∑rank(blocks of J)\mathrm{rank}(A) = \mathrm{rank}(J) = \sum \mathrm{rank}(\text{blocks of } J).

Dual Sequence

Any non-increasing sequence of natural numbers has a dual sequence, created by transposing its Ferrers diagram. The ℓ\ellth entry of the dual counts the entries that are ≥ℓ\ge \ell.

E.g. 5,3,3,2→4,4,3,1,15, 3, 3, 2 \to 4, 4, 3, 1, 1. The diagram has rows of lengths 5,3,3,25, 3, 3, 2, so its columns have lengths 4,4,3,1,14, 4, 3, 1, 1. Taking the dual twice gives back the original sequence.

Block Sizes from Ranks

First let A∼⨁iJni(0)A \sim \bigoplus_i J_{n_i}(0), e.g. J=J5(0)⊕J3(0)⊕J3(0)⊕J2(0)J = J_5(0) \oplus J_3(0) \oplus J_3(0) \oplus J_2(0). Since rank Aℓ=∑imax⁡(ni−ℓ,0)\mathrm{rank}\, A^\ell = \sum_i \max(n_i - \ell, 0),

sℓ:=rank Aℓ−1−rank Aℓ=#{i:ni≥ℓ}s_\ell := \mathrm{rank}\, A^{\ell - 1} - \mathrm{rank}\, A^\ell = \#\{i : n_i \ge \ell\}

which is the ℓ\ellth entry of the dual sequence of the block sizes. So the block sizes are the dual of the rank differences,

(t1,t2,…)=(s1,s2,…)∗  ⟹  A∼⨁iJti(0)(t_1, t_2, \ldots) = (s_1, s_2, \ldots)^* \implies A \sim \bigoplus_i J_{t_i}(0)

In general let A∼⨁i=1kJni(λi)A \sim \bigoplus_{i=1}^{k} J_{n_i}(\lambda_i) and γ∈C\gamma \in \mathbb{C}. Then

A−γI=SJS−1−γI=S(J−γI)S−1=S[⨁iJni(λi−γ)]S−1A - \gamma I = SJS^{-1} - \gamma I = S(J - \gamma I)S^{-1} = S\Big[\bigoplus_i J_{n_i}(\lambda_i - \gamma)\Big]S^{-1}

The blocks with λi≠γ\lambda_i \neq \gamma are invertible and keep full rank under powers, so they are dead weight in the rank differences. For γ=λj\gamma = \lambda_j only the blocks of λj\lambda_j lose rank, so

sℓ=rank(A−λjI)ℓ−1−rank(A−λjI)ℓs_\ell = \mathrm{rank}(A - \lambda_j I)^{\ell - 1} - \mathrm{rank}(A - \lambda_j I)^\ell

and taking the dual gives the Jordan block sizes of λj\lambda_j.

  • The first rank difference is s1=gemu(λj)s_1 = \mathrm{gemu}(\lambda_j), the number of Jordan blocks of λj\lambda_j.
  • The dual of the rank differences gives the Jordan blocks.
  • The size of the largest Jordan block gives the minimal polynomial.

AA is Similar to ATA^T

For all A∈MnA \in M_n, A∼ATA \sim A^T.

Proof

  1. pA(t)=det⁡(tI−A)=det⁡(tI−AT)=pAT(t)p_A(t) = \det(tI - A) = \det(tI - A^T) = p_{A^T}(t), so they have the same eigenvalues.
  2. rank(A−λI)ℓ=rank[(A−λI)ℓ]T=rank(AT−λI)ℓ\mathrm{rank}(A - \lambda I)^\ell = \mathrm{rank}\big[(A - \lambda I)^\ell\big]^T = \mathrm{rank}(A^T - \lambda I)^\ell for all λ∈σ(A)\lambda \in \sigma(A) and all ℓ\ell.
  3. So they have the same Jordan matrix.