上一節的零空間判準
det(A−λI)=0
已經告訴你如何測試一個候選數 λ 是否為特徵值。下一步,是把這個測試包成一個單一的多項式,而它的根正是矩陣的特徵值。
特徵多項式限制特徵值的數目,記錄重複特徵值的情形,並給出可對角化的實用判準。
本節固定以 F=R 為基礎域。除非明確討論複數擴張,矩陣、向量、純量與特徵值都在 F 中;所有行列式均來自 n≥1 的方陣。多項式未定元統一使用 t,並約定 pA(t)=det(A−tIn)。
特徵多項式
設 A∈Mn(F),其中 F=R 且 n≥1。多項式
pA(t)=det(A−tIn)稱為 A 的 特徵多項式
變數 t 在這裏只是未定元。一旦你把某個純量 λ 代入進去,
pA(λ) 便會變成 A−λI 的行列式。
設 A 是 n×n 矩陣,而 λ 是純量。下列兩件事等價:
- λ 是 A 的特徵值;
- λ 是 pA(t) 的根
這其實只是把上一節的行列式判準整理成較有系統的寫法:
λ 是特徵值⟺det(A−λI)=0⟺pA(λ)=0.
根確定特徵值,特徵空間提供方向
pA(t)=det(A−tI) 的根確定特徵值;相應特徵空間的維數則決定能否找到足夠多的線性無關特徵向量來對角化。
由特徵根走到對角化測試觀看 det(A-xI) 如何變成特徵多項式,重根如何分成代數與幾何重數,以及特徵空間維數如何決定可否對角化。
行列式流程
特徵值測試 det(A-lambda I)=0 會被打包成一個多項式 p_A(x)=det(A-xI),而它的根正正就是特徵值。
因式分解讀根
對 A=[[3,2],[3,-2]],行列式化簡為 x^2-x-12=(x-4)(x+3),所以特徵值是 4 與 -3。
代數重數
對 B=[[2,3],[0,2]],p_B(x)=(x-2)^2。重根表示特徵值 2 的代數重數是 2。
幾何重數
解 (B-2I)v=0 得到 E_B(2)=span{(1,0)^T},所以幾何重數只有 1。
維數總和測試
對角化成功,當且僅當所有特徵空間的維數總和等於矩陣大小 n。
互異特徵值捷徑
如果 n 乘 n 矩陣有 n 個互異特徵值,對應特徵向量會自動線性無關,所以矩陣可對角化。
特徵多項式找出可能的對角線項;但對角化仍要問特徵空間能否提供完整特徵基底。
若 A 是 n×n 矩陣,則 pA(t) 是一個 n 次多項式,其最高次項系數是 (−1)n,而常數項是 det(A)
因此,特徵多項式絕不是任意的表達式。它的次數由矩陣大小固定,而其常數項已經記住了矩陣的行列式。
把 A 的欄寫成 a1,…,an。A−tIn 的各欄是 aj−tej。由欄的多重線性,行列式展開為每一欄選擇 aj 或 −tej 的和。所有欄都選擇 −tej
時,得到
(−t)ndet(In)=(−1)ntn其餘選擇至少含有一個 aj,所以至多含有 n−1 個 t 因子。因此次數正好是
n,最高次項系數是 (−1)n。不含 t 的項是 det(A),所以 pA(0)=det(A)。這裏堅持使用 det(A−tIn);另一約定 det(tIn−A)=(−1)npA(t),根與重數相同,但奇數次數的系數符號會不同
第一批例子
令
A=[332−2].則
pA(t)=det[3−t32−2−t]=(3−t)(−2−t)−6.化簡後得
pA(t)=t2−t−12=(t−4)(t+3).這個係數計算值得直接檢查,而不應只記住因式分解的形式:
(3−t)(−2−t)−6=t2−t−6−6=t2−t−12按本節約定 det(A−tI),常數項是 −12,與 det(A)=−12 一致;最高次項係數是 +1,正好符合二階矩陣的
(−1)n 符號。於是兩個因子給出根 4 與 −3,根定理再把這些純量根解釋為特徵值。
令
B=[2032].則
pB(t)=det[2−t032−t]=(2−t)2=(t−2)2.所以 2 是唯一特徵值,但它以重數 2 出現
代數重數與幾何重數
一旦根可以重複,便需要更精細的語言。
設 λ∈F 是 pA(t) 的根。它的 代數重數 ma(λ) 是使
(t−λ)m 整除 pA(t) 的最大正整數 m。等價地,
pA(t)=(t−λ)mq(t),其中 q(λ)=0。它的 幾何重數 mg(λ)
是 dimFEA(λ)=dimFN(A−λIn)。
若 pA 在 F 上完全分裂,可以把互異的根列作 λ1,…,λs,並有
pA(t)=(−1)n(t−λ1)m1⋯(t−λs)ms,i∑mi=n這個分解與重數總和要求所選域上完全分裂;局部的重數定義本身不需要這個條件
代數重數來自多項式本身;幾何重數來自零空間 N(A−λI)。它們衡量的不是同一件事,而且未必相等。
若 λ 是 n×n 矩陣 A 的特徵值,則
1≤mg(λ)≤ma(λ)≤n.
因此,每個特徵空間至少有一維,但它的維數永遠不會超過相應特徵值的代數重數。
固定特徵值 λ,令 g=dimEA(λ)。取該特徵空間的一組基,並把它擴展為
Fn 的一組基。令 P 的欄組成這組擴展基。矩陣 B=P−1AP 的前 g 欄是
λe1,…,λeg,因此
B=[λIg0CD]相似變換保持特徵多項式。若 g<n,B−tI 的前 g 欄各自至多有一個可能非零的元素,即 λ−t;當 t=λ 時該元素也為零,連續沿這些欄展開這個多項式恆等式,得到
pA(t)=det(B−tI)=(λ−t)gdet(D−tI)所以 (t−λ)g 整除 pA(t),從而 g≤ma(λ)。若 g=n,則
B=λIn,且 pA(t)=(λ−t)n。這個證明使用前面建立的擴展基與相似變換,不假設特徵多項式完全分裂
對
B=[2032],特徵多項式是 (t−2)2,所以特徵值 2 的代數重數是 2。
現在解 (B−2I)x=0:
B−2I=[0030].因此 x2=0,x1 自由。特徵空間為
span{[10]},只有一維。也就是說
ma(2)=2,mg(2)=1.這個不相等正是它不可對角化的原因。
互異特徵值會強迫線性無關
定理 8.3.9
(互異特徵值所對應的特徵向量必線性無關)
若 v1,…,vk 是 A 的特徵向量,而它們分別對應於互不相同的特徵值 λ1,…,λk,則
v1,…,vk必定線性無關
設 λ1,…,λs 兩兩不同,且 wj∈EA(λj)。我們用 s 歸納證明:若 j=1∑swj=0,則每個 wj=0,允許其中某些 wj 本來就是零向量。
s=1 時結論顯然。歸納步驟中,對關係式應用 A−λsI,最後一項消失,得到
j=1∑s−1(λj−λs)wj=0每個數乘後的向量仍在自己的特徵空間中。歸納假設使所有
(λj−λs)wj 為零;互異性使這些純量非零,所以 j<s 時 wj=0,再由原關係得到 ws=0。把按特徵值分組的任意關係應用這個結論,再利用每個特徵空間內部基的獨立性,就得到所有系數都為零;因此各特徵空間的基的並集線性無關
這條定理立即帶來兩個重要推論。
一個 n×n 矩陣最多只可擁有 n 個互異特徵值
若一個 n×n 矩陣恰好有 n 個互異特徵值,則它必定可對角化
前面的獨立性結果同時給出兩項推論:互異特徵值最多有 n 個;若恰好有 n 個,從每個特徵值選取一個非零特徵向量,就得到長度為 n 的線性無關列表,因而是 Fn 的一組基。
其逆命題是假的。可對角化矩陣仍然可以有重複特徵值。最簡單的例子就是單位矩陣:每一個非零向量都是特徵值 1 的特徵向量,而矩陣本身早已是對角矩陣。
更精確的可對角化判準
互異特徵值是一個充分條件,但不是唯一條件。
把所選域 F 中 A 的所有互異特徵值列作 λ1,…,λs;這個列表可以為空。則 A 在 F 上可對角化,當且僅當
dimFEA(λ1)+⋯+dimFEA(λs)=n
這項判準是一個基底計數,並不要求特徵多項式完全分裂。特徵空間獨立性證明說明,各特徵空間基的並集線性無關;若維數總和是 n,它就是一組基,因此可對角化。反過來,對角化基底可按對角元素的特徵值分組,每個向量落在相應特徵空間中。各特徵空間基的並集線性無關,維數至多是 n;它的張成空間包含對角化基底中的每個向量,因此張成整個 Fn,維數至少是 n。所以總和正好是 n。若 F 中沒有特徵值,總和是 0,而 n≥1 的空間不可能在 F 上對角化。
更精細的重數判準需要完全分裂:
若 pA 在 F 上完全分裂,則 A 在 F 上可對角化,當且僅當對 F 中的每個特徵值
λ 都有 mg(λ)=ma(λ)
完全分裂時,i∑ma(λi)=n,而重數不等式逐項給出
mg(λi)≤ma(λi)。因此每一對重數都相等,正好等價於維數總和為
n,再應用前面的判準即可。不完全分裂時,這個條件不能安全替代維數總和判準。在實數上,
R=[11−11],pR(t)=(t−1)2+1
沒有實特徵值。「每個實特徵值都滿足 mg=ma」會變成空命題,但 R 不能在 R 上對角化。在 C 上,根 1+i 與 1−i 互異,帶有域條件的判準才適用。
令
A=211121112令 d=2−t。沿第一行展開時,第一個餘子式為
det[d11d]=d2−1,符號為正;第二個餘子式為
det[111d]=d−1,餘子式符號為負;第三個餘子式為
det[11d1]=1−d,符號為正。因此三個帶符號的項分別是
d(d2−1)、−(d−1) 和 +(1−d),所以
d(d2−1)−(d−1)+(1−d)=d3−3d+2=(d−1)2(d+2)代回 d=2−t 得
pA(t)=(1−t)2(4−t)=−(t−1)2(t−4)這個因式分解顯示根 1 出現兩次、根 4 出現一次;三次多項式的奇次數也解釋了負的最高次項係數。由於多項式次數為三,這些就是它在 R 中的全部根,後面的特徵空間計數因此涵蓋了這個矩陣的每個實特徵值。
當 λ=1 時,解 (A−I)x=0 得 x+y+z=0,所以
EA(1)=span⎩⎨⎧1−10,10−1⎭⎬⎫當 λ=4 時,方程給出 x=y=z,故 EA(4)=span{(1,1,1)T}。
EA(1) 中的兩個向量由第二、第三座標立即看出線性無關。再用互異特徵空間引理,它們與
EA(4) 中的向量合起來仍然線性無關,並且
dimEA(1)+dimEA(4)=2+1=3這裏也可以直接檢查零空間:A−I 的三行均給出 x+y+z=0,所以維數是 2;而 A−4I
的方程等價於 x=y=z,所以維數是 1。這與多項式中根 1 的重數 2 和根 4 的重數 1
完全一致;三個方向合起來構成 F3 的基。
因此 A 雖有重根,仍然可對角化
Cayley-Hamilton 定理給出一條多項式恆等式
本節再向前一步,把特徵多項式重新代回矩陣本身。
定理 8.3.15
(Cayley-Hamilton 定理)
若
pA(t)=c0+c1t+c2t2+⋯+cntn,則
c0I+c1A+c2A2+⋯+cnAn=0.等價寫法是 pA(A)=0。這個定理是一般性的;下面只證明可對角化情形,對包括不可對角化矩陣在內的一般情形,見 MATH1550 的係數比較證明。其先修是該證明前已建立的交換環 F[t] 上的行列式與伴隨恆等式。那裏採用 χA(t)=det(tI−A)=(−1)npA(t),所以 χA(A)=0 與本節約定等價。方法是在 (tI−A)adj(tI−A)=χA(t)I 中比較係數,再左乘 A 的冪並求和相消;不能直接在含矩陣係數的恆等式中代入 t=A。讀完後返回這裏,比較可對角化情形的初等證明和 Jordan 例的直接驗證。
對可對角化矩陣來說,這條定理很好理解:若 A=SDS−1,那麼把多項式 pA 套到 A 上,等於把它逐項套到 D 的對角線各項上,而每個對角線項本身都是 A 的特徵值,因此也是 pA 的根。
證明
(可對角化情形的 Cayley-Hamilton 證明)
假設 A=SDS−1,其中 D=diag(λ1,…,λn)。對純量多項式
q(t)=c0+⋯+crtr,定義 q(A)=c0I+c1A+⋯+crAr;這是對矩陣冪作標量代入,不是在行列式內部作形式代換。因為冪保持相似關係,q(A)=Sq(D)S−1,而
q(D) 的對角元素是 q(λ1),…,q(λn)。取 q=pA 後,根據特徵值—
根定理這些對角元素全為零。常數項使用 I=SIS−1,所以 pA(A)=0
Cayley–Hamilton 恆等式把高次冪化為低次冪。若
pA(t)=c0+c1t+⋯+cn−1tn−1+cntn,
則 cn=(−1)n=0,所以可以解出 An,再依次乘以 A 並化簡,得到更高次冪。當 c0=det(A)=0 時,A 可逆;把恆等式乘以 A−1 並整理,得到
A−1=−c0c1I+c2A+⋯+cnAn−1.
這裏必須先由 c0=0 確定可逆性,才能使用 A−1。
例題 8.3.16
(一條小型 Cayley-Hamilton 恆等式)
令
J=[1041]其特徵多項式為
pJ(t)=(1−t)2=t2−2t+1因此 Cayley-Hamilton 告訴你
J2−2J+I=0直接計算
J2=[1081],J2−2J+I=0把這個已驗證的恆等式乘以 J,得到
J3=2J2−J=2(2J−I)−J=3J−2I更具體地,將 J2−2J+I=0 乘以 Jk,對每個 k≥0 都得到遞推式
Jk+2=2Jk+1−Jk。因此 J4=2J3−J2=4J−3I,無需從頭計算較大的矩陣乘積。這正是多項式恆等式提供的具體化簡機制。事實上,對每個正整數 m,
Jm=mJ−(m−1)I=[104m1]m=1 時這是 J 的定義。若公式對 m 成立,則
Jm+1=mJ2−(m−1)J=m(2J−I)−(m−1)J=(m+1)J−mI,於是下一個情形也成立。這說明即使 J 不可對角化,化簡仍然可以具體進行。
現在可以具體應用前面的係數公式:這裏 c0=1、c1=−2、c2=1,所以
J−1=−1−2I+J=2I−J=[10−41]兩個消去次序都可以直接檢查:
J(2I−J)=2J−J2=I(2I−J)J=2J−J2=I
常見錯誤
反例模式
重複特徵值仍可能可對角化
錯誤說法:任何重複特徵值都會阻止對角化。取 A=2I2,則
pA(t)=(2−t)2,所以 2 的代數重數是 2;但 EA(2)=R2 的基是
e1,e2,矩陣本身已經是對角矩陣。前面的 B 具有相同的重根模式,卻只有一維特徵空間,所以不可對角化。修正方法是在已聲明的分裂條件下比較 mg(λ) 與
ma(λ),或在所選域上使用不附帶分裂條件的維數總和判準
快速檢查
快速檢查
([2005] 的特徵多項式是甚麼?)
用 pA(t)=det(A−tI)
解答 · 答案
pA(t)=2−t005−t=(2−t)(5−t).
快速檢查
(若一個 n×n 矩陣有 n 個互異特徵值,你可以立刻得出甚麼結論?)
解答 · 答案
它必定可對角化,因為互異特徵值所對應的特徵向量必線性無關
解答 · 答案
不會。對每個特徵值 λ,都有
mg(λ)≤ma(λ).
練習
快速檢查
(求 [0−213] 的特徵多項式,並列出其特徵值。)
先算 det(A−tI),再因式分解
解答 · 引導解答
pA(t)=−t−213−t=t2−3t+2=(t−1)(t−2).所以特徵值是 1 與 2
快速檢查
(一個實 4×4 矩陣有四個互異的實特徵值。若從每個特徵值各選一個特徵向量,這四個向量所張成空間的維數是多少?)
解答 · 引導解答
這四個特徵向量線性無關,因此它們張成的空間維數是 4。所以它們本身已構成 R4 的一組基底
快速檢查
(若 pA(t)=−t3+6t2−11t+6,Cayley-Hamilton 會給出甚麼矩陣恆等式?)
解答 · 引導解答
−A3+6A2−11A+6I=0或者兩邊乘以 −1:
A3−6A2+11A−6I=0這就是對應的 Cayley-Hamilton 恆等式
相關筆記
請把
8.1 特徵值、特徵向量與特徵空間
打開,因為特徵值與零空間的關係是這一節的基礎。
同時也應保留
8.2 對角化與相似
在旁,因為本節提供的多項式工具正是用來判斷那一節的對角化是否成立。
而 pA(t)=det(A−tI) 的行列式部分,則依賴
7.1 行列式與餘因子展開。