在上一节中,我们介绍了 Jordan 标准型的存在性以及通用的构造方法(寻找广义特征向量链)。
然而,直接求解高阶矩阵的广义特征空间往往计算量巨大。
Cayley-Hamilton 定理告诉我们 FA(A)=O,即特征多项式是一个可以将矩阵 “零化” 的多项式。
但特征多项式往往显得 “太大” 了,是否存在一个次数更低的多项式也能让矩阵归零?这个多项式是否蕴含了关于 Jordan 块结构的秘密?
# 极小多项式
定义
令 A 为 F 上的 n 阶方阵。
在所有满足 P(A)=O 的非零多项式 P(λ)∈F[λ] 中,次数最低且首项系数为 1 的多项式,称为 A 的 极小多项式 (Minimal Polynomial)「最小多項式」,记作 mA(λ)。
极小多项式具有极其优良的性质,它是矩阵代数结构的 “基因”。
命题 极小多项式的性质
- 唯一性:矩阵 A 的极小多项式是唯一的。
- 整除性:若多项式 g(λ) 满足 g(A)=O,则 mA(λ) 整除 g(λ)。特别地,mA(λ) 整除特征多项式 FA(λ)。
- 根的同一性:mA(λ) 与 FA(λ) 拥有完全相同的根(特征值),只是重数可能不同。
证明
(1) & (2)
设 g(λ) 满足 g(A)=O。做带余除法:g(λ)=mA(λ)q(λ)+r(λ),其中 degr<degmA 或 r=0。
代入 A 得 g(A)=mA(A)q(A)+r(A)⟹O=O⋅q(A)+r(A)⟹r(A)=O。
由于 mA 是次数最低的零化多项式,若 r 非零且次数更低,则矛盾。故 r(λ)=0,即 mA∣g。
若有两个首一的极小多项式 m1,m2,则互为整除且首一,故 m1=m2。
(3)
由 (2) 可知 FA(λ)=mA(λ)k(λ),故 mA 的根必为 FA 的根。
反之,设 λ0 为 FA 的根(特征值),对应特征向量 x=0。
有 Ax=λ0x⟹Akx=λ0kx⟹mA(A)x=mA(λ0)x。
因 mA(A)=O,故 mA(λ0)x=0。因 x=0,故 mA(λ0)=0。
即特征值必为极小多项式的根。
□
利用极小多项式,我们可以给出一个判定对角化的最强充要条件。
定理
n 阶方阵 A 可对角化的充分必要条件是:
A 的极小多项式 mA(λ) 可以分解为互不相同的线性因式的乘积(即无重根)。
mA(λ)=(λ−λ1)(λ−λ2)…(λ−λk),λi=λj
- 推论:若 A 是幂等矩阵 (A2=A),则 mA(λ)∣λ(λ−1),无重根,故幂等矩阵必可对角化。
# 不变因子与初等因子
为了精确确定 Jordan 标准型,仅知道特征多项式是不够的,通常还需要结合秩的信息。但从多项式理论的角度,我们可以引入 “不变因子” 来彻底解构矩阵。
考虑特征矩阵 λE−A。这是一个元素为 λ 的多项式的矩阵。
通过行列变换(允许乘以多项式系数,互换,加减),可以将 λE−A 化为对角型,称为 Smith 标准型:
λE−A∼⎝⎜⎜⎜⎛d1(λ)d2(λ)⋱dn(λ)⎠⎟⎟⎟⎞
满足 d1(λ)∣d2(λ)∣⋯∣dn(λ)(整除关系)。
这些非零多项式 di(λ) 称为矩阵 A 的 不变因子 (Invariant Factor)「不変因子」。
命题 不变因子与多项式的关系
- FA(λ)=d1(λ)d2(λ)⋯dn(λ) (不计符号差)
- mA(λ)=dn(λ) (即最大的那个不变因子)
将每个不变因子 di(λ) 在复数域上彻底分解为一次因式幂积的乘积:
di(λ)=(λ−λ1)ei1(λ−λ2)ei2⋯
所有这些分解出来的形如 (λ−λj)k 的幂因子(不含常数项),汇集在一起,统称为矩阵 A 的 初等因子 (Elementary Divisor)「単因子」。
初等因子与 Jordan 块的一一对应关系:
每一个初等因子 (λ−c)k,都严格对应一个特征值为 c,阶数为 k 的 Jordan 块 Jk(c)。
# 极小多项式与 Jordan 型的速解
在实际计算中(尤其是低阶矩阵),我们不需要去求 Smith 标准型。
结合特征多项式 FA(λ)、极小多项式 mA(λ) 以及秩的信息,通常足以快速确定 Jordan 标准型。
以下是判定规则的总结:
- 特征值:FA(λ) 的根决定了对角线上的元素。
- 代数重数:FA(λ)=⋯(λ−c)n⋯,指数 n 决定了特征值 c 对应的所有 Jordan 块的总阶数之和(即 Jordan 矩阵中 c 出现的总次数)。
- 几何重数:dimVc=n−rank(A−cE) 决定了特征值 c 对应的 Jordan 块的个数。
- 极小多项式:mA(λ)=⋯(λ−c)m⋯,指数 m 决定了特征值 c 对应的 最大的那个 Jordan 块的阶数。
示例
设 4 阶方阵 A 的特征多项式为 FA(λ)=(λ−2)4。
此时特征值只有 2。这就意味着 Jordan 矩阵对角线上全为 2。
我们需要确定 Jordan 块的结构(是 4 还是 3+1 还是 2+2 等)。
讨论极小多项式 mA(λ) 的可能情况:
-
若 mA(λ)=(λ−2):
最大块阶数为 1。即所有块都是 J1(2)。
此时 J=diag(2,2,2,2)。且 A 必须等于 2E。
-
若 mA(λ)=(λ−2)2:
最大块阶数为 2。可能的划分:4=2+2 或 4=2+1+1。
- 如何区分?看几何重数(特征向量个数)。
- 若 rank(A−2E)=2,则几何重数 4−2=2,有两个块,故为 J2(2)⊕J2(2)。
- 若 rank(A−2E)=1,则几何重数 4−1=3,有三个块,故为 J2(2)⊕J1(2)⊕J1(2)。
-
若 mA(λ)=(λ−2)3:
最大块阶数为 3。划分只能是 4=3+1。
此时 J=J3(2)⊕J1(2)。
-
若 mA(λ)=(λ−2)4:
最大块阶数为 4。划分只能是 4=4。
此时 J=J4(2)。
如何快速求极小多项式?
- 先写出 FA(λ)。
- 设 mA(λ) 为包含所有不同特征值因子的多项式,指数先设为 1 或未知。
- 验证 (A−λiE)k 是否为零矩阵。从低次幂开始尝试。
示例
求矩阵
A=⎝⎛210021002⎠⎞
的 Jordan 标准型。
解
方法一:利用秩与链
特征多项式显然为 FA(λ)=(λ−2)3。特征值 2,代数重数 3。
计算 A−2E:
A−2E=⎝⎛010001000⎠⎞
rank(A−2E)=2。
几何重数 = 3−2=1。
说明只有一个 Jordan 块。
既然总阶数是 3,又只有一个块,那必然是 J3(2)。
J=⎝⎛200120012⎠⎞
方法二:利用极小多项式
FA(λ)=(λ−2)3。
mA(λ) 可能是 (λ−2),(λ−2)2,(λ−2)3。
计算 (A−2E):非零。
计算 (A−2E)2:
(A−2E)2=⎝⎛010001000⎠⎞⎝⎛010001000⎠⎞=⎝⎛001000000⎠⎞=O
计算 (A−2E)3:
(A−2E)3=(A−2E)2(A−2E)=O
所以 mA(λ)=(λ−2)3。
最大 Jordan 块的阶数是 3。因为矩阵是 3 阶的,所以只能有一个块 J3(2)。
□
示例
求矩阵
A=⎝⎛3−441−1−800−2⎠⎞
的 Jordan 标准型。
解
计算特征多项式
FA(λ)=det(λE−A)=(λ+2)[(λ−3)(λ+1)+4]=(λ+2)(λ2−2λ+1)=(λ+2)(λ−1)2
特征值:−2 (重数 1),1 (重数 2)。
对于 λ=−2:代数重数 1,对应的 Jordan 块必然是 J1(−2)(即 [−2])。
对于 λ=1:代数重数 2。可能的块结构是 J2(1) 或 J1(1)⊕J1(1)。
这取决于 A−1E 的秩,或者极小多项式中 (λ−1) 的次数。
计算 A−E:
A−E=⎝⎛2−441−2−800−3⎠⎞
第一行乘以 -2 等于第二行,显然前两行线性相关。
但第三行显然无法由第一行消去(第三列非零)。
实际上 rank(A−E)=2。
几何重数 = 3−2=1。
说明对应 λ=1 只有一个 Jordan 块。
因为代数重数为 2,所以这个块必然是 J2(1)。
综上,Jordan 标准型为
J=⎝⎛−200010011⎠⎞
(注:块的顺序不唯一)
□
通过结合特征多项式(决定总成分)、极小多项式(决定最大块)、秩与几何重数(决定块的数量),我们可以像拼积木一样快速确定 Jordan 标准型,而无需进行繁琐的基底构造过程。
内容已经过 Gemini 3.0 Pro 审查