LinearAlgebra-Chapter5:特征值和特征向量/矩阵对角化
这一章的核心问题是:
能否通过更换一组基,把一个复杂矩阵化成对角矩阵,从而把矩阵运算转化成对角元素的运算?
若存在可逆矩阵 P,使
P−1AP=Λ=diag(λ1,λ2,…,λn),则称 A 可对角化。此时:
- λ1,λ2,…,λn 是 A 的特征值;
- P 的各列是与这些特征值依次对应的特征向量;
- Ak、矩阵多项式 f(A) 等均可借助对角矩阵快速计算;
- 当 A 为实对称矩阵时,可以进一步选取正交矩阵 U 完成对角化。
本章的主线为:
- 求特征值;
- 求各特征值对应的特征子空间;
- 判断能否选出 n 个线性无关的特征向量;
- 构造可逆矩阵 P 或正交矩阵 U;
- 利用对角化计算矩阵的方幂、矩阵多项式和行列式。
NOTE教学范围
两次课程完整讲授了教材中的:
- 5.1 特征值和特征向量
- 5.2 矩阵对角化
- 5.3 矩阵相似的理论和应用
- 5.4 实对称矩阵的对角化
教材 5.5 为全章概要与小结,课堂未将其作为新的知识板块单独讲授,因此本笔记不另设 5.5。
课堂明确略去的证明:
- 特征子空间维数不超过特征根重数的证明;
- 任意实对称矩阵均可正交对角化的完整证明。
这两个结论仍需掌握并会使用。
5.1 特征值和特征向量#
从矩阵对角化引出特征值#
设 A 为 n 阶方阵。若 A 可对角化,则存在可逆矩阵 P,使
P−1AP=Λ,Λ=diag(λ1,λ2,…,λn).两边左乘 P,得到
AP=PΛ.将 P 按列分块:
P=[ξ1,ξ2,…,ξn].由于 P 可逆,其列向量 ξ1,…,ξn 线性无关。另一方面,
AP=[Aξ1,Aξ2,…,Aξn],而
PΛ=[λ1ξ1,λ2ξ2,…,λnξn].因此
Aξi=λiξi,i=1,2,…,n.这说明矩阵能否对角化,最终取决于能否找到足够多的、满足 Aξ=λξ 的线性无关向量。
特征值与特征向量的定义#
设 A 为 n 阶方阵。若存在数 λ 和非零向量 ξ,使
Aξ=λξ,ξ=0,则称:
- λ 为 A 的一个 特征值;
- ξ 为 A 的属于特征值 λ 的一个 特征向量。
WARNING定义中必须有 ξ=0。
若允许 ξ=0,则对任意 λ 都有 A0=λ0,定义将失去意义。
特征值与特征向量的基本关系#
一个特征向量只能属于一个特征值#
若同一非零向量 ξ 同时满足
Aξ=λ1ξ,Aξ=λ2ξ,则
(λ1−λ2)ξ=0.因为 ξ=0,所以
λ1=λ2.一个特征值有无穷多个特征向量#
若
Aξ=λξ,ξ=0,则对任意非零数 k,
A(kξ)=kAξ=kλξ=λ(kξ).因此 kξ 仍是属于 λ 的特征向量。
这也说明:特征向量通常成族出现,真正需要确定的是其所在的特征子空间。
特征多项式与特征方程#
由
Aξ=λξ可得
(λE−A)ξ=0.要使该齐次线性方程组有非零解,必须且只需
∣λE−A∣=0.定义
fA(λ)=∣λE−A∣为 A 的 特征多项式,方程
fA(λ)=0称为 A 的 特征方程。
由于 fA(λ) 是首项系数为 1 的 n 次多项式,根据代数基本定理,在复数域中共有 n 个根,按重数计算。这些根正是 A 的全部特征值。
两种常见表述:
- λ1,λ2,…,λn 表示全部特征值,允许相等;
- λ1,λ2,…,λs 表示全部互异特征值,其重数分别为 n1,n2,…,ns,满足
n1+n2+⋯+ns=n.此时
fA(λ)=(λ−λ1)n1(λ−λ2)n2⋯(λ−λs)ns.特征子空间#
设 λi 为 A 的特征值。齐次线性方程组
(λiE−A)X=0的解空间记为
Wλi=ker(λiE−A),称为属于 λi 的 特征子空间。
其维数为
dimWλi=n−r(λiE−A).若 Wλi 的一组基为
ξi1,ξi2,…,ξiri,则属于 λi 的全部特征向量为
ξ=t1ξi1+t2ξi2+⋯+triξiri,其中 t1,t2,…,tri 不全为零。
TIP特征子空间中包含零向量;特征向量不包含零向量。
因此:
{属于 λi 的特征向量}=Wλi∖{0}.
求特征值和特征向量的标准步骤#
- 计算特征多项式
fA(λ)=∣λE−A∣.
- 解特征方程
fA(λ)=0,求出全部互异特征值 λ1,…,λs 及其重数。
- 对每个 λi,解
(λiE−A)X=0,求出特征子空间的一组基。
- 将基础解系作不全为零的线性组合,写出全部特征向量。
例 5.1.1:对角元相同、非对角元相同的矩阵#
设
A=ab⋮bba⋮bbb⋮b⋯⋯⋯bb⋮a.取
ξ=11⋮1,则每一行与 ξ 相乘的结果均为
a+(n−1)b.所以
Aξ=[a+(n−1)b]ξ.因此
λ=a+(n−1)b是 A 的一个特征值,ξ 是其对应的特征向量。
进一步可得
fA(λ)=[λ−a−(n−1)b](λ−a+b)n−1.故全部特征值为
λ1=a+(n−1)b,λ2=⋯=λn=a−b.当 a=0 时,特征值为
(n−1)b,−b(n−1 重).例 5.1.3:完整求特征值和特征向量#
设
A=−1−41130002.先求特征多项式:
fA(λ)=∣λE−A∣=λ+14−1−1λ−3000λ−2=(λ−2)(λ−1)2.因此
λ1=2,λ2=λ3=1.属于 λ=2 的特征向量#
解
(2E−A)X=0,即
34−1−1−10000x1x2x3=0.其基础解系可取
ξ1=001.所以属于 2 的全部特征向量为
kξ1,k=0.属于 λ=1 的特征向量#
解
(E−A)X=0,即
24−1−1−2000−1x1x2x3=0.其基础解系可取
ξ2=12−1.所以属于 1 的全部特征向量为
kξ2,k=0.实矩阵可能具有复特征值#
设
A=0−13102001.则
fA(λ)=(λ−1)(λ2+1).
- 若只在实数域中讨论,特征值只有 1;
- 若在复数域中讨论,全部特征值为
1,i,−i.因此,讨论特征值时必须明确所处数域。
5.2 矩阵对角化#
属于不同特征值的特征向量线性无关#
设 λ1,λ2,…,λs 是 A 的互异特征值,ξi 是属于 λi 的特征向量,则
ξ1,ξ2,…,ξs线性无关。
证明思路:数学归纳法#
当 s=1 时,单个非零向量线性无关。
假设 ξ1,…,ξk 线性无关。对 k+1 个向量,设
t1ξ1+t2ξ2+⋯+tkξk+tk+1ξk+1=0.(1)两边左乘 A,得
t1λ1ξ1+⋯+tkλkξk+tk+1λk+1ξk+1=0.(2)用 λk+1×(1) 减去 (2),消去最后一项:
i=1∑kti(λk+1−λi)ξi=0.由归纳假设,ξ1,…,ξk 线性无关,且
λk+1−λi=0,所以
t1=t2=⋯=tk=0.代回式 (1),由 ξk+1=0 得
tk+1=0.因此结论成立。
TIP若 n 阶矩阵有 n 个互异特征值,则必有 n 个线性无关的特征向量,因此必可对角化。
特征子空间基的合并仍线性无关#
设 λ1,…,λs 为 A 的互异特征值,Wλi 的一组基为
ξi1,ξi2,…,ξiri.则合并后的向量组
ξ11,…,ξ1r1,ξ21,…,ξ2r2,…,ξs1,…,ξsrs线性无关。
证明仍采用归纳法:把最后一个特征子空间对应的整组向量看作一组,分别对线性组合左乘 A 与乘以最后一个特征值,再作差,先消去最后一组;利用归纳假设推出前面各组系数为零,再利用最后一组本身线性无关推出剩余系数为零。
这个性质把“每个特征值取一个特征向量”的结论推广为“每个特征子空间取一组基”。
可对角化的三个充要条件#
设 A 为 n 阶方阵,互异特征值为
λ1,λ2,…,λs.对应特征子空间为 Wλi,记
ri=dimWλi.特征值 λi 的代数重数记为 ni。
以下三个条件彼此等价。
条件一:有 n 个线性无关的特征向量#
A 可对角化⟺A 有 n 个线性无关的特征向量.若这些特征向量为 ξ1,…,ξn,其对应特征值为 λ1,…,λn,则
P=[ξ1,ξ2,…,ξn]可逆,且
P−1AP=diag(λ1,λ2,…,λn).这个条件在理论证明中最直接。
条件二:全部特征子空间维数之和为 n#
A 可对角化⟺r1+r2+⋯+rs=n.因为各特征子空间的基合并后线性无关;当向量总数达到 n 时,恰好得到 n 个线性无关的特征向量。
这个条件最适合计算判定。
条件三:几何重数等于代数重数#
A 可对角化⟺ri=ni,i=1,2,…,s.由于
n1+n2+⋯+ns=n,所以条件二和条件三等价。
代数重数与几何重数#
对特征值 λi:
- ni:λi 作为特征方程根的重数,称为 代数重数;
- ri=dimWλi:对应特征子空间维数,称为 几何重数。
有基本不等式
1≤ri≤ni.课堂未证明上界 ri≤ni,但要求会用。
TIP若 λi 是单根,则 ni=1。又因为 1≤ri≤ni,所以必有
ri=1.因此判断可对角化时,只需重点检查重根。
判断可对角化的计算流程#
- 求全部特征值及其代数重数;
- 对每个重根 λi 计算
ri=n−r(λiE−A);
- 检查是否对每个 i 都有 ri=ni,或检查
r1+r2+⋯+rs=n;
- 若只要求判断,不必把基础解系完整求出,计算矩阵的秩即可;
- 若还要求求出 P,则需要解各齐次方程组,取各特征子空间的基作为 P 的列。
例 5.2.1:互异特征值保证可对角化#
设
B=−101130002.其特征多项式为
fB(λ)=(λ+1)(λ−2)(λ−3).所以特征值为
−1,2,3.它们互不相同,故 B 可对角化。
分别求得特征向量可取
ξ1=30−1,ξ2=001,ξ3=141.令
P=[ξ1,ξ2,ξ3],则
P−1BP=diag(−1,2,3).例 5.2.2:重根对应的特征向量不足#
仍取例 5.1.3 中的矩阵
A=−1−41130002.其特征值为
2,1(2 重).前面已经求得
dimW2=1,dimW1=1.所以
dimW2+dimW1=2<3.因此 A 不能对角化。
关键原因:特征值 1 的代数重数为 2,几何重数只有 1。
例 5.2.3:重根仍可能可对角化#
设
A=122212221.特征多项式为
fA(λ)=(λ−5)(λ+1)2.故特征值为
5,−1(2 重).属于 5 的特征子空间可取基
ξ1=111.属于 −1 的特征子空间可取基
ξ2=10−1,ξ3=01−1.因此
dimW5=1,dimW−1=2,且
1+2=3.所以 A 可对角化。令
P=[ξ1,ξ2,ξ3]=11110−101−1,则
P−1AP=diag(5,−1,−1).NOTE有重根并不能直接推出矩阵不可对角化。真正要检查的是重根对应的特征子空间维数是否达到其代数重数。
例 5.2.4:利用重数直接求秩#
已知三阶矩阵 A 的特征值为
1,2,2,且 A 可对角化,求
r(E−A),r(2E−A).特征值 1 为单根,因此
dimW1=1.又
dimW1=3−r(E−A),所以
r(E−A)=2.特征值 2 为二重根。由于 A 可对角化,几何重数等于代数重数,因此
dimW2=2.于是
3−r(2E−A)=2,故
r(2E−A)=1.例 5.2.5:含参数的秩一矩阵#
设
α=1234,β=3−2−1a,A=αβT.问 a 取何值时 A 可对角化。
先计算
fA(λ)=λ3(λ−4a+4).所以特征值为
0(3 重),4(a−1).又因 A=αβT 且 α,β=0,所以
r(A)=1.当 a=1 时#
特征值为
0(3 重),4(a−1)=0(1 重).对特征值 0,
dimW0=4−r(A)=3,恰好等于其代数重数。另一个特征值为单根,故 A 可对角化。
当 a=1 时#
唯一特征值为
0(4 重).但
dimW0=4−r(A)=3<4.因此 A 不能对角化。
结论:
a=1时 A 可对角化。
例 5.2.6:非零幂零矩阵不能对角化#
设 A=0,且存在正整数 k≥2,使
Ak=0.证明 A 不能对角化。
设 λ 是 A 的任一特征值,ξ 是对应的非零特征向量,则
Aξ=λξ.连续作用 A,得到
Akξ=λkξ.由 Ak=0,有
λkξ=0.因为 ξ=0,所以
λ=0.即 A 的全部特征值均为 0。
若 A 可对角化,则存在可逆矩阵 P,使
P−1AP=diag(0,0,…,0)=0.于是
A=0,与题设 A=0 矛盾。因此 A 不能对角化。
5.3 矩阵相似的理论和应用#
特征多项式的两个重要系数#
设
A=[aij]n×n,全部特征值为 λ1,…,λn,按重数计算。
特征多项式可写成
fA(λ)=∣λE−A∣=λn−(a11+a22+⋯+ann)λn−1+⋯+(−1)n∣A∣.另一方面,
fA(λ)=(λ−λ1)(λ−λ2)⋯(λ−λn).比较 λn−1 的系数和常数项,得到
λ1+λ2+⋯+λn=trA以及
λ1λ2⋯λn=∣A∣.其中
trA=a11+a22+⋯+ann称为矩阵 A 的迹。
WARNING特征值之和、乘积均按重数计算。
由乘积公式立即得到
A 可逆⟺∣A∣=0⟺0 不是 A 的特征值.矩阵方幂与矩阵多项式的特征值#
若
Aξ=λξ,ξ=0,则
A2ξ=A(Aξ)=A(λξ)=λAξ=λ2ξ.逐次递推可得
Akξ=λkξ.因此:
- λk 是 Ak 的特征值;
- ξ 仍是对应的特征向量。
设
g(x)=akxk+⋯+a1x+a0,则
g(A)=akAk+⋯+a1A+a0E.于是
g(A)ξ=akAkξ+⋯+a1Aξ+a0ξ=[akλk+⋯+a1λ+a0]ξ=g(λ)ξ.所以 g(λ) 是 g(A) 的特征值,ξ 仍是相应特征向量。
在复数域中有谱映射关系
σ(g(A))=g(σ(A)).也就是说,g(A) 的全部特征值由 A 的特征值代入 g 得到。若不同的 λi 代入后得到同一个数,对应特征值会合并。
为什么 A2 的特征值不会遗漏#
设 μ 是 A2 的特征值,则
∣μE−A2∣=0.在复数域中取 μ,有
μE−A2=(μE−A)(μE+A).因此
∣μE−A∣⋅∣μE+A∣=0.所以 μ 或 −μ 是 A 的特征值,从而
μ=λ2可由 A 的某个特征值 λ 平方得到。
逆矩阵的特征值#
设 A 可逆,且
Aξ=λξ.因为可逆矩阵的特征值均非零,可在两边左乘 A−1:
ξ=λA−1ξ.因此
A−1ξ=λ−1ξ.所以:
σ(A−1)={λ1:λ∈σ(A)}且特征向量保持不变。
相似矩阵的共同性质#
若 A 与 B 相似,则存在可逆矩阵 P,使
B=P−1AP.相同的特征多项式#
∣λE−B∣=∣λE−P−1AP∣=∣P−1(λE−A)P∣=∣P−1∣∣λE−A∣∣P∣=∣λE−A∣.因此相似矩阵有:
- 相同的特征多项式;
- 相同的全部特征值及其代数重数;
- 相同的迹;
- 相同的行列式;
- 相同的秩。
逆命题通常不成立#
例如
E=[1001],B=[1011].二者的特征多项式均为
(λ−1)2,但 E 只与自身相似,而 B=E,所以 E 与 B 不相似。
因此,相同特征多项式只是矩阵相似的必要条件。
相同特征值何时能推出相似#
若 A,B 都可对角化,则
A∼B⟺A,B 有相同的特征值及重数.理由:若二者都可对角化,并具有相同的特征值,则它们都相似于同一个对角矩阵,利用相似关系的对称性与传递性可得 A∼B。
特别地,两个实对称矩阵一定可对角化,所以判断两个实对称矩阵是否相似,只需比较其特征值及重数。
相似关系与矩阵多项式#
若
B=P−1AP,则
B2=P−1A2P,B3=P−1A3P,…,Bk=P−1AkP.故
Ak∼Bk.对任意多项式
g(x)=akxk+⋯+a1x+a0,有
g(B)=P−1g(A)P.因此
g(A)∼g(B).而且使用的是同一个相似变换矩阵 P。
利用对角化计算方幂与多项式#
若
P−1AP=Λ=diag(λ1,…,λn),则
A=PΛP−1.所以
Ak=PΛkP−1=Pdiag(λ1k,…,λnk)P−1.对多项式 g,
g(A)=Pg(Λ)P−1=Pdiag[g(λ1),…,g(λn)]P−1.进而
∣g(A)∣=g(λ1)g(λ2)⋯g(λn).TIP计算高次方时,本章提供了第二种常用情形:
- 若 r(A)=1,可利用 A=αβT 的结构;
- 若 A 可对角化,可利用 Ak=PΛkP−1。
例 5.3.1:矩阵多项式的特征值、行列式与矩阵值#
设
A=12−221−22−21.第一步:求 A 的特征值和特征向量#
fA(λ)=(λ+1)(λ−1)(λ−3).所以特征值为
−1,1,3.对应特征向量可取
ξ1=1−10,ξ2=1−11,ξ3=01−1.由于三个特征值互异,A 可对角化。令
P=[ξ1,ξ2,ξ3],则
P−1AP=diag(−1,1,3).第二步:求矩阵多项式的特征值#
令
g(x)=2x3−3x2−4x+2.则
g(A)=2A3−3A2−4A+2E.其特征值为
g(−1)=1,g(1)=−3,g(3)=17.相应特征向量仍分别为 ξ1,ξ2,ξ3。
第三步:求行列式#
∣g(A)∣=g(−1)g(1)g(3)=1×(−3)×17=−51.第四步:求矩阵 g(A)#
g(A)=Pdiag(1,−3,17)P−1.计算得
g(A)=−320−20−421−20−44−3.例 5.3.2:伴随矩阵与矩阵多项式#
设三阶方阵 A 满足
∣A∣=6,且 A 有一个特征值 −2。
求 A∗ 的一个特征值#
由 ∣A∣=0,A 可逆,且
A∗=∣A∣A−1=6A−1.因为 −2 是 A 的特征值,所以
−21是 A−1 的特征值。因此
6(−21)=−3是 A∗ 的一个特征值。
求 A3+4A2+8A+8E 的一个特征值#
取
g(x)=x3+4x2+8x+8.则
g(−2)=0.所以 0 是
A3+4A2+8A+8E的一个特征值。
因此
∣A3+4A2+8A+8E∣=0.例 5.3.5:计算高次方#
设三阶方阵 A,已知
ξ1=120,ξ2=230,ξ3=002,并满足
Aξ1=ξ1,Aξ2=−ξ2,Aξ3=2ξ3.因此特征值为
1,−1,2,互不相同,A 可对角化。
令
P=[ξ1,ξ2,ξ3]=120230002,Λ=diag(1,−1,2).则
A=PΛP−1=−7−120470002.同时
A1000=Pdiag(11000,(−1)1000,21000)P−1.因为
11000=(−1)1000=1,计算得
A1000=1000100021000.
5.4 实对称矩阵的对角化#
设 A 为实矩阵且
AT=A,则称 A 为 实对称矩阵。
实对称矩阵是一类结构很强的矩阵:
- 特征值全部为实数;
- 不同特征值对应的特征向量彼此正交;
- 一定可以通过正交矩阵对角化。
实对称矩阵的特征值均为实数#
设 λ 是实对称矩阵 A 的特征值,ξ=0 为相应特征向量:
Aξ=λξ.即使 ξ 暂时取复向量,也有
ξ∗Aξ=λξ∗ξ.(1)另一方面,由 A∗=A,
ξ∗Aξ=(Aξ)∗ξ=(λξ)∗ξ=λξ∗ξ.(2)比较式 (1) 与式 (2):
(λ−λ)ξ∗ξ=0.因为 ξ=0,所以
ξ∗ξ>0.从而
λ=λ.因此 λ 为实数。
又因为 λE−A 是实矩阵,其齐次方程组可在实数域中求解,所以相应特征向量可以取为实向量。
不同特征值对应的特征向量正交#
设
Aξ1=λ1ξ1,Aξ2=λ2ξ2,λ1=λ2.则
(Aξ1)Tξ2=λ1ξ1Tξ2.另一方面,利用 AT=A:
(Aξ1)Tξ2=ξ1TATξ2=ξ1TAξ2=λ2ξ1Tξ2.所以
(λ1−λ2)ξ1Tξ2=0.由于 λ1=λ2,故
ξ1Tξ2=0.即
ξ1⊥ξ2.相比一般矩阵“不同特征值对应的特征向量线性无关”,这里的结论更强:它们直接正交。
实对称矩阵的正交对角化定理#
若 A 为 n 阶实对称矩阵,则存在正交矩阵 U,使
U−1AU=UTAU=diag(λ1,λ2,…,λn).其中 λ1,…,λn 是 A 的全部特征值,按重数排列。
课堂未证明该定理,只要求理解其含义并掌握构造 U 的方法。
TIP对一般可对角化矩阵,通常只能保证存在可逆矩阵 P。
对实对称矩阵,可以把 P 进一步选成正交矩阵 U,于是
U−1=UT.
求正交矩阵的标准步骤#
设 A 为实对称矩阵。
第一步:求全部互异特征值#
求解
∣λE−A∣=0,得到
λ1,λ2,…,λs及其代数重数 n1,…,ns。
第二步:求各特征子空间的基#
对每个 λi 解
(λiE−A)X=0,得到基础解系
ξi1,ξi2,…,ξini.因为实对称矩阵一定可对角化,所以每个特征子空间维数必等于对应特征值的重数。若计算结果不相等,说明前面的特征值、秩或方程组求解出现错误。
第三步:在每个特征子空间内部正交化、单位化#
- 若 ni=1,只需将唯一基向量单位化;
- 若 ni>1,先对该特征子空间的一组基作施密特正交化,再分别单位化。
得到标准正交基
ηi1,ηi2,…,ηini.不同特征子空间之间本来就彼此正交,因此只需处理同一重根内部的向量。
第四步:按特征值顺序组成正交矩阵#
U=[η11,…,η1n1,η21,…,ηsns].则
UTAU=diag(n1λ1,…,λ1,…,nsλs,…,λs).NOTE若题目只要求求一个可逆矩阵 P 使 P−1AP 为对角矩阵,到第二步后直接将各特征子空间的基按列排入 P 即可,无须正交化和单位化。
例 5.4.1:三个互异特征值#
设
A=110110001.特征多项式为
fA(λ)=λ(λ−1)(λ−2).所以特征值为
0,1,2.相应特征向量可取
ξ1=1−10,ξ2=001,ξ3=110.三个特征值互异,三个特征向量已经两两正交,只需单位化:
η1=211−10,η2=001,η3=21110.令
U=[η1,η2,η3],则 U 为正交矩阵,且
UTAU=diag(0,1,2).例 5.4.2:重根情形需要施密特正交化#
设
A=122212221.由例 5.2.3,特征值为
5,−1(2 重).可取
ξ1=111作为 W5 的基;取
ξ2=10−1,ξ3=01−1作为 W−1 的一组基。
ξ2 与 ξ3 属于同一特征值,未必正交,因此对它们施密特正交化。
取
ξ2′=ξ2.再令
ξ3′=ξ3−(ξ2′,ξ2′)(ξ3,ξ2′)ξ2′=01−1−2110−1=21−12−1.单位化后得到
η1=31111,η2=2110−1,η3=61−12−1.因此
U=[η1,η2,η3]=313131210−21−6162−61.于是
UTAU=diag(5,−1,−1).
课堂综合例题#
正交矩阵的行列式为负时必有特征值 -1#
设 U 为实正交矩阵,且
∣U∣=−1.证明 −1 是 U 的特征值。
只需证明
∣−E−U∣=0.由正交性
UTU=E,有
−E−U=(−UT−E)U.因此
∣−E−U∣=∣−UT−E∣∣U∣=−∣−UT−E∣.又因为行列式在转置下不变,
∣−UT−E∣=∣(−U−E)T∣=∣−U−E∣=∣−E−U∣.所以
∣−E−U∣=−∣−E−U∣,从而
∣−E−U∣=0.故 −1 是 U 的特征值。
幂等矩阵的行列式计算#
设 A 为 n 阶方阵,满足
A2=A,且
r(A)=r,0<r≤n.求
∣5E+A∣.第一步:确定 A 的可能特征值#
若 λ 是 A 的特征值,则由
g(x)=x2−x及
g(A)=A2−A=0可得
g(λ)=λ2−λ=0.所以
λ=0或λ=1.第二步:证明 A 可对角化#
属于 0 的特征子空间维数为
dimW0=n−r(A).属于 1 的特征子空间维数为
dimW1=n−r(E−A).需要证明
dimW0+dimW1=n,等价于
r(A)+r(E−A)=n.一方面,由
A+(E−A)=E及秩的次可加性,
n=r(E)≤r(A)+r(E−A).另一方面,由 Sylvester 秩不等式,
r(A)+r(E−A)≤n+r[A(E−A)].而
A(E−A)=A−A2=0,所以
r(A)+r(E−A)≤n.综合可得
r(A)+r(E−A)=n.因此 A 可对角化。
第三步:确定重数#
因为
dimW0=n−r(A)=n−r,且 A 可对角化,所以特征值 0 的代数重数为 n−r;特征值 1 的代数重数为 r。
第四步:计算目标行列式#
矩阵
5E+A的特征值为
5+0=5(n−r 重),以及
5+1=6(r 重).因此
∣5E+A∣=5n−r6r.
全章解题框架#
类型一:求全部特征值和特征向量#
- 写出 λE−A;
- 计算 ∣λE−A∣;
- 因式分解并写清重数;
- 对每个互异特征值解 (λiE−A)X=0;
- 用基础解系的不全为零线性组合表示全部特征向量。
类型二:判断能否对角化#
优先使用:
i∑dimWλi=n或
dimWλi=ni.计算时:
dimWλi=n−r(λiE−A).
- 单根自动满足;
- 只检查重根;
- 只需判断时可只算秩;
- 还要求 P 时再求基础解系。
类型三:构造相似变换矩阵 P#
将各特征子空间的一组基按列排列:
P=[线性无关的特征向量].对角矩阵中各特征值的顺序必须与 P 中相应特征向量的列顺序一致。
类型四:计算 Ak 或 g(A)#
若
P−1AP=Λ,则
Ak=PΛkP−1,g(A)=Pg(Λ)P−1.计算行列式时可直接用
∣g(A)∣=i=1∏ng(λi).类型五:实对称矩阵正交对角化#
- 求特征值;
- 求各特征子空间的基;
- 同一重根对应的基作施密特正交化;
- 全部向量单位化;
- 按列组成 U;
- 写出
UTAU=Λ.
易错点#
-
特征向量必须非零。
-
特征子空间包含零向量,特征向量集合不包含零向量。
-
有重根不代表一定不能对角化。 应比较几何重数和代数重数。
-
有 n 个特征值不代表有 n 个互异特征值。 特征值总数按重数计算始终为 n。
-
求 P 时,列向量顺序要和对角矩阵中特征值的顺序对应。
-
相同特征多项式不能直接推出矩阵相似。 若两矩阵都可对角化,才可通过相同特征值推出相似。
-
g(A) 的特征向量可沿用 A 的特征向量,但不同特征值经 g 映射后可能合并。
-
实矩阵的特征值可能为复数。 实对称矩阵的特征值才保证全为实数。
-
实对称矩阵中,不同特征值对应的特征向量已经正交。 施密特正交化只需在同一重根对应的特征子空间内部进行。
-
计算特征多项式时,∣λE−A∣ 与 ∣A−λE∣ 相差 (−1)n。 求根相同,但后续公式与教材符号应保持一致。