广义特征空间分解 #
这讲讲什么 :当一个算子的特征向量不够凑成一组基时,用「被 ( T − λ I ) (T-\lambda I) ( T − λ I ) 反复作用后归零」的向量把空间补满,得到按特征值划分的直和分解;每一块上算子只是「一个标量 + 一个幂零抖动」。前置 :17 讲(不变子空间与特征值)、20 讲(可对角化)、19 讲(上三角矩阵)、08 讲(零空间与值域)、06 讲(维数与直和)、14 讲(带余除法,只在证明互素分解时用一次)一句话 :G ( λ , T ) = null ( T − λ I ) dim V G(\lambda,T)=\operatorname{null}(T-\lambda I)^{\dim V} G ( λ , T ) = null ( T − λ I ) d i m V 永远是子空间;只要特征值全在数域里,V V V 就是各 G ( λ , T ) G(\lambda,T) G ( λ , T ) 的直和,而 dim G ( λ , T ) \dim G(\lambda,T) dim G ( λ , T ) 就是 λ \lambda λ 的重数。
一、为什么需要它 #
1.1 一个特征空间只有一维的算子 #
看这个最小的例子。取
A = ( 3 1 0 3 ) , A=\begin{pmatrix}3&1\\0&3\end{pmatrix}, A = ( 3 0 1 3 ) ,
T T T 是它在 C 2 \mathbb{C}^2 C 2 的标准基下对应的算子。求特征值:A − 3 I = ( 0 1 0 0 ) A-3I=\begin{pmatrix}0&1\\0&0\end{pmatrix} A − 3 I = ( 0 0 1 0 ) ,它的零空间是 span ( e 1 ) \operatorname{span}(e_1) span ( e 1 ) ,只有一维。
所以 T T T 的特征空间加起来只有一维,C 2 \mathbb{C}^2 C 2 里凑不出由特征向量组成的基。按 20 讲的判据,T T T 不可对角化。
可是 T T T 真的「什么都没管住」吗?把 A − 3 I A-3I A − 3 I 平方一下:
( A − 3 I ) 2 = ( 0 1 0 0 ) 2 = ( 0 0 0 0 ) = 0. (A-3I)^2=\begin{pmatrix}0&1\\0&0\end{pmatrix}^2=\begin{pmatrix}0&0\\0&0\end{pmatrix}=0. ( A − 3 I ) 2 = ( 0 0 1 0 ) 2 = ( 0 0 0 0 ) = 0.
也就是说,C 2 \mathbb{C}^2 C 2 里每一个向量 v v v 都满足 ( T − 3 I ) 2 v = 0 (T-3I)^2v=0 ( T − 3 I ) 2 v = 0 。第一次作用把 e 2 e_2 e 2 压成 e 1 e_1 e 1 那一层,第二次压成 0 0 0 。T T T 并不是抛弃了 e 2 e_2 e 2 ,而是用「3 e 2 3e_2 3 e 2 再叠一个 e 1 e_1 e 1 的偏移」这种带抖动 的方式管住它。
把这种「抖动了几下之后归零」的向量全部收集起来,就得到本讲的主角。它是特征空间的补足:特征空间太小的时候,它往往正好够大。
1.2 上三角化只做到了一半 #
19 讲证明了:有限维复向量空间上的任何算子,都有一组基使矩阵成为上三角。这已经把所有对角元以下的位置清零,但上三角矩阵平均还有大约 n 2 / 2 n^2/2 n 2 /2 个位置是自由的,T 100 T^{100} T 100 依旧不好算。
能不能再多要一点?能。 只要把基按特征值分组——同一特征值的基向量排在一起——不同组之间矩阵元全是零,矩阵成为分块对角 :
( A 1 0 ⋱ 0 A m ) , A k = ( λ k ∗ ⋱ 0 λ k ) . \begin{pmatrix}A_1&&0\\&\ddots&\\0&&A_m\end{pmatrix},\qquad A_k=\begin{pmatrix}\lambda_k&&*\\&\ddots&\\0&&\lambda_k\end{pmatrix}. A 1 0 ⋱ 0 A m , A k = λ k 0 ⋱ ∗ λ k .
跨块的零不是碰巧,而是因为每个 G ( λ k , T ) G(\lambda_k,T) G ( λ k , T ) 在 T T T 下不变(第三节证)。于是「对角化」这个完全失败的方案,被「分块对角化」接住了:块内不是对角,但块内只有一种特征值 ,抖动全被关在块里。求幂、求逆、代多项式,都可以逐块进行。
对角化是奢侈要求,分块对角化是普遍可得的。 这句话是本讲的动机,也是整章的目标。
1.3 读完你应该能回答 #
给定 λ \lambda λ 与 T T T ,G ( λ , T ) G(\lambda,T) G ( λ , T ) 怎么算?为什么指数取 dim V \dim V dim V 就够,取更大不会有新东西?
为什么不同特征值的广义特征空间加起来一定是直和?「直」这件事到底靠什么保证?
「重数」为什么可以用一个零空间的维数来定义,而不必先有行列式?
二、定义与最小例子 #
2.1 定义 #
设 T ∈ L ( V ) T\in L(V) T ∈ L ( V ) ,λ ∈ F \lambda\in\mathbb{F} λ ∈ F 。T T T 对应于 λ \lambda λ 的广义特征空间 定义为
G ( λ , T ) = { v ∈ V : ( T − λ I ) k v = 0 , k 为某个正整数 } . G(\lambda,T)=\{\,v\in V:\ (T-\lambda I)^kv=0,\ k\ \text{为某个正整数}\,\}. G ( λ , T ) = { v ∈ V : ( T − λ I ) k v = 0 , k 为某个正整数 } .
其中的非零向量称为 T T T 对应于 λ \lambda λ 的广义特征向量 。
读这个定义要抓住三处。
第一,量词是 ∀ v ∃ k \forall v\,\exists k ∀ v ∃ k 。 每个 v v v 各自有一个 k k k ,不同的 v v v 可以用不同的 k k k 。定义里没有 说存在一个统一的 k k k 通吃全部向量——但第三节会证明,统一取 k = dim V k=\dim V k = dim V 恰好就够。这个「事后加强」是本讲第一个技术点。
第二,0 0 0 永远在 G ( λ , T ) G(\lambda,T) G ( λ , T ) 里 (取 k = 1 k=1 k = 1 ,( T − λ I ) 0 = 0 (T-\lambda I)0=0 ( T − λ I ) 0 = 0 )。所以 G ( λ , T ) G(\lambda,T) G ( λ , T ) 从不空,问题只在于它是不是只有 0 0 0 。
第三,E ( λ , T ) ⊆ G ( λ , T ) E(\lambda,T)\subseteq G(\lambda,T) E ( λ , T ) ⊆ G ( λ , T ) 。 特征向量满足 ( T − λ I ) v = 0 (T-\lambda I)v=0 ( T − λ I ) v = 0 ,取 k = 1 k=1 k = 1 即符合定义。反过来的包含一般不成立,1.1 的例子就是:E ( 3 , T ) = span ( e 1 ) E(3,T)=\operatorname{span}(e_1) E ( 3 , T ) = span ( e 1 ) ,而 G ( 3 , T ) = C 2 G(3,T)=\mathbb{C}^2 G ( 3 , T ) = C 2 。
2.2 最小例子:把 1.1 算完整 #
仍取 A = ( 3 1 0 3 ) A=\begin{pmatrix}3&1\\0&3\end{pmatrix} A = ( 3 0 1 3 ) ,N = T − 3 I N=T-3I N = T − 3 I 。上面已经算出 N 2 = 0 N^2=0 N 2 = 0 ,所以对任何 v v v 都有 ( T − 3 I ) 2 v = 0 (T-3I)^2v=0 ( T − 3 I ) 2 v = 0 ,于是
G ( 3 , T ) = C 2 , dim G ( 3 , T ) = 2. G(3,T)=\mathbb{C}^2,\qquad \dim G(3,T)=2. G ( 3 , T ) = C 2 , dim G ( 3 , T ) = 2.
而
E ( 3 , T ) = null ( T − 3 I ) = span ( e 1 ) , dim E ( 3 , T ) = 1. E(3,T)=\operatorname{null}(T-3I)=\operatorname{span}(e_1),\qquad \dim E(3,T)=1. E ( 3 , T ) = null ( T − 3 I ) = span ( e 1 ) , dim E ( 3 , T ) = 1.
同一个特征值 3 3 3 ,两个维数不一样大。这个差距就是「不可对角化」的量化形式:E E E 少了 1 1 1 维,G G G 补上了它。
同时注意 G ( 3 , T ) G(3,T) G ( 3 , T ) 的基可以取 e 1 , e 2 e_1,e_2 e 1 , e 2 ,而 T T T 关于这组基的矩阵还是 A A A 本身——上三角、对角元全是 3 3 3 。这正是 1.2 里那种「块」的最简形态。
2.3 一个真正拆开的例子 #
只看一个特征值看不出分解。取 C 3 \mathbb{C}^3 C 3 上的算子
T ( z 1 , z 2 , z 3 ) = ( 4 z 2 , 0 , 5 z 3 ) . T(z_1,z_2,z_3)=(4z_2,\ 0,\ 5z_3). T ( z 1 , z 2 , z 3 ) = ( 4 z 2 , 0 , 5 z 3 ) .
它关于标准基的矩阵是
A = ( 0 4 0 0 0 0 0 0 5 ) , A=\begin{pmatrix}0&4&0\\0&0&0\\0&0&5\end{pmatrix}, A = 0 0 0 4 0 0 0 0 5 ,
上三角,对角元是 0 , 0 , 5 0,0,5 0 , 0 , 5 ,所以特征值只有 0 0 0 与 5 5 5 (上三角矩阵的特征值就是对角元,见 19 讲)。
先算 T 2 T^2 T 2 :T ( e 1 ) = 0 T(e_1)=0 T ( e 1 ) = 0 ,T ( e 2 ) = 4 e 1 T(e_2)=4e_1 T ( e 2 ) = 4 e 1 而 T ( e 1 ) = 0 T(e_1)=0 T ( e 1 ) = 0 ,T ( e 3 ) = 5 e 3 T(e_3)=5e_3 T ( e 3 ) = 5 e 3 而 T ( 5 e 3 ) = 25 e 3 T(5e_3)=25e_3 T ( 5 e 3 ) = 25 e 3 。故 T 2 ( e 1 ) = T 2 ( e 2 ) = 0 T^2(e_1)=T^2(e_2)=0 T 2 ( e 1 ) = T 2 ( e 2 ) = 0 ,T 2 ( e 3 ) = 25 e 3 T^2(e_3)=25e_3 T 2 ( e 3 ) = 25 e 3 。于是
G ( 0 , T ) = null T 3 = null T 2 = span ( e 1 , e 2 ) , dim = 2. G(0,T)=\operatorname{null}T^3=\operatorname{null}T^2=\operatorname{span}(e_1,e_2),\qquad \dim=2. G ( 0 , T ) = null T 3 = null T 2 = span ( e 1 , e 2 ) , dim = 2.
再算 T − 5 I T-5I T − 5 I :A − 5 I = ( − 5 4 0 0 − 5 0 0 0 0 ) A-5I=\begin{pmatrix}-5&4&0\\0&-5&0\\0&0&0\end{pmatrix} A − 5 I = − 5 0 0 4 − 5 0 0 0 0 ,零空间由 − 5 z 1 + 4 z 2 = 0 -5z_1+4z_2=0 − 5 z 1 + 4 z 2 = 0 、− 5 z 2 = 0 -5z_2=0 − 5 z 2 = 0 给出,即 z 2 = z 1 = 0 z_2=z_1=0 z 2 = z 1 = 0 、z 3 z_3 z 3 自由,所以
G ( 5 , T ) = span ( e 3 ) , dim = 1. G(5,T)=\operatorname{span}(e_3),\qquad \dim=1. G ( 5 , T ) = span ( e 3 ) , dim = 1.
两者的交只有 0 0 0 (span ( e 1 , e 2 ) ∩ span ( e 3 ) = { 0 } \operatorname{span}(e_1,e_2)\cap\operatorname{span}(e_3)=\{0\} span ( e 1 , e 2 ) ∩ span ( e 3 ) = { 0 } ),维数相加 2 + 1 = 3 = dim C 3 2+1=3=\dim\mathbb{C}^3 2 + 1 = 3 = dim C 3 ,于是
C 3 = G ( 0 , T ) ⊕ G ( 5 , T ) . \mathbb{C}^3=G(0,T)\oplus G(5,T). C 3 = G ( 0 , T ) ⊕ G ( 5 , T ) .
这就是本讲主定理在具体例子里的样子:空间被完整地切成两块,每块对应一个特征值。
2.4 回看定义:三个立刻能问的问题 #
问题一:λ \lambda λ 不是特征值时,G ( λ , T ) G(\lambda,T) G ( λ , T ) 是什么?
答案是 { 0 } \{0\} { 0 } 。若存在 v ≠ 0 v\neq0 v = 0 与正整数 k k k 使 ( T − λ I ) k v = 0 (T-\lambda I)^kv=0 ( T − λ I ) k v = 0 ,取最小 的这样的 k k k 。由最小性,w = ( T − λ I ) k − 1 v ≠ 0 w=(T-\lambda I)^{k-1}v\neq0 w = ( T − λ I ) k − 1 v = 0 ,而 ( T − λ I ) w = 0 (T-\lambda I)w=0 ( T − λ I ) w = 0 ,所以 w w w 是特征向量,λ \lambda λ 是特征值。逆命题也不用证:特征值当然有非零的 G G G 。
于是「G ( λ , T ) ≠ { 0 } G(\lambda,T)\neq\{0\} G ( λ , T ) = { 0 } 」与「λ \lambda λ 是特征值」是同一件事。这条小结论在自测里会用到。
问题二:G ( λ , T ) G(\lambda,T) G ( λ , T ) 为什么一定是子空间?
定义里它是若干个零空间的并 ⋃ k null ( T − λ I ) k \bigcup_k\operatorname{null}(T-\lambda I)^k ⋃ k null ( T − λ I ) k ,而零空间的并一般不是 子空间。它能成为子空间,全靠第三节那条链稳定:这个并实际上等于其中某一项,而零空间是子空间。
先有「并是否等于某一项」这个技术事实,才有「G G G 是子空间」。直接说「G G G 是子空间,因为它是零空间」是循环论证——它本来不是单个零空间。
问题三:k k k 到底要试到多大?
试到不再增长为止。第三节证明这一步最多走 dim V \dim V dim V 次。所以在实际计算中,永远只需要算到 ( T − λ I ) dim V (T-\lambda I)^{\dim V} ( T − λ I ) d i m V 。
三、几个必须自己算一遍的结论 #
3.1 链稳定:为什么指数取 dim V \dim V dim V 就够 #
记 N = T − λ I N=T-\lambda I N = T − λ I ,U k = null N k U_k=\operatorname{null}N^k U k = null N k 。首先有包含关系:
U 1 ⊆ U 2 ⊆ U 3 ⊆ ⋯ U_1\subseteq U_2\subseteq U_3\subseteq\cdots U 1 ⊆ U 2 ⊆ U 3 ⊆ ⋯
因为 N k v = 0 N^kv=0 N k v = 0 推出 N k + 1 v = N ( N k v ) = 0 N^{k+1}v=N(N^kv)=0 N k + 1 v = N ( N k v ) = 0 。
若某个 k k k 使 U k = U k + 1 U_k=U_{k+1} U k = U k + 1 ,则对一切 j ≥ k j\ge k j ≥ k 都有 U j = U k U_j=U_k U j = U k 。
证明只有两行。设 U k = U k + 1 U_k=U_{k+1} U k = U k + 1 ,要证 U k + 1 = U k + 2 U_{k+1}=U_{k+2} U k + 1 = U k + 2 :取 v ∈ U k + 2 v\in U_{k+2} v ∈ U k + 2 ,即 N k + 2 v = 0 N^{k+2}v=0 N k + 2 v = 0 。把 N k + 2 v = N k + 1 ( N v ) N^{k+2}v=N^{k+1}(Nv) N k + 2 v = N k + 1 ( N v ) 读一遍:N v ∈ U k + 1 = U k Nv\in U_{k+1}=U_k N v ∈ U k + 1 = U k ,也就是 N k ( N v ) = 0 N^k(Nv)=0 N k ( N v ) = 0 ,即 N k + 1 v = 0 N^{k+1}v=0 N k + 1 v = 0 ,故 v ∈ U k + 1 v\in U_{k+1} v ∈ U k + 1 。反包含已在上面。于是逐步右移,全链从第 k k k 项起恒定。
于是一条只升不降的维数序列 dim U 1 ≤ dim U 2 ≤ ⋯ \dim U_1\le\dim U_2\le\cdots dim U 1 ≤ dim U 2 ≤ ⋯ ,只要还没稳定,每步至少涨 1 1 1 ;而它不能超过 dim V \dim V dim V 。所以最多走 dim V \dim V dim V 步必停 ,并且
G ( λ , T ) = ⋃ k ≥ 1 U k = U dim V = null ( T − λ I ) dim V . G(\lambda,T)=\bigcup_{k\ge1}U_k=U_{\dim V}=\operatorname{null}(T-\lambda I)^{\dim V}. G ( λ , T ) = k ≥ 1 ⋃ U k = U d i m V = null ( T − λ I ) d i m V .
这一条同时解决三件事:G G G 是子空间(它是 U dim V U_{\dim V} U d i m V )、G G G 可以有限步算出、以及定义里那个「某个 k k k 」可以统一换成 dim V \dim V dim V 。
因为那时你还不知道 T − λ I T-\lambda I T − λ I 在 G G G 上的幂零指数是多少——那正是要算的东西。dim V \dim V dim V 是一个不需要任何信息 就能写下的上界,这正是它的价值。
3.2 手算一个 3 × 3 3\times3 3 × 3 :两种重数分道扬镳 #
取
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 6 与 7 7 7 。
先算 λ = 6 \lambda=6 λ = 6 。 令 N = A − 6 I = ( 0 3 4 0 0 2 0 0 1 ) N=A-6I=\begin{pmatrix}0&3&4\\0&0&2\\0&0&1\end{pmatrix} N = A − 6 I = 0 0 0 3 0 0 4 2 1 。它对 ( z 1 , z 2 , z 3 ) (z_1,z_2,z_3) ( z 1 , z 2 , z 3 ) 的作用是
N ( z 1 , z 2 , z 3 ) = ( 3 z 2 + 4 z 3 , 2 z 3 , z 3 ) . N(z_1,z_2,z_3)=(3z_2+4z_3,\ 2z_3,\ z_3). N ( z 1 , z 2 , z 3 ) = ( 3 z 2 + 4 z 3 , 2 z 3 , z 3 ) .
令其为零:z 3 = 0 z_3=0 z 3 = 0 ,再由第一式 z 2 = 0 z_2=0 z 2 = 0 ,z 1 z_1 z 1 自由。故 null N = span ( e 1 ) \operatorname{null}N=\operatorname{span}(e_1) null N = span ( e 1 ) ,维数 1 1 1 。
再算 N 2 N^2 N 2 :
N 2 = ( 0 0 10 0 0 2 0 0 1 ) , N 2 ( z 1 , z 2 , z 3 ) = ( 10 z 3 , 2 z 3 , z 3 ) . N^2=\begin{pmatrix}0&0&10\\0&0&2\\0&0&1\end{pmatrix},\qquad N^2(z_1,z_2,z_3)=(10z_3,\ 2z_3,\ z_3). N 2 = 0 0 0 0 0 0 10 2 1 , N 2 ( z 1 , z 2 , z 3 ) = ( 10 z 3 , 2 z 3 , z 3 ) .
零空间是 z 3 = 0 z_3=0 z 3 = 0 ,即 span ( e 1 , e 2 ) \operatorname{span}(e_1,e_2) span ( e 1 , e 2 ) ,维数 2 2 2 。再乘一次:N 3 = N 2 N^3=N^2 N 3 = N 2 (N 2 N^2 N 2 的各行都是第三行的倍数,再乘 N N N 不变),故零空间停在维数 2 2 2 。于是
G ( 6 , T ) = span ( e 1 , e 2 ) , dim G ( 6 , T ) = 2. G(6,T)=\operatorname{span}(e_1,e_2),\qquad \dim G(6,T)=2. G ( 6 , T ) = span ( e 1 , e 2 ) , dim G ( 6 , T ) = 2.
再算 λ = 7 \lambda=7 λ = 7 。 M = A − 7 I = ( − 1 3 4 0 − 1 2 0 0 0 ) M=A-7I=\begin{pmatrix}-1&3&4\\0&-1&2\\0&0&0\end{pmatrix} M = A − 7 I = − 1 0 0 3 − 1 0 4 2 0 。令 M ( z ) = 0 M(z)=0 M ( z ) = 0 :第三行自动成立,第二行给 z 2 = 2 z 3 z_2=2z_3 z 2 = 2 z 3 ,第一行给 z 1 = 3 z 2 + 4 z 3 = 10 z 3 z_1=3z_2+4z_3=10z_3 z 1 = 3 z 2 + 4 z 3 = 10 z 3 。于是
G ( 7 , T ) = null M = span ( ( 10 , 2 , 1 ) ) , dim G ( 7 , T ) = 1. G(7,T)=\operatorname{null}M=\operatorname{span}\big((10,2,1)\big),\qquad \dim G(7,T)=1. G ( 7 , T ) = null M = span ( ( 10 , 2 , 1 ) ) , dim G ( 7 , T ) = 1.
这一块一步就到底,没有再往上加维数。
核对。 两维数相加 2 + 1 = 3 = dim C 3 2+1=3=\dim\mathbb{C}^3 2 + 1 = 3 = dim C 3 ,两组基向量 ( 1 , 0 , 0 ) , ( 0 , 1 , 0 ) , ( 10 , 2 , 1 ) (1,0,0),(0,1,0),(10,2,1) ( 1 , 0 , 0 ) , ( 0 , 1 , 0 ) , ( 10 , 2 , 1 ) 线性无关,构成 C 3 \mathbb{C}^3 C 3 的基。T T T 关于这组基的矩阵是
( 6 3 0 0 6 0 0 0 7 ) , \begin{pmatrix}6&3&0\\0&6&0\\0&0&7\end{pmatrix}, 6 0 0 3 6 0 0 0 7 ,
分块对角,两块分别是 2 × 2 2\times2 2 × 2 与 1 × 1 1\times1 1 × 1 。验证一下交叉项:T ( 10 , 2 , 1 ) = ( 60 + 6 + 4 , 12 + 2 , 7 ) = ( 70 , 14 , 7 ) = 7 ( 10 , 2 , 1 ) T(10,2,1)=(60+6+4,\ 12+2,\ 7)=(70,14,7)=7(10,2,1) T ( 10 , 2 , 1 ) = ( 60 + 6 + 4 , 12 + 2 , 7 ) = ( 70 , 14 , 7 ) = 7 ( 10 , 2 , 1 ) ,确实被 7 7 7 特住;T ( e 2 ) = 3 e 1 + 6 e 2 T(e_2)=3e_1+6e_2 T ( e 2 ) = 3 e 1 + 6 e 2 ,块内的那个 3 3 3 就是抖动的痕迹。
注意这里最容易被忽略的一点 :dim E ( 6 , T ) = 1 \dim E(6,T)=1 dim E ( 6 , T ) = 1 而 dim G ( 6 , T ) = 2 \dim G(6,T)=2 dim G ( 6 , T ) = 2 。特征向量只有 e 1 e_1 e 1 一个方向,但广义特征向量有两个方向。所以这组基里有 2 2 2 个来自 G ( 6 ) G(6) G ( 6 ) ,其中只有 1 1 1 个是特征向量——C 3 \mathbb{C}^3 C 3 没有由特征向量组成的基 ,这正是 20 讲判据说它不可对角化的原因。
3.3 主分解定理 #
设 F = C \mathbb{F}=\mathbb{C} F = C ,T ∈ L ( V ) T\in L(V) T ∈ L ( V ) ,λ 1 , … , λ m \lambda_1,\dots,\lambda_m λ 1 , … , λ m 是 T T T 的全部互异特征值。则
(a) 每个 G ( λ k , T ) G(\lambda_k,T) G ( λ k , T ) 在 T T T 下不变;(b) 每个 ( T − λ k I ) ∣ G ( λ k , T ) (T-\lambda_kI)|_{G(\lambda_k,T)} ( T − λ k I ) ∣ G ( λ k , T ) 是幂零的;(c) 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 ) 。
(a) 的证明。 由 3.1,G ( λ k , T ) = null p ( T ) G(\lambda_k,T)=\operatorname{null}p(T) G ( λ k , T ) = null p ( T ) ,其中 p ( z ) = ( z − λ k ) dim V p(z)=(z-\lambda_k)^{\dim V} p ( z ) = ( z − λ k ) d i m V 。多项式算子与 T T T 可交换:p ( T ) T = T p ( T ) p(T)T=Tp(T) p ( T ) T = T p ( T ) 。于是若 p ( T ) v = 0 p(T)v=0 p ( T ) v = 0 ,则
p ( T ) ( T v ) = T ( p ( T ) v ) = T 0 = 0 , p(T)(Tv)=T\big(p(T)v\big)=T0=0, p ( T ) ( T v ) = T ( p ( T ) v ) = T 0 = 0 ,
即 T v Tv T v 也在零空间里。这就是 17 讲那条「多项式零空间是不变子空间」的直接应用。
(b) 的证明。 定义就摆在那里:对 v ∈ G ( λ k , T ) v\in G(\lambda_k,T) v ∈ G ( λ k , T ) 有 ( T − λ k I ) dim V v = 0 (T-\lambda_kI)^{\dim V}v=0 ( T − λ k I ) d i m V v = 0 ,所以限制映射的 dim V \dim V dim V 次幂是零算子。幂零指数还更小,见第六节。
(c) 的证明分两半:直和、张满。 直和这一半靠一条引理。
若 λ 1 , … , λ m \lambda_1,\dots,\lambda_m λ 1 , … , λ m 互异,v k ∈ G ( λ k , T ) v_k\in G(\lambda_k,T) v k ∈ G ( λ k , T ) ,且 v 1 + ⋯ + v m = 0 v_1+\cdots+v_m=0 v 1 + ⋯ + v m = 0 ,则每个 v k = 0 v_k=0 v k = 0 。
对项数归纳。m = 1 m=1 m = 1 时显然。设 m ≥ 2 m\ge2 m ≥ 2 ,取 p = dim V p=\dim V p = dim V ,令 N = ( T − λ m I ) p N=(T-\lambda_mI)^p N = ( T − λ m I ) p ,则 N v m = 0 Nv_m=0 N v m = 0 。对 k < m k<m k < m ,N N N 限制在 G ( λ k , T ) G(\lambda_k,T) G ( λ k , T ) 上可逆 :在该子空间上
T − λ m I = ( λ k − λ m ) I + ( T − λ k I ) , T-\lambda_mI=(\lambda_k-\lambda_m)I+\big(T-\lambda_kI\big), T − λ m I = ( λ k − λ m ) I + ( T − λ k I ) ,
后半是幂零的,而 λ k − λ m ≠ 0 \lambda_k-\lambda_m\neq0 λ k − λ m = 0 ,故它是「非零标量 + 幂零」,可逆;可逆算子的幂也可逆。
对 v 1 + ⋯ + v m = 0 v_1+\cdots+v_m=0 v 1 + ⋯ + v m = 0 两边作用 N N N :
N v 1 + ⋯ + N v m − 1 = 0. Nv_1+\cdots+Nv_{m-1}=0. N v 1 + ⋯ + N v m − 1 = 0.
各项仍分别在 G ( λ 1 , T ) , … , G ( λ m − 1 , T ) G(\lambda_1,T),\dots,G(\lambda_{m-1},T) G ( λ 1 , T ) , … , G ( λ m − 1 , T ) 中(因为 N N N 是 T T T 的多项式,这些子空间在 N N N 下不变),由归纳假设(只有 m − 1 m-1 m − 1 项)得 N v k = 0 Nv_k=0 N v k = 0 ,再由 N N N 可逆得 v k = 0 v_k=0 v k = 0 (k < m k<m k < m )。最后 v m = 0 v_m=0 v m = 0 。
有了这条引理,直和立刻成立:若 ∑ k v k = 0 \sum_k v_k=0 ∑ k v k = 0 且 v k ∈ G ( λ k , T ) v_k\in G(\lambda_k,T) v k ∈ G ( λ k , T ) ,则所有 v k = 0 v_k=0 v k = 0 ——这正是直和的定义(06 讲)。
「几个子空间的维数加起来等于 dim V \dim V dim V 」并不能保证是直和:R 3 \mathbb{R}^3 R 3 里两个过原点的平面,维数 2 + 2 = 4 > 3 2+2=4>3 2 + 2 = 4 > 3 ,它们必然交于一条线。直和是要证 的,靠的是上面的线性无关性,不是维数。
3.4 张满:为什么这些块能拼回整个 V V V #
直和只说明「不重叠」,还得说明「不缺」。这一步是本讲技术含量最高的地方,但路线很短:上三角化 + 一条恒等式 + 一次互素分解。
第一步:上三角化带来一条零化恒等式。 由 19 讲,取基 v 1 , … , v n v_1,\dots,v_n v 1 , … , v n 使 T T T 的矩阵 A A A 上三角,对角元记作 λ 1 , … , λ n \lambda_1,\dots,\lambda_n λ 1 , … , λ n (按位置记,可以重复)。则
∏ j = 1 n ( T − λ j I ) = 0. \prod_{j=1}^{n}(T-\lambda_jI)=0. j = 1 ∏ n ( T − λ j I ) = 0.
理由是一个三行归纳。记 M k = ∏ j = 1 k ( T − λ j I ) M_k=\prod_{j=1}^{k}(T-\lambda_jI) M k = ∏ j = 1 k ( T − λ j I ) ,断言 M k M_k M k 杀死 v 1 , … , v k v_1,\dots,v_k v 1 , … , v k 。k = 1 k=1 k = 1 时,上三角矩阵 A − λ 1 I A-\lambda_1I A − λ 1 I 的第一列全为零,故 M 1 v 1 = 0 M_1v_1=0 M 1 v 1 = 0 。归纳步:M k + 1 v i = M k ( ( A − λ k + 1 I ) v i ) M_{k+1}v_i=M_k\big((A-\lambda_{k+1}I)v_i\big) M k + 1 v i = M k ( ( A − λ k + 1 I ) v i ) ;当 i ≤ k i\le k i ≤ k 时,( A − λ k + 1 I ) v i ∈ span ( v 1 , … , v i ) (A-\lambda_{k+1}I)v_i\in\operatorname{span}(v_1,\dots,v_i) ( A − λ k + 1 I ) v i ∈ span ( v 1 , … , v i ) ,被 M k M_k M k 杀死;当 i = k + 1 i=k+1 i = k + 1 时,该向量属于 span ( v 1 , … , v k + 1 ) \operatorname{span}(v_1,\dots,v_{k+1}) span ( v 1 , … , v k + 1 ) ,而其中 v k + 1 v_{k+1} v k + 1 的系数是 λ k + 1 − λ k + 1 = 0 \lambda_{k+1}-\lambda_{k+1}=0 λ k + 1 − λ k + 1 = 0 ,所以它其实落在 span ( v 1 , … , v k ) \operatorname{span}(v_1,\dots,v_k) span ( v 1 , … , v k ) 里,同样被杀死。取 k = n k=n k = n ,M n M_n M n 杀死全部基向量,故 M n = 0 M_n=0 M n = 0 。
第二步:把重复的因式合并。 设互异特征值为 μ 1 , … , μ r \mu_1,\dots,\mu_r μ 1 , … , μ r ,μ k \mu_k μ k 在对角线上出现 n k n_k n k 次(n 1 + ⋯ + n r = n n_1+\cdots+n_r=n n 1 + ⋯ + n r = n )。则
p ( z ) = ∏ k = 1 r ( z − μ k ) n k , p ( T ) = 0. p(z)=\prod_{k=1}^{r}(z-\mu_k)^{n_k},\qquad p(T)=0. p ( z ) = k = 1 ∏ r ( z − μ k ) n k , p ( T ) = 0.
第三步:互素分解,把 v v v 拆开。 令 p k = p / ( z − μ k ) n k = ∏ j ≠ k ( z − μ j ) n j p_k=p/(z-\mu_k)^{n_k}=\prod_{j\neq k}(z-\mu_j)^{n_j} p k = p / ( z − μ k ) n k = ∏ j = k ( z − μ j ) n j 。诸 p k p_k p k 没有公共零点,故两两互素,于是存在多项式 a 1 , … , a r a_1,\dots,a_r a 1 , … , a r 使
a 1 p 1 + ⋯ + a r p r = 1 a_1p_1+\cdots+a_rp_r=1 a 1 p 1 + ⋯ + a r p r = 1
(14 讲的带余除法给出这个裴蜀型恒等式)。把 T T T 代进去,对任意 v ∈ V v\in V v ∈ V :
v = ∑ k = 1 r a k ( T ) p k ( T ) v . v=\sum_{k=1}^{r}a_k(T)p_k(T)v. v = k = 1 ∑ r a k ( T ) p k ( T ) v .
而第 k k k 项确实落在 G ( μ k , T ) G(\mu_k,T) G ( μ k , T ) 里,因为
( T − μ k I ) n k ( a k ( T ) p k ( T ) v ) = a k ( T ) ( T − μ k I ) n k p k ( T ) ⏟ = p ( T ) v = a k ( T ) 0 = 0. (T-\mu_kI)^{n_k}\Big(a_k(T)p_k(T)v\Big)=a_k(T)\underbrace{(T-\mu_kI)^{n_k}p_k(T)}_{=\,p(T)}v=a_k(T)\,0=0. ( T − μ k I ) n k ( a k ( T ) p k ( T ) v ) = a k ( T ) = p ( T ) ( T − μ k I ) n k p k ( T ) v = a k ( T ) 0 = 0.
所以每个 v v v 都是若干个广义特征向量之和:
V = G ( μ 1 , T ) + ⋯ + G ( μ r , T ) . V=G(\mu_1,T)+\cdots+G(\mu_r,T). V = G ( μ 1 , T ) + ⋯ + G ( μ r , T ) .
配上 3.3 的直和,(c) 证完。
第一步用上三角化(19 讲),第二步是用词,第三步用多项式互素(14 讲)。没有一步用到行列式 ,也没有一步用到「基由广义特征向量组成」这个待证的结论——这正是这条路线的好处:它不循环。
3.5 重数:dim G ( λ , T ) \dim G(\lambda,T) dim G ( λ , T ) #
设 T ∈ L ( V ) T\in L(V) T ∈ L ( V ) ,λ \lambda λ 是 T T T 的特征值。λ \lambda λ 的重数 定义为其广义特征空间的维数:
mult ( λ ) = dim G ( λ , T ) = dim null ( T − λ I ) dim V . \operatorname{mult}(\lambda)=\dim G(\lambda,T)=\dim\operatorname{null}(T-\lambda I)^{\dim V}. mult ( λ ) = dim G ( λ , T ) = dim null ( T − λ I ) d i m V .
由 3.3(c) 和直和的维数相加(06 讲),立刻得到本讲的第二个主结论:
∑ k = 1 m dim G ( λ k , T ) = dim V . \sum_{k=1}^{m}\dim G(\lambda_k,T)=\dim V. k = 1 ∑ m dim G ( λ k , T ) = dim V .
这个等式看起来平常,实际很强:它说明 G G G 的维数不是随意的,几个特征值之间必须恰好 把 dim V \dim V dim V 分完。
又因为 E ( λ , T ) ⊆ G ( λ , T ) E(\lambda,T)\subseteq G(\lambda,T) E ( λ , T ) ⊆ G ( λ , T ) ,总有
1 ≤ dim E ( λ , T ) ≤ mult ( λ ) ( λ 是特征值 ) . 1\le\dim E(\lambda,T)\le\operatorname{mult}(\lambda)\qquad(\lambda\ \text{是特征值}). 1 ≤ dim E ( λ , T ) ≤ mult ( λ ) ( λ 是特征值 ) .
左边取 1 1 1 是因为特征值一定有非零特征向量。这个双重不等式把「可对角化」量化了:
T T T 可对角化 ⟺ \iff ⟺ 对每个 特征值 λ \lambda λ 都有 dim E ( λ , T ) = mult ( λ ) \dim E(\lambda,T)=\operatorname{mult}(\lambda) dim E ( λ , T ) = mult ( λ ) 。
证明两句话。若处处相等,则 V = ⨁ k G ( λ k , T ) = ⨁ k E ( λ k , T ) V=\bigoplus_k G(\lambda_k,T)=\bigoplus_k E(\lambda_k,T) V = ⨁ k G ( λ k , T ) = ⨁ k E ( λ k , T ) (维数相同、后者含于前者,必相等),于是各 E E E 的基拼起来就是 V V V 的基,全是特征向量,故可对角化。反之若可对角化,V V V 有特征向量基,即 V = ⨁ k E ( λ k , T ) V=\bigoplus_k E(\lambda_k,T) V = ⨁ k E ( λ k , T ) 且维数和为 dim V \dim V dim V ;而 V = ⨁ k G ( λ k , T ) V=\bigoplus_k G(\lambda_k,T) V = ⨁ k G ( λ k , T ) 的维数和也是 dim V \dim V dim V ,逐项比较即得相等。
拿 3.2 的例子对一下:dim E ( 6 , T ) = 1 < 2 = mult ( 6 ) \dim E(6,T)=1<2=\operatorname{mult}(6) dim E ( 6 , T ) = 1 < 2 = mult ( 6 ) ,所以不可对角化。拿 1.1 的例子:dim E ( 3 , T ) = 1 < 2 = mult ( 3 ) \dim E(3,T)=1<2=\operatorname{mult}(3) dim E ( 3 , T ) = 1 < 2 = mult ( 3 ) ,同样不可对角化。判据用起来只需算两个零空间。
3.6 矩阵形态:分块对角 #
把 3.3 的三条翻成矩阵语言。对每个 k k k ,取 G ( λ k , T ) G(\lambda_k,T) G ( λ k , T ) 的一组基,使幂零算子 ( T − λ k I ) ∣ G ( λ k , T ) (T-\lambda_kI)|_{G(\lambda_k,T)} ( T − λ k I ) ∣ G ( λ k , T ) 的矩阵为上三角且对角全零(19 讲对幂零算子给的就是这个形态)。因为
T ∣ G ( λ k , T ) = ( T − λ k I ) ∣ G ( λ k , T ) + λ k I , T|_{G(\lambda_k,T)}=(T-\lambda_kI)|_{G(\lambda_k,T)}+\lambda_kI, T ∣ G ( λ k , T ) = ( T − λ k I ) ∣ G ( λ k , T ) + λ k I ,
T T T 在那一块上的矩阵就是「对角 λ k \lambda_k λ k 的上三角 d k × d k d_k\times d_k d k × d k 矩阵」,d k = mult ( λ k ) d_k=\operatorname{mult}(\lambda_k) d k = mult ( λ k ) 。把各块的基按顺序拼起来(3.3(c) 保证拼成 V V V 的基),得到
M = ( A 1 0 ⋱ 0 A m ) , A k = ( λ k ∗ ⋱ 0 λ k ) , A k 为 d k × d k . M=\begin{pmatrix}A_1&&0\\&\ddots&\\0&&A_m\end{pmatrix},\qquad A_k=\begin{pmatrix}\lambda_k&&*\\&\ddots&\\0&&\lambda_k\end{pmatrix},\ \ A_k\ \text{为}\ d_k\times d_k. M = A 1 0 ⋱ 0 A m , A k = λ k 0 ⋱ ∗ λ k , A k 为 d k × d k .
这就是 1.2 里承诺的形态:块数 = 互异特征值的个数,块的大小 = 各重数,块内只有一种对角元。 3.2 末尾那个矩阵就是它的一个实例。
四、几何直观 / 记忆钩子 #
本讲没有配图,用三个画面代替。
4.1 漏斗:U 1 ⊆ U 2 ⊆ ⋯ U_1\subseteq U_2\subseteq\cdots U 1 ⊆ U 2 ⊆ ⋯ #
把 ( T − λ I ) (T-\lambda I) ( T − λ I ) 想成一台往下吸的泵 。U 1 U_1 U 1 是「一吸就到底」的向量,U 2 U_2 U 2 是「吸两次到底」的,依此类推。链
{ 0 } = U 0 ⊆ U 1 ⊆ U 2 ⊆ ⋯ ⊆ U dim V = G ( λ , T ) \{0\}=U_0\subseteq U_1\subseteq U_2\subseteq\cdots\subseteq U_{\dim V}=G(\lambda,T) { 0 } = U 0 ⊆ U 1 ⊆ U 2 ⊆ ⋯ ⊆ U d i m V = G ( λ , T )
像漏斗的一层层:越往下口径越大,一旦某一层不再变大,下面的层就全都一样了。G ( λ , T ) G(\lambda,T) G ( λ , T ) 就是最终会落到漏斗底部的所有向量 。
E ( λ , T ) = U 1 E(\lambda,T)=U_1 E ( λ , T ) = U 1 :一层到底的,最陡。
mult ( λ ) = dim U dim V \operatorname{mult}(\lambda)=\dim U_{\dim V} mult ( λ ) = dim U d i m V :漏斗最宽处的口径。
所以「重数 ≥ 几何重数」在这个画面里就是「漏斗是往上变宽的」。
4.2 局部图像:T = λ I + N T=\lambda I+N T = λ I + N ,N N N 幂零 #
在 G ( λ , T ) G(\lambda,T) G ( λ , T ) 内部,把 T T T 写成
T = λ I + N , N = ( T − λ I ) ∣ G ( λ , T ) 幂零 . T=\lambda I+N,\qquad N=(T-\lambda I)|_{G(\lambda,T)}\ \text{幂零}. T = λ I + N , N = ( T − λ I ) ∣ G ( λ , T ) 幂零 .
这是本讲最有用的一句话:在一个广义特征空间里,算子只有两种动作——把所有方向均匀拉 λ \lambda λ 倍,再沿某个方向擦掉一层高阶信息。 没有第三种。
好处立刻体现:因为 I I I 与 N N N 可交换,二项式定理照用,N d = 0 N^d=0 N d = 0 让展开截断:
T m = ( λ I + N ) m = ∑ j = 0 d − 1 ( m j ) λ m − j N j ( d = 幂零指数 ) . T^m=(\lambda I+N)^m=\sum_{j=0}^{d-1}\binom{m}{j}\lambda^{m-j}N^j\qquad(d=\text{幂零指数}). T m = ( λ I + N ) m = j = 0 ∑ d − 1 ( j m ) λ m − j N j ( d = 幂零指数 ) .
若 N = 0 N=0 N = 0 (即这一块是特征空间),右边只剩 λ m I \lambda^mI λ m I ,回到对角化的情形。若 N ≠ 0 N\neq0 N = 0 ,多出来的项是「m m m 的低次多项式 × 幂零算子」——它增长得比 λ m \lambda^m λ m 慢,但在 λ = 0 \lambda=0 λ = 0 时它就是全部。这就是为什么幂零块求高次幂仍然便宜 :只需算 d d d 项,d d d 不超过块的边长。
4.3 房间图:按特征值分房 #
V = ⨁ k G ( λ k , T ) V=\bigoplus_k G(\lambda_k,T) V = ⨁ k G ( λ k , T ) 的画面是:把空间隔成 m m m 个房间,第 k k k 间房门上写着 λ k \lambda_k λ k 。
跨房间没有路 :T T T 不会把一个房间的向量送到另一个房间(不变性)。矩阵里的那些零就是墙。
房间大小是固定的 :dim G ( λ k , T ) \dim G(\lambda_k,T) dim G ( λ k , T ) 由 T T T 决定,不随你挑基而变。房间总面积正好是 dim V \dim V dim V 。
可对角化 = 每个房间都「一层到底」 :每间房里没有抖动,N = 0 N=0 N = 0 ,于是每个房间的基向量全是特征向量。
4.4 四句口诀 #
口诀
含义
E E E 是停住的,G G G 是停稳的
E v = 0 Ev=0 E v = 0 是一步停;G G G 是若干步后停
E E E 管能不能对角化,G G G 管能不能分块
判据看 E E E ,分解看 G G G
重数是房间面积,几何重数是门口宽度
mult = dim G \operatorname{mult}=\dim G mult = dim G ,dim E \dim E dim E 是入口
房间里只有「缩放 + 幂零」
T = λ I + N T=\lambda I+N T = λ I + N ,没有第三种动作
4.5 一个「没有房间」的对照 #
不是每个算子都能分房。R 2 \mathbb{R}^2 R 2 上逆时针转 90 ∘ 90^\circ 9 0 ∘ 的算子 T T T ,矩阵是 ( 0 − 1 1 0 ) \begin{pmatrix}0&-1\\1&0\end{pmatrix} ( 0 1 − 1 0 ) ,特征多项式 z 2 + 1 z^2+1 z 2 + 1 在 R \mathbb{R} R 上没有零点,没有任何实特征值 。于是每个 G ( λ , T ) G(\lambda,T) G ( λ , T ) 都是 { 0 } \{0\} { 0 } ,R 2 \mathbb{R}^2 R 2 分不出房间。
但把数域换成 C \mathbb{C} C ,情况立刻不同:特征值是 ± i \pm\mathrm{i} ± i ,而 T T T 可对角化(两个一维特征空间),C 2 = G ( i , T ) ⊕ G ( − i , T ) \mathbb{C}^2=G(\mathrm{i},T)\oplus G(-\mathrm{i},T) C 2 = G ( i , T ) ⊕ G ( − i , T ) 。同一个算子,能不能分解取决于在哪个数域上看。 这一点第五节还会踩一次。
五、常见误区 #
误区 1:把 G ( λ , T ) G(\lambda,T) G ( λ , T ) 当成 null ( T − λ I ) \operatorname{null}(T-\lambda I) null ( T − λ I ) 这是最常见的错,因为它「看起来」和定义只差一个指数。反例就是 1.1 的 A = ( 3 1 0 3 ) A=\begin{pmatrix}3&1\\0&3\end{pmatrix} A = ( 3 0 1 3 ) :null ( T − 3 I ) = span ( e 1 ) \operatorname{null}(T-3I)=\operatorname{span}(e_1) null ( T − 3 I ) = span ( e 1 ) ,只有一维;而 G ( 3 , T ) = C 2 G(3,T)=\mathbb{C}^2 G ( 3 , T ) = C 2 ,有两维。漏掉幂次会把整个分解定理读成一句假话(V V V 根本拼不起来)。
误区 2:以为「不同特征值的广义特征空间互不相交」是显然的 「λ \lambda λ 不同,所以 G G G 不同,所以不交」——这个推理跳了步。子空间不交(交为 { 0 } \{0\} { 0 } )是要证的,而且不能靠维数:R 3 \mathbb{R}^3 R 3 里两个过原点的平面维数都是 2 2 2 ,必然交于一条线。正确的依据是 3.3 的线性无关性引理,它用到的是「( T − λ m I ) (T-\lambda_mI) ( T − λ m I ) 在别的 G G G 上可逆」这个具体的代数事实 ,不是「特征值不同」这四个字。
误区 3:在 R \mathbb{R} R 上无条件使用分解定理 主定理要求 F = C \mathbb{F}=\mathbb{C} F = C ,不能省。4.5 的旋转算子就是反例:在 R \mathbb{R} R 上它没有特征值,V ≠ ⨁ G V\neq\bigoplus G V = ⨁ G (左边二维,右边零维)。C \mathbb{C} C 上之所以总成立,靠的是每个复多项式都能分解成一次因式——这是 15 讲的内容,不是线性代数自动给的。 在 R \mathbb{R} R 上的正确条件是「极小多项式在 R \mathbb{R} R 上完全分解」,见第六节。
有人写「G ( λ , T ) = null ( T − λ I ) d G(\lambda,T)=\operatorname{null}(T-\lambda I)^{d} G ( λ , T ) = null ( T − λ I ) d ,其中 d d d 是重数」——这句话本身是对的(第六节会证),但不能拿来当初次计算的公式 :重数被定义成 dim G ( λ , T ) \dim G(\lambda,T) dim G ( λ , T ) ,而 G G G 又要用幂次算,这是转圈。第一次算永远取 k = dim V k=\dim V k = dim V ;等算出 d d d 之后,才知道取 d d d 也够。
代数重数(= 本讲的重数)是 dim G \dim G dim G ,管的是「块有多大」;几何重数是 dim E \dim E dim E ,管的是「特征向量有几个方向」。两者永远满足 dim E ≤ dim G \dim E\le\dim G dim E ≤ dim G ,等号成立才可对角化。若记反,就会得出「dim E ≤ dim G \dim E\le\dim G dim E ≤ dim G 意味着特征向量比广义特征向量多」这种荒谬结论。记法:几何的多,是「看得见的」;代数的大,是「算得出来的」。 这里「几何重数小、代数重数大」才是常态。
误区 6:以为 G ( λ , T ) G(\lambda,T) G ( λ , T ) 只有在 λ \lambda λ 是特征值时才需要算 反过来更要注意:对非 特征值的 λ \lambda λ ,G ( λ , T ) = { 0 } G(\lambda,T)=\{0\} G ( λ , T ) = { 0 } (2.4 已证)。如果计算中你算出某个「广义特征空间」有正维数,那 λ \lambda λ 一定是特征值——这是一个很好的自查手段。算错了幂次或矩阵时,常常会在这里露馅。
六、拓展(可选,可跳过) #
6.1 R \mathbb{R} R 上的正确条件 #
主定理在 R \mathbb{R} R 上失败,是因为复数特征值「无处安放」。补上这一块,条件可以写得很干净:
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 ) ⟺ \iff ⟺ T T T 的极小多项式 在 F \mathbb{F} F 上完全分解成一次因式之积。
C \mathbb{C} C 上「完全分解」自动成立(15 讲),所以主定理是它的特例。R \mathbb{R} R 上的旋转算子极小多项式是 z 2 + 1 z^2+1 z 2 + 1 ,不可分解,故分解失败——与 4.5 的观察一致。
为什么条件恰好是这一条?回看 3.4 的证明:那里只用了「存在一个在 F \mathbb{F} F 上分解成一次因式、且零化 T T T 的多项式」。极小多项式正是次数最低的那个零化多项式(18 讲),它若能分解,整条证明原样照搬;它若不能分解,就有一个不可约因式「没有特征值可对应」,分解必然差一块。
6.2 幂次可以取到重数,而不是 dim V \dim V dim V #
由 (b),( T − λ I ) (T-\lambda I) ( T − λ I ) 在 G ( λ , T ) G(\lambda,T) G ( λ , T ) 上幂零,其幂零指数不超过 dim G ( λ , T ) = mult ( λ ) \dim G(\lambda,T)=\operatorname{mult}(\lambda) dim G ( λ , T ) = mult ( λ ) 。于是
G ( λ , T ) = null ( T − λ I ) mult ( λ ) . G(\lambda,T)=\operatorname{null}(T-\lambda I)^{\operatorname{mult}(\lambda)}. G ( λ , T ) = null ( T − λ I ) mult ( λ ) .
这比 dim V \dim V dim V 紧,但要注意上面的循环问题:它适合验算 ,不适合第一遍计算。例如 3.2 里 mult ( 6 ) = 2 \operatorname{mult}(6)=2 mult ( 6 ) = 2 ,而 null ( T − 6 I ) 2 \operatorname{null}(T-6I)^2 null ( T − 6 I ) 2 确实已经稳定,不必算到 3 3 3 。
6.3 凯莱–哈密顿定理的三行证明 #
有了分块结构,这条经典定理几乎是免费的。设 q ( z ) = ∏ k ( z − λ k ) d k q(z)=\prod_k(z-\lambda_k)^{d_k} q ( z ) = ∏ k ( z − λ k ) d k ,d k = mult ( λ k ) d_k=\operatorname{mult}(\lambda_k) d k = mult ( λ k ) 。由 6.2,( T − λ k I ) d k (T-\lambda_kI)^{d_k} ( T − λ k I ) d k 在 G ( λ k , T ) G(\lambda_k,T) G ( λ k , T ) 上是零算子。又 V V V 是这些 G G G 的直和,故只需验证 q ( T ) q(T) q ( T ) 在每一块上为零。在 G ( λ k , T ) G(\lambda_k,T) G ( λ k , T ) 上,q ( T ) q(T) q ( T ) 的各因式可交换,把 ( T − λ k I ) d k (T-\lambda_kI)^{d_k} ( T − λ k I ) d k 挪到最右边,整串乘积就是零:
q ( T ) ∣ G ( λ k , T ) = ( T − λ 1 I ) d 1 ⋯ ( T − λ k I ) d k ⋯ ( T − λ m I ) d m ∣ G ( λ k , T ) = 0. q(T)|_{G(\lambda_k,T)}=(T-\lambda_1I)^{d_1}\cdots(T-\lambda_kI)^{d_k}\cdots(T-\lambda_mI)^{d_m}\Big|_{G(\lambda_k,T)}=0. q ( T ) ∣ G ( λ k , T ) = ( T − λ 1 I ) d 1 ⋯ ( T − λ k I ) d k ⋯ ( T − λ m I ) d m G ( λ k , T ) = 0.
于是 q ( T ) = 0 q(T)=0 q ( T ) = 0 。关键就一句:算子多项式彼此可交换,所以能把「致命的那一项」挪到最后。 这也顺带说明 q q q 是零化多项式,即特征多项式是最小多项式的倍式(18 讲)。
6.4 重数 = 上三角矩阵对角线上出现的次数 #
把 T T T 上三角化后,每个特征值在对角线上出现若干次。这个次数等于重数,思路是:上三角矩阵对角线上零的个数 d d d 限制了零空间的维数(dim null ≤ d \dim\operatorname{null}\le d dim null ≤ d ),把这条估计用到 ( T − λ I ) (T-\lambda I) ( T − λ I ) 与它的各次幂上,得到「重数 ≤ \le ≤ 出现次数」;而两边对所有特征值求和都等于 dim V \dim V dim V ,于是逐项只能是等号。这类「先证不等号、再靠总和相等逼出等号」的手法很值得记。
6.5 再往前一步 #
Jordan 形 :把每个 A k A_k A k 再细分,可以做到块内只有紧邻上对角线上有 1 1 1 、其余为零。本讲的分块对角是它的粗版本。
D + N D+N D + N 分解 :T T T 可以写成一个可对角化算子 D D D 与一个幂零算子 N N N 之和,且 D N = N D DN=ND D N = N D 。构造很直接:在每个 G ( λ k , T ) G(\lambda_k,T) G ( λ k , T ) 上取 D = λ k I D=\lambda_kI D = λ k I 、N = T − λ k I N=T-\lambda_kI N = T − λ k I ,再拼起来。这条分解在解线性微分方程组时是标准工具——e t T = e t D e t N e^{tT}=e^{tD}e^{tN} e tT = e t D e tN ,而后一项是有限和。
七、自测(先做,再看答案) #
1. 设 T T T 在标准基下的矩阵是 A = ( 2 1 0 0 2 0 0 0 3 ) A=\begin{pmatrix}2&1&0\\0&2&0\\0&0&3\end{pmatrix} A = 2 0 0 1 2 0 0 0 3 。求 G ( 2 , T ) G(2,T) G ( 2 , T ) 与 G ( 3 , T ) G(3,T) G ( 3 , T ) 及其重数,并判断 T T T 是否可对角化。
2. 证明:G ( λ , T ) ≠ { 0 } G(\lambda,T)\neq\{0\} G ( λ , T ) = { 0 } 当且仅当 λ \lambda λ 是 T T T 的特征值。
3. 证明 dim E ( λ , T ) ≤ dim G ( λ , T ) \dim E(\lambda,T)\le\dim G(\lambda,T) dim E ( λ , T ) ≤ dim G ( λ , T ) ,并说明「对每个特征值都取等号」等价于可对角化。
4. 设 T T T 可逆,λ ≠ 0 \lambda\neq0 λ = 0 。证明 G ( λ , T ) = G ( 1 / λ , T − 1 ) G(\lambda,T)=G(1/\lambda,T^{-1}) G ( λ , T ) = G ( 1/ λ , T − 1 ) 。
5. 构造一个 C 4 \mathbb{C}^4 C 4 上的算子,使 mult ( 2 ) = 3 \operatorname{mult}(2)=3 mult ( 2 ) = 3 、mult ( 5 ) = 1 \operatorname{mult}(5)=1 mult ( 5 ) = 1 ,并写出它的一组由广义特征向量组成的基。
答案(点开) 1. 令 N = A − 2 I = ( 0 1 0 0 0 0 0 0 1 ) N=A-2I=\begin{pmatrix}0&1&0\\0&0&0\\0&0&1\end{pmatrix} N = A − 2 I = 0 0 0 1 0 0 0 0 1 。null N \operatorname{null}N null N :z 2 = 0 , z 3 = 0 z_2=0,z_3=0 z 2 = 0 , z 3 = 0 ,得 span ( e 1 ) \operatorname{span}(e_1) span ( e 1 ) ,维数 1 1 1 。N 2 = ( 0 0 0 0 0 0 0 0 1 ) N^2=\begin{pmatrix}0&0&0\\0&0&0\\0&0&1\end{pmatrix} N 2 = 0 0 0 0 0 0 0 0 1 ,零空间是 z 3 = 0 z_3=0 z 3 = 0 ,即 span ( e 1 , e 2 ) \operatorname{span}(e_1,e_2) span ( e 1 , e 2 ) ,维数 2 2 2 ;N 3 = N 2 N^3=N^2 N 3 = N 2 ,故稳定。 所以 G ( 2 , T ) = span ( e 1 , e 2 ) G(2,T)=\operatorname{span}(e_1,e_2) G ( 2 , T ) = span ( e 1 , e 2 ) ,mult ( 2 ) = 2 \operatorname{mult}(2)=2 mult ( 2 ) = 2 ,而 dim E ( 2 , T ) = 1 \dim E(2,T)=1 dim E ( 2 , T ) = 1 。 对 λ = 3 \lambda=3 λ = 3 :A − 3 I = ( − 1 1 0 0 − 1 0 0 0 0 ) A-3I=\begin{pmatrix}-1&1&0\\0&-1&0\\0&0&0\end{pmatrix} A − 3 I = − 1 0 0 1 − 1 0 0 0 0 ,零空间为 span ( e 3 ) \operatorname{span}(e_3) span ( e 3 ) ,故 G ( 3 , T ) = span ( e 3 ) G(3,T)=\operatorname{span}(e_3) G ( 3 , T ) = span ( e 3 ) ,mult ( 3 ) = 1 \operatorname{mult}(3)=1 mult ( 3 ) = 1 。 维数和 2 + 1 = 3 2+1=3 2 + 1 = 3 ✓。因 dim E ( 2 , T ) = 1 < 2 = mult ( 2 ) \dim E(2,T)=1<2=\operatorname{mult}(2) dim E ( 2 , T ) = 1 < 2 = mult ( 2 ) ,不可对角化 。
2. (⇐ \Leftarrow ⇐ )若 λ \lambda λ 是特征值,取特征向量 v ≠ 0 v\neq0 v = 0 ,则 ( T − λ I ) v = 0 (T-\lambda I)v=0 ( T − λ I ) v = 0 ,即 v ∈ G ( λ , T ) v\in G(\lambda,T) v ∈ G ( λ , T ) 。 (⇒ \Rightarrow ⇒ )设 0 ≠ v ∈ G ( λ , T ) 0\neq v\in G(\lambda,T) 0 = v ∈ G ( λ , T ) ,取使 ( T − λ I ) k v = 0 (T-\lambda I)^kv=0 ( T − λ I ) k v = 0 的最小 正整数 k k k 。由最小性 w = ( T − λ I ) k − 1 v ≠ 0 w=(T-\lambda I)^{k-1}v\neq0 w = ( T − λ I ) k − 1 v = 0 ,且 ( T − λ I ) w = 0 (T-\lambda I)w=0 ( T − λ I ) w = 0 ,故 w w w 是特征向量,λ \lambda λ 是特征值。
3. E ( λ , T ) = null ( T − λ I ) ⊆ null ( T − λ I ) dim V = G ( λ , T ) E(\lambda,T)=\operatorname{null}(T-\lambda I)\subseteq\operatorname{null}(T-\lambda I)^{\dim V}=G(\lambda,T) E ( λ , T ) = null ( T − λ I ) ⊆ null ( T − λ I ) d i m V = G ( λ , T ) ,维数不等式随之成立。 若每个特征值都取等号:由主分解定理 V = ⨁ k G ( λ k , T ) = ⨁ k E ( λ k , T ) V=\bigoplus_kG(\lambda_k,T)=\bigoplus_kE(\lambda_k,T) V = ⨁ k G ( λ k , T ) = ⨁ k E ( λ k , T ) ,各 E E E 的基拼成 V V V 的基,全是特征向量,故可对角化。 若可对角化:V V V 有特征向量基,即 V = ⨁ k E ( λ k , T ) V=\bigoplus_kE(\lambda_k,T) V = ⨁ k E ( λ k , T ) 且 ∑ dim E = dim V \sum\dim E=\dim V ∑ dim E = dim V ;又 V = ⨁ k G ( λ k , T ) V=\bigoplus_kG(\lambda_k,T) V = ⨁ k G ( λ k , T ) 且 ∑ dim G = dim V \sum\dim G=\dim V ∑ dim G = dim V 。逐项 dim E ≤ dim G \dim E\le\dim G dim E ≤ dim G 而总和相等,故处处相等。
4. 注意 ( T − λ I ) = T ( I − λ T − 1 ) = − λ T ( T − 1 − λ − 1 I ) (T-\lambda I)=T(I-\lambda T^{-1})=-\lambda\,T\,(T^{-1}-\lambda^{-1}I) ( T − λ I ) = T ( I − λ T − 1 ) = − λ T ( T − 1 − λ − 1 I ) (最后一个等号:I − λ T − 1 = − λ ( T − 1 − λ − 1 I ) I-\lambda T^{-1}=-\lambda(T^{-1}-\lambda^{-1}I) I − λ T − 1 = − λ ( T − 1 − λ − 1 I ) ,两边左乘 T T T 即得)。于是
( T − λ I ) k = ( − λ ) k T k ( T − 1 − λ − 1 I ) k . (T-\lambda I)^k=(-\lambda)^kT^k\big(T^{-1}-\lambda^{-1}I\big)^k. ( T − λ I ) k = ( − λ ) k T k ( T − 1 − λ − 1 I ) k .
T k T^k T k 可逆,故 ( T − λ I ) k v = 0 ⟺ ( T − 1 − λ − 1 I ) k v = 0 (T-\lambda I)^kv=0\iff(T^{-1}-\lambda^{-1}I)^kv=0 ( T − λ I ) k v = 0 ⟺ ( T − 1 − λ − 1 I ) k v = 0 。对任意 k k k 成立,取 k = dim V k=\dim V k = dim V 即得两个广义特征空间相等。
5. 取 A = ( 2 1 0 0 0 2 1 0 0 0 2 0 0 0 0 5 ) A=\begin{pmatrix}2&1&0&0\\0&2&1&0\\0&0&2&0\\0&0&0&5\end{pmatrix} A = 2 0 0 0 1 2 0 0 0 1 2 0 0 0 0 5 (一个 3 × 3 3\times3 3 × 3 的幂零块配 2 2 2 ,加一个 1 × 1 1\times1 1 × 1 的块配 5 5 5 )。 令 N = A − 2 I N=A-2I N = A − 2 I ,则 N ( z 1 , z 2 , z 3 , z 4 ) = ( z 2 , z 3 , 0 , 3 z 4 ) N(z_1,z_2,z_3,z_4)=(z_2,\ z_3,\ 0,\ 3z_4) N ( z 1 , z 2 , z 3 , z 4 ) = ( z 2 , z 3 , 0 , 3 z 4 ) 。于是
N 2 ( z ) = ( z 3 , 0 , 0 , 9 z 4 ) , N 3 ( z ) = ( 0 , 0 , 0 , 27 z 4 ) , N 4 ( z ) = ( 0 , 0 , 0 , 81 z 4 ) . N^2(z)=(z_3,\ 0,\ 0,\ 9z_4),\qquad N^3(z)=(0,\ 0,\ 0,\ 27z_4),\qquad N^4(z)=(0,0,0,81z_4). N 2 ( z ) = ( z 3 , 0 , 0 , 9 z 4 ) , N 3 ( z ) = ( 0 , 0 , 0 , 27 z 4 ) , N 4 ( z ) = ( 0 , 0 , 0 , 81 z 4 ) .
零空间依次是 span ( e 1 ) \operatorname{span}(e_1) span ( e 1 ) (维数 1 1 1 )、span ( e 1 , e 2 ) \operatorname{span}(e_1,e_2) span ( e 1 , e 2 ) (2 2 2 )、span ( e 1 , e 2 , e 3 ) \operatorname{span}(e_1,e_2,e_3) span ( e 1 , e 2 , e 3 ) (3 3 3 ),到 N 4 N^4 N 4 不再增长。故 G ( 2 , T ) = span ( e 1 , e 2 , e 3 ) G(2,T)=\operatorname{span}(e_1,e_2,e_3) G ( 2 , T ) = span ( e 1 , e 2 , e 3 ) ,重数 3 3 3 。 对 λ = 5 \lambda=5 λ = 5 :A − 5 I A-5I A − 5 I 的第 4 4 4 行全为零,而前 3 3 3 行依次给出 z 3 = 0 z_3=0 z 3 = 0 、z 2 = 0 z_2=0 z 2 = 0 、z 1 = 0 z_1=0 z 1 = 0 ,故零空间为 span ( e 4 ) \operatorname{span}(e_4) span ( e 4 ) ,重数 1 1 1 。 维数和 3 + 1 = 4 3+1=4 3 + 1 = 4 ✓。一组基:e 1 , e 2 , e 3 , e 4 e_1,e_2,e_3,e_4 e 1 , e 2 , e 3 , e 4 ——注意 e 2 , e 3 e_2,e_3 e 2 , e 3 不是特征向量(T e 2 = 2 e 2 + e 1 Te_2=2e_2+e_1 T e 2 = 2 e 2 + e 1 ),但它们是广义特征向量。
八、复习卡 #
项
内容
定义
G ( λ , T ) = { v : ( T − λ I ) k v = 0 , k 为某正整数 } G(\lambda,T)=\{v:(T-\lambda I)^kv=0,\ k\ \text{为某正整数}\} G ( λ , T ) = { v : ( T − λ I ) k v = 0 , k 为某正整数 }
可算形式
G ( λ , T ) = null ( T − λ I ) dim V G(\lambda,T)=\operatorname{null}(T-\lambda I)^{\dim V} G ( λ , T ) = null ( T − λ I ) d i m V (链最多 dim V \dim V dim V 步稳定)
链稳定
U k = U k + 1 ⇒ U j = U k ( j ≥ k ) U_k=U_{k+1}\Rightarrow U_j=U_k\ (j\ge k) U k = U k + 1 ⇒ U j = U k ( j ≥ k ) ;维数每步至少 + 1 +1 + 1
子空间性
因为并等于某一项 U dim V U_{\dim V} U d i m V ;不能直接说「零空间的并」
包含关系
E ( λ , T ) ⊆ G ( λ , T ) E(\lambda,T)\subseteq G(\lambda,T) E ( λ , T ) ⊆ G ( λ , T ) ,等号一般不成立
何时为零
G ( λ , T ) = { 0 } ⟺ λ G(\lambda,T)=\{0\}\iff\lambda G ( λ , T ) = { 0 } ⟺ λ 不是特征值
不变性
G ( λ , T ) = null p ( T ) G(\lambda,T)=\operatorname{null}p(T) G ( λ , T ) = null p ( T ) ,p ( z ) = ( z − λ ) dim V p(z)=(z-\lambda)^{\dim V} p ( z ) = ( z − λ ) d i m V ,故 T T T -不变
幂零性
( T − λ I ) ∣ G ( λ , T ) (T-\lambda I)\big|_{G(\lambda,T)} ( T − λ I ) G ( λ , T ) 幂零,指数 ≤ mult ( λ ) \le\operatorname{mult}(\lambda) ≤ mult ( λ )
主分解
C \mathbb{C} C 上: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 )
直和靠什么
不同特征值的广义特征向量线性无关(不能用维数糊)
张满靠什么
上三角化 ⇒ \Rightarrow ⇒ 零化多项式 ⇒ \Rightarrow ⇒ 互素分解 v = ∑ a k ( T ) p k ( T ) v v=\sum a_k(T)p_k(T)v v = ∑ a k ( T ) p k ( T ) v
重数
mult ( λ ) = dim G ( λ , T ) \operatorname{mult}(\lambda)=\dim G(\lambda,T) mult ( λ ) = dim G ( λ , T ) ;∑ k mult ( λ k ) = dim V \sum_k\operatorname{mult}(\lambda_k)=\dim V ∑ k mult ( λ k ) = dim V
代数 vs 几何
代数重数 = dim G =\dim G = dim G ;几何重数 = dim E =\dim E = dim E ;总有 dim E ≤ dim G \dim E\le\dim G dim E ≤ dim G
可对角化判据
∀ λ : dim E ( λ , T ) = mult ( λ ) \forall\lambda:\ \dim E(\lambda,T)=\operatorname{mult}(\lambda) ∀ λ : dim E ( λ , T ) = mult ( λ )
矩阵形态
分块对角,第 k k k 块为对角元全 λ k \lambda_k λ k 的 d k × d k d_k\times d_k d k × d k 上三角
块内图像
T = λ I + N T=\lambda I+N T = λ I + N ,N N N 幂零;T m = ∑ j < d ( m j ) λ m − j N j T^m=\sum_{j<d}\binom{m}{j}\lambda^{m-j}N^j T m = ∑ j < d ( j m ) λ m − j N j
R \mathbb{R} R 上条件
极小多项式完全分解 ⟺ V = ⨁ G \iff V=\bigoplus G ⟺ V = ⨁ G
反面例子
R 2 \mathbb{R}^2 R 2 转 90 ∘ 90^\circ 9 0 ∘ :无实特征值,G G G 全为零,分解失败
九、出处与延伸(可选,不影响阅读) #
主要出处 :Sheldon Axler《线性代数应该这样学》(第四版,吴俊达、何阳译)第 8 章 8B 节「广义特征空间分解」,印刷页 257–265(PDF p271–279)。本节的定义、分解定理、重数、特征多项式与分块对角矩阵都在这个范围内。
术语别名 :「广义特征空间」在部分文献里叫「根空间」(root space),「广义特征向量」叫「根向量」;「重数」在旧文献里叫「代数重数」,与之并列的「几何重数」是 dim E ( λ , T ) \dim E(\lambda,T) dim E ( λ , T ) 。不同书对「重数」默认指哪一个并不统一,读别的材料时要看它的定义式。
重数的行列式定义 :多数教材用「特征多项式中 ( z − λ ) (z-\lambda) ( z − λ ) 的指数」定义重数,需要先有行列式;本讲的路线用 dim null ( T − λ I ) dim V \dim\operatorname{null}(T-\lambda I)^{\dim V} dim null ( T − λ I ) d i m V ,不需要行列式,两者最终等价。用行列式算重数在数值上手快,但它不告诉你 G G G 里长什么样。
R \mathbb{R} R 上的版本 :8B 节习题 8 给出「V = ⨁ G V=\bigoplus G V = ⨁ G 当且仅当极小多项式完全分解」的等价刻画,正文只提了一句「R \mathbb{R} R 情形的类似结论见习题」。第六节 6.1 把它补全了。
想深挖 :Jordan 形与本讲的直和分解是同一件事的两种精度;若关心数值计算,可查「广义特征向量」在矩阵函数(e t A e^{tA} e t A 、A 1 / 2 A^{1/2} A 1/2 )计算中的用法,那里 4.2 的二项式展开是核心公式。
一处提醒 :本讲所有「上三角化」「极小多项式」的引用都只借结论、不重述证明,需要回看时按正文标出的讲次编号去查。