Jordan 形 #
这讲讲什么 :给复向量空间上的每个算子挑一组基,把矩阵压到只剩对角线和紧挨其上的一条线——对角线填特征值,紧邻上方填 0 0 0 或 1 1 1 。这就是 Jordan 形。前置 :32 讲(广义特征空间分解);31 讲(幂零算子、广义特征向量);19 讲(上三角化);20 讲(可对角化);08 讲的线性映射基本定理(第三节的维数估计要用一次)一句话 :先在 V V V 上按特征值分堆,再把每一堆里的幂零部分拆成若干条「链」,每条链正好给一个 Jordan 块。
一、为什么需要它 #
1.1 三步走,每一步都留下一点没收拾干净的东西 #
20 讲已经说清:最理想的基是特征向量基,使 T T T 的矩阵是对角阵。可惜不是每个算子都有这么多特征向量。
19 讲退一步:复向量空间上每个算子都有上三角矩阵,对角线上按重数排着全部特征值。代价是上三角部分可以塞任意数字。取
A = ( 6 3 4 0 6 2 0 0 7 ) A=\begin{pmatrix}6&3&4\\0&6&2\\0&0&7\end{pmatrix} A = 6 0 0 3 6 0 4 2 7
对角线上读出特征值 6 , 6 , 7 6,6,7 6 , 6 , 7 。可右上角那三个数 3 , 4 , 2 3,4,2 3 , 4 , 2 有什么规律?没有。换一组基它们就变。
32 讲再进一步:V V V 能按特征值拆成广义特征空间的直和
V = G ( λ 1 , T ) ⊕ ⋯ ⊕ G ( λ m , T ) , V=G(\lambda_1,T)\oplus\cdots\oplus G(\lambda_m,T), V = G ( λ 1 , T ) ⊕ ⋯ ⊕ G ( λ m , T ) ,
每一块在 T T T 下不变,而且 T T T 在 G ( λ k , T ) G(\lambda_k,T) G ( λ k , T ) 上等于 λ k I \lambda_kI λ k I 加上一个幂零算子。于是矩阵成了分块对角,每块是 λ k \lambda_k λ k 在对角线的上三角。可块内还有一串 ∗ * ∗ 没定。
这一讲要做的,就是把最后那串 ∗ * ∗ 也收拾掉。 收拾的结果是:块内除了对角线上的 λ \lambda λ ,只剩紧挨对角线上方那条线上的 0 0 0 和 1 1 1 。剩下的 1 1 1 不是垃圾——它精确记录了「这个块有多长」。
1.2 三个层次,一句话各说一遍 #
同一个算子,换不同的基,能得到的矩阵整齐程度不同:
能保证的
靠什么
没收拾掉的东西
上三角
19 讲
上三角部分任意
分块对角,块内上三角、对角元相同
32 讲
块内上三角部分任意
分块对角,块内只有超对角线的 0 0 0 与 1 1 1
本讲
只剩块长与 λ \lambda λ
一句话:上三角是「把垃圾推到角落」,分块对角是「按特征值分堆」,Jordan 形是「每堆里只留下链长的记录」。
1.3 读完这一讲你应该能回答 #
一个 n × n n\times n n × n 矩阵的 Jordan 形里,到底有几个块?块长怎么定?
为什么 λ \lambda λ 相同的两个块不能随随便便合并成一个?
每个复矩阵都有 Jordan 形——那它唯一吗?本节能证到哪一步?
二、定义与最小例子 #
2.1 先把幂零算子摆出来 #
31 讲定义过:T T T 幂零,指存在正整数 m m m 使 T m = 0 T^m=0 T m = 0 。使 T m = 0 T^m=0 T m = 0 成立的最小 m m m 叫 T T T 的幂零指数(T = 0 T=0 T = 0 时指数为 1 1 1 )。
幂零算子的全部特征值都是 0 0 0 :若 T v = λ v Tv=\lambda v T v = λ v 且 v ≠ 0 v\neq0 v = 0 ,则
0 = T m v = λ m v , 0=T^mv=\lambda^mv, 0 = T m v = λ m v ,
故 λ = 0 \lambda=0 λ = 0 。所以下面几个例子里,矩阵对角线上全是 0 0 0 。
2.2 例 A:一条链,一个块 #
在 C 4 \mathbb{C}^4 C 4 上定义
T ( z 1 , z 2 , z 3 , z 4 ) = ( 0 , z 1 , z 2 , z 3 ) . T(z_1,z_2,z_3,z_4)=(0,\ z_1,\ z_2,\ z_3). T ( z 1 , z 2 , z 3 , z 4 ) = ( 0 , z 1 , z 2 , z 3 ) .
容易验 T 4 = 0 T^4=0 T 4 = 0 (每作用一次,非零分量往右挪一格,挪出四格就没了),而 T 3 ≠ 0 T^3\neq0 T 3 = 0 。取 v = ( 1 , 0 , 0 , 0 ) v=(1,0,0,0) v = ( 1 , 0 , 0 , 0 ) ,算出
v = ( 1 , 0 , 0 , 0 ) , T v = ( 0 , 1 , 0 , 0 ) , T 2 v = ( 0 , 0 , 1 , 0 ) , T 3 v = ( 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). v = ( 1 , 0 , 0 , 0 ) , T v = ( 0 , 1 , 0 , 0 ) , T 2 v = ( 0 , 0 , 1 , 0 ) , T 3 v = ( 0 , 0 , 0 , 1 ) .
这四个向量分别是 e 1 , e 2 , e 3 , e 4 e_1,e_2,e_3,e_4 e 1 , e 2 , e 3 , e 4 ,当然是 C 4 \mathbb{C}^4 C 4 的基。把它们倒过来 排:
B = ( T 3 v , T 2 v , T v , v ) = ( e 4 , e 3 , e 2 , e 1 ) . \mathcal{B}=(T^3v,\ T^2v,\ Tv,\ v)=(e_4,e_3,e_2,e_1). B = ( T 3 v , T 2 v , T v , v ) = ( e 4 , e 3 , e 2 , e 1 ) .
逐列算 T T T 关于 B \mathcal{B} B 的矩阵。第 j j j 列是 T T T 作用在第 j j j 个基向量上的像,用 B \mathcal{B} B 表示出来的坐标:
T ( e 4 ) = 0 T(e_4)=0 T ( e 4 ) = 0 ,第 1 列是 ( 0 , 0 , 0 , 0 ) T (0,0,0,0)^{\mathsf T} ( 0 , 0 , 0 , 0 ) T ;
T ( e 3 ) = e 4 T(e_3)=e_4 T ( e 3 ) = e 4 ,第 2 列是 ( 1 , 0 , 0 , 0 ) T (1,0,0,0)^{\mathsf T} ( 1 , 0 , 0 , 0 ) T ;
T ( e 2 ) = e 3 T(e_2)=e_3 T ( e 2 ) = e 3 ,第 3 列是 ( 0 , 1 , 0 , 0 ) T (0,1,0,0)^{\mathsf T} ( 0 , 1 , 0 , 0 ) T ;
T ( e 1 ) = e 2 T(e_1)=e_2 T ( e 1 ) = e 2 ,第 4 列是 ( 0 , 0 , 1 , 0 ) T (0,0,1,0)^{\mathsf T} ( 0 , 0 , 1 , 0 ) T 。
于是
M ( T , B ) = ( 0 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 ) . M(T,\mathcal{B})=\begin{pmatrix}0&1&0&0\\0&0&1&0\\0&0&0&1\\0&0&0&0\end{pmatrix}. M ( T , B ) = 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 0 .
对角线上全 0 0 0 ,紧挨对角线上方四个位置全是 1 1 1 ,其余全 0 0 0 。这就是一个 4 × 4 4\times4 4 × 4 的 Jordan 块 。
因为 T T T 把链上的向量往前推 :v ↦ T v ↦ T 2 v ↦ ⋯ v\mapsto Tv\mapsto T^2v\mapsto\cdots v ↦ T v ↦ T 2 v ↦ ⋯ 。 要让矩阵里出现「上一行的 1」,就得让「T T T 把第 j j j 个基向量送到第 j − 1 j-1 j − 1 个」。 顺排 B ′ = ( v , T v , T 2 v , T 3 v ) \mathcal{B}'=(v,Tv,T^2v,T^3v) B ′ = ( v , T v , T 2 v , T 3 v ) 得到的矩阵,1 1 1 全在次对角线 上——同一个算子,只是写法不合约定。
2.3 例 B:三条链,三个块 #
再看一个「一个循环向量不够用」的算子。在 C 6 \mathbb{C}^6 C 6 上定义
T ( z 1 , z 2 , z 3 , z 4 , z 5 , z 6 ) = ( 0 , z 1 , z 2 , 0 , z 4 , 0 ) . T(z_1,z_2,z_3,z_4,z_5,z_6)=(0,\ z_1,\ z_2,\ 0,\ z_4,\ 0). T ( z 1 , z 2 , z 3 , z 4 , z 5 , z 6 ) = ( 0 , z 1 , z 2 , 0 , z 4 , 0 ) .
T 3 = 0 T^3=0 T 3 = 0 。注意 z 4 , z 6 z_4,z_6 z 4 , z 6 这两个坐标从来没被送到别处,所以这个算子处理不了 6 6 6 维——不存在 u u u 使 u , T u , … , T 5 u u,Tu,\dots,T^5u u , T u , … , T 5 u 成为基(那要求 T 2 u ≠ 0 T^2u\neq0 T 2 u = 0 且这六个向量撑满空间,可 T T T 的像只有 span ( e 2 , e 3 , e 5 ) \operatorname{span}(e_2,e_3,e_5) span ( e 2 , e 3 , e 5 ) 三个方向)。
改成分头找链。取
v 1 = e 1 , v 2 = e 4 , v 3 = e 6 . v_1=e_1,\qquad v_2=e_4,\qquad v_3=e_6. v 1 = e 1 , v 2 = e 4 , v 3 = e 6 .
三条链分别是
v 1 ↦ e 2 ↦ e 3 ↦ 0 , v 2 ↦ e 5 ↦ 0 , v 3 ↦ 0. v_1\mapsto e_2\mapsto e_3\mapsto 0,\qquad v_2\mapsto e_5\mapsto 0,\qquad v_3\mapsto 0. v 1 ↦ e 2 ↦ e 3 ↦ 0 , v 2 ↦ e 5 ↦ 0 , v 3 ↦ 0.
长度依次是 3 , 2 , 1 3,2,1 3 , 2 , 1 。把每条链倒着 排,依次拼起来:
B = ( e 3 , e 2 , e 1 ⏟ 3 , e 5 , e 4 ⏟ 2 , e 6 ⏟ 1 ) . \mathcal{B}=(\underbrace{e_3,e_2,e_1}_{3},\ \underbrace{e_5,e_4}_{2},\ \underbrace{e_6}_{1}). B = ( 3 e 3 , e 2 , e 1 , 2 e 5 , e 4 , 1 e 6 ) .
逐列验:T ( e 3 ) = 0 T(e_3)=0 T ( e 3 ) = 0 、T ( e 2 ) = e 3 T(e_2)=e_3 T ( e 2 ) = e 3 、T ( e 1 ) = e 2 T(e_1)=e_2 T ( e 1 ) = e 2 、T ( e 5 ) = 0 T(e_5)=0 T ( e 5 ) = 0 、T ( e 4 ) = e 5 T(e_4)=e_5 T ( e 4 ) = e 5 、T ( e 6 ) = 0 T(e_6)=0 T ( e 6 ) = 0 。矩阵是
M ( T , B ) = ( 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 ) . 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}. M ( T , B ) = 0 0 0 0 0 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 .
它是一个 3 × 3 3\times3 3 × 3 块、一个 2 × 2 2\times2 2 × 2 块、一个 1 × 1 1\times1 1 × 1 块拼成的分块对角矩阵。
例 A 会让人误以为「n n n 维幂零算子的 Jordan 形就是一个 n × n n\times n n × n 块」。 例 B 一句话推翻它:同一个算子可以有好几条链,块可以不止一个。 链的长短由算子自己决定,选基的人没有自由。
2.4 一般情形的定义 #
现在把 λ I \lambda I λ I 加上去。
定义(Jordan 块) 设 λ ∈ F \lambda\in\mathbb{F} λ ∈ F 、d ≥ 1 d\ge1 d ≥ 1 。d × d d\times d d × d 的 Jordan 块 指
J d ( λ ) = ( λ 1 0 λ 1 ⋱ ⋱ λ 1 0 λ ) , J_d(\lambda)=\begin{pmatrix}\lambda&1&&&0\\&\lambda&1&&\\&&\ddots&\ddots&\\&&&\lambda&1\\0&&&&\lambda\end{pmatrix}, J d ( λ ) = λ 0 1 λ 1 ⋱ ⋱ λ 0 1 λ ,
即:对角线上全是 λ \lambda λ ,紧挨对角线上方那条线(超对角线 )上全是 1 1 1 ,其余全 0 0 0 。
改写一下:
J d ( λ ) = λ I d + N d , J_d(\lambda)=\lambda I_d+N_d, J d ( λ ) = λ I d + N d ,
其中 N d N_d N d 就是把例 A 那个幂零块缩到 d d d 维,满足 N d d = 0 N_d^d=0 N d d = 0 、N d d − 1 ≠ 0 N_d^{d-1}\neq0 N d d − 1 = 0 。
定义(Jordan 基、Jordan 形) 设 T ∈ L ( V ) T\in\mathcal{L}(V) T ∈ L ( V ) 。若 V V V 有一个基,使 T T T 关于该基的矩阵是分块对角矩阵
( A 1 0 ⋱ 0 A p ) , A k = J d k ( λ k ) , \begin{pmatrix}A_1&&0\\&\ddots&\\0&&A_p\end{pmatrix},\qquad A_k=J_{d_k}(\lambda_k), A 1 0 ⋱ 0 A p , A k = J d k ( λ k ) ,
就说这个基是 T T T 的 Jordan 基 ,这个矩阵是 T T T 的 Jordan 形 。
注意定义里没有任何要求说 λ 1 , … , λ p \lambda_1,\dots,\lambda_p λ 1 , … , λ p 互不相同 。同一个 λ \lambda λ 可以对应好几个块(例 B 就是 λ = 0 \lambda=0 λ = 0 配三个块)。
2.5 回看定义,三个条件各管什么 #
分块对角 :各块作用在互不相交的基向量组上,块与块之间没有耦合。这是 32 讲直和分解的矩阵版本。
每块对角元相同 :一个块只属于一个 λ \lambda λ ,也就是只落在一个 G ( λ , T ) G(\lambda,T) G ( λ , T ) 里。不同特征值的向量不可能进同一个块。
超对角线全是 1 1 1 :这是规范化的结果。1 1 1 的含义是「T T T 把这个基向量送向下一个,同时只加一个 λ \lambda λ 」,也就是 ( T − λ I ) (T-\lambda I) ( T − λ I ) 沿着链往下推一格。换成 2 2 2 行不行?行,但那是没把基选好(第三节 3.4 会亲手把 3 3 3 改成 1 1 1 )。
三、几个必须自己算一遍的结论 #
3.1 一条链的矩阵必然是超对角线上的一串 1 1 1 #
这一条第二节已经算过两遍,这里给出一般论证,因为它解释了「块」到底从哪来。
设 N N N 幂零,u ∈ V u\in V u ∈ V 满足 N m u = 0 N^mu=0 N m u = 0 、N m − 1 u ≠ 0 N^{m-1}u\neq0 N m − 1 u = 0 。则
u , N u , N 2 u , … , N m − 1 u u,\ Nu,\ N^2u,\ \dots,\ N^{m-1}u u , N u , N 2 u , … , N m − 1 u
线性无关。证明:若 ∑ i = 0 m − 1 c i N i u = 0 \sum_{i=0}^{m-1}c_iN^iu=0 ∑ i = 0 m − 1 c i N i u = 0 ,设 j j j 是最小的使 c j ≠ 0 c_j\neq0 c j = 0 的下标,两边作用 N m − 1 − j N^{m-1-j} N m − 1 − j ,得 c j N m − 1 u = 0 c_jN^{m-1}u=0 c j N m − 1 u = 0 ,与 N m − 1 u ≠ 0 N^{m-1}u\neq0 N m − 1 u = 0 矛盾。
记 U = span ( u , N u , … , N m − 1 u ) U=\operatorname{span}(u,Nu,\dots,N^{m-1}u) U = span ( u , N u , … , N m − 1 u ) ,则 dim U = m \dim U=m dim U = m ,且 U U U 在 N N N 下不变。取基
B U = ( N m − 1 u , … , N u , u ) , \mathcal{B}_U=(N^{m-1}u,\ \dots,\ Nu,\ u), B U = ( N m − 1 u , … , N u , u ) ,
即倒序。第 j j j 个基向量是 N m − j u N^{m-j}u N m − j u ,而
N ( N m − j u ) = { N m − j + 1 u , 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} N ( N m − j u ) = { N m − j + 1 u , 0 , j ≥ 2 , j = 1.
于是 j = 1 j=1 j = 1 时第 1 列是零列(N ( N m − 1 u ) = N m u = 0 N(N^{m-1}u)=N^mu=0 N ( N m − 1 u ) = N m u = 0 );j ≥ 2 j\ge2 j ≥ 2 时 N m − j + 1 u N^{m-j+1}u N m − j + 1 u 正是 B U \mathcal{B}_U B U 的第 j − 1 j-1 j − 1 个向量,第 j j j 列就是 e j − 1 e_{j-1} e j − 1 。矩阵恰是 J m ( 0 ) J_m(0) J m ( 0 ) 。
全部内容就是一句:N N N 把链上每个向量送到它的前一个,把最后一个送到 0 0 0 。 一旦这么理解,块长、1 1 1 的位置、为什么倒序,都不必背。
3.2 每个幂零算子都有 Jordan 基:证明的三步骨架 #
这是本讲的第一个核心结论。
定理 A 设 T ∈ L ( V ) T\in\mathcal{L}(V) T ∈ L ( V ) 幂零。则 V V V 有一个基是 T T T 的 Jordan 基。
对 dim V \dim V dim V 归纳。dim V = 1 \dim V=1 dim V = 1 时唯一的幂零算子是 0 0 0 ,随便取一个非零向量即可。
设 dim V > 1 \dim V>1 dim V > 1 ,且结论对所有维数更小的空间成立。令 m m m 是 T T T 的幂零指数,取 u u u 使 T m − 1 u ≠ 0 T^{m-1}u\neq0 T m − 1 u = 0 ,令
U = span ( u , T u , … , T m − 1 u ) . U=\operatorname{span}(u,Tu,\dots,T^{m-1}u). U = span ( u , T u , … , T m − 1 u ) .
U U U 在 T T T 下不变,dim U = m \dim U=m dim U = m 。若 U = V U=V U = V ,由 3.1 就得到 Jordan 基,证完。
难点在 U ≠ V U\neq V U = V 的情形。思路 :找到另一个不变子空间 W W W ,使 V = U ⊕ W V=U\oplus W V = U ⊕ W ;然后对 U U U 、W W W 分别用归纳假设,两组基拼起来就是 V V V 的 Jordan 基。
怎么找 W W W ?关键一步是用线性泛函切一刀 。取 φ ∈ V ′ \varphi\in V' φ ∈ V ′ 使 φ ( T m − 1 u ) ≠ 0 \varphi(T^{m-1}u)\neq0 φ ( T m − 1 u ) = 0 (存在:T m − 1 u ≠ 0 T^{m-1}u\neq0 T m − 1 u = 0 ,可以把它扩充成基,再用对偶基取泛函)。令
W = { v ∈ V : φ ( T k v ) = 0 对每个 k = 0 , 1 , … , m − 1 } . W=\big\{\,v\in V:\ \varphi(T^kv)=0\ \text{对每个}\ k=0,1,\dots,m-1\ \big\}. W = { v ∈ V : φ ( T k v ) = 0 对每个 k = 0 , 1 , … , m − 1 } .
W W W 在 T T T 下不变 :若 v ∈ W v\in W v ∈ W ,则对 k = 0 , … , m − 2 k=0,\dots,m-2 k = 0 , … , m − 2 有 φ ( T k ( T v ) ) = φ ( T k + 1 v ) = 0 \varphi\big(T^k(Tv)\big)=\varphi(T^{k+1}v)=0 φ ( T k ( T v ) ) = φ ( T k + 1 v ) = 0 ;对 k = m − 1 k=m-1 k = m − 1 有 φ ( T m − 1 ( T v ) ) = φ ( T m v ) = 0 \varphi\big(T^{m-1}(Tv)\big)=\varphi(T^mv)=0 φ ( T m − 1 ( T v ) ) = φ ( T m v ) = 0 (因为 T m = 0 T^m=0 T m = 0 )。所以 T v ∈ W Tv\in W T v ∈ W 。
U ∩ W = { 0 } U\cap W=\{0\} U ∩ W = { 0 } :设 v ∈ U ∩ W v\in U\cap W v ∈ U ∩ W 且 v ≠ 0 v\neq0 v = 0 。写
v = c 0 u + c 1 T u + ⋯ + c m − 1 T m − 1 u , v=c_0u+c_1Tu+\cdots+c_{m-1}T^{m-1}u, v = c 0 u + c 1 T u + ⋯ + c m − 1 T m − 1 u ,
令 j j j 是最小的使 c j ≠ 0 c_j\neq0 c j = 0 的下标。两边作用 T m − j − 1 T^{m-j-1} T m − j − 1 :
T m − j − 1 v = c j T m − 1 u , T^{m-j-1}v=c_jT^{m-1}u, T m − j − 1 v = c j T m − 1 u ,
因为次数更高的项都被 T m = 0 T^m=0 T m = 0 杀掉了(T m − j − 1 ⋅ T i u = T m − 1 + ( i − j ) u T^{m-j-1}\cdot T^iu=T^{m-1+(i-j)}u T m − j − 1 ⋅ T i u = T m − 1 + ( i − j ) u ,当 i > j i>j i > j 时次数 ≥ m \ge m ≥ m ,为零)。再作用 φ \varphi φ :
φ ( T m − j − 1 v ) = c j φ ( T m − 1 u ) ≠ 0 , \varphi\big(T^{m-j-1}v\big)=c_j\,\varphi\big(T^{m-1}u\big)\neq0, φ ( T m − j − 1 v ) = c j φ ( T m − 1 u ) = 0 ,
所以 v ∉ W v\notin W v ∈ / W ,矛盾。故 U ∩ W = { 0 } U\cap W=\{0\} U ∩ W = { 0 } 。
U + W = V U+W=V U + W = V :定义 S : V → F m S:V\to\mathbb{F}^m S : V → F m ,
S v = ( φ ( v ) , φ ( T v ) , … , φ ( T m − 1 v ) ) . Sv=\big(\varphi(v),\ \varphi(Tv),\ \dots,\ \varphi(T^{m-1}v)\big). S v = ( φ ( v ) , φ ( T v ) , … , φ ( T m − 1 v ) ) .
S S S 线性,null S = W \operatorname{null}S=W 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, dim W = dim null S = dim V − dim range S ≥ dim V − m ,
因为 range S ⊆ F m \operatorname{range}S\subseteq\mathbb{F}^m range S ⊆ F m ,维数至多 m m m 。于是
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, dim ( U + W ) = dim U + dim W ≥ m + ( dim V − m ) = dim V ,
故 U + W = V U+W=V U + W = V ;与 U ∩ W = { 0 } U\cap W=\{0\} U ∩ W = { 0 } 合起来得 V = U ⊕ W V=U\oplus W V = U ⊕ W 。
U U U 在 T T T 下不变、dim U = m < dim V \dim U=m<\dim V dim U = m < dim V ,由归纳假设 U U U 有 T ∣ U T|_U T ∣ U 的 Jordan 基;W W W 在 T T T 下不变、dim W < dim V \dim W<\dim V dim W < dim V ,同理也有。两组基合并就是 V V V 的基,T T T 关于它的矩阵是两块分块对角,每块都是 Jordan 块。证毕。
归纳法的难点从来不是「怎么归纳」,而是「补空间怎么造」。 这里的构造值得记住:用泛函把链的最高次项「测」出来,再把测不出东西的向量收进 W W W 。 W W W 之所以够大,靠的是「S S S 的值域最多 m m m 维」这个上界。
3.3 从幂零推广到一般 #
定理 B(存在性) 设 F = C \mathbb{F}=\mathbb{C} F = C 、T ∈ L ( V ) T\in\mathcal{L}(V) T ∈ L ( V ) 。则 V V V 有 T T T 的 Jordan 基。
由 32 讲的广义特征空间分解,
V = G ( λ 1 , T ) ⊕ ⋯ ⊕ G ( λ m , T ) , V=G(\lambda_1,T)\oplus\cdots\oplus G(\lambda_m,T), V = G ( λ 1 , T ) ⊕ ⋯ ⊕ G ( λ m , T ) ,
每块在 T T T 下不变,且每个 ( T − λ k I ) ∣ G ( λ k , T ) (T-\lambda_kI)|_{G(\lambda_k,T)} ( T − λ k I ) ∣ G ( λ k , T ) 幂零。在 G ( λ k , T ) G(\lambda_k,T) G ( λ k , T ) 上写成
T ∣ G ( λ k , T ) = λ k I + N k , T\big|_{G(\lambda_k,T)}=\lambda_kI+N_k, T G ( λ k , T ) = λ k I + N k ,
由定理 A,存在 G ( λ k , T ) G(\lambda_k,T) G ( λ k , T ) 的基使 N k N_k N k 的矩阵是若干块 J d ( 0 ) J_d(0) J d ( 0 ) 。同一组基下 λ k I \lambda_kI λ k I 的矩阵就是 λ k I \lambda_kI λ k I ,两者相加(同一个基,逐项相加)把每个块的 λ \lambda λ 补上:
J d ( 0 ) + λ k I = J d ( λ k ) . J_d(0)+\lambda_kI=J_d(\lambda_k). J d ( 0 ) + λ k I = J d ( λ k ) .
把各 G ( λ k , T ) G(\lambda_k,T) G ( λ k , T ) 的基合并,就得到 V V V 的 Jordan 基。
注意定理 B 里 F = C \mathbb{F}=\mathbb{C} F = C 是硬条件。 定理 A 对实空间也成立(证明只用了 T m = 0 T^m=0 T m = 0 、线性泛函、维数估计),但定理 B 不行:R \mathbb{R} R 上不保证特征值存在,广义特征空间分解就无从谈起。R 2 \mathbb{R}^2 R 2 上转 90 ° 90° 90° 就是反例(见第五节误区 1)。
3.4 亲手把一个 3 3 3 改成 1 1 1 #
回到第一节那个矩阵
A = ( 6 3 4 0 6 2 0 0 7 ) , A=\begin{pmatrix}6&3&4\\0&6&2\\0&0&7\end{pmatrix}, A = 6 0 0 3 6 0 4 2 7 ,
它定义的算子是
T ( z 1 , z 2 , z 3 ) = ( 6 z 1 + 3 z 2 + 4 z 3 , 6 z 2 + 2 z 3 , 7 z 3 ) , T(z_1,z_2,z_3)=(6z_1+3z_2+4z_3,\ 6z_2+2z_3,\ 7z_3), T ( z 1 , z 2 , z 3 ) = ( 6 z 1 + 3 z 2 + 4 z 3 , 6 z 2 + 2 z 3 , 7 z 3 ) ,
特征值 6 , 6 , 7 6,6,7 6 , 6 , 7 。
先看 G ( 7 , T ) 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), T ( 10 , 2 , 1 ) = ( 70 , 14 , 7 ) = 7 ( 10 , 2 , 1 ) ,
而 7 7 7 的重数是 1 1 1 ,所以 G ( 7 , T ) = span ( 10 , 2 , 1 ) G(7,T)=\operatorname{span}(10,2,1) G ( 7 , T ) = span ( 10 , 2 , 1 ) 。
再看 G ( 6 , T ) G(6,T) G ( 6 , T ) 。令 z 3 = 0 z_3=0 z 3 = 0 :
T ( z 1 , z 2 , 0 ) = ( 6 z 1 + 3 z 2 , 6 z 2 , 0 ) , T(z_1,z_2,0)=(6z_1+3z_2,\ 6z_2,\ 0), T ( z 1 , z 2 , 0 ) = ( 6 z 1 + 3 z 2 , 6 z 2 , 0 ) ,
可见 { ( z 1 , z 2 , 0 ) } \{(z_1,z_2,0)\} {( z 1 , z 2 , 0 )} 在 T T T 下不变、维数 2 2 2 ,就是 G ( 6 , T ) G(6,T) G ( 6 , T ) 。
在这个二维块上,T T T 关于 ( e 1 , e 2 ) (e_1,e_2) ( e 1 , e 2 ) 的矩阵是
( 6 3 0 6 ) , \begin{pmatrix}6&3\\0&6\end{pmatrix}, ( 6 0 3 6 ) ,
也就是
T ( e 1 ) = 6 e 1 , T ( e 2 ) = 3 e 1 + 6 e 2 . T(e_1)=6e_1,\qquad T(e_2)=3e_1+6e_2. T ( e 1 ) = 6 e 1 , T ( e 2 ) = 3 e 1 + 6 e 2 .
令 N = T − 6 I N=T-6I N = T − 6 I ,则 N ( e 1 ) = 0 N(e_1)=0 N ( e 1 ) = 0 、N ( e 2 ) = 3 e 1 N(e_2)=3e_1 N ( e 2 ) = 3 e 1 、N 2 = 0 N^2=0 N 2 = 0 。要按 3.1 的套路造链,不能从 e 1 e_1 e 1 起(它已经在 N N N 的零空间里,链长只有 1 1 1 ),要从 N N N 作用下非零 的向量起。取 u = e 2 u=e_2 u = e 2 :
N u = 3 e 1 ≠ 0 , N 2 u = 0. Nu=3e_1\neq0,\qquad N^2u=0. N u = 3 e 1 = 0 , N 2 u = 0.
链是 e 2 ↦ 3 e 1 ↦ 0 e_2\mapsto3e_1\mapsto0 e 2 ↦ 3 e 1 ↦ 0 ,长 2 2 2 。倒着排得基
( N u , u ) = ( 3 e 1 , e 2 ) = ( ( 3 , 0 , 0 ) , ( 0 , 1 , 0 ) ) . \big(Nu,\ u\big)=\big(3e_1,\ e_2\big)=\big((3,0,0),\ (0,1,0)\big). ( N u , u ) = ( 3 e 1 , e 2 ) = ( ( 3 , 0 , 0 ) , ( 0 , 1 , 0 ) ) .
验算:
T ( 3 e 1 ) = 18 e 1 = 6 ⋅ ( 3 e 1 ) , T ( e 2 ) = 3 e 1 + 6 e 2 = ( 3 e 1 ) + 6 e 2 . T(3e_1)=18e_1=6\cdot(3e_1),\qquad T(e_2)=3e_1+6e_2=(3e_1)+6\,e_2. T ( 3 e 1 ) = 18 e 1 = 6 ⋅ ( 3 e 1 ) , T ( e 2 ) = 3 e 1 + 6 e 2 = ( 3 e 1 ) + 6 e 2 .
关于基 ( ( 3 , 0 , 0 ) , ( 0 , 1 , 0 ) , ( 10 , 2 , 1 ) ) \big((3,0,0),\ (0,1,0),\ (10,2,1)\big) ( ( 3 , 0 , 0 ) , ( 0 , 1 , 0 ) , ( 10 , 2 , 1 ) ) ,
M ( T ) = ( 6 1 0 0 6 0 0 0 7 ) = diag ( J 2 ( 6 ) , J 1 ( 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). M ( T ) = 6 0 0 1 6 0 0 0 7 = diag ( J 2 ( 6 ) , J 1 ( 7 ) ) .
那个讨厌的 3 3 3 变成了 1 1 1 。它不是被算掉的,是被基 3 e 1 3e_1 3 e 1 吃掉的 ——把第一个基向量放缩 3 3 3 倍,非对角元就跟着除以 3 3 3 。
块内超对角线上的非零数字 都只是标度,换基就能改成 1 1 1 。真正消不掉的是能不能连成链 : 若某个位置上的耦合是 0 0 0 ,那就是块与块的接缝,再怎么换基也不会变出 1 1 1 来。所以 Jordan 形要回答的问题不是「数字是多少」,而是「块怎么切」。
3.5 存在之后:块长分布是不是被算子唯一决定 #
必须坦白说清楚:上面证的是「存在」,不是「唯一」。
一个算子的 Jordan 形当然不可能字面上唯一——把块换个顺序排,矩阵变了,但仍然分块对角、每块仍是 Jordan 块,还是同一个算子关于另一组基的矩阵。
真正的问题是:块长分布是不是被算子唯一决定? 也就是「6 6 6 对应一个 2 × 2 2\times2 2 × 2 块加一个 1 × 1 1\times1 1 × 1 块」这件事,会不会因为换基变成「三个 1 × 1 1\times1 1 × 1 块」?
答案是「不会」(否则这套理论就没意义)。要证这一点,得再加一段关于 dim null ( T − λ I ) k \dim\operatorname{null}(T-\lambda I)^k dim null ( T − λ I ) k 的分析——下面 3.5 末就把它补完 。
不过,从已有的证明里能免费读出三个必然不变 的量:
量
等于什么
为什么
λ \lambda λ 的块数
dim E ( λ , T ) \dim E(\lambda,T) dim E ( λ , T ) (几何重数)
每条链的末端都落在 null ( T − λ I ) \operatorname{null}(T-\lambda I) null ( T − λ I ) 里,且这些末端恰好构成它的基
该 λ \lambda λ 所有块长之和
dim G ( λ , T ) \dim G(\lambda,T) dim G ( λ , T ) (代数重数)
这些链上的向量全体构成 G ( λ , T ) G(\lambda,T) G ( λ , T ) 的基
最长块的长度
极小多项式中 ( z − λ ) (z-\lambda) ( z − λ ) 的指数
最长链要走到第 d d d 步才归零,d d d 就是使 ( T − λ I ) d (T-\lambda I)^d ( T − λ I ) d 在 G G G 上为零的最小 d d d
在幂零情形下,这三条退化成:块数 = dim null T =\dim\operatorname{null}T = dim null T ;块长之和 = dim V =\dim V = dim V ;最长块长 = = = 幂零指数。
举例:dim V = 8 \dim V=8 dim V = 8 、块数 3 3 3 、最长块 4 4 4 。 分布可以是 4 + 3 + 1 4+3+1 4 + 3 + 1 ,也可以是 4 + 2 + 2 4+2+2 4 + 2 + 2 ,两者都满足上面三条。 它们确实是不同的算子——一个可亲手验证的区分是 dim null T 2 \dim\operatorname{null}T^2 dim null T 2 : 对 4 + 3 + 1 4+3+1 4 + 3 + 1 它是 2 + 2 + 1 = 5 2+2+1=5 2 + 2 + 1 = 5 ,对 4 + 2 + 2 4+2+2 4 + 2 + 2 它是 2 + 2 + 2 = 6 2+2+2=6 2 + 2 + 2 = 6 。所以「块数 + 总长 + 最长块」只是必要条件,不是完整的不变量。
补全:唯一性的完整判据 #
上面那个反例反而指了路——要区分块长分布,得看更高次的零空间 。事实上只要把所有 k k k 次的零空间维数都看一遍,分布就被完全钉死了。
记 N = T − λ I N=T-\lambda I N = T − λ I (把它限制在 G ( λ , T ) G(\lambda,T) G ( λ , T ) 上看)。设长度为 s s s 的链(即对应 s × s s\times s s × s 的 Jordan 块)有 c s c_s c s 条。核心观察是:
dim null N k = ∑ s ≥ 1 min ( k , s ) c s \dim\operatorname{null}N^k=\sum_{s\ge 1}\min(k,s)\,c_s dim null N k = s ≥ 1 ∑ min ( k , s ) c s
为什么 :一条长度为 s s s 的链上,N k N^k N k 把它打掉前 min ( k , s ) \min(k,s) min ( k , s ) 个向量(链走 s s s 步归零,走到第 k k k 步就归零了 k k k 个)。各条链独立,直接相加。
于是相邻两次的差
dim null N k − dim null N k − 1 = # { s ≥ k } = 「长度 ≥ k 的链有多少条」 \dim\operatorname{null}N^{k}-\dim\operatorname{null}N^{k-1}=\#\{s\ge k\}=\text{「长度}\ge k\text{ 的链有多少条」} dim null N k − dim null N k − 1 = # { s ≥ k } = 「长度 ≥ k 的链有多少条」
这就把分布完全确定了 :设 b k : = # { s ≥ k } b_k:=\#\{s\ge k\} b k := # { s ≥ k } ,则长度恰好为 k k k 的链数为 c k = b k − b k + 1 c_k=b_k-b_{k+1} c k = b k − b k + 1 。b k b_k b k 全部由算子的零空间维数算出,与基的选择无关。所以块长分布是算子的不变量 。
回到 4 + 3 + 1 4+3+1 4 + 3 + 1 与 4 + 2 + 2 4+2+2 4 + 2 + 2 (dim V = 8 \dim V=8 dim V = 8 ):
k k k
b k = # { s ≥ k } b_k=\#\{s\ge k\} b k = # { s ≥ k } (4+3+1)
dim null N k \dim\operatorname{null}N^k dim null N k
b k b_k b k (4+2+2)
dim null N k \dim\operatorname{null}N^k dim null N k
1
3
3
3
3
2
2
5
3
6
3
2
7
1
7
4
1
8
1
8
两列在 k = 2 k=2 k = 2 处首次分岔(5 5 5 vs 6 6 6 ),正是前面提到的 dim null T 2 \dim\operatorname{null}T^2 dim null T 2 。 而把每行的差算出来(b k − b k + 1 b_k-b_{k+1} b k − b k + 1 ),分别得 c 1 , c 2 , c 3 , c 4 c_1,c_2,c_3,c_4 c 1 , c 2 , c 3 , c 4 :
4 + 3 + 1 4+3+1 4 + 3 + 1 :c = ( 1 , 0 , 1 , 1 ) c=(1,0,1,1) c = ( 1 , 0 , 1 , 1 ) → 链长 4 , 3 , 1 4,3,1 4 , 3 , 1 ✓
4 + 2 + 2 4+2+2 4 + 2 + 2 :c = ( 0 , 2 , 0 , 1 ) c=(0,2,0,1) c = ( 0 , 2 , 0 , 1 ) → 链长 4 , 2 , 2 4,2,2 4 , 2 , 2 ✓
分布因此被唯一还原。
在不计块的排列顺序 的意义下,Jordan 形由算子唯一决定。 换个说法:两个算子相似 ⟺ \iff ⟺ 它们有相同的 Jordan 形(每个特征值的块长分布都相同)。 这条判据把「相似」这个看似无从下手的等价关系,变成了一次可以算完的核对。
注意 :本讲前面(3.2/3.3)证的存在性、加上这里证的唯一性,合起来才是完整的 Jordan 形定理。本节之前只做到存在性,现在补齐了。
有些教材不用「块长分布」这个说法,而说「初等因子」:特征值 λ \lambda λ 配上一个长度 s s s 的块,就记一个因子 ( z − λ ) s (z-\lambda)^s ( z − λ ) s 。 于是「Jordan 形唯一」等价于「初等因子组唯一」。两种说法是同一件事的两种记账方式。下面把这条记账法做成可操作的工具。
3.6 初等因子:从特征矩阵直接算出来 #
本节约定 F = C \mathbb{F}=\mathbb{C} F = C (或至少特征值全在 F \mathbb{F} F 内)。这样才能把初等因子写成 ( λ − λ 0 ) s (\lambda-\lambda_0)^s ( λ − λ 0 ) s 并对应到 Jordan 块;一般域上的情形见 3.7。
前面判断相似,要逐个 k k k 算 dim null N k \dim\operatorname{null}N^k dim null N k 。实际做题时有一条纯代数 的路:不必碰基,只算多项式。
特征矩阵 。设 A ∈ F n , n A\in\mathbb{F}^{n,n} A ∈ F n , n ,把 λ \lambda λ 当不定元,A A A 的特征矩阵 是 λ I − A \lambda I-A λ I − A (元素都是 F [ λ ] \mathbb{F}[\lambda] F [ λ ] 里的多项式)。
行列式因子 。对 k = 1 , … , n k=1,\dots,n k = 1 , … , n ,定义 D k ( λ ) D_k(\lambda) D k ( λ ) 为 λ I − A \lambda I-A λ I − A 的全部 k k k 阶子式的最大公因式 (首一化)。约定 D 0 = 1 D_0=1 D 0 = 1 。
不变因子 。令
d k ( λ ) = D k ( λ ) D k − 1 ( λ ) ( k = 1 , … , n ) d_k(\lambda)=\frac{D_k(\lambda)}{D_{k-1}(\lambda)}\qquad(k=1,\dots,n) d k ( λ ) = D k − 1 ( λ ) D k ( λ ) ( k = 1 , … , n )
则 d 1 ∣ d 2 ∣ ⋯ ∣ d n d_1\mid d_2\mid\cdots\mid d_n d 1 ∣ d 2 ∣ ⋯ ∣ d n ,且 d n d_n d n 就是极小多项式 ,而 ∏ k = 1 n d k = det ( λ I − A ) \prod_{k=1}^{n}d_k=\det(\lambda I-A) ∏ k = 1 n d k = det ( λ I − A ) 是特征多项式 。
[!note] 这三条不是「显然」 整除链 d 1 ∣ ⋯ ∣ d n d_1\mid\cdots\mid d_n d 1 ∣ ⋯ ∣ d n 、d n = d_n= d n = 极小多项式、∏ d k = det \prod d_k=\det ∏ d k = det ——这三条都是 Smith 标准形定理 的内容, 严格证明要用 λ \lambda λ -矩阵的初等变换理论。本节把它当已知结论用,只演示「怎么算」。
初等因子 。把每个 d k d_k d k 分解成不可约因式的幂,把所有非平凡的幂次因子收集起来——这就是初等因子组。
A ∼ B ⟺ λ I − A 与 λ I − B 有相同的初等因子组 A\sim B\iff \lambda I-A\ \text{与}\ \lambda I-B\ \text{有相同的初等因子组} A ∼ B ⟺ λ I − A 与 λ I − B 有相同的初等因子组
更简洁的版本:λ I − A \lambda I-A λ I − A 与 λ I − B \lambda I-B λ I − B 等价 (即能通过初等行/列变换互变)。 这条路不需要求特征向量、不需要算零空间维数,只要会算行列式与最大公因式。
初等因子 ↔ Jordan 块 ,对应关系很简单:
初等因子
Jordan 块
( λ − λ 0 ) s (\lambda-\lambda_0)^s ( λ − λ 0 ) s
一个 s × s s\times s s × s 的 J s ( λ 0 ) J_s(\lambda_0) J s ( λ 0 )
λ − λ 0 \lambda-\lambda_0 λ − λ 0 (即 s = 1 s=1 s = 1 )
一个 1 × 1 1\times1 1 × 1 块
一个完整的手算例 。设
A = ( 2 1 0 0 2 1 0 0 2 ) A=\begin{pmatrix}2&1&0\\0&2&1\\0&0&2\end{pmatrix} A = 2 0 0 1 2 0 0 1 2
第一步 ,写特征矩阵:
λ I − A = ( λ − 2 − 1 0 0 λ − 2 − 1 0 0 λ − 2 ) \lambda I-A=\begin{pmatrix}\lambda-2&-1&0\\0&\lambda-2&-1\\0&0&\lambda-2\end{pmatrix} λ I − A = λ − 2 0 0 − 1 λ − 2 0 0 − 1 λ − 2
第二步 ,算行列式因子。D 3 = det ( λ I − A ) = ( λ − 2 ) 3 D_3=\det(\lambda I-A)=(\lambda-2)^3 D 3 = det ( λ I − A ) = ( λ − 2 ) 3 。D 2 D_2 D 2 是全部二阶子式的最大公因式。3 × 3 3\times3 3 × 3 矩阵的二阶子式有 ( 3 2 ) ( 3 2 ) = 9 \binom32\binom32=9 ( 2 3 ) ( 2 3 ) = 9 个,但只需抽出三个就能定案(三个都取自第 1 , 2 1,2 1 , 2 行 ):
∣ λ − 2 − 1 0 λ − 2 ∣ = ( λ − 2 ) 2 ( 第 1 , 2 列 ) , ∣ λ − 2 0 0 − 1 ∣ = − ( λ − 2 ) ( 第 1 , 3 列 ) , ∣ − 1 0 λ − 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列). λ − 2 0 − 1 λ − 2 = ( λ − 2 ) 2 ( 第 1 , 2 列 ) , λ − 2 0 0 − 1 = − ( λ − 2 ) ( 第 1 , 3 列 ) , − 1 λ − 2 0 − 1 = 1 ( 第 2 , 3 列 ) .
最大公因式是 1 1 1 (因为出现了 1 1 1 ),所以 D 2 = 1 D_2=1 D 2 = 1 ;同理 D 1 = 1 D_1=1 D 1 = 1 (有元素 1 1 1 )。
第三步 ,不变因子:d 1 = D 1 / D 0 = 1 d_1=D_1/D_0=1 d 1 = D 1 / D 0 = 1 ,d 2 = D 2 / D 1 = 1 d_2=D_2/D_1=1 d 2 = D 2 / D 1 = 1 ,d 3 = D 3 / D 2 = ( λ − 2 ) 3 d_3=D_3/D_2=(\lambda-2)^3 d 3 = D 3 / D 2 = ( λ − 2 ) 3 。
第四步 ,初等因子组:只有 { ( λ − 2 ) 3 } \{(\lambda-2)^3\} {( λ − 2 ) 3 } 一个。
结论 :A A A 的 Jordan 形就是一个 3 × 3 3\times3 3 × 3 块 J 3 ( 2 ) J_3(2) J 3 ( 2 ) ——它本身就已是 Jordan 形,与自己相似,不稀奇。
换一个更有意思的例子 。设
B = ( 2 1 0 0 2 0 0 0 2 ) B=\begin{pmatrix}2&1&0\\0&2&0\\0&0&2\end{pmatrix} B = 2 0 0 1 2 0 0 0 2
D 3 = ( λ − 2 ) 3 D_3=(\lambda-2)^3 D 3 = ( λ − 2 ) 3 不变。D 2 D_2 D 2 是全部二阶子式的最大公因式。3 × 3 3\times3 3 × 3 矩阵有 ( 3 2 ) ( 3 2 ) = 9 \binom32\binom32=9 ( 2 3 ) ( 2 3 ) = 9 个二阶子式,逐个枚举后,非零的只有四个 :
取的行
取的列
子式
值
1 , 2 1,2 1 , 2
1 , 2 1,2 1 , 2
∣ λ − 2 − 1 0 λ − 2 ∣ \begin{vmatrix}\lambda-2&-1\\0&\lambda-2\end{vmatrix} λ − 2 0 − 1 λ − 2
( λ − 2 ) 2 (\lambda-2)^2 ( λ − 2 ) 2
1 , 3 1,3 1 , 3
1 , 3 1,3 1 , 3
∣ λ − 2 0 0 λ − 2 ∣ \begin{vmatrix}\lambda-2&0\\0&\lambda-2\end{vmatrix} λ − 2 0 0 λ − 2
( λ − 2 ) 2 (\lambda-2)^2 ( λ − 2 ) 2
2 , 3 2,3 2 , 3
2 , 3 2,3 2 , 3
∣ λ − 2 0 0 λ − 2 ∣ \begin{vmatrix}\lambda-2&0\\0&\lambda-2\end{vmatrix} λ − 2 0 0 λ − 2
( λ − 2 ) 2 (\lambda-2)^2 ( λ − 2 ) 2
1 , 3 1,3 1 , 3
2 , 3 2,3 2 , 3
∣ − 1 0 0 λ − 2 ∣ \begin{vmatrix}-1&0\\0&\lambda-2\end{vmatrix} − 1 0 0 λ − 2
− ( λ − 2 ) -(\lambda-2) − ( λ − 2 )
其余 5 5 5 个全为 0 0 0 (例如第 1 , 2 1,2 1 , 2 行第 1 , 3 1,3 1 , 3 列 ∣ λ − 2 0 0 0 ∣ = 0 \begin{vmatrix}\lambda-2&0\\0&0\end{vmatrix}=0 λ − 2 0 0 0 = 0 )。
于是非零项恰是 ( λ − 2 ) 2 (\lambda-2)^2 ( λ − 2 ) 2 (三个)与 − ( λ − 2 ) -(\lambda-2) − ( λ − 2 ) (一个),取公因式:
D 2 = λ − 2. D_2=\lambda-2 . D 2 = λ − 2.
这里最容易错的一步是 :因为出现了 − ( λ − 2 ) -(\lambda-2) − ( λ − 2 ) 这个一次因子,而不是常数,所以公因式是 λ − 2 \lambda-2 λ − 2 ,不是 1 1 1 。正是这个 D 2 D_2 D 2 把 B B B 与上面的 A A A 区分开了——A A A 的 D 2 D_2 D 2 是 1 1 1 。
于是
d 1 = 1 , d 2 = λ − 2 1 = λ − 2 , d 3 = ( λ − 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 . d 1 = 1 , d 2 = 1 λ − 2 = λ − 2 , d 3 = λ − 2 ( λ − 2 ) 3 = ( λ − 2 ) 2 .
初等因子组:{ ( λ − 2 ) , ( λ − 2 ) 2 } \{(\lambda-2),\ (\lambda-2)^2\} {( λ − 2 ) , ( λ − 2 ) 2 } ——两个块 ,尺寸 1 1 1 与 2 2 2 。
与 A A A 对比:A A A 与 B B B 的特征多项式都是 ( λ − 2 ) 3 (\lambda-2)^3 ( λ − 2 ) 3 、特征值都是 2 2 2 ,但两者不相似 (A A A 是一个 3 × 3 3\times3 3 × 3 块,B B B 是 2 + 1 2+1 2 + 1 两个块)。注意几何重数已经不同 :dim null ( A − 2 I ) = 1 \dim\operatorname{null}(A-2I)=1 dim null ( A − 2 I ) = 1 ,而 dim null ( B − 2 I ) = 2 \dim\operatorname{null}(B-2I)=2 dim null ( B − 2 I ) = 2 (因为 rank ( B − 2 I ) = 1 \operatorname{rank}(B-2I)=1 rank ( B − 2 I ) = 1 )——这一对其实用几何重数就能分开。差别也体现在极小多项式:A A A 的是 ( λ − 2 ) 3 (\lambda-2)^3 ( λ − 2 ) 3 ,B B B 的是 ( λ − 2 ) 2 (\lambda-2)^2 ( λ − 2 ) 2 ;以及 D 2 D_2 D 2 :A A A 的 D 2 = 1 D_2=1 D 2 = 1 ,B B B 的 D 2 = λ − 2 D_2=\lambda-2 D 2 = λ − 2 。
上面这一对用极小多项式分开了,容易让人以为「特征多项式 + + + 极小多项式」就够。 反例见 3.8:J 2 ( 2 ) ⊕ J 2 ( 2 ) J_2(2)\oplus J_2(2) J 2 ( 2 ) ⊕ J 2 ( 2 ) 与 J 2 ( 2 ) ⊕ J 1 ( 2 ) ⊕ J 1 ( 2 ) J_2(2)\oplus J_1(2)\oplus J_1(2) J 2 ( 2 ) ⊕ J 1 ( 2 ) ⊕ J 1 ( 2 ) 的特征多项式、极小多项式、最长块长全都一样 ,仍不相似。 差别只在更有量的东西上:前者一个 2 2 2 阶块,后者一个 2 2 2 阶块和两个 1 1 1 阶块。这个差别在不变因子里看得见(见 3.6 末),在「dim null ( B − 2 I ) \dim\operatorname{null}(B-2I) dim null ( B − 2 I ) 」里也看得见(2 2 2 vs 3 3 3 )。
一、D k D_k D k 是最大公因式,不是子式的乘积。 B B B 的非零二阶子式有三个都是 ( λ − 2 ) 2 (\lambda-2)^2 ( λ − 2 ) 2 、一个 − ( λ − 2 ) -(\lambda-2) − ( λ − 2 ) ;取公因式得 λ − 2 \lambda-2 λ − 2 ,不能把 ( λ − 2 ) 2 (\lambda-2)^2 ( λ − 2 ) 2 当答案。二、D k D_k D k 永远不会是零多项式。 因为 det ( λ I − A ) ≠ 0 \det(\lambda I-A)\neq0 det ( λ I − A ) = 0 ,λ I − A \lambda I-A λ I − A 在 F ( λ ) \mathbb{F}(\lambda) F ( λ ) 上可逆,秩恒为 n n n ,故至少有一个 k k k 阶子式非零。所以「D k = 0 D_k=0 D k = 0 」这种情形根本不会发生——真算出 0 0 0 一定是算错了。三、λ I − A \lambda I-A λ I − A 的秩不必等于 A A A 的秩。 把 λ \lambda λ 取成某个特征值时会掉秩,但那只是某一个 λ \lambda λ 处的取值,与 D k D_k D k 作为多项式的公因式不是一回事。B B B 在 λ = 2 \lambda=2 λ = 2 处 λ I − B = 2 I − B \lambda I-B=2I-B λ I − B = 2 I − B 的秩是 1 1 1 (零空间 2 2 2 维),可 D 2 D_2 D 2 仍是 λ − 2 \lambda-2 λ − 2 而不是 0 0 0 。
3.7 有理标准形:不依赖代数闭域的那种标准形 #
Jordan 形有个前提:特征值必须在 F \mathbb{F} F 里 。实矩阵若有复特征值(比如 ( 0 − 1 1 0 ) \begin{pmatrix}0&-1\\1&0\end{pmatrix} ( 0 1 − 1 0 ) ),在 R \mathbb{R} R 上就没有 Jordan 形。
要一个「任何域上都有」的标准形,就得把块从 Jordan 块换成伴侣矩阵 ,并允许素因子是不可约多项式 (而不只是一次因子)。这就是有理标准形 (也叫 Frobenius 标准形)。
伴侣矩阵 。给定首一多项式 p ( z ) = a 0 + a 1 z + ⋯ + a m − 1 z m − 1 + z m p(z)=a_0+a_1z+\cdots+a_{m-1}z^{m-1}+z^m p ( z ) = a 0 + a 1 z + ⋯ + a m − 1 z m − 1 + z m ,它的伴侣矩阵 是
C ( p ) = ( 0 0 ⋯ 0 − a 0 1 0 ⋯ 0 − a 1 0 1 ⋯ 0 − a 2 ⋮ ⋱ ⋮ 0 0 ⋯ 1 − a m − 1 ) m × m C(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 ) = 0 1 0 ⋮ 0 0 0 1 0 ⋯ ⋯ ⋯ ⋱ ⋯ 0 0 0 1 − a 0 − a 1 − a 2 ⋮ − a m − 1 m × m
关键性质 :C ( p ) C(p) C ( p ) 的特征多项式与极小多项式都等于 p p p 。这正好反过来印证 18 讲那句话——「每个首一多项式都是某个算子的极小多项式」,那个算子就是伴侣矩阵。
有理标准形定理 。任何 A ∈ F n , n A\in\mathbb{F}^{n,n} A ∈ F n , n 都相似于若干个伴侣矩阵的块对角和:
A ∼ diag ( C ( p 1 ) , C ( p 2 ) , … , C ( p r ) ) . A\ \sim\ \operatorname{diag}\bigl(C(p_1),\,C(p_2),\,\dots,\,C(p_r)\bigr). A ∼ diag ( C ( p 1 ) , C ( p 2 ) , … , C ( p r ) ) .
这里有一个必须说清的岔路:p i p_i p i 有两种取法,对应两种不同 的写法。
初等因子形式 :p i p_i p i 取全部初等因子(每个素幂配一个块,同一素多项式可以有多个块)。这与 3.6 的初等因子组一一对应。
不变因子形式 :p i p_i p i 取全部非常数不变因子 ,此时 p 1 ∣ p 2 ∣ ⋯ ∣ p r p_1\mid p_2\mid\cdots\mid p_r p 1 ∣ p 2 ∣ ⋯ ∣ p r ,块数等于非平凡不变因子的个数。
两种写法只差把块重排、合并(用中国剩余定理:C ( p q ) ∼ diag ( C ( p ) , C ( q ) ) C(pq)\sim\operatorname{diag}(C(p),C(q)) C ( pq ) ∼ diag ( C ( p ) , C ( q )) 当 p , q p,q p , q 互素),相似类完全一样。但「p i p_i p i 恰是非常数不变因子」这句只对第二种写法成立 ,不能同时套在第一种上。
两者都是「不变量 + 标准形」的思路,区别在于:
Jordan 形
有理标准形
需要特征值在 F \mathbb{F} F 内
是
否
块最简单的样子
对角 + 超对角 1 1 1
伴侣矩阵
在 C \mathbb{C} C 上
判别力最强,块形最简单
判别力相同 (都是完全不变量),但块是伴侣块、不好一眼读出块长
在 R \mathbb{R} R 上
可能不存在
总存在
所以:代数闭域上用 Jordan 形,一般域上用有理标准形。 有理标准形是把 Jordan 形的初等因子换成了不可约多项式因子。
回到前面的例子 。A = ( 2 1 0 0 2 1 0 0 2 ) A=\begin{pmatrix}2&1&0\\0&2&1\\0&0&2\end{pmatrix} A = 2 0 0 1 2 0 0 1 2 :不变因子是 1 , 1 , ( λ − 2 ) 3 1,1,(\lambda-2)^3 1 , 1 , ( λ − 2 ) 3 。非常数不变因子只有 ( λ − 2 ) 3 (\lambda-2)^3 ( λ − 2 ) 3 ,所以有理标准形就是一个 伴侣矩阵 C ( ( λ − 2 ) 3 ) C\bigl((\lambda-2)^3\bigr) C ( ( λ − 2 ) 3 ) :
( λ − 2 ) 3 = λ 3 − 6 λ 2 + 12 λ − 8 , (\lambda-2)^3=\lambda^3-6\lambda^2+12\lambda-8, ( λ − 2 ) 3 = λ 3 − 6 λ 2 + 12 λ − 8 ,
即 a 0 = − 8 a_0=-8 a 0 = − 8 、a 1 = 12 a_1=12 a 1 = 12 、a 2 = − 6 a_2=-6 a 2 = − 6 。代入定义(最后一列放 − a i -a_i − a i ):
C = ( 0 0 − a 0 1 0 − a 1 0 1 − a 2 ) = ( 0 0 8 1 0 − 12 0 1 6 ) . 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}. C = 0 1 0 0 0 1 − a 0 − a 1 − a 2 = 0 1 0 0 0 1 8 − 12 6 .
验一下:tr C = 6 \operatorname{tr}C=6 tr C = 6 ,A A A 的迹也是 6 6 6 ✓;det C = 8 \det C=8 det C = 8 ,而三个特征值 2 , 2 , 2 2,2,2 2 , 2 , 2 之积是 8 8 8 ✓。
再看 B = ( 2 1 0 0 2 0 0 0 2 ) B=\begin{pmatrix}2&1&0\\0&2&0\\0&0&2\end{pmatrix} B = 2 0 0 1 2 0 0 0 2 ,它的不变因子是 1 , λ − 2 , ( λ − 2 ) 2 1,\lambda-2,(\lambda-2)^2 1 , λ − 2 , ( λ − 2 ) 2 。非常数不变因子有两个,所以有理标准形有两个块 :
C ( λ − 2 ) = ( 2 ) , C ( ( λ − 2 ) 2 ) = ( 0 − 4 1 4 ) , C(\lambda-2)=(\,2\,),\qquad C\bigl((\lambda-2)^2\bigr)=\begin{pmatrix}0&-4\\1&4\end{pmatrix}, C ( λ − 2 ) = ( 2 ) , C ( ( λ − 2 ) 2 ) = ( 0 1 − 4 4 ) ,
于是 B ∼ diag ( ( 2 ) , ( 0 − 4 1 4 ) ) B\sim\operatorname{diag}\bigl((2),\ \begin{pmatrix}0&-4\\1&4\end{pmatrix}\bigr) B ∼ diag ( ( 2 ) , ( 0 1 − 4 4 ) ) 。
A A A 与 B B B 的不变因子不同(A A A 是 1 , 1 , ( λ − 2 ) 3 1,1,(\lambda-2)^3 1 , 1 , ( λ − 2 ) 3 ,B B B 是 1 , λ − 2 , ( λ − 2 ) 2 1,\lambda-2,(\lambda-2)^2 1 , λ − 2 , ( λ − 2 ) 2 ),所以不相似 ——这点两个标准形都能判定。 但注意:A A A 的有理标准形只有一个 块、B B B 有两个 块,这已经是两者的全部差别。有理标准形把「块怎么切」记录在不变因子的个数与次数 里,而 Jordan 形记录在块的尺寸 里——信息量相同,摆法不同。 在 C \mathbb{C} C 上二者的对应是一对一的:C ( ( λ − λ 0 ) s ) ∼ J s ( λ 0 ) C\bigl((\lambda-\lambda_0)^s\bigr)\sim J_s(\lambda_0) C ( ( λ − λ 0 ) s ) ∼ J s ( λ 0 ) 。所以有理标准形不是「更粗」,只是换了一套允许的块 (伴侣矩阵而不是 Jordan 块),代价是块长得不直观,好处是不要求特征值落在域内。
一个真正需要有理标准形的例子 。R 2 \mathbb{R}^2 R 2 上 A = ( 0 − 1 1 0 ) A=\begin{pmatrix}0&-1\\1&0\end{pmatrix} A = ( 0 1 − 1 0 ) (转 90 ° 90° 90° )。λ I − A = ( λ 1 − 1 λ ) \lambda I-A=\begin{pmatrix}\lambda&1\\-1&\lambda\end{pmatrix} λ I − A = ( λ − 1 1 λ ) ,D 1 = 1 D_1=1 D 1 = 1 (有常数 1 1 1 ),D 2 = λ 2 + 1 D_2=\lambda^2+1 D 2 = λ 2 + 1 ,故不变因子是 ( 1 , λ 2 + 1 ) (1,\ \lambda^2+1) ( 1 , λ 2 + 1 ) 。λ 2 + 1 \lambda^2+1 λ 2 + 1 在 R \mathbb{R} R 上不可约,初等因子就是它自己,有理标准形是单个伴侣块
C ( λ 2 + 1 ) = ( 0 − 1 1 0 ) = A , C(\lambda^2+1)=\begin{pmatrix}0&-1\\1&0\end{pmatrix}=A, C ( λ 2 + 1 ) = ( 0 1 − 1 0 ) = A ,
即 A A A 本身已经是有理标准形。它在 R \mathbb{R} R 上没有 Jordan 形 (第五节误区 1),但有有理标准形——这正是这套工具存在的理由。同一个矩阵搬到 C \mathbb{C} C 上,λ 2 + 1 = ( λ − i ) ( λ + i ) \lambda^2+1=(\lambda-\mathrm{i})(\lambda+\mathrm{i}) λ 2 + 1 = ( λ − i ) ( λ + i ) ,初等因子裂成两个,Jordan 形成了 diag ( i , − i ) \operatorname{diag}(\mathrm{i},-\mathrm{i}) diag ( i , − i ) 。
上例里不变因子 ( 1 , λ 2 + 1 ) (1,\lambda^2+1) ( 1 , λ 2 + 1 ) 在 R \mathbb{R} R 与 C \mathbb{C} C 上一字不变 ,初等因子却从 { λ 2 + 1 } \{\lambda^2+1\} { λ 2 + 1 } 变成 { λ − i , λ + i } \{\lambda-\mathrm{i},\lambda+\mathrm{i}\} { λ − i , λ + i } 。 这就是「不变因子与基域无关,初等因子依赖基域 」的实例,也是「不变」这个形容词的来由。 实用推论:判断两个实矩阵是否在 R \mathbb{R} R 上相似,看不变因子;它们是否在 C \mathbb{C} C 上相似,也看不变因子。 两者答案永远一致。
「有理」在这里不是指有理数,而是指有理标准形只用到域 F \mathbb{F} F 自身的运算 (有理运算),不需要扩域到代数闭包。 这正是它在实矩阵、有理数矩阵上都能用的原因;每个块的元素都是 A A A 的系数的有理表达式。
3.8 相似:把「无从下手」变成「算一串多项式」 #
λ I − A \lambda I-A λ I − A 与 λ I − B \lambda I-B λ I − B 之间能通过初等行、列变换互相化过去,就叫等价 。判据是:
设 A , B A,B A , B 是同阶方阵。下列四条等价:
A ∼ B A\sim B A ∼ B ;
λ I − A \lambda I-A λ I − A 与 λ I − B \lambda I-B λ I − B 等价 ;
λ I − A \lambda I-A λ I − A 与 λ I − B \lambda I-B λ I − B 有相同的 Smith 标准形 ,即不变因子相同;
初等因子组相同。
在代数闭域上还可再添一条:Jordan 形相同(不计块序);在任何域上:有理标准形相同。
「1 ⇒ 2 1\Rightarrow2 1 ⇒ 2 」几乎是白送的:若 A = P B P − 1 A=PBP^{-1} A = P B P − 1 且 P P P 是常数 矩阵,则
λ I − A = λ I − P B P − 1 = P ( λ I − B ) P − 1 , \lambda I-A=\lambda I-PBP^{-1}=P\,(\lambda I-B)\,P^{-1}, λ I − A = λ I − P B P − 1 = P ( λ I − B ) P − 1 ,
右端正是对 λ I − B \lambda I-B λ I − B 左右乘可逆的 λ \lambda λ -矩阵,故两者等价。
「2 ⇒ 1 2\Rightarrow1 2 ⇒ 1 」才是实质内容,靠 Smith 形是不变量、且能从它读回标准形(有理标准形或 Jordan 形)——见 3.7。这两个方向合起来,把「相似」这个关于基变换 的问题,翻译成了关于多项式整除 的问题。
下面三对矩阵,特征多项式全都一样,却两两不相似,差别一层比一层深:
一对矩阵
特征多项式
极小多项式
几何重数
不相似的原因
J 2 ( 2 ) ⊕ J 2 ( 2 ) J_2(2)\oplus J_2(2) J 2 ( 2 ) ⊕ J 2 ( 2 ) 与 J 2 ( 2 ) ⊕ J 1 ( 2 ) ⊕ J 1 ( 2 ) J_2(2)\oplus J_1(2)\oplus J_1(2) J 2 ( 2 ) ⊕ J 1 ( 2 ) ⊕ J 1 ( 2 )
都是 ( λ − 2 ) 4 (\lambda-2)^4 ( λ − 2 ) 4
都是 ( λ − 2 ) 2 (\lambda-2)^2 ( λ − 2 ) 2
2 2 2 vs 3 3 3
块数不同
J 3 ( 2 ) ⊕ J 2 ( 2 ) J_3(2)\oplus J_2(2) J 3 ( 2 ) ⊕ J 2 ( 2 ) 与 J 3 ( 2 ) ⊕ J 1 ( 2 ) ⊕ J 1 ( 2 ) J_3(2)\oplus J_1(2)\oplus J_1(2) J 3 ( 2 ) ⊕ J 1 ( 2 ) ⊕ J 1 ( 2 )
都是 ( λ − 2 ) 5 (\lambda-2)^5 ( λ − 2 ) 5
都是 ( λ − 2 ) 3 (\lambda-2)^3 ( λ − 2 ) 3
2 2 2 vs 3 3 3
块数不同
J 3 ( 2 ) J_3(2) J 3 ( 2 ) 与 J 2 ( 2 ) ⊕ J 1 ( 2 ) J_2(2)\oplus J_1(2) J 2 ( 2 ) ⊕ J 1 ( 2 )
都是 ( λ − 2 ) 3 (\lambda-2)^3 ( λ − 2 ) 3
( λ − 2 ) 3 (\lambda-2)^3 ( λ − 2 ) 3 vs ( λ − 2 ) 2 (\lambda-2)^2 ( λ − 2 ) 2
1 1 1 vs 2 2 2
极小多项式不同
第一、二行是最阴的:特征多项式相同、极小多项式也相同 、最长块也相同,仍不相似。 能一次分开它们的是不变因子:J 2 ( 2 ) ⊕ J 2 ( 2 ) J_2(2)\oplus J_2(2) J 2 ( 2 ) ⊕ J 2 ( 2 ) 的不变因子是 ( 1 , 1 , ( λ − 2 ) 2 , ( λ − 2 ) 2 ) \bigl(1,1,(\lambda-2)^2,(\lambda-2)^2\bigr) ( 1 , 1 , ( λ − 2 ) 2 , ( λ − 2 ) 2 ) ,J 2 ( 2 ) ⊕ J 1 ( 2 ) ⊕ J 1 ( 2 ) J_2(2)\oplus J_1(2)\oplus J_1(2) J 2 ( 2 ) ⊕ J 1 ( 2 ) ⊕ J 1 ( 2 ) 的是 ( 1 , λ − 2 , λ − 2 , ( λ − 2 ) 2 ) \bigl(1,\lambda-2,\lambda-2,(\lambda-2)^2\bigr) ( 1 , λ − 2 , λ − 2 , ( λ − 2 ) 2 ) 。必须比到不变因子(等价地:初等因子)这一层,判断才有保证。
3.9 给定 A A A ,一步步算出 Jordan 形 #
流程四步:
写出 λ I − A \lambda I-A λ I − A ,求 D 1 , … , D n D_1,\dots,D_n D 1 , … , D n ——D k D_k D k 是全部 k k k 阶子式的首一最大公因式。
相除得不变因子 d k = D k / D k − 1 d_k=D_k/D_{k-1} d k = D k / D k − 1 ,等于 1 1 1 的直接丢掉。
把每个非平凡 d k d_k d k 分解成素幂,收集全部素幂,这是初等因子。
每个 ( λ − λ 0 ) s (\lambda-\lambda_0)^s ( λ − λ 0 ) s 写一个 J s ( λ 0 ) J_s(\lambda_0) J s ( λ 0 ) ,按任意顺序拼成分块对角。
第 1 步是唯一费力的。三处能省力:D n D_n D n 就是特征多项式(不必另算);D 1 D_1 D 1 只要矩阵里有一个非零常数元素就等于 1 1 1 ;D k D_k D k 只需算到「公因式稳定」为止。
**完整手算例。**设
A = ( 1 1 0 0 0 1 1 0 0 0 2 0 0 0 0 2 ) , λ I − A = ( λ − 1 − 1 0 0 0 λ − 1 − 1 0 0 0 λ − 2 0 0 0 0 λ − 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}. A = 1 0 0 0 1 1 0 0 0 1 2 0 0 0 0 2 , λ I − A = λ − 1 0 0 0 − 1 λ − 1 0 0 0 − 1 λ − 2 0 0 0 0 λ − 2 .
**第一步:**D 1 = 1 D_1=1 D 1 = 1 ,因为一阶子式里有常数 − 1 -1 − 1 。
**第二步:**D 2 D_2 D 2 。取第 1 , 2 1,2 1 , 2 行与第 2 , 3 2,3 2 , 3 列,子式是 ∣ − 1 0 λ − 1 − 1 ∣ = 1 \begin{vmatrix}-1&0\\\lambda-1&-1\end{vmatrix}=1 − 1 λ − 1 0 − 1 = 1 。故 D 2 = 1 D_2=1 D 2 = 1 。
**第三步:**D 3 D_3 D 3 。三阶子式共 ( 4 3 ) ( 4 3 ) = 16 \binom43\binom43=16 ( 3 4 ) ( 3 4 ) = 16 个。先看重一个事实:第 3 3 3 行只有第 3 3 3 列那个 λ − 2 \lambda-2 λ − 2 非零,第 4 4 4 行只有第 4 4 4 列那个 λ − 2 \lambda-2 λ − 2 非零(λ I − A \lambda I-A λ I − A 是上三角的)。
任取一个三阶子式,它要用三行。若不用第 3 , 4 3,4 3 , 4 行,就只剩第 1 , 2 1,2 1 , 2 两行,凑不满三行——所以任何三阶子式都至少用上第 3 3 3 行或第 4 4 4 行 。沿那一行展开,它就等于 ( λ − 2 ) (\lambda-2) ( λ − 2 ) 乘上某个二阶子式,因此每个三阶子式都被 λ − 2 \lambda-2 λ − 2 整除 。
再算三个具体的,确认公因式正好是 λ − 2 \lambda-2 λ − 2 、没有更高次:
∣ λ − 1 − 1 0 0 λ − 1 − 1 0 0 λ − 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{ 行与列}), λ − 1 0 0 − 1 λ − 1 0 0 − 1 λ − 2 = ( λ − 1 ) 2 ( λ − 2 ) ( 第 1 , 2 , 3 行与列 ) ,
∣ − 1 0 0 λ − 1 − 1 0 0 0 λ − 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{ 列}), − 1 λ − 1 0 0 − 1 0 0 0 λ − 2 = λ − 2 ( 第 1 , 2 , 4 行,第 2 , 3 , 4 列 ) ,
∣ − 1 0 0 0 λ − 2 0 0 0 λ − 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{ 列}). − 1 0 0 0 λ − 2 0 0 0 λ − 2 = − ( λ − 2 ) 2 ( 第 1 , 3 , 4 行,第 2 , 3 , 4 列 ) .
第二个恰等于 λ − 2 \lambda-2 λ − 2 本身,说明公因式的次数不可能更高。故
D 3 = λ − 2. D_3=\lambda-2 . D 3 = λ − 2.
**第四步:**D 4 = det ( λ I − A ) = ( λ − 1 ) 2 ( λ − 2 ) 2 D_4=\det(\lambda I-A)=(\lambda-1)^2(\lambda-2)^2 D 4 = det ( λ I − A ) = ( λ − 1 ) 2 ( λ − 2 ) 2 (上三角,直接读对角线)。
于是
d 1 = 1 , d 2 = 1 , d 3 = λ − 2 , d 4 = ( λ − 1 ) 2 ( λ − 2 ) . d_1=1,\quad d_2=1,\quad d_3=\lambda-2,\quad d_4=(\lambda-1)^2(\lambda-2). d 1 = 1 , d 2 = 1 , d 3 = λ − 2 , d 4 = ( λ − 1 ) 2 ( λ − 2 ) .
**第五步:**分解。d 3 = λ − 2 d_3=\lambda-2 d 3 = λ − 2 与 d 4 = ( λ − 1 ) 2 ( λ − 2 ) d_4=(\lambda-1)^2(\lambda-2) d 4 = ( λ − 1 ) 2 ( λ − 2 ) ,素幂是
( λ − 1 ) 2 , ( λ − 2 ) , ( λ − 2 ) . (\lambda-1)^2,\qquad (\lambda-2),\qquad (\lambda-2). ( λ − 1 ) 2 , ( λ − 2 ) , ( λ − 2 ) .
结论:
Jordan 形 = diag ( J 2 ( 1 ) , J 1 ( 2 ) , J 1 ( 2 ) ) . \text{Jordan 形}=\operatorname{diag}\bigl(J_2(1),\ J_1(2),\ J_1(2)\bigr). Jordan 形 = diag ( J 2 ( 1 ) , J 1 ( 2 ) , J 1 ( 2 ) ) .
核对:块长之和 2 + 1 + 1 = 4 = dim V 2+1+1=4=\dim V 2 + 1 + 1 = 4 = dim V ✓;最长块 2 2 2 ,而极小多项式 d 4 = ( λ − 1 ) 2 ( λ − 2 ) d_4=(\lambda-1)^2(\lambda-2) d 4 = ( λ − 1 ) 2 ( λ − 2 ) 中 ( λ − 1 ) (\lambda-1) ( λ − 1 ) 的指数确实是 2 2 2 ✓;λ = 2 \lambda=2 λ = 2 配两个 1 × 1 1\times1 1 × 1 块,故 dim null ( A − 2 I ) = 2 \dim\operatorname{null}(A-2I)=2 dim null ( A − 2 I ) = 2 ,与直接数「A − 2 I A-2I A − 2 I 的零空间维数」相符 ✓。
值得,因为它是唯一 一条不看基、不求特征向量、纯靠行列式就能走完的路。 手算时若矩阵不大(n ≤ 4 n\le 4 n ≤ 4 ),把子式按「结构」筛一遍比硬列 16 16 16 个可靠;n n n 大时任一手算都会崩,那时应该交给计算机,而判据本身(不变因子相同 ⟺ \iff ⟺ 相似)仍然是那句话 。 还有一个加速技巧:先算特征多项式并因式分解,再对每个特征值 λ 0 \lambda_0 λ 0 只看 dim null ( A − λ 0 I ) k \dim\operatorname{null}(A-\lambda_0I)^k dim null ( A − λ 0 I ) k (3.5 的公式),块长直接读出。两条路可以互相验算——本讲的例题都用两种方法对过。
3.1–3.5 走的是「算子—链—零空间维数」的路线,证明的是存在性 与块长分布的唯一性 ,用的是线性映射的语言。 3.6–3.9 走的是「矩阵—多项式—不变因子」的路线,给的是算法 与判据 ,用的是 F [ λ ] \mathbb{F}[\lambda] F [ λ ] 的语言。 两条路结论完全一致:3.9 的 4 × 4 4\times4 4 × 4 例题算出的 diag ( J 2 ( 1 ) , J 1 ( 2 ) , J 1 ( 2 ) ) \operatorname{diag}(J_2(1),J_1(2),J_1(2)) diag ( J 2 ( 1 ) , J 1 ( 2 ) , J 1 ( 2 )) ,用 3.5 的零空间公式复核也得同样的块长(λ = 2 \lambda=2 λ = 2 的 dim null \dim\operatorname{null} dim null 是 2 2 2 而 dim null 2 \dim\operatorname{null}^2 dim null 2 是 2 2 2 ,故两个一阶块;λ = 1 \lambda=1 λ = 1 的 dim null \dim\operatorname{null} dim null 是 1 1 1 而 dim null 2 \dim\operatorname{null}^2 dim null 2 是 2 2 2 ,故一个二阶块)。 一条解释「为什么」,一条回答「怎么算」。别把它们看成两套矛盾的理论。
四、几何直观 / 记忆钩子 #
4.1 一张图:链 #
每个 Jordan 块背后是一条链。幂零情形(λ = 0 \lambda=0 λ = 0 )画成
u → T T u → T T 2 u → T ⋯ → T T m − 1 u → T 0 u\ \xrightarrow{T}\ Tu\ \xrightarrow{T}\ T^2u\ \xrightarrow{T}\ \cdots\ \xrightarrow{T}\ T^{m-1}u\ \xrightarrow{T}\ 0 u T T u T T 2 u T ⋯ T T m − 1 u T 0
一般情形(λ ≠ 0 \lambda\neq0 λ = 0 ),链的末端不是零,而是一个特征向量:
v → T T v → T ⋯ → T ( T 不再把它推到新方向 ) v\ \xrightarrow{T}\ Tv\ \xrightarrow{T}\ \cdots\ \xrightarrow{T}\ \big(T\text{ 不再把它推到新方向}\big) v T T v T ⋯ T ( T 不再把它推到新方向 )
矩阵里的 1 1 1 就是这些箭头的记录:每个 1 1 1 都是一次「被推着走」;一旦不再被推,块就到头了。
4.2 读矩阵的三句口诀 #
拿到一个 Jordan 形矩阵,不看数字只看形状:
数块 :块与块的分界处上方是 0 0 0 ;块数 = = = 该 λ \lambda λ 的几何重数。
量长 :从左上一路往右下数连续 1 1 1 的个数,加 1 1 1 就是块长;最长的块长 = = = 幂零指数。
定 λ \lambda λ :块内对角元的共同值。
4.3 三个极端,帮你定位 #
情形
Jordan 形长什么样
与前面讲的对应
T T T 可对角化
全是 1 × 1 1\times1 1 × 1 块,即对角阵
20 讲
T T T 幂零
所有块都是 J d ( 0 ) J_d(0) J d ( 0 ) ,对角线上全 0 0 0
31 讲
T = λ I + N T=\lambda I+N T = λ I + N ,N N N 幂零
一个(或几个)J d ( λ ) J_d(\lambda) J d ( λ )
本讲
可对角化就是「每条链都只有一步」 ——这句话把 20 讲和本讲接上了。
五、常见误区 #
定理 B 的前提是 F = C \mathbb{F}=\mathbb{C} F = C 。 反例:R 2 \mathbb{R}^2 R 2 上 T ( x , y ) = ( − y , x ) T(x,y)=(-y,x) T ( x , y ) = ( − y , x ) (逆时针转 90 ° 90° 90° )。特征多项式是 z 2 + 1 z^2+1 z 2 + 1 ,没有实根,连一个特征值都拿不出来。 它不可能有 Jordan 形——Jordan 形的对角线上必须填特征值。 (在 C 2 \mathbb{C}^2 C 2 上看它反而简单:特征值 ± i \pm\mathrm{i} ± i ,两个 1 × 1 1\times1 1 × 1 块。)实空间上能做到 32 讲的分块上三角,做不到 Jordan 形。
32 讲给的是「每块 λ k \lambda_k λ k 在对角线、上面任意」;Jordan 形要求超对角线上只能是 0 0 0 或 1 1 1 。 第一节那个矩阵是上三角,也「差不多」是分块形式,但超对角线上是 3 3 3 而不是 1 1 1 ——它不是 Jordan 形。 用 3.4 的基换过之后才成为 ( 6 1 0 0 6 0 0 0 7 ) \begin{pmatrix}6&1&0\\0&6&0\\0&0&7\end{pmatrix} 6 0 0 1 6 0 0 0 7 。
重数(代数重数)= dim G ( λ , T ) =\dim G(\lambda,T) = dim G ( λ , T ) ,是该 λ \lambda λ 所有块的尺寸之和 。 块数(几何重数)= dim E ( λ , T ) =\dim E(\lambda,T) = dim E ( λ , T ) ,是块的个数 。 例:( 2 1 0 0 2 1 0 0 2 ) \begin{pmatrix}2&1&0\\0&2&1\\0&0&2\end{pmatrix} 2 0 0 1 2 0 0 1 2 中 λ = 2 \lambda=2 λ = 2 的重数是 3 3 3 ,块数却只有 1 1 1 (就一个 J 3 ( 2 ) J_3(2) J 3 ( 2 ) )。 换个写法:( 2 1 0 0 2 0 0 0 2 ) \begin{pmatrix}2&1&0\\0&2&0\\0&0&2\end{pmatrix} 2 0 0 1 2 0 0 0 2 重数还是 3 3 3 ,但块数是 2 2 2 (J 2 ( 2 ) ⊕ J 1 ( 2 ) J_2(2)\oplus J_1(2) J 2 ( 2 ) ⊕ J 1 ( 2 ) )。 两者只在可对角化时相等——那时块数 = = = 重数。
误区 4:以为 λ \lambda λ 相同的块该合并 两个 J 2 ( 6 ) J_2(6) J 2 ( 6 ) 块不能合并成 J 4 ( 6 ) J_4(6) J 4 ( 6 ) 。用 dim null ( T − 6 I ) \dim\operatorname{null}(T-6I) dim null ( T − 6 I ) 就能看出区别:J 2 ( 6 ) ⊕ J 2 ( 6 ) J_2(6)\oplus J_2(6) J 2 ( 6 ) ⊕ J 2 ( 6 ) 的零空间是 2 2 2 维,J 4 ( 6 ) J_4(6) J 4 ( 6 ) 的只有 1 1 1 维。块的个数是算子自己的属性,不是排版习惯。
误区 5:以为 Jordan 块的对角线一定是 0 0 0 0 0 0 只是幂零情形的特例。λ ≠ 0 \lambda\neq0 λ = 0 时,超对角线上的 1 1 1 表示「T T T 沿着链推一格,同时加一个 λ \lambda λ 」, 而不是 「T T T 把向量变成 0 0 0 」。别把 J d ( 0 ) J_d(0) J d ( 0 ) 的直觉套到 J d ( λ ) J_d(\lambda) J d ( λ ) 上。
误区 6:把 Jordan 形当成「像对角形那样的唯一标准形」 它只在「块的顺序可以任意」这个意义下唯一;3.2/3.3 证的是存在性,唯一性见 3.5 末 。 想拿它当判据(例如判断两个矩阵是否相似),必须先把唯一性那块补上。
六、拓展(可跳过) #
6.1 用 Jordan 结构算函数:幂零为什么这么好用 #
设 N N N 幂零、N m = 0 N^m=0 N m = 0 。那么对任意多项式 p p p ,
p ( λ I + N ) = ∑ k ≥ 0 p ( k ) ( λ ) k ! N k p(\lambda I+N)=\sum_{k\ge0}\frac{p^{(k)}(\lambda)}{k!}\,N^k p ( λ I + N ) = k ≥ 0 ∑ k ! p ( k ) ( λ ) N k
这个展开里只有前 m m m 项非零。更一般地,任何形式幂级数在 λ I + N \lambda I+N λ I + N 上都自动截断成有限和 ,收敛性问题消失。这是 Jordan 形最实用的地方。
举一个本讲内容直接支持的例子:I + T I+T I + T 的平方根(T T T 幂零)。考虑
1 + x = 1 + a 1 x + a 2 x 2 + a 3 x 3 + ⋯ , a 1 = 1 2 , a 2 = − 1 8 , a 3 = 1 16 . \sqrt{1+x}=1+a_1x+a_2x^2+a_3x^3+\cdots,\qquad a_1=\tfrac12,\ a_2=-\tfrac18,\ a_3=\tfrac1{16}. 1 + x = 1 + a 1 x + a 2 x 2 + a 3 x 3 + ⋯ , a 1 = 2 1 , a 2 = − 8 1 , a 3 = 16 1 .
把 x x x 换成 T T T 、1 1 1 换成 I I I ,令
R = I + a 1 T + a 2 T 2 + ⋯ + a m − 1 T m − 1 . R=I+a_1T+a_2T^2+\cdots+a_{m-1}T^{m-1}. R = I + a 1 T + a 2 T 2 + ⋯ + a m − 1 T m − 1 .
要 R 2 = I + T R^2=I+T R 2 = I + T ,只需逐次让 T k T^k T k 的系数对上。前几项展开:
R 2 = I + 2 a 1 T + ( 2 a 2 + a 1 2 ) T 2 + ( 2 a 3 + 2 a 1 a 2 ) T 3 + ⋯ R^2=I+2a_1T+\big(2a_2+a_1^2\big)T^2+\big(2a_3+2a_1a_2\big)T^3+\cdots R 2 = I + 2 a 1 T + ( 2 a 2 + a 1 2 ) T 2 + ( 2 a 3 + 2 a 1 a 2 ) T 3 + ⋯
令 T T T 的系数为 1 1 1 :2 a 1 = 1 2a_1=1 2 a 1 = 1 ,得 a 1 = 1 2 a_1=\tfrac12 a 1 = 2 1 。 令 T 2 T^2 T 2 的系数为 0 0 0 :2 a 2 + a 1 2 = 0 2a_2+a_1^2=0 2 a 2 + a 1 2 = 0 ,得 a 2 = − 1 8 a_2=-\tfrac18 a 2 = − 8 1 。 令 T 3 T^3 T 3 的系数为 0 0 0 :2 a 3 + 2 a 1 a 2 = 0 2a_3+2a_1a_2=0 2 a 3 + 2 a 1 a 2 = 0 ,得 a 3 = 1 16 a_3=\tfrac1{16} a 3 = 16 1 。
第 k k k 步的方程总形如 2 a k + ( 只含 a 1 , … , a k − 1 的项 ) = 0 2a_k+(\text{只含 }a_1,\dots,a_{k-1}\text{ 的项})=0 2 a k + ( 只含 a 1 , … , a k − 1 的项 ) = 0 ,因为 2 ≠ 0 2\neq0 2 = 0 ,总能解出 a k a_k a k 。于是 R R R 存在。这个构造对实、复空间都成立 ,只用到 T T T 幂零。
再往上一层:T T T 可逆时,在每个 G ( λ k , T ) G(\lambda_k,T) G ( λ k , T ) 上写成
T ∣ G ( λ k , T ) = λ k ( I + N k λ k ) , T\big|_{G(\lambda_k,T)}=\lambda_k\Big(I+\frac{N_k}{\lambda_k}\Big), T G ( λ k , T ) = λ k ( I + λ k N k ) ,
其中 N k / λ k N_k/\lambda_k N k / λ k 幂零,所以 I + N k / λ k I+N_k/\lambda_k I + N k / λ k 有平方根;把复数 λ k \lambda_k λ k 的平方根与它相乘,就得到该块上的平方根;再用直和拼起来,得到整个 T T T 的平方根。这就是「分块 + 幂零」这台机器的一次真实运转。
C 3 \mathbb{C}^3 C 3 上 T ( z 1 , z 2 , z 3 ) = ( z 2 , z 3 , 0 ) T(z_1,z_2,z_3)=(z_2,z_3,0) T ( z 1 , z 2 , z 3 ) = ( z 2 , z 3 , 0 ) 恰是一个 J 3 ( 0 ) J_3(0) J 3 ( 0 ) ,它没有 平方根。 原因不只是「不可逆」:更精确地说,是链太长。 设 R 2 = T R^2=T R 2 = T ,则 R R R 与 T T T 交换,于是 R R R 必须把 null T \operatorname{null}T null T 送进 null T \operatorname{null}T null T …… 顺着推几步就会撞上矛盾(第七节自测第 5 题)。方根存在与否,取决于链的结构,而不只是可逆性。
6.2 可对角化 = Jordan 形全是一阶块 #
若 T T T 可对角化,Jordan 形就是对角阵(全部 1 × 1 1\times1 1 × 1 块),反之亦然。这条把 20 讲与本节并成一句:可对角化是「每个特征值只有块长为 1 的链」的极端情形。
七、自测(先做,再看答案) #
1. 设 T ∈ L ( C 3 ) T\in\mathcal{L}(\mathbb{C}^3) T ∈ L ( C 3 ) 为 T ( z 1 , z 2 , z 3 ) = ( 0 , z 1 , 0 ) T(z_1,z_2,z_3)=(0,z_1,0) T ( z 1 , z 2 , z 3 ) = ( 0 , z 1 , 0 ) 。求 T T T 的 Jordan 形,并说出块数与各块长。
2. 判断:矩阵 ( 6 3 4 0 6 2 0 0 7 ) \begin{pmatrix}6&3&4\\0&6&2\\0&0&7\end{pmatrix} 6 0 0 3 6 0 4 2 7 是不是 Jordan 形?若不是,该往哪个方向改?
3. 矩阵 B = ( 2 1 0 0 0 2 0 0 0 0 2 1 0 0 0 2 ) B=\begin{pmatrix}2&1&0&0\\0&2&0&0\\0&0&2&1\\0&0&0&2\end{pmatrix} B = 2 0 0 0 1 2 0 0 0 0 2 0 0 0 1 2 中 λ = 2 \lambda=2 λ = 2 的重数与块数各是多少?它能对角化吗?
4. 用本讲的方法,写出 T ( z 1 , z 2 , z 3 ) = ( 6 z 1 + 3 z 2 + 4 z 3 , 6 z 2 + 2 z 3 , 7 z 3 ) T(z_1,z_2,z_3)=(6z_1+3z_2+4z_3,\ 6z_2+2z_3,\ 7z_3) T ( z 1 , z 2 , z 3 ) = ( 6 z 1 + 3 z 2 + 4 z 3 , 6 z 2 + 2 z 3 , 7 z 3 ) 的一组 Jordan 基,并验证矩阵。
5.(难) 证明 C 3 \mathbb{C}^3 C 3 上 T ( z 1 , z 2 , z 3 ) = ( z 2 , z 3 , 0 ) T(z_1,z_2,z_3)=(z_2,z_3,0) T ( z 1 , z 2 , z 3 ) = ( z 2 , z 3 , 0 ) 没有平方根。
答案(点开) 1. T e 1 = e 2 Te_1=e_2 T e 1 = e 2 、T e 2 = 0 Te_2=0 T e 2 = 0 、T e 3 = 0 Te_3=0 T e 3 = 0 。 链:e 1 ↦ e 2 ↦ 0 e_1\mapsto e_2\mapsto0 e 1 ↦ e 2 ↦ 0 (长 2 2 2 ),e 3 ↦ 0 e_3\mapsto0 e 3 ↦ 0 (长 1 1 1 )。 取基 ( e 2 , e 1 , e 3 ) (e_2,e_1,e_3) ( e 2 , e 1 , e 3 ) (每条链倒排、依次拼接):
M = ( 0 1 0 0 0 0 0 0 0 ) = J 2 ( 0 ) ⊕ J 1 ( 0 ) . M=\begin{pmatrix}0&1&0\\0&0&0\\0&0&0\end{pmatrix}=J_2(0)\oplus J_1(0). M = 0 0 0 1 0 0 0 0 0 = J 2 ( 0 ) ⊕ J 1 ( 0 ) .
块数 2 2 2 ,块长 2 2 2 与 1 1 1 。核对:null T = span ( e 2 , e 3 ) \operatorname{null}T=\operatorname{span}(e_2,e_3) null T = span ( e 2 , e 3 ) 是 2 2 2 维,正好等于块数。
2. 不是。它是上三角,超对角线上是 3 3 3 (不是 1 1 1 ),还多出右上角两个数。 方向:在 G ( 6 , T ) = span ( e 1 , e 2 ) G(6,T)=\operatorname{span}(e_1,e_2) G ( 6 , T ) = span ( e 1 , e 2 ) 里把基由 ( e 1 , e 2 ) (e_1,e_2) ( e 1 , e 2 ) 换成 ( 3 e 1 , e 2 ) (3e_1,e_2) ( 3 e 1 , e 2 ) , 第三个基向量换成 ( 10 , 2 , 1 ) (10,2,1) ( 10 , 2 , 1 ) ,得到
diag ( J 2 ( 6 ) , J 1 ( 7 ) ) = ( 6 1 0 0 6 0 0 0 7 ) . \operatorname{diag}\big(J_2(6),J_1(7)\big)=\begin{pmatrix}6&1&0\\0&6&0\\0&0&7\end{pmatrix}. diag ( J 2 ( 6 ) , J 1 ( 7 ) ) = 6 0 0 1 6 0 0 0 7 .
3. 两个 2 × 2 2\times2 2 × 2 块。重数 = dim G ( 2 , B ) = 4 =\dim G(2,B)=4 = dim G ( 2 , B ) = 4 ;块数 = 2 =2 = 2 (把 B − 2 I B-2I B − 2 I 写出来,它只有两个非零元,零空间 2 2 2 维)。 不能对角化:块数 2 < 4 2<4 2 < 4 ,特征向量不够 4 4 4 个。
4. λ = 6 \lambda=6 λ = 6 :N = T − 6 I N=T-6I N = T − 6 I 满足 N ( e 2 ) = 3 e 1 N(e_2)=3e_1 N ( e 2 ) = 3 e 1 、N 2 = 0 N^2=0 N 2 = 0 。 链 e 2 ↦ 3 e 1 ↦ 0 e_2\mapsto3e_1\mapsto0 e 2 ↦ 3 e 1 ↦ 0 ,取基 ( ( 3 , 0 , 0 ) , ( 0 , 1 , 0 ) ) \big((3,0,0),(0,1,0)\big) ( ( 3 , 0 , 0 ) , ( 0 , 1 , 0 ) ) 。λ = 7 \lambda=7 λ = 7 :T ( 10 , 2 , 1 ) = 7 ( 10 , 2 , 1 ) T(10,2,1)=7(10,2,1) T ( 10 , 2 , 1 ) = 7 ( 10 , 2 , 1 ) ,取基 ( 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) ( ( 3 , 0 , 0 ) , ( 0 , 1 , 0 ) , ( 10 , 2 , 1 ) ) ,矩阵 ( 6 1 0 0 6 0 0 0 7 ) \begin{pmatrix}6&1&0\\0&6&0\\0&0&7\end{pmatrix} 6 0 0 1 6 0 0 0 7 。 验算:T ( 3 , 0 , 0 ) = ( 18 , 0 , 0 ) = 6 ( 3 , 0 , 0 ) 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 ( 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) T ( 10 , 2 , 1 ) = ( 70 , 14 , 7 ) = 7 ( 10 , 2 , 1 ) ✓。
5. 记 e 1 = ( 1 , 0 , 0 ) e_1=(1,0,0) e 1 = ( 1 , 0 , 0 ) 等。按题面 T ( z 1 , z 2 , z 3 ) = ( z 2 , z 3 , 0 ) T(z_1,z_2,z_3)=(z_2,z_3,0) T ( z 1 , z 2 , z 3 ) = ( z 2 , z 3 , 0 ) ,所以
T e 1 = 0 , T e 2 = e 1 , T e 3 = e 2 , Te_1=0,\qquad Te_2=e_1,\qquad Te_3=e_2, T e 1 = 0 , T e 2 = e 1 , T e 3 = e 2 ,
即 T T T 的矩阵是 ( 0 1 0 0 0 1 0 0 0 ) = J 3 ( 0 ) \begin{pmatrix}0&1&0\\0&0&1\\0&0&0\end{pmatrix}=J_3(0) 0 0 0 1 0 0 0 1 0 = J 3 ( 0 ) (超对角线上全 1 1 1 ),null T = span ( e 1 ) \operatorname{null}T=\operatorname{span}(e_1) null T = span ( e 1 ) 。 设 R 2 = T R^2=T R 2 = T 。由 R 4 = T 2 R^4=T^2 R 4 = T 2 、R 6 = T 3 = 0 R^6=T^3=0 R 6 = T 3 = 0 。
路线一(顺着交换性) :R T = R ⋅ R 2 = R 3 = R 2 ⋅ R = T R RT=R\cdot R^2=R^3=R^2\cdot R=TR R T = R ⋅ R 2 = R 3 = R 2 ⋅ R = T R ,所以 R R R 与 T T T 交换。于是 R R R 把 null T \operatorname{null}T null T 送进 null T \operatorname{null}T null T :若 T v = 0 Tv=0 T v = 0 ,则 T ( R v ) = R ( T v ) = 0 T(Rv)=R(Tv)=0 T ( R v ) = R ( T v ) = 0 。故 R e 1 = c e 1 Re_1=ce_1 R e 1 = c e 1 。 由 R 2 = T R^2=T R 2 = T 得 R 2 e 1 = T e 1 = 0 R^2e_1=Te_1=\mathbf 0 R 2 e 1 = T e 1 = 0 ,而 R 2 e 1 = R ( c e 1 ) = c 2 e 1 R^2e_1=R(ce_1)=c^2e_1 R 2 e 1 = R ( c e 1 ) = c 2 e 1 ,于是 c 2 = 0 c^2=0 c 2 = 0 ,即 c = 0 c=0 c = 0 、R e 1 = 0 Re_1=0 R e 1 = 0 。 再看 e 2 e_2 e 2 :T e 2 = e 1 ≠ 0 Te_2=e_1\ne0 T e 2 = e 1 = 0 ,故 e 2 ∈ ker T 2 e_2\in\ker T^2 e 2 ∈ ker T 2 (因 T 2 e 2 = T e 1 = 0 T^2e_2=Te_1=0 T 2 e 2 = T e 1 = 0 )。而 R R R 保持 ker T 2 \ker T^2 ker T 2 不变(与 T T T 交换),且 ker T 2 = span ( e 1 , e 2 ) \ker T^2=\operatorname{span}(e_1,e_2) ker T 2 = span ( e 1 , e 2 ) 。设 R e 2 = α e 1 + β e 2 Re_2=\alpha e_1+\beta e_2 R e 2 = α e 1 + β e 2 ,则
R 2 e 2 = R e 1 ⋅ α + R ( β e 2 ) = α R e 1 + β ( α e 1 + β e 2 ) = β α e 1 + β 2 e 2 , 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, R 2 e 2 = R e 1 ⋅ α + R ( β e 2 ) = α R e 1 + β ( α e 1 + β e 2 ) = β α e 1 + β 2 e 2 ,
它必须等于 T e 2 = e 1 Te_2=e_1 T e 2 = e 1 ,故 β 2 = 1 \beta^2=1 β 2 = 1 且 α β = 1 \alpha\beta=1 α β = 1 ,即 α = β = ± 1 \alpha=\beta=\pm1 α = β = ± 1 。代入验 R 2 e 2 = e 1 R^2e_2=e_1 R 2 e 2 = e 1 无矛盾——这条路要再往下才出矛盾,写起来偏长。
路线二(走幂零指数,最干净) :由 R 2 = T R^2=T R 2 = T 得 R 4 = T 2 R^4=T^2 R 4 = T 2 。T 2 T^2 T 2 的矩阵是 T T T 的平方,直接算得 T 2 = ( 0 0 1 0 0 0 0 0 0 ) ≠ 0 T^2=\begin{pmatrix}0&0&1\\0&0&0\\0&0&0\end{pmatrix}\ne0 T 2 = 0 0 0 0 0 0 1 0 0 = 0 ,故 R 4 ≠ 0 R^4\ne0 R 4 = 0 。 另一方面 R 6 = T 3 = 0 R^6=T^3=0 R 6 = T 3 = 0 ,即 R R R 是 R 3 \mathbb{R}^3 R 3 上的幂零算子,幂零指数 ≤ 3 \le3 ≤ 3 (3 维空间上幂零算子的指数不超过维数),所以 R 3 = 0 R^3=0 R 3 = 0 。 于是 T 2 = R 4 = R ⋅ R 3 = 0 T^2=R^4=R\cdot R^3=0 T 2 = R 4 = R ⋅ R 3 = 0 ,与 R 4 ≠ 0 R^4\ne0 R 4 = 0 矛盾。 所以 T T T 没有平方根。(这个论证只用到「T T T 是一条长 3 3 3 的链」,正好印证 6.1 的说法。)
八、复习卡 #
项
内容
Jordan 块 J d ( λ ) J_d(\lambda) J d ( λ )
d × d d\times d d × d :对角全 λ \lambda λ ,超对角线全 1 1 1 ,其余 0 0 0
等价写法
J d ( λ ) = λ I d + N d J_d(\lambda)=\lambda I_d+N_d J d ( λ ) = λ I d + N d ,N d d = 0 N_d^d=0 N d d = 0 ,N d d − 1 ≠ 0 N_d^{d-1}\neq0 N d d − 1 = 0
Jordan 基
使 T T T 的矩阵为各 J d k ( λ k ) J_{d_k}(\lambda_k) J d k ( λ k ) 的分块对角的基
Jordan 形
该矩阵;λ k \lambda_k λ k 不必互异
链(幂零)
u ↦ T u ↦ ⋯ ↦ T m − 1 u ↦ 0 u\mapsto Tu\mapsto\cdots\mapsto T^{m-1}u\mapsto0 u ↦ T u ↦ ⋯ ↦ T m − 1 u ↦ 0
块从哪来
一条链一个块;基倒序排即得超对角线上的 1 1 1
块数
= dim E ( λ , T ) = dim null ( T − λ I ) =\dim E(\lambda,T)=\dim\operatorname{null}(T-\lambda I) = dim E ( λ , T ) = dim null ( T − λ I ) ,即几何重数
块长之和
= dim G ( λ , T ) =\dim G(\lambda,T) = dim G ( λ , T ) ,即代数重数
最长块长
极小多项式中 ( z − λ ) (z-\lambda) ( z − λ ) 的指数;幂零时即幂零指数
存在性(幂零)
定理 A:任何域都成立;归纳 + 泛函造补空间 W W W
存在性(一般)
定理 B:需 F = C \mathbb{F}=\mathbb{C} F = C ;广义特征空间分解 + 定理 A
唯一性
块长分布由 dim null N k \dim\operatorname{null}N^k dim null N k 逐层确定(b k − b k + 1 = b_k-b_{k+1}= b k − b k + 1 = 长度恰为 k k k 的链数);不计块顺序时 Jordan 形唯一
数域限制
R \mathbb{R} R 上一般没有;反例:90 ° 90° 90° 旋转
与对角化
可对角化 ⟺ \iff ⟺ 全为 1 × 1 1\times1 1 × 1 块(每条链只有一步)
规范化
块内超对角线的非零元可靠换基改成 1 1 1 ;0 0 0 表示块的接缝
行列式因子 D k D_k D k
λ I − A \lambda I-A λ I − A 的全部 k k k 阶子式的首一最大公因式;D n = D_n= D n = 特征多项式
不变因子 d k d_k d k
d k = D k / D k − 1 d_k=D_k/D_{k-1} d k = D k / D k − 1 ,满足 d 1 ∣ ⋯ ∣ d n d_1\mid\cdots\mid d_n d 1 ∣ ⋯ ∣ d n ;d n = d_n= d n = 极小多项式
初等因子
各 d k d_k d k 分解成素幂后,收齐全部 ( λ − λ 0 ) s (\lambda-\lambda_0)^s ( λ − λ 0 ) s (按重数重复)
主对应
一个 ( λ − λ 0 ) s (\lambda-\lambda_0)^s ( λ − λ 0 ) s ⟷ \longleftrightarrow ⟷ 一个 J s ( λ 0 ) J_s(\lambda_0) J s ( λ 0 ) ;次数之和 = n =n = n
块数/块长怎么读
λ 0 \lambda_0 λ 0 配几个块 = = = 含 λ − λ 0 \lambda-\lambda_0 λ − λ 0 的初等因子个数;块长 = = = 各指数
相似判据
A ∼ B ⟺ λ I − A A\sim B\iff\lambda I-A A ∼ B ⟺ λ I − A 与 λ I − B \lambda I-B λ I − B 等价 ⟺ \iff ⟺ 不变因子相同 ⟺ \iff ⟺ 初等因子组相同
最常见的误判
特征多项式相同(甚至极小多项式也相同)不足 以保证相似
有理标准形
素幂 p e p^{e} p e 配伴侣块 C [ p e ] C[p^{e}] C [ p e ] ;任何域上都存在,不需特征值
伴侣矩阵
次对角线全 1 1 1 ,最后一列放 − a i -a_i − a i ;det ( λ I − C [ p ] ) = p ( λ ) \det(\lambda I-C[p])=p(\lambda) det ( λ I − C [ p ]) = p ( λ )
两种挂法
初等因子形式(每个素幂一块)/不变因子形式(每个非常数 d k d_k d 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 dim null ( T − λ I ) k (k = 1 , 2 , 3 , … k=1,2,3,\dots k = 1 , 2 , 3 , … )这一族数刻画块长分布,可查任何一本标准教材的「Jordan 标准形」章节,或「有理标准形 / Frobenius 标准形」相关材料。
术语别名 :「Jordan 形」=「若尔当型」=「Jordan 标准形」=「Jordan canonical form」;「广义特征空间」=「根空间(root space)」;「超对角线」=「superdiagonal」;「幂零指数」=「nilpotency index」。
初等因子 / 不变因子 / 有理标准形 :Steven Roman《高等线性代数》(Advanced Linear Algebra , 3rd ed., GTM 135)第 7 章「有理标准形与 Jordan 标准形」。初等因子与不变因子的定义(D k D_k D k 与 d k d_k d 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 的「D k D_k D k 是 k k k 阶子式的最大公因式、Smith 形主对角就是 d k d_k d k 」与 3.9 的四步流程都与这套做法一致。
伴侣矩阵与 Frobenius 标准形 :上述两书均有;「伴侣矩阵的极小多项式与特征多项式都等于 p p p 」是该理论的标准入口,本讲 3.7 用它说明 18 讲那句「每个首一多项式都是某个算子的极小多项式」。
历史 :Camille Jordan(1838–1922)于 1870 年发表该分解的证明;Frobenius 的有理标准形(1887 前后)比 Jordan 形更早给出「任意域上的标准形」。
关于平方根 :a 1 = 1 2 , a 2 = − 1 8 , a 3 = 1 16 a_1=\tfrac12,\ a_2=-\tfrac18,\ a_3=\tfrac1{16} a 1 = 2 1 , a 2 = − 8 1 , a 3 = 16 1 与二项级数 ( 1 + x ) 1 / 2 = ∑ k ( 1 / 2 k ) x k (1+x)^{1/2}=\sum_k\binom{1/2}{k}x^k ( 1 + x ) 1/2 = ∑ k ( k 1/2 ) x k 一致;本讲只用到「能逐步解出」这一点,不依赖级数收敛。