syntax.garden

33 Jordan-形

33 Jordan-形

Jordan 形

这讲讲什么:给复向量空间上的每个算子挑一组基,把矩阵压到只剩对角线和紧挨其上的一条线——对角线填特征值,紧邻上方填 00 或 11。这就是 Jordan 形。
前置:32 讲(广义特征空间分解);31 讲(幂零算子、广义特征向量);19 讲(上三角化);20 讲(可对角化);08 讲的线性映射基本定理(第三节的维数估计要用一次)
一句话:先在 VV 上按特征值分堆,再把每一堆里的幂零部分拆成若干条「链」,每条链正好给一个 Jordan 块。


一、为什么需要它

1.1 三步走,每一步都留下一点没收拾干净的东西

20 讲已经说清:最理想的基是特征向量基,使 TT 的矩阵是对角阵。可惜不是每个算子都有这么多特征向量。

19 讲退一步:复向量空间上每个算子都有上三角矩阵,对角线上按重数排着全部特征值。代价是上三角部分可以塞任意数字。取

A=(634062007)A=\begin{pmatrix}6&3&4\\0&6&2\\0&0&7\end{pmatrix}

对角线上读出特征值 6,6,76,6,7。可右上角那三个数 3,4,23,4,2 有什么规律?没有。换一组基它们就变。

32 讲再进一步:VV 能按特征值拆成广义特征空间的直和

V=G(λ1,T)⊕⋯⊕G(λm,T),V=G(\lambda_1,T)\oplus\cdots\oplus G(\lambda_m,T),

每一块在 TT 下不变,而且 TT 在 G(λk,T)G(\lambda_k,T) 上等于 λkI\lambda_kI 加上一个幂零算子。于是矩阵成了分块对角,每块是 λk\lambda_k 在对角线的上三角。可块内还有一串 ∗* 没定。

这一讲要做的,就是把最后那串 ∗* 也收拾掉。 收拾的结果是:块内除了对角线上的 λ\lambda,只剩紧挨对角线上方那条线上的 00 和 11。剩下的 11 不是垃圾——它精确记录了「这个块有多长」。

1.2 三个层次,一句话各说一遍

同一个算子,换不同的基,能得到的矩阵整齐程度不同:

能保证的 靠什么 没收拾掉的东西
上三角 19 讲 上三角部分任意
分块对角,块内上三角、对角元相同 32 讲 块内上三角部分任意
分块对角,块内只有超对角线的 00 与 11 本讲 只剩块长与 λ\lambda

一句话:上三角是「把垃圾推到角落」,分块对角是「按特征值分堆」,Jordan 形是「每堆里只留下链长的记录」。

1.3 读完这一讲你应该能回答

  1. 一个 n×nn\times n 矩阵的 Jordan 形里,到底有几个块?块长怎么定?
  2. 为什么 λ\lambda 相同的两个块不能随随便便合并成一个?
  3. 每个复矩阵都有 Jordan 形——那它唯一吗?本节能证到哪一步?

二、定义与最小例子

2.1 先把幂零算子摆出来

31 讲定义过:TT 幂零,指存在正整数 mm 使 Tm=0T^m=0。使 Tm=0T^m=0 成立的最小 mm 叫 TT 的幂零指数(T=0T=0 时指数为 11)。

幂零算子的全部特征值都是 00:若 Tv=λvTv=\lambda v 且 v≠0v\neq0,则

0=Tmv=λmv,0=T^mv=\lambda^mv,

故 λ=0\lambda=0。所以下面几个例子里,矩阵对角线上全是 00。

2.2 例 A:一条链,一个块

在 C4\mathbb{C}^4 上定义

T(z1,z2,z3,z4)=(0, z1, z2, z3).T(z_1,z_2,z_3,z_4)=(0,\ z_1,\ z_2,\ z_3).

容易验 T4=0T^4=0(每作用一次,非零分量往右挪一格,挪出四格就没了),而 T3≠0T^3\neq0。取 v=(1,0,0,0)v=(1,0,0,0),算出

v=(1,0,0,0),Tv=(0,1,0,0),T2v=(0,0,1,0),T3v=(0,0,0,1).v=(1,0,0,0),\quad Tv=(0,1,0,0),\quad T^2v=(0,0,1,0),\quad T^3v=(0,0,0,1).

这四个向量分别是 e1,e2,e3,e4e_1,e_2,e_3,e_4,当然是 C4\mathbb{C}^4 的基。把它们倒过来排:

B=(T3v, T2v, Tv, v)=(e4,e3,e2,e1).\mathcal{B}=(T^3v,\ T^2v,\ Tv,\ v)=(e_4,e_3,e_2,e_1).

逐列算 TT 关于 B\mathcal{B} 的矩阵。第 jj 列是 TT 作用在第 jj 个基向量上的像,用 B\mathcal{B} 表示出来的坐标:

  • T(e4)=0T(e_4)=0,第 1 列是 (0,0,0,0)T(0,0,0,0)^{\mathsf T};
  • T(e3)=e4T(e_3)=e_4,第 2 列是 (1,0,0,0)T(1,0,0,0)^{\mathsf T};
  • T(e2)=e3T(e_2)=e_3,第 3 列是 (0,1,0,0)T(0,1,0,0)^{\mathsf T};
  • T(e1)=e2T(e_1)=e_2,第 4 列是 (0,0,1,0)T(0,0,1,0)^{\mathsf T}。

于是

M(T,B)=(0100001000010000).M(T,\mathcal{B})=\begin{pmatrix}0&1&0&0\\0&0&1&0\\0&0&0&1\\0&0&0&0\end{pmatrix}.

对角线上全 00,紧挨对角线上方四个位置全是 11,其余全 00。这就是一个 4×44\times4 的 Jordan 块。

为什么倒着排

因为 TT 把链上的向量往前推:v↦Tv↦T2v↦⋯v\mapsto Tv\mapsto T^2v\mapsto\cdots。
要让矩阵里出现「上一行的 1」,就得让「TT 把第 jj 个基向量送到第 j−1j-1 个」。
顺排 B′=(v,Tv,T2v,T3v)\mathcal{B}'=(v,Tv,T^2v,T^3v) 得到的矩阵,11 全在次对角线上——同一个算子,只是写法不合约定。

2.3 例 B:三条链,三个块

再看一个「一个循环向量不够用」的算子。在 C6\mathbb{C}^6 上定义

T(z1,z2,z3,z4,z5,z6)=(0, z1, z2, 0, z4, 0).T(z_1,z_2,z_3,z_4,z_5,z_6)=(0,\ z_1,\ z_2,\ 0,\ z_4,\ 0).

T3=0T^3=0。注意 z4,z6z_4,z_6 这两个坐标从来没被送到别处,所以这个算子处理不了 66 维——不存在 uu 使 u,Tu,…,T5uu,Tu,\dots,T^5u 成为基(那要求 T2u≠0T^2u\neq0 且这六个向量撑满空间,可 TT 的像只有 span⁡(e2,e3,e5)\operatorname{span}(e_2,e_3,e_5) 三个方向)。

改成分头找链。取

v1=e1,v2=e4,v3=e6.v_1=e_1,\qquad v_2=e_4,\qquad v_3=e_6.

三条链分别是

v1↦e2↦e3↦0,v2↦e5↦0,v3↦0.v_1\mapsto e_2\mapsto e_3\mapsto 0,\qquad v_2\mapsto e_5\mapsto 0,\qquad v_3\mapsto 0.

长度依次是 3,2,13,2,1。把每条链倒着排,依次拼起来:

B=(e3,e2,e1⏟3, e5,e4⏟2, e6⏟1).\mathcal{B}=(\underbrace{e_3,e_2,e_1}_{3},\ \underbrace{e_5,e_4}_{2},\ \underbrace{e_6}_{1}).

逐列验:T(e3)=0T(e_3)=0、T(e2)=e3T(e_2)=e_3、T(e1)=e2T(e_1)=e_2、T(e5)=0T(e_5)=0、T(e4)=e5T(e_4)=e_5、T(e6)=0T(e_6)=0。矩阵是

M(T,B)=(010000001000000000000010000000000000).M(T,\mathcal{B})=\begin{pmatrix}0&1&0&0&0&0\\0&0&1&0&0&0\\0&0&0&0&0&0\\0&0&0&0&1&0\\0&0&0&0&0&0\\0&0&0&0&0&0\end{pmatrix}.

它是一个 3×33\times3 块、一个 2×22\times2 块、一个 1×11\times1 块拼成的分块对角矩阵。

这个例子为什么关键

例 A 会让人误以为「nn 维幂零算子的 Jordan 形就是一个 n×nn\times n 块」。
例 B 一句话推翻它:同一个算子可以有好几条链,块可以不止一个。
链的长短由算子自己决定,选基的人没有自由。

2.4 一般情形的定义

现在把 λI\lambda I 加上去。

定义(Jordan 块) 设 λ∈F\lambda\in\mathbb{F}、d≥1d\ge1。d×dd\times d 的 Jordan 块指

Jd(λ)=(λ10λ1⋱⋱λ10λ),J_d(\lambda)=\begin{pmatrix}\lambda&1&&&0\\&\lambda&1&&\\&&\ddots&\ddots&\\&&&\lambda&1\\0&&&&\lambda\end{pmatrix},

即:对角线上全是 λ\lambda,紧挨对角线上方那条线(超对角线)上全是 11,其余全 00。

改写一下:

Jd(λ)=λId+Nd,J_d(\lambda)=\lambda I_d+N_d,

其中 NdN_d 就是把例 A 那个幂零块缩到 dd 维,满足 Ndd=0N_d^d=0、Ndd−1≠0N_d^{d-1}\neq0。

定义(Jordan 基、Jordan 形) 设 T∈L(V)T\in\mathcal{L}(V)。若 VV 有一个基,使 TT 关于该基的矩阵是分块对角矩阵

(A10⋱0Ap),Ak=Jdk(λk),\begin{pmatrix}A_1&&0\\&\ddots&\\0&&A_p\end{pmatrix},\qquad A_k=J_{d_k}(\lambda_k),

就说这个基是 TT 的 Jordan 基,这个矩阵是 TT 的 Jordan 形。

注意定义里没有任何要求说 λ1,…,λp\lambda_1,\dots,\lambda_p 互不相同。同一个 λ\lambda 可以对应好几个块(例 B 就是 λ=0\lambda=0 配三个块)。

2.5 回看定义,三个条件各管什么

  • 分块对角:各块作用在互不相交的基向量组上,块与块之间没有耦合。这是 32 讲直和分解的矩阵版本。
  • 每块对角元相同:一个块只属于一个 λ\lambda,也就是只落在一个 G(λ,T)G(\lambda,T) 里。不同特征值的向量不可能进同一个块。
  • 超对角线全是 11:这是规范化的结果。11 的含义是「TT 把这个基向量送向下一个,同时只加一个 λ\lambda」,也就是 (T−λI)(T-\lambda I) 沿着链往下推一格。换成 22 行不行?行,但那是没把基选好(第三节 3.4 会亲手把 33 改成 11)。

三、几个必须自己算一遍的结论

3.1 一条链的矩阵必然是超对角线上的一串 11

这一条第二节已经算过两遍,这里给出一般论证,因为它解释了「块」到底从哪来。

设 NN 幂零,u∈Vu\in V 满足 Nmu=0N^mu=0、Nm−1u≠0N^{m-1}u\neq0。则

u, Nu, N2u, …, Nm−1uu,\ Nu,\ N^2u,\ \dots,\ N^{m-1}u

线性无关。证明:若 ∑i=0m−1ciNiu=0\sum_{i=0}^{m-1}c_iN^iu=0,设 jj 是最小的使 cj≠0c_j\neq0 的下标,两边作用 Nm−1−jN^{m-1-j},得 cjNm−1u=0c_jN^{m-1}u=0,与 Nm−1u≠0N^{m-1}u\neq0 矛盾。

记 U=span⁡(u,Nu,…,Nm−1u)U=\operatorname{span}(u,Nu,\dots,N^{m-1}u),则 dim⁡U=m\dim U=m,且 UU 在 NN 下不变。取基

BU=(Nm−1u, …, Nu, u),\mathcal{B}_U=(N^{m-1}u,\ \dots,\ Nu,\ u),

即倒序。第 jj 个基向量是 Nm−juN^{m-j}u,而

N(Nm−ju)={Nm−j+1u,j≥2,0,j=1.N\big(N^{m-j}u\big)=\begin{cases}N^{m-j+1}u,&j\ge2,\\ 0,&j=1.\end{cases}

于是 j=1j=1 时第 1 列是零列(N(Nm−1u)=Nmu=0N(N^{m-1}u)=N^mu=0);j≥2j\ge2 时 Nm−j+1uN^{m-j+1}u 正是 BU\mathcal{B}_U 的第 j−1j-1 个向量,第 jj 列就是 ej−1e_{j-1}。矩阵恰是 Jm(0)J_m(0)。

这里没有任何「技巧」

全部内容就是一句:NN 把链上每个向量送到它的前一个,把最后一个送到 00。
一旦这么理解,块长、11 的位置、为什么倒序,都不必背。

3.2 每个幂零算子都有 Jordan 基:证明的三步骨架

这是本讲的第一个核心结论。

定理 A 设 T∈L(V)T\in\mathcal{L}(V) 幂零。则 VV 有一个基是 TT 的 Jordan 基。

对 dim⁡V\dim V 归纳。dim⁡V=1\dim V=1 时唯一的幂零算子是 00,随便取一个非零向量即可。

设 dim⁡V>1\dim V>1,且结论对所有维数更小的空间成立。令 mm 是 TT 的幂零指数,取 uu 使 Tm−1u≠0T^{m-1}u\neq0,令

U=span⁡(u,Tu,…,Tm−1u).U=\operatorname{span}(u,Tu,\dots,T^{m-1}u).

UU 在 TT 下不变,dim⁡U=m\dim U=m。若 U=VU=V,由 3.1 就得到 Jordan 基,证完。

难点在 U≠VU\neq V 的情形。思路:找到另一个不变子空间 WW,使 V=U⊕WV=U\oplus W;然后对 UU、WW 分别用归纳假设,两组基拼起来就是 VV 的 Jordan 基。

怎么找 WW?关键一步是用线性泛函切一刀。取 φ∈V′\varphi\in V' 使 φ(Tm−1u)≠0\varphi(T^{m-1}u)\neq0(存在:Tm−1u≠0T^{m-1}u\neq0,可以把它扩充成基,再用对偶基取泛函)。令

W={ v∈V: φ(Tkv)=0 对每个 k=0,1,…,m−1 }.W=\big\{\,v\in V:\ \varphi(T^kv)=0\ \text{对每个}\ k=0,1,\dots,m-1\ \big\}.

WW 在 TT 下不变:若 v∈Wv\in W,则对 k=0,…,m−2k=0,\dots,m-2 有 φ(Tk(Tv))=φ(Tk+1v)=0\varphi\big(T^k(Tv)\big)=\varphi(T^{k+1}v)=0;对 k=m−1k=m-1 有 φ(Tm−1(Tv))=φ(Tmv)=0\varphi\big(T^{m-1}(Tv)\big)=\varphi(T^mv)=0(因为 Tm=0T^m=0)。所以 Tv∈WTv\in W。

U∩W={0}U\cap W=\{0\}:设 v∈U∩Wv\in U\cap W 且 v≠0v\neq0。写

v=c0u+c1Tu+⋯+cm−1Tm−1u,v=c_0u+c_1Tu+\cdots+c_{m-1}T^{m-1}u,

令 jj 是最小的使 cj≠0c_j\neq0 的下标。两边作用 Tm−j−1T^{m-j-1}:

Tm−j−1v=cjTm−1u,T^{m-j-1}v=c_jT^{m-1}u,

因为次数更高的项都被 Tm=0T^m=0 杀掉了(Tm−j−1⋅Tiu=Tm−1+(i−j)uT^{m-j-1}\cdot T^iu=T^{m-1+(i-j)}u,当 i>ji>j 时次数 ≥m\ge m,为零)。再作用 φ\varphi:

φ(Tm−j−1v)=cj φ(Tm−1u)≠0,\varphi\big(T^{m-j-1}v\big)=c_j\,\varphi\big(T^{m-1}u\big)\neq0,

所以 v∉Wv\notin W,矛盾。故 U∩W={0}U\cap W=\{0\}。

U+W=VU+W=V:定义 S:V→FmS:V\to\mathbb{F}^m,

Sv=(φ(v), φ(Tv), …, φ(Tm−1v)).Sv=\big(\varphi(v),\ \varphi(Tv),\ \dots,\ \varphi(T^{m-1}v)\big).

SS 线性,null⁡S=W\operatorname{null}S=W。由 08 讲的线性映射基本定理,

dim⁡W=dim⁡null⁡S=dim⁡V−dim⁡range⁡S ≥ dim⁡V−m,\dim W=\dim\operatorname{null}S=\dim V-\dim\operatorname{range}S\ \ge\ \dim V-m,

因为 range⁡S⊆Fm\operatorname{range}S\subseteq\mathbb{F}^m,维数至多 mm。于是

dim⁡(U+W)=dim⁡U+dim⁡W ≥ m+(dim⁡V−m)=dim⁡V,\dim(U+W)=\dim U+\dim W\ \ge\ m+(\dim V-m)=\dim V,

故 U+W=VU+W=V;与 U∩W={0}U\cap W=\{0\} 合起来得 V=U⊕WV=U\oplus W。

UU 在 TT 下不变、dim⁡U=m<dim⁡V\dim U=m<\dim V,由归纳假设 UU 有 T∣UT|_U 的 Jordan 基;WW 在 TT 下不变、dim⁡W<dim⁡V\dim W<\dim V,同理也有。两组基合并就是 VV 的基,TT 关于它的矩阵是两块分块对角,每块都是 Jordan 块。证毕。

这一步的「为什么」

归纳法的难点从来不是「怎么归纳」,而是「补空间怎么造」。
这里的构造值得记住:用泛函把链的最高次项「测」出来,再把测不出东西的向量收进 WW。
WW 之所以够大,靠的是「SS 的值域最多 mm 维」这个上界。

3.3 从幂零推广到一般

定理 B(存在性) 设 F=C\mathbb{F}=\mathbb{C}、T∈L(V)T\in\mathcal{L}(V)。则 VV 有 TT 的 Jordan 基。

由 32 讲的广义特征空间分解,

V=G(λ1,T)⊕⋯⊕G(λm,T),V=G(\lambda_1,T)\oplus\cdots\oplus G(\lambda_m,T),

每块在 TT 下不变,且每个 (T−λkI)∣G(λk,T)(T-\lambda_kI)|_{G(\lambda_k,T)} 幂零。在 G(λk,T)G(\lambda_k,T) 上写成

T∣G(λk,T)=λkI+Nk,T\big|_{G(\lambda_k,T)}=\lambda_kI+N_k,

由定理 A,存在 G(λk,T)G(\lambda_k,T) 的基使 NkN_k 的矩阵是若干块 Jd(0)J_d(0)。同一组基下 λkI\lambda_kI 的矩阵就是 λkI\lambda_kI,两者相加(同一个基,逐项相加)把每个块的 λ\lambda 补上:

Jd(0)+λkI=Jd(λk).J_d(0)+\lambda_kI=J_d(\lambda_k).

把各 G(λk,T)G(\lambda_k,T) 的基合并,就得到 VV 的 Jordan 基。

注意定理 B 里 F=C\mathbb{F}=\mathbb{C} 是硬条件。 定理 A 对实空间也成立(证明只用了 Tm=0T^m=0、线性泛函、维数估计),但定理 B 不行:R\mathbb{R} 上不保证特征值存在,广义特征空间分解就无从谈起。R2\mathbb{R}^2 上转 90°90° 就是反例(见第五节误区 1)。

3.4 亲手把一个 33 改成 11

回到第一节那个矩阵

A=(634062007),A=\begin{pmatrix}6&3&4\\0&6&2\\0&0&7\end{pmatrix},

它定义的算子是

T(z1,z2,z3)=(6z1+3z2+4z3, 6z2+2z3, 7z3),T(z_1,z_2,z_3)=(6z_1+3z_2+4z_3,\ 6z_2+2z_3,\ 7z_3),

特征值 6,6,76,6,7。

先看 G(7,T)G(7,T)。算一下

T(10,2,1)=(70,14,7)=7 (10,2,1),T(10,2,1)=(70,14,7)=7\,(10,2,1),

而 77 的重数是 11,所以 G(7,T)=span⁡(10,2,1)G(7,T)=\operatorname{span}(10,2,1)。

再看 G(6,T)G(6,T)。令 z3=0z_3=0:

T(z1,z2,0)=(6z1+3z2, 6z2, 0),T(z_1,z_2,0)=(6z_1+3z_2,\ 6z_2,\ 0),

可见 {(z1,z2,0)}\{(z_1,z_2,0)\} 在 TT 下不变、维数 22,就是 G(6,T)G(6,T)。

在这个二维块上,TT 关于 (e1,e2)(e_1,e_2) 的矩阵是

(6306),\begin{pmatrix}6&3\\0&6\end{pmatrix},

也就是

T(e1)=6e1,T(e2)=3e1+6e2.T(e_1)=6e_1,\qquad T(e_2)=3e_1+6e_2.

令 N=T−6IN=T-6I,则 N(e1)=0N(e_1)=0、N(e2)=3e1N(e_2)=3e_1、N2=0N^2=0。要按 3.1 的套路造链,不能从 e1e_1 起(它已经在 NN 的零空间里,链长只有 11),要从 NN 作用下非零的向量起。取 u=e2u=e_2:

Nu=3e1≠0,N2u=0.Nu=3e_1\neq0,\qquad N^2u=0.

链是 e2↦3e1↦0e_2\mapsto3e_1\mapsto0,长 22。倒着排得基

(Nu, u)=(3e1, e2)=((3,0,0), (0,1,0)).\big(Nu,\ u\big)=\big(3e_1,\ e_2\big)=\big((3,0,0),\ (0,1,0)\big).

验算:

T(3e1)=18e1=6⋅(3e1),T(e2)=3e1+6e2=(3e1)+6 e2.T(3e_1)=18e_1=6\cdot(3e_1),\qquad T(e_2)=3e_1+6e_2=(3e_1)+6\,e_2.

关于基 ((3,0,0), (0,1,0), (10,2,1))\big((3,0,0),\ (0,1,0),\ (10,2,1)\big),

M(T)=(610060007)=diag⁡(J2(6), J1(7)).M(T)=\begin{pmatrix}6&1&0\\0&6&0\\0&0&7\end{pmatrix}=\operatorname{diag}\big(J_2(6),\ J_1(7)\big).

那个讨厌的 33 变成了 11。它不是被算掉的,是被基 3e13e_1 吃掉的——把第一个基向量放缩 33 倍,非对角元就跟着除以 33。

由这个例子能抽出的一般规律

块内超对角线上的非零数字都只是标度,换基就能改成 11。真正消不掉的是能不能连成链:
若某个位置上的耦合是 00,那就是块与块的接缝,再怎么换基也不会变出 11 来。
所以 Jordan 形要回答的问题不是「数字是多少」,而是「块怎么切」。

3.5 存在之后:块长分布是不是被算子唯一决定

必须坦白说清楚:上面证的是「存在」,不是「唯一」。

一个算子的 Jordan 形当然不可能字面上唯一——把块换个顺序排,矩阵变了,但仍然分块对角、每块仍是 Jordan 块,还是同一个算子关于另一组基的矩阵。

真正的问题是:块长分布是不是被算子唯一决定? 也就是「66 对应一个 2×22\times2 块加一个 1×11\times1 块」这件事,会不会因为换基变成「三个 1×11\times1 块」?

答案是「不会」(否则这套理论就没意义)。要证这一点,得再加一段关于 dim⁡null⁡(T−λI)k\dim\operatorname{null}(T-\lambda I)^k 的分析——下面 3.5 末就把它补完。

不过,从已有的证明里能免费读出三个必然不变的量:

量 等于什么 为什么
λ\lambda 的块数 dim⁡E(λ,T)\dim E(\lambda,T)(几何重数) 每条链的末端都落在 null⁡(T−λI)\operatorname{null}(T-\lambda I) 里,且这些末端恰好构成它的基
该 λ\lambda 所有块长之和 dim⁡G(λ,T)\dim G(\lambda,T)(代数重数) 这些链上的向量全体构成 G(λ,T)G(\lambda,T) 的基
最长块的长度 极小多项式中 (z−λ)(z-\lambda) 的指数 最长链要走到第 dd 步才归零,dd 就是使 (T−λI)d(T-\lambda I)^d 在 GG 上为零的最小 dd

在幂零情形下,这三条退化成:块数 =dim⁡null⁡T=\dim\operatorname{null}T;块长之和 =dim⁡V=\dim V;最长块长 == 幂零指数。

这三条不足以唯一确定块长分布

举例:dim⁡V=8\dim V=8、块数 33、最长块 44。
分布可以是 4+3+14+3+1,也可以是 4+2+24+2+2,两者都满足上面三条。
它们确实是不同的算子——一个可亲手验证的区分是 dim⁡null⁡T2\dim\operatorname{null}T^2:
对 4+3+14+3+1 它是 2+2+1=52+2+1=5,对 4+2+24+2+2 它是 2+2+2=62+2+2=6。
所以「块数 + 总长 + 最长块」只是必要条件,不是完整的不变量。

补全:唯一性的完整判据

上面那个反例反而指了路——要区分块长分布,得看更高次的零空间。事实上只要把所有 kk 次的零空间维数都看一遍,分布就被完全钉死了。

记 N=T−λIN=T-\lambda I(把它限制在 G(λ,T)G(\lambda,T) 上看)。设长度为 ss 的链(即对应 s×ss\times s 的 Jordan 块)有 csc_s 条。核心观察是:

dim⁡null⁡Nk=∑s≥1min⁡(k,s) cs\dim\operatorname{null}N^k=\sum_{s\ge 1}\min(k,s)\,c_s

为什么:一条长度为 ss 的链上,NkN^k 把它打掉前 min⁡(k,s)\min(k,s) 个向量(链走 ss 步归零,走到第 kk 步就归零了 kk 个)。各条链独立,直接相加。

于是相邻两次的差

dim⁡null⁡Nk−dim⁡null⁡Nk−1=#{s≥k}=「长度≥k 的链有多少条」\dim\operatorname{null}N^{k}-\dim\operatorname{null}N^{k-1}=\#\{s\ge k\}=\text{「长度}\ge k\text{ 的链有多少条」}

这就把分布完全确定了:设 bk:=#{s≥k}b_k:=\#\{s\ge k\},则长度恰好为 kk 的链数为 ck=bk−bk+1c_k=b_k-b_{k+1}。bkb_k 全部由算子的零空间维数算出,与基的选择无关。所以块长分布是算子的不变量。

用反例验一遍

回到 4+3+14+3+1 与 4+2+24+2+2(dim⁡V=8\dim V=8):

kk bk=#{s≥k}b_k=\#\{s\ge k\}(4+3+1) dim⁡null⁡Nk\dim\operatorname{null}N^k bkb_k(4+2+2) dim⁡null⁡Nk\dim\operatorname{null}N^k
1 3 3 3 3
2 2 5 3 6
3 2 7 1 7
4 1 8 1 8

两列在 k=2k=2 处首次分岔(55 vs 66),正是前面提到的 dim⁡null⁡T2\dim\operatorname{null}T^2。
而把每行的差算出来(bk−bk+1b_k-b_{k+1}),分别得 c1,c2,c3,c4c_1,c_2,c_3,c_4:

  • 4+3+14+3+1:c=(1,0,1,1)c=(1,0,1,1) → 链长 4,3,14,3,1 ✓
  • 4+2+24+2+2:c=(0,2,0,1)c=(0,2,0,1) → 链长 4,2,24,2,2 ✓

分布因此被唯一还原。

唯一性的准确陈述

在不计块的排列顺序的意义下,Jordan 形由算子唯一决定。
换个说法:两个算子相似   ⟺  \iff 它们有相同的 Jordan 形(每个特征值的块长分布都相同)。
这条判据把「相似」这个看似无从下手的等价关系,变成了一次可以算完的核对。

注意:本讲前面(3.2/3.3)证的存在性、加上这里证的唯一性,合起来才是完整的 Jordan 形定理。本节之前只做到存在性,现在补齐了。

术语:初等因子

有些教材不用「块长分布」这个说法,而说「初等因子」:特征值 λ\lambda 配上一个长度 ss 的块,就记一个因子 (z−λ)s(z-\lambda)^s。
于是「Jordan 形唯一」等价于「初等因子组唯一」。两种说法是同一件事的两种记账方式。
下面把这条记账法做成可操作的工具。

3.6 初等因子:从特征矩阵直接算出来

本节约定 F=C\mathbb{F}=\mathbb{C}(或至少特征值全在 F\mathbb{F} 内)。这样才能把初等因子写成 (λ−λ0)s(\lambda-\lambda_0)^s 并对应到 Jordan 块;一般域上的情形见 3.7。

前面判断相似,要逐个 kk 算 dim⁡null⁡Nk\dim\operatorname{null}N^k。实际做题时有一条纯代数的路:不必碰基,只算多项式。

特征矩阵。设 A∈Fn,nA\in\mathbb{F}^{n,n},把 λ\lambda 当不定元,AA 的特征矩阵是 λI−A\lambda I-A(元素都是 F[λ]\mathbb{F}[\lambda] 里的多项式)。

行列式因子。对 k=1,…,nk=1,\dots,n,定义 Dk(λ)D_k(\lambda) 为 λI−A\lambda I-A 的全部 kk 阶子式的最大公因式(首一化)。约定 D0=1D_0=1。

不变因子。令

dk(λ)=Dk(λ)Dk−1(λ)(k=1,…,n)d_k(\lambda)=\frac{D_k(\lambda)}{D_{k-1}(\lambda)}\qquad(k=1,\dots,n)

则 d1∣d2∣⋯∣dnd_1\mid d_2\mid\cdots\mid d_n,且 dnd_n 就是极小多项式,而 ∏k=1ndk=det⁡(λI−A)\prod_{k=1}^{n}d_k=\det(\lambda I-A) 是特征多项式。

[!note] 这三条不是「显然」
整除链 d1∣⋯∣dnd_1\mid\cdots\mid d_n、dn=d_n= 极小多项式、∏dk=det⁡\prod d_k=\det——这三条都是 Smith 标准形定理的内容,
严格证明要用 λ\lambda-矩阵的初等变换理论。本节把它当已知结论用,只演示「怎么算」。

初等因子。把每个 dkd_k 分解成不可约因式的幂,把所有非平凡的幂次因子收集起来——这就是初等因子组。

为什么这条路线值得记

它把「判断两个矩阵是否相似」变成了纯多项式计算:

A∼B  ⟺  λI−A 与 λI−B 有相同的初等因子组A\sim B\iff \lambda I-A\ \text{与}\ \lambda I-B\ \text{有相同的初等因子组}

更简洁的版本:λI−A\lambda I-A 与 λI−B\lambda I-B 等价(即能通过初等行/列变换互变)。
这条路不需要求特征向量、不需要算零空间维数,只要会算行列式与最大公因式。

初等因子 ↔ Jordan 块,对应关系很简单:

初等因子 Jordan 块
(λ−λ0)s(\lambda-\lambda_0)^s 一个 s×ss\times s 的 Js(λ0)J_s(\lambda_0)
λ−λ0\lambda-\lambda_0(即 s=1s=1) 一个 1×11\times1 块

一个完整的手算例。设

A=(210021002)A=\begin{pmatrix}2&1&0\\0&2&1\\0&0&2\end{pmatrix}

第一步,写特征矩阵:

λI−A=(λ−2−100λ−2−100λ−2)\lambda I-A=\begin{pmatrix}\lambda-2&-1&0\\0&\lambda-2&-1\\0&0&\lambda-2\end{pmatrix}

第二步,算行列式因子。D3=det⁡(λI−A)=(λ−2)3D_3=\det(\lambda I-A)=(\lambda-2)^3。D2D_2 是全部二阶子式的最大公因式。3×33\times3 矩阵的二阶子式有 (32)(32)=9\binom32\binom32=9 个,但只需抽出三个就能定案(三个都取自第 1,21,2 行):

∣λ−2−10λ−2∣=(λ−2)2  (第1,2列),∣λ−200−1∣=−(λ−2)  (第1,3列),∣−10λ−2−1∣=1  (第2,3列).\left|\begin{matrix}\lambda-2&-1\\0&\lambda-2\end{matrix}\right|=(\lambda-2)^2\ \ (第1,2列),\qquad \left|\begin{matrix}\lambda-2&0\\0&-1\end{matrix}\right|=-(\lambda-2)\ \ (第1,3列),\qquad \left|\begin{matrix}-1&0\\\lambda-2&-1\end{matrix}\right|=1\ \ (第2,3列).

最大公因式是 11(因为出现了 11),所以 D2=1D_2=1;同理 D1=1D_1=1(有元素 11)。

第三步,不变因子:d1=D1/D0=1d_1=D_1/D_0=1,d2=D2/D1=1d_2=D_2/D_1=1,d3=D3/D2=(λ−2)3d_3=D_3/D_2=(\lambda-2)^3。

第四步,初等因子组:只有 {(λ−2)3}\{(\lambda-2)^3\} 一个。

结论:AA 的 Jordan 形就是一个 3×33\times3 块 J3(2)J_3(2)——它本身就已是 Jordan 形,与自己相似,不稀奇。

换一个更有意思的例子。设

B=(210020002)B=\begin{pmatrix}2&1&0\\0&2&0\\0&0&2\end{pmatrix}

D3=(λ−2)3D_3=(\lambda-2)^3 不变。D2D_2 是全部二阶子式的最大公因式。3×33\times3 矩阵有 (32)(32)=9\binom32\binom32=9 个二阶子式,逐个枚举后,非零的只有四个:

取的行 取的列 子式 值
1,21,2 1,21,2 ∣λ−2−10λ−2∣\begin{vmatrix}\lambda-2&-1\\0&\lambda-2\end{vmatrix} (λ−2)2(\lambda-2)^2
1,31,3 1,31,3 ∣λ−200λ−2∣\begin{vmatrix}\lambda-2&0\\0&\lambda-2\end{vmatrix} (λ−2)2(\lambda-2)^2
2,32,3 2,32,3 ∣λ−200λ−2∣\begin{vmatrix}\lambda-2&0\\0&\lambda-2\end{vmatrix} (λ−2)2(\lambda-2)^2
1,31,3 2,32,3 ∣−100λ−2∣\begin{vmatrix}-1&0\\0&\lambda-2\end{vmatrix} −(λ−2)-(\lambda-2)

其余 55 个全为 00(例如第 1,21,2 行第 1,31,3 列 ∣λ−2000∣=0\begin{vmatrix}\lambda-2&0\\0&0\end{vmatrix}=0)。

于是非零项恰是 (λ−2)2(\lambda-2)^2(三个)与 −(λ−2)-(\lambda-2)(一个),取公因式:

D2=λ−2.D_2=\lambda-2 .

这里最容易错的一步是:因为出现了 −(λ−2)-(\lambda-2) 这个一次因子,而不是常数,所以公因式是 λ−2\lambda-2,不是 11。正是这个 D2D_2 把 BB 与上面的 AA 区分开了——AA 的 D2D_2 是 11。

于是

d1=1,d2=λ−21=λ−2,d3=(λ−2)3λ−2=(λ−2)2.d_1=1,\quad d_2=\frac{\lambda-2}{1}=\lambda-2,\quad d_3=\frac{(\lambda-2)^3}{\lambda-2}=(\lambda-2)^2 .

初等因子组:{(λ−2), (λ−2)2}\{(\lambda-2),\ (\lambda-2)^2\}——两个块,尺寸 11 与 22。

与 AA 对比:AA 与 BB 的特征多项式都是 (λ−2)3(\lambda-2)^3、特征值都是 22,但两者不相似(AA 是一个 3×33\times3 块,BB 是 2+12+1 两个块)。注意几何重数已经不同:dim⁡null⁡(A−2I)=1\dim\operatorname{null}(A-2I)=1,而 dim⁡null⁡(B−2I)=2\dim\operatorname{null}(B-2I)=2(因为 rank⁡(B−2I)=1\operatorname{rank}(B-2I)=1)——这一对其实用几何重数就能分开。差别也体现在极小多项式:AA 的是 (λ−2)3(\lambda-2)^3,BB 的是 (λ−2)2(\lambda-2)^2;以及 D2D_2:AA 的 D2=1D_2=1,BB 的 D2=λ−2D_2=\lambda-2。

光靠极小多项式其实也不够

上面这一对用极小多项式分开了,容易让人以为「特征多项式 ++ 极小多项式」就够。
反例见 3.8:J2(2)⊕J2(2)J_2(2)\oplus J_2(2) 与 J2(2)⊕J1(2)⊕J1(2)J_2(2)\oplus J_1(2)\oplus J_1(2) 的特征多项式、极小多项式、最长块长全都一样,仍不相似。
差别只在更有量的东西上:前者一个 22 阶块,后者一个 22 阶块和两个 11 阶块。这个差别在不变因子里看得见(见 3.6 末),在「dim⁡null⁡(B−2I)\dim\operatorname{null}(B-2I)」里也看得见(22 vs 33)。

三个容易踩的坑

一、DkD_k 是最大公因式,不是子式的乘积。 BB 的非零二阶子式有三个都是 (λ−2)2(\lambda-2)^2、一个 −(λ−2)-(\lambda-2);取公因式得 λ−2\lambda-2,不能把 (λ−2)2(\lambda-2)^2 当答案。
二、DkD_k 永远不会是零多项式。 因为 det⁡(λI−A)≠0\det(\lambda I-A)\neq0,λI−A\lambda I-A 在 F(λ)\mathbb{F}(\lambda) 上可逆,秩恒为 nn,故至少有一个 kk 阶子式非零。所以「Dk=0D_k=0」这种情形根本不会发生——真算出 00 一定是算错了。
三、λI−A\lambda I-A 的秩不必等于 AA 的秩。 把 λ\lambda 取成某个特征值时会掉秩,但那只是某一个 λ\lambda 处的取值,与 DkD_k 作为多项式的公因式不是一回事。BB 在 λ=2\lambda=2 处 λI−B=2I−B\lambda I-B=2I-B 的秩是 11(零空间 22 维),可 D2D_2 仍是 λ−2\lambda-2 而不是 00。

3.7 有理标准形:不依赖代数闭域的那种标准形

Jordan 形有个前提:特征值必须在 F\mathbb{F} 里。实矩阵若有复特征值(比如 (0−110)\begin{pmatrix}0&-1\\1&0\end{pmatrix}),在 R\mathbb{R} 上就没有 Jordan 形。

要一个「任何域上都有」的标准形,就得把块从 Jordan 块换成伴侣矩阵,并允许素因子是不可约多项式(而不只是一次因子)。这就是有理标准形(也叫 Frobenius 标准形)。

伴侣矩阵。给定首一多项式 p(z)=a0+a1z+⋯+am−1zm−1+zmp(z)=a_0+a_1z+\cdots+a_{m-1}z^{m-1}+z^m,它的伴侣矩阵是

C(p)=(00⋯0−a010⋯0−a101⋯0−a2⋮⋱⋮00⋯1−am−1)m×mC(p)=\begin{pmatrix}0&0&\cdots&0&-a_0\\ 1&0&\cdots&0&-a_1\\ 0&1&\cdots&0&-a_2\\ \vdots&&\ddots&&\vdots\\ 0&0&\cdots&1&-a_{m-1}\end{pmatrix}_{m\times m}

关键性质:C(p)C(p) 的特征多项式与极小多项式都等于 pp。这正好反过来印证 18 讲那句话——「每个首一多项式都是某个算子的极小多项式」,那个算子就是伴侣矩阵。

有理标准形定理。任何 A∈Fn,nA\in\mathbb{F}^{n,n} 都相似于若干个伴侣矩阵的块对角和:

A ∼ diag⁡(C(p1), C(p2), …, C(pr)).A\ \sim\ \operatorname{diag}\bigl(C(p_1),\,C(p_2),\,\dots,\,C(p_r)\bigr).

这里有一个必须说清的岔路:pip_i 有两种取法,对应两种不同的写法。

  • 初等因子形式:pip_i 取全部初等因子(每个素幂配一个块,同一素多项式可以有多个块)。这与 3.6 的初等因子组一一对应。
  • 不变因子形式:pip_i 取全部非常数不变因子,此时 p1∣p2∣⋯∣prp_1\mid p_2\mid\cdots\mid p_r,块数等于非平凡不变因子的个数。

两种写法只差把块重排、合并(用中国剩余定理:C(pq)∼diag⁡(C(p),C(q))C(pq)\sim\operatorname{diag}(C(p),C(q)) 当 p,qp,q 互素),相似类完全一样。但「pip_i 恰是非常数不变因子」这句只对第二种写法成立,不能同时套在第一种上。

Jordan 形与有理标准形的关系

两者都是「不变量 + 标准形」的思路,区别在于:

Jordan 形 有理标准形
需要特征值在 F\mathbb{F} 内 是 否
块最简单的样子 对角 + 超对角 11 伴侣矩阵
在 C\mathbb{C} 上 判别力最强,块形最简单 判别力相同(都是完全不变量),但块是伴侣块、不好一眼读出块长
在 R\mathbb{R} 上 可能不存在 总存在

所以:代数闭域上用 Jordan 形,一般域上用有理标准形。 有理标准形是把 Jordan 形的初等因子换成了不可约多项式因子。

回到前面的例子。A=(210021002)A=\begin{pmatrix}2&1&0\\0&2&1\\0&0&2\end{pmatrix}:不变因子是 1,1,(λ−2)31,1,(\lambda-2)^3。非常数不变因子只有 (λ−2)3(\lambda-2)^3,所以有理标准形就是一个伴侣矩阵 C((λ−2)3)C\bigl((\lambda-2)^3\bigr):

(λ−2)3=λ3−6λ2+12λ−8,(\lambda-2)^3=\lambda^3-6\lambda^2+12\lambda-8,

即 a0=−8a_0=-8、a1=12a_1=12、a2=−6a_2=-6。代入定义(最后一列放 −ai-a_i):

C=(00−a010−a101−a2)=(00810−12016).C=\begin{pmatrix}0&0&-a_0\\ 1&0&-a_1\\ 0&1&-a_2\end{pmatrix}=\begin{pmatrix}0&0&8\\ 1&0&-12\\ 0&1&6\end{pmatrix}.

验一下:tr⁡C=6\operatorname{tr}C=6,AA 的迹也是 66 ✓;det⁡C=8\det C=8,而三个特征值 2,2,22,2,2 之积是 88 ✓。

再看 B=(210020002)B=\begin{pmatrix}2&1&0\\0&2&0\\0&0&2\end{pmatrix},它的不变因子是 1,λ−2,(λ−2)21,\lambda-2,(\lambda-2)^2。非常数不变因子有两个,所以有理标准形有两个块:

C(λ−2)=( 2 ),C((λ−2)2)=(0−414),C(\lambda-2)=(\,2\,),\qquad C\bigl((\lambda-2)^2\bigr)=\begin{pmatrix}0&-4\\1&4\end{pmatrix},

于是 B∼diag⁡((2), (0−414))B\sim\operatorname{diag}\bigl((2),\ \begin{pmatrix}0&-4\\1&4\end{pmatrix}\bigr)。

这里能看出有理标准形与 Jordan 形的分工

AA 与 BB 的不变因子不同(AA 是 1,1,(λ−2)31,1,(\lambda-2)^3,BB 是 1,λ−2,(λ−2)21,\lambda-2,(\lambda-2)^2),所以不相似——这点两个标准形都能判定。
但注意:AA 的有理标准形只有一个块、BB 有两个块,这已经是两者的全部差别。有理标准形把「块怎么切」记录在不变因子的个数与次数里,而 Jordan 形记录在块的尺寸里——信息量相同,摆法不同。
在 C\mathbb{C} 上二者的对应是一对一的:C((λ−λ0)s)∼Js(λ0)C\bigl((\lambda-\lambda_0)^s\bigr)\sim J_s(\lambda_0)。所以有理标准形不是「更粗」,只是换了一套允许的块(伴侣矩阵而不是 Jordan 块),代价是块长得不直观,好处是不要求特征值落在域内。

一个真正需要有理标准形的例子。R2\mathbb{R}^2 上 A=(0−110)A=\begin{pmatrix}0&-1\\1&0\end{pmatrix}(转 90°90°)。λI−A=(λ1−1λ)\lambda I-A=\begin{pmatrix}\lambda&1\\-1&\lambda\end{pmatrix},D1=1D_1=1(有常数 11),D2=λ2+1D_2=\lambda^2+1,故不变因子是 (1, λ2+1)(1,\ \lambda^2+1)。λ2+1\lambda^2+1 在 R\mathbb{R} 上不可约,初等因子就是它自己,有理标准形是单个伴侣块

C(λ2+1)=(0−110)=A,C(\lambda^2+1)=\begin{pmatrix}0&-1\\1&0\end{pmatrix}=A,

即 AA 本身已经是有理标准形。它在 R\mathbb{R} 上没有 Jordan 形(第五节误区 1),但有有理标准形——这正是这套工具存在的理由。同一个矩阵搬到 C\mathbb{C} 上,λ2+1=(λ−i)(λ+i)\lambda^2+1=(\lambda-\mathrm{i})(\lambda+\mathrm{i}),初等因子裂成两个,Jordan 形成了 diag⁡(i,−i)\operatorname{diag}(\mathrm{i},-\mathrm{i})。

换域时什么变了、什么没变

上例里不变因子 (1,λ2+1)(1,\lambda^2+1) 在 R\mathbb{R} 与 C\mathbb{C} 上一字不变,初等因子却从 {λ2+1}\{\lambda^2+1\} 变成 {λ−i,λ+i}\{\lambda-\mathrm{i},\lambda+\mathrm{i}\}。
这就是「不变因子与基域无关,初等因子依赖基域」的实例,也是「不变」这个形容词的来由。
实用推论:判断两个实矩阵是否在 R\mathbb{R} 上相似,看不变因子;它们是否在 C\mathbb{C} 上相似,也看不变因子。 两者答案永远一致。

为什么叫「有理」

「有理」在这里不是指有理数,而是指有理标准形只用到域 F\mathbb{F} 自身的运算(有理运算),不需要扩域到代数闭包。
这正是它在实矩阵、有理数矩阵上都能用的原因;每个块的元素都是 AA 的系数的有理表达式。


3.8 相似:把「无从下手」变成「算一串多项式」

λI−A\lambda I-A 与 λI−B\lambda I-B 之间能通过初等行、列变换互相化过去,就叫等价。判据是:

相似判据(两个版本)

设 A,BA,B 是同阶方阵。下列四条等价:

  1. A∼BA\sim B;
  2. λI−A\lambda I-A 与 λI−B\lambda I-B 等价;
  3. λI−A\lambda I-A 与 λI−B\lambda I-B 有相同的 Smith 标准形,即不变因子相同;
  4. 初等因子组相同。

在代数闭域上还可再添一条:Jordan 形相同(不计块序);在任何域上:有理标准形相同。

「1⇒21\Rightarrow2」几乎是白送的:若 A=PBP−1A=PBP^{-1} 且 PP 是常数矩阵,则

λI−A=λI−PBP−1=P (λI−B) P−1,\lambda I-A=\lambda I-PBP^{-1}=P\,(\lambda I-B)\,P^{-1},

右端正是对 λI−B\lambda I-B 左右乘可逆的 λ\lambda-矩阵,故两者等价。

「2⇒12\Rightarrow1」才是实质内容,靠 Smith 形是不变量、且能从它读回标准形(有理标准形或 Jordan 形)——见 3.7。这两个方向合起来,把「相似」这个关于基变换的问题,翻译成了关于多项式整除的问题。

别把它退化成「特征多项式相同」

下面三对矩阵,特征多项式全都一样,却两两不相似,差别一层比一层深:

一对矩阵 特征多项式 极小多项式 几何重数 不相似的原因
J2(2)⊕J2(2)J_2(2)\oplus J_2(2) 与 J2(2)⊕J1(2)⊕J1(2)J_2(2)\oplus J_1(2)\oplus J_1(2) 都是 (λ−2)4(\lambda-2)^4 都是 (λ−2)2(\lambda-2)^2 22 vs 33 块数不同
J3(2)⊕J2(2)J_3(2)\oplus J_2(2) 与 J3(2)⊕J1(2)⊕J1(2)J_3(2)\oplus J_1(2)\oplus J_1(2) 都是 (λ−2)5(\lambda-2)^5 都是 (λ−2)3(\lambda-2)^3 22 vs 33 块数不同
J3(2)J_3(2) 与 J2(2)⊕J1(2)J_2(2)\oplus J_1(2) 都是 (λ−2)3(\lambda-2)^3 (λ−2)3(\lambda-2)^3 vs (λ−2)2(\lambda-2)^2 11 vs 22 极小多项式不同

第一、二行是最阴的:特征多项式相同、极小多项式也相同、最长块也相同,仍不相似。
能一次分开它们的是不变因子:J2(2)⊕J2(2)J_2(2)\oplus J_2(2) 的不变因子是 (1,1,(λ−2)2,(λ−2)2)\bigl(1,1,(\lambda-2)^2,(\lambda-2)^2\bigr),
J2(2)⊕J1(2)⊕J1(2)J_2(2)\oplus J_1(2)\oplus J_1(2) 的是 (1,λ−2,λ−2,(λ−2)2)\bigl(1,\lambda-2,\lambda-2,(\lambda-2)^2\bigr)。
必须比到不变因子(等价地:初等因子)这一层,判断才有保证。

3.9 给定 AA,一步步算出 Jordan 形

流程四步:

  1. 写出 λI−A\lambda I-A,求 D1,…,DnD_1,\dots,D_n——DkD_k 是全部 kk 阶子式的首一最大公因式。
  2. 相除得不变因子 dk=Dk/Dk−1d_k=D_k/D_{k-1},等于 11 的直接丢掉。
  3. 把每个非平凡 dkd_k 分解成素幂,收集全部素幂,这是初等因子。
  4. 每个 (λ−λ0)s(\lambda-\lambda_0)^s 写一个 Js(λ0)J_s(\lambda_0),按任意顺序拼成分块对角。

第 1 步是唯一费力的。三处能省力:DnD_n 就是特征多项式(不必另算);D1D_1 只要矩阵里有一个非零常数元素就等于 11;DkD_k 只需算到「公因式稳定」为止。

**完整手算例。**设

A=(1100011000200002),λI−A=(λ−1−1000λ−1−1000λ−20000λ−2).A=\begin{pmatrix}1&1&0&0\\0&1&1&0\\0&0&2&0\\0&0&0&2\end{pmatrix},\qquad \lambda I-A=\begin{pmatrix}\lambda-1&-1&0&0\\0&\lambda-1&-1&0\\0&0&\lambda-2&0\\0&0&0&\lambda-2\end{pmatrix}.

**第一步:**D1=1D_1=1,因为一阶子式里有常数 −1-1。

**第二步:**D2D_2。取第 1,21,2 行与第 2,32,3 列,子式是 ∣−10λ−1−1∣=1\begin{vmatrix}-1&0\\\lambda-1&-1\end{vmatrix}=1。故 D2=1D_2=1。

**第三步:**D3D_3。三阶子式共 (43)(43)=16\binom43\binom43=16 个。先看重一个事实:第 33 行只有第 33 列那个 λ−2\lambda-2 非零,第 44 行只有第 44 列那个 λ−2\lambda-2 非零(λI−A\lambda I-A 是上三角的)。

任取一个三阶子式,它要用三行。若不用第 3,43,4 行,就只剩第 1,21,2 两行,凑不满三行——所以任何三阶子式都至少用上第 33 行或第 44 行。沿那一行展开,它就等于 (λ−2)(\lambda-2) 乘上某个二阶子式,因此每个三阶子式都被 λ−2\lambda-2 整除。

再算三个具体的,确认公因式正好是 λ−2\lambda-2、没有更高次:

∣λ−1−100λ−1−100λ−2∣=(λ−1)2(λ−2)(第 1,2,3 行与列),\begin{vmatrix}\lambda-1&-1&0\\0&\lambda-1&-1\\0&0&\lambda-2\end{vmatrix}=(\lambda-1)^2(\lambda-2)\quad(\text{第 }1,2,3\text{ 行与列}),

∣−100λ−1−1000λ−2∣=λ−2(第 1,2,4 行,第 2,3,4 列),\begin{vmatrix}-1&0&0\\\lambda-1&-1&0\\0&0&\lambda-2\end{vmatrix}=\lambda-2\quad(\text{第 }1,2,4\text{ 行,第 }2,3,4\text{ 列}),

∣−1000λ−2000λ−2∣=−(λ−2)2(第 1,3,4 行,第 2,3,4 列).\begin{vmatrix}-1&0&0\\0&\lambda-2&0\\0&0&\lambda-2\end{vmatrix}=-(\lambda-2)^2\quad(\text{第 }1,3,4\text{ 行,第 }2,3,4\text{ 列}).

第二个恰等于 λ−2\lambda-2 本身,说明公因式的次数不可能更高。故

D3=λ−2.D_3=\lambda-2 .

**第四步:**D4=det⁡(λI−A)=(λ−1)2(λ−2)2D_4=\det(\lambda I-A)=(\lambda-1)^2(\lambda-2)^2(上三角,直接读对角线)。

于是

d1=1,d2=1,d3=λ−2,d4=(λ−1)2(λ−2).d_1=1,\quad d_2=1,\quad d_3=\lambda-2,\quad d_4=(\lambda-1)^2(\lambda-2).

**第五步:**分解。d3=λ−2d_3=\lambda-2 与 d4=(λ−1)2(λ−2)d_4=(\lambda-1)^2(\lambda-2),素幂是

(λ−1)2,(λ−2),(λ−2).(\lambda-1)^2,\qquad (\lambda-2),\qquad (\lambda-2).

结论:

Jordan 形=diag⁡(J2(1), J1(2), J1(2)).\text{Jordan 形}=\operatorname{diag}\bigl(J_2(1),\ J_1(2),\ J_1(2)\bigr).

核对:块长之和 2+1+1=4=dim⁡V2+1+1=4=\dim V ✓;最长块 22,而极小多项式 d4=(λ−1)2(λ−2)d_4=(\lambda-1)^2(\lambda-2) 中 (λ−1)(\lambda-1) 的指数确实是 22 ✓;λ=2\lambda=2 配两个 1×11\times1 块,故 dim⁡null⁡(A−2I)=2\dim\operatorname{null}(A-2I)=2,与直接数「A−2IA-2I 的零空间维数」相符 ✓。

这套流程值不值得背

值得,因为它是唯一一条不看基、不求特征向量、纯靠行列式就能走完的路。
手算时若矩阵不大(n≤4n\le 4),把子式按「结构」筛一遍比硬列 1616 个可靠;nn 大时任一手算都会崩,那时应该交给计算机,而判据本身(不变因子相同   ⟺  \iff 相似)仍然是那句话。
还有一个加速技巧:先算特征多项式并因式分解,再对每个特征值 λ0\lambda_0 只看 dim⁡null⁡(A−λ0I)k\dim\operatorname{null}(A-\lambda_0I)^k(3.5 的公式),块长直接读出。两条路可以互相验算——本讲的例题都用两种方法对过。

与 3.1–3.5 的关系

3.1–3.5 走的是「算子—链—零空间维数」的路线,证明的是存在性与块长分布的唯一性,用的是线性映射的语言。
3.6–3.9 走的是「矩阵—多项式—不变因子」的路线,给的是算法与判据,用的是 F[λ]\mathbb{F}[\lambda] 的语言。
两条路结论完全一致:3.9 的 4×44\times4 例题算出的 diag⁡(J2(1),J1(2),J1(2))\operatorname{diag}(J_2(1),J_1(2),J_1(2)),用 3.5 的零空间公式复核也得同样的块长(λ=2\lambda=2 的 dim⁡null⁡\dim\operatorname{null} 是 22 而 dim⁡null⁡2\dim\operatorname{null}^2 是 22,故两个一阶块;λ=1\lambda=1 的 dim⁡null⁡\dim\operatorname{null} 是 11 而 dim⁡null⁡2\dim\operatorname{null}^2 是 22,故一个二阶块)。
一条解释「为什么」,一条回答「怎么算」。别把它们看成两套矛盾的理论。


四、几何直观 / 记忆钩子

4.1 一张图:链

每个 Jordan 块背后是一条链。幂零情形(λ=0\lambda=0)画成

u →T Tu →T T2u →T ⋯ →T Tm−1u →T 0u\ \xrightarrow{T}\ Tu\ \xrightarrow{T}\ T^2u\ \xrightarrow{T}\ \cdots\ \xrightarrow{T}\ T^{m-1}u\ \xrightarrow{T}\ 0

一般情形(λ≠0\lambda\neq0),链的末端不是零,而是一个特征向量:

v →T Tv →T ⋯ →T (T 不再把它推到新方向)v\ \xrightarrow{T}\ Tv\ \xrightarrow{T}\ \cdots\ \xrightarrow{T}\ \big(T\text{ 不再把它推到新方向}\big)

矩阵里的 11 就是这些箭头的记录:每个 11 都是一次「被推着走」;一旦不再被推,块就到头了。

4.2 读矩阵的三句口诀

拿到一个 Jordan 形矩阵,不看数字只看形状:

  1. 数块:块与块的分界处上方是 00;块数 == 该 λ\lambda 的几何重数。
  2. 量长:从左上一路往右下数连续 11 的个数,加 11 就是块长;最长的块长 == 幂零指数。
  3. 定 λ\lambda:块内对角元的共同值。

4.3 三个极端,帮你定位

情形 Jordan 形长什么样 与前面讲的对应
TT 可对角化 全是 1×11\times1 块,即对角阵 20 讲
TT 幂零 所有块都是 Jd(0)J_d(0),对角线上全 00 31 讲
T=λI+NT=\lambda I+N,NN 幂零 一个(或几个)Jd(λ)J_d(\lambda) 本讲

可对角化就是「每条链都只有一步」——这句话把 20 讲和本讲接上了。


五、常见误区

误区 1:以为实数域上也有 Jordan 形

定理 B 的前提是 F=C\mathbb{F}=\mathbb{C}。
反例:R2\mathbb{R}^2 上 T(x,y)=(−y,x)T(x,y)=(-y,x)(逆时针转 90°90°)。特征多项式是 z2+1z^2+1,没有实根,连一个特征值都拿不出来。
它不可能有 Jordan 形——Jordan 形的对角线上必须填特征值。
(在 C2\mathbb{C}^2 上看它反而简单:特征值 ±i\pm\mathrm{i},两个 1×11\times1 块。)
实空间上能做到 32 讲的分块上三角,做不到 Jordan 形。

误区 2:把分块上三角当成 Jordan 形

32 讲给的是「每块 λk\lambda_k 在对角线、上面任意」;Jordan 形要求超对角线上只能是 00 或 11。
第一节那个矩阵是上三角,也「差不多」是分块形式,但超对角线上是 33 而不是 11——它不是 Jordan 形。
用 3.4 的基换过之后才成为 (610060007)\begin{pmatrix}6&1&0\\0&6&0\\0&0&7\end{pmatrix}。

误区 3:以为块数等于重数

重数(代数重数)=dim⁡G(λ,T)=\dim G(\lambda,T),是该 λ\lambda 所有块的尺寸之和。
块数(几何重数)=dim⁡E(λ,T)=\dim E(\lambda,T),是块的个数。
例:(210021002)\begin{pmatrix}2&1&0\\0&2&1\\0&0&2\end{pmatrix} 中 λ=2\lambda=2 的重数是 33,块数却只有 11(就一个 J3(2)J_3(2))。
换个写法:(210020002)\begin{pmatrix}2&1&0\\0&2&0\\0&0&2\end{pmatrix} 重数还是 33,但块数是 22(J2(2)⊕J1(2)J_2(2)\oplus J_1(2))。
两者只在可对角化时相等——那时块数 == 重数。

误区 4:以为 λ\lambda 相同的块该合并

两个 J2(6)J_2(6) 块不能合并成 J4(6)J_4(6)。用 dim⁡null⁡(T−6I)\dim\operatorname{null}(T-6I) 就能看出区别:
J2(6)⊕J2(6)J_2(6)\oplus J_2(6) 的零空间是 22 维,J4(6)J_4(6) 的只有 11 维。
块的个数是算子自己的属性,不是排版习惯。

误区 5:以为 Jordan 块的对角线一定是 00

00 只是幂零情形的特例。λ≠0\lambda\neq0 时,超对角线上的 11 表示「TT 沿着链推一格,同时加一个 λ\lambda」,
而不是「TT 把向量变成 00」。别把 Jd(0)J_d(0) 的直觉套到 Jd(λ)J_d(\lambda) 上。

误区 6:把 Jordan 形当成「像对角形那样的唯一标准形」

它只在「块的顺序可以任意」这个意义下唯一;3.2/3.3 证的是存在性,唯一性见 3.5 末。
想拿它当判据(例如判断两个矩阵是否相似),必须先把唯一性那块补上。


六、拓展(可跳过)

6.1 用 Jordan 结构算函数:幂零为什么这么好用

设 NN 幂零、Nm=0N^m=0。那么对任意多项式 pp,

p(λI+N)=∑k≥0p(k)(λ)k! Nkp(\lambda I+N)=\sum_{k\ge0}\frac{p^{(k)}(\lambda)}{k!}\,N^k

这个展开里只有前 mm 项非零。更一般地,任何形式幂级数在 λI+N\lambda I+N 上都自动截断成有限和,收敛性问题消失。这是 Jordan 形最实用的地方。

举一个本讲内容直接支持的例子:I+TI+T 的平方根(TT 幂零)。考虑

1+x=1+a1x+a2x2+a3x3+⋯ ,a1=12, a2=−18, a3=116.\sqrt{1+x}=1+a_1x+a_2x^2+a_3x^3+\cdots,\qquad a_1=\tfrac12,\ a_2=-\tfrac18,\ a_3=\tfrac1{16}.

把 xx 换成 TT、11 换成 II,令

R=I+a1T+a2T2+⋯+am−1Tm−1.R=I+a_1T+a_2T^2+\cdots+a_{m-1}T^{m-1}.

要 R2=I+TR^2=I+T,只需逐次让 TkT^k 的系数对上。前几项展开:

R2=I+2a1T+(2a2+a12)T2+(2a3+2a1a2)T3+⋯R^2=I+2a_1T+\big(2a_2+a_1^2\big)T^2+\big(2a_3+2a_1a_2\big)T^3+\cdots

令 TT 的系数为 11:2a1=12a_1=1,得 a1=12a_1=\tfrac12。
令 T2T^2 的系数为 00:2a2+a12=02a_2+a_1^2=0,得 a2=−18a_2=-\tfrac18。
令 T3T^3 的系数为 00:2a3+2a1a2=02a_3+2a_1a_2=0,得 a3=116a_3=\tfrac1{16}。

第 kk 步的方程总形如 2ak+(只含 a1,…,ak−1 的项)=02a_k+(\text{只含 }a_1,\dots,a_{k-1}\text{ 的项})=0,因为 2≠02\neq0,总能解出 aka_k。于是 RR 存在。这个构造对实、复空间都成立,只用到 TT 幂零。

再往上一层:TT 可逆时,在每个 G(λk,T)G(\lambda_k,T) 上写成

T∣G(λk,T)=λk(I+Nkλk),T\big|_{G(\lambda_k,T)}=\lambda_k\Big(I+\frac{N_k}{\lambda_k}\Big),

其中 Nk/λkN_k/\lambda_k 幂零,所以 I+Nk/λkI+N_k/\lambda_k 有平方根;把复数 λk\lambda_k 的平方根与它相乘,就得到该块上的平方根;再用直和拼起来,得到整个 TT 的平方根。这就是「分块 + 幂零」这台机器的一次真实运转。

反过来也有意思

C3\mathbb{C}^3 上 T(z1,z2,z3)=(z2,z3,0)T(z_1,z_2,z_3)=(z_2,z_3,0) 恰是一个 J3(0)J_3(0),它没有平方根。
原因不只是「不可逆」:更精确地说,是链太长。
设 R2=TR^2=T,则 RR 与 TT 交换,于是 RR 必须把 null⁡T\operatorname{null}T 送进 null⁡T\operatorname{null}T……
顺着推几步就会撞上矛盾(第七节自测第 5 题)。
方根存在与否,取决于链的结构,而不只是可逆性。

6.2 可对角化 = Jordan 形全是一阶块

若 TT 可对角化,Jordan 形就是对角阵(全部 1×11\times1 块),反之亦然。这条把 20 讲与本节并成一句:可对角化是「每个特征值只有块长为 1 的链」的极端情形。


七、自测(先做,再看答案)

1. 设 T∈L(C3)T\in\mathcal{L}(\mathbb{C}^3) 为 T(z1,z2,z3)=(0,z1,0)T(z_1,z_2,z_3)=(0,z_1,0)。求 TT 的 Jordan 形,并说出块数与各块长。

2. 判断:矩阵 (634062007)\begin{pmatrix}6&3&4\\0&6&2\\0&0&7\end{pmatrix} 是不是 Jordan 形?若不是,该往哪个方向改?

3. 矩阵 B=(2100020000210002)B=\begin{pmatrix}2&1&0&0\\0&2&0&0\\0&0&2&1\\0&0&0&2\end{pmatrix} 中 λ=2\lambda=2 的重数与块数各是多少?它能对角化吗?

4. 用本讲的方法,写出 T(z1,z2,z3)=(6z1+3z2+4z3, 6z2+2z3, 7z3)T(z_1,z_2,z_3)=(6z_1+3z_2+4z_3,\ 6z_2+2z_3,\ 7z_3) 的一组 Jordan 基,并验证矩阵。

5.(难) 证明 C3\mathbb{C}^3 上 T(z1,z2,z3)=(z2,z3,0)T(z_1,z_2,z_3)=(z_2,z_3,0) 没有平方根。

答案(点开)

1. Te1=e2Te_1=e_2、Te2=0Te_2=0、Te3=0Te_3=0。
链:e1↦e2↦0e_1\mapsto e_2\mapsto0(长 22),e3↦0e_3\mapsto0(长 11)。
取基 (e2,e1,e3)(e_2,e_1,e_3)(每条链倒排、依次拼接):

M=(010000000)=J2(0)⊕J1(0).M=\begin{pmatrix}0&1&0\\0&0&0\\0&0&0\end{pmatrix}=J_2(0)\oplus J_1(0).

块数 22,块长 22 与 11。核对:null⁡T=span⁡(e2,e3)\operatorname{null}T=\operatorname{span}(e_2,e_3) 是 22 维,正好等于块数。

2. 不是。它是上三角,超对角线上是 33(不是 11),还多出右上角两个数。
方向:在 G(6,T)=span⁡(e1,e2)G(6,T)=\operatorname{span}(e_1,e_2) 里把基由 (e1,e2)(e_1,e_2) 换成 (3e1,e2)(3e_1,e_2),
第三个基向量换成 (10,2,1)(10,2,1),得到

diag⁡(J2(6),J1(7))=(610060007).\operatorname{diag}\big(J_2(6),J_1(7)\big)=\begin{pmatrix}6&1&0\\0&6&0\\0&0&7\end{pmatrix}.

3. 两个 2×22\times2 块。重数 =dim⁡G(2,B)=4=\dim G(2,B)=4;块数 =2=2(把 B−2IB-2I 写出来,它只有两个非零元,零空间 22 维)。
不能对角化:块数 2<42<4,特征向量不够 44 个。

4. λ=6\lambda=6:N=T−6IN=T-6I 满足 N(e2)=3e1N(e_2)=3e_1、N2=0N^2=0。
链 e2↦3e1↦0e_2\mapsto3e_1\mapsto0,取基 ((3,0,0),(0,1,0))\big((3,0,0),(0,1,0)\big)。
λ=7\lambda=7:T(10,2,1)=7(10,2,1)T(10,2,1)=7(10,2,1),取基 (10,2,1)(10,2,1)。
合并成一族:((3,0,0), (0,1,0), (10,2,1))\big((3,0,0),\ (0,1,0),\ (10,2,1)\big),矩阵 (610060007)\begin{pmatrix}6&1&0\\0&6&0\\0&0&7\end{pmatrix}。
验算:T(3,0,0)=(18,0,0)=6(3,0,0)T(3,0,0)=(18,0,0)=6(3,0,0) ✓;
T(0,1,0)=(3,6,0)=(3,0,0)+6(0,1,0)T(0,1,0)=(3,6,0)=(3,0,0)+6(0,1,0) ✓;
T(10,2,1)=(70,14,7)=7(10,2,1)T(10,2,1)=(70,14,7)=7(10,2,1) ✓。

5. 记 e1=(1,0,0)e_1=(1,0,0) 等。按题面 T(z1,z2,z3)=(z2,z3,0)T(z_1,z_2,z_3)=(z_2,z_3,0),所以

Te1=0,Te2=e1,Te3=e2,Te_1=0,\qquad Te_2=e_1,\qquad Te_3=e_2,

即 TT 的矩阵是 (010001000)=J3(0)\begin{pmatrix}0&1&0\\0&0&1\\0&0&0\end{pmatrix}=J_3(0)(超对角线上全 11),null⁡T=span⁡(e1)\operatorname{null}T=\operatorname{span}(e_1)。
设 R2=TR^2=T。由 R4=T2R^4=T^2、R6=T3=0R^6=T^3=0。

路线一(顺着交换性):RT=R⋅R2=R3=R2⋅R=TRRT=R\cdot R^2=R^3=R^2\cdot R=TR,所以 RR 与 TT 交换。于是 RR 把 null⁡T\operatorname{null}T 送进 null⁡T\operatorname{null}T:若 Tv=0Tv=0,则 T(Rv)=R(Tv)=0T(Rv)=R(Tv)=0。故 Re1=ce1Re_1=ce_1。
由 R2=TR^2=T 得 R2e1=Te1=0R^2e_1=Te_1=\mathbf 0,而 R2e1=R(ce1)=c2e1R^2e_1=R(ce_1)=c^2e_1,于是 c2=0c^2=0,即 c=0c=0、Re1=0Re_1=0。
再看 e2e_2:Te2=e1≠0Te_2=e_1\ne0,故 e2∈ker⁡T2e_2\in\ker T^2(因 T2e2=Te1=0T^2e_2=Te_1=0)。而 RR 保持 ker⁡T2\ker T^2 不变(与 TT 交换),且 ker⁡T2=span⁡(e1,e2)\ker T^2=\operatorname{span}(e_1,e_2)。设 Re2=αe1+βe2Re_2=\alpha e_1+\beta e_2,则

R2e2=Re1⋅α+R(βe2)=αRe1+β(αe1+βe2)=βαe1+β2e2,R^2e_2=Re_1\cdot\alpha+R(\beta e_2)=\alpha Re_1+\beta(\alpha e_1+\beta e_2)=\beta\alpha e_1+\beta^2e_2,

它必须等于 Te2=e1Te_2=e_1,故 β2=1\beta^2=1 且 αβ=1\alpha\beta=1,即 α=β=±1\alpha=\beta=\pm1。代入验 R2e2=e1R^2e_2=e_1 无矛盾——这条路要再往下才出矛盾,写起来偏长。

路线二(走幂零指数,最干净):由 R2=TR^2=T 得 R4=T2R^4=T^2。T2T^2 的矩阵是 TT 的平方,直接算得 T2=(001000000)≠0T^2=\begin{pmatrix}0&0&1\\0&0&0\\0&0&0\end{pmatrix}\ne0,故 R4≠0R^4\ne0。
另一方面 R6=T3=0R^6=T^3=0,即 RR 是 R3\mathbb{R}^3 上的幂零算子,幂零指数 ≤3\le3(3 维空间上幂零算子的指数不超过维数),所以 R3=0R^3=0。
于是 T2=R4=R⋅R3=0T^2=R^4=R\cdot R^3=0,与 R4≠0R^4\ne0 矛盾。
所以 TT 没有平方根。(这个论证只用到「TT 是一条长 33 的链」,正好印证 6.1 的说法。)


八、复习卡

项 内容
Jordan 块 Jd(λ)J_d(\lambda) d×dd\times d:对角全 λ\lambda,超对角线全 11,其余 00
等价写法 Jd(λ)=λId+NdJ_d(\lambda)=\lambda I_d+N_d,Ndd=0N_d^d=0,Ndd−1≠0N_d^{d-1}\neq0
Jordan 基 使 TT 的矩阵为各 Jdk(λk)J_{d_k}(\lambda_k) 的分块对角的基
Jordan 形 该矩阵;λk\lambda_k 不必互异
链(幂零) u↦Tu↦⋯↦Tm−1u↦0u\mapsto Tu\mapsto\cdots\mapsto T^{m-1}u\mapsto0
块从哪来 一条链一个块;基倒序排即得超对角线上的 11
块数 =dim⁡E(λ,T)=dim⁡null⁡(T−λI)=\dim E(\lambda,T)=\dim\operatorname{null}(T-\lambda I),即几何重数
块长之和 =dim⁡G(λ,T)=\dim G(\lambda,T),即代数重数
最长块长 极小多项式中 (z−λ)(z-\lambda) 的指数;幂零时即幂零指数
存在性(幂零) 定理 A:任何域都成立;归纳 + 泛函造补空间 WW
存在性(一般) 定理 B:需 F=C\mathbb{F}=\mathbb{C};广义特征空间分解 + 定理 A
唯一性 块长分布由 dim⁡null⁡Nk\dim\operatorname{null}N^k 逐层确定(bk−bk+1=b_k-b_{k+1}= 长度恰为 kk 的链数);不计块顺序时 Jordan 形唯一
数域限制 R\mathbb{R} 上一般没有;反例:90°90° 旋转
与对角化 可对角化   ⟺  \iff 全为 1×11\times1 块(每条链只有一步)
规范化 块内超对角线的非零元可靠换基改成 11;00 表示块的接缝
行列式因子 DkD_k λI−A\lambda I-A 的全部 kk 阶子式的首一最大公因式;Dn=D_n= 特征多项式
不变因子 dkd_k dk=Dk/Dk−1d_k=D_k/D_{k-1},满足 d1∣⋯∣dnd_1\mid\cdots\mid d_n;dn=d_n= 极小多项式
初等因子 各 dkd_k 分解成素幂后,收齐全部 (λ−λ0)s(\lambda-\lambda_0)^s(按重数重复)
主对应 一个 (λ−λ0)s(\lambda-\lambda_0)^s ⟷\longleftrightarrow 一个 Js(λ0)J_s(\lambda_0);次数之和 =n=n
块数/块长怎么读 λ0\lambda_0 配几个块 == 含 λ−λ0\lambda-\lambda_0 的初等因子个数;块长 == 各指数
相似判据 A∼B  ⟺  λI−AA\sim B\iff\lambda I-A 与 λI−B\lambda I-B 等价   ⟺  \iff 不变因子相同   ⟺  \iff 初等因子组相同
最常见的误判 特征多项式相同(甚至极小多项式也相同)不足以保证相似
有理标准形 素幂 pep^{e} 配伴侣块 C[pe]C[p^{e}];任何域上都存在,不需特征值
伴侣矩阵 次对角线全 11,最后一列放 −ai-a_i;det⁡(λI−C[p])=p(λ)\det(\lambda I-C[p])=p(\lambda)
两种挂法 初等因子形式(每个素幂一块)/不变因子形式(每个非常数 dkd_k 一块)
与基域的关系 不变因子与基域无关;初等因子随基域变化

九、出处与延伸(可选,不影响阅读)

  • 主线索:Sheldon Axler《线性代数应该这样学》第 4 版(吴俊达、何阳 译)§8C「广义特征空间分解的推论」,印刷页码 p266–271(PDF p280–285)。本讲的存在性定理、两个幂零算子的例子、Jordan 基的定义都出自这一节。
  • 上一层的准备:广义特征空间与它的直和分解在同一章 §8B(印刷页码 p257 起),本讲只在 3.3 用它来拼接,细节见 32 讲。
  • 幂零算子的基本性质:31 讲;上三角化:19 讲;可对角化:20 讲。
  • 唯一性的证明:本讲 3.5 末已给完整路线——用 dim⁡null⁡(T−λI)k\dim\operatorname{null}(T-\lambda I)^k(k=1,2,3,…k=1,2,3,\dots)这一族数刻画块长分布,可查任何一本标准教材的「Jordan 标准形」章节,或「有理标准形 / Frobenius 标准形」相关材料。
  • 术语别名:「Jordan 形」=「若尔当型」=「Jordan 标准形」=「Jordan canonical form」;「广义特征空间」=「根空间(root space)」;「超对角线」=「superdiagonal」;「幂零指数」=「nilpotency index」。
  • 初等因子 / 不变因子 / 有理标准形:Steven Roman《高等线性代数》(Advanced Linear Algebra, 3rd ed., GTM 135)第 7 章「有理标准形与 Jordan 标准形」。初等因子与不变因子的定义(DkD_k 与 dkd_k)、有理标准形的两种形式、Jordan 标准形作为代数闭域上的特例,均见该章(PDF p148–p171)。本讲 3.7 的「不变因子形式 / 初等因子形式」之分就是那一章的术语;「不变因子与基域无关、初等因子随基域变化」是该章末的定理,3.7 换域的例子正是在演示它。
  • λ\lambda-矩阵与初等因子的算法:Siegfried Bosch《线性代数》(Lineare Algebra, Springer 2014)第 6 章「初等因子理论」。那里从主理想环上的矩阵消元讲起,给出「逐次用带余除法压低次数、把主对角线摆成整除链」的构造(PDF p188–p199),并用一个手算例子走完全程(PDF p218–p219)。本讲 3.6 的「DkD_k 是 kk 阶子式的最大公因式、Smith 形主对角就是 dkd_k」与 3.9 的四步流程都与这套做法一致。
  • 伴侣矩阵与 Frobenius 标准形:上述两书均有;「伴侣矩阵的极小多项式与特征多项式都等于 pp」是该理论的标准入口,本讲 3.7 用它说明 18 讲那句「每个首一多项式都是某个算子的极小多项式」。
  • 历史:Camille Jordan(1838–1922)于 1870 年发表该分解的证明;Frobenius 的有理标准形(1887 前后)比 Jordan 形更早给出「任意域上的标准形」。
  • 关于平方根:a1=12, a2=−18, a3=116a_1=\tfrac12,\ a_2=-\tfrac18,\ a_3=\tfrac1{16} 与二项级数 (1+x)1/2=∑k(1/2k)xk(1+x)^{1/2}=\sum_k\binom{1/2}{k}x^k 一致;本讲只用到「能逐步解出」这一点,不依赖级数收敛。