行列式会告诉你一个方阵是否可逆。特征值问的是另一个问题:哪些非零向量
在矩阵作用之下仍然保持它所张成的直线?若
Av=λv,则直线 span{v} 会被映入自身。当
λ>0 时,直线上的方向保持;当 λ<0 时,方向反转;当
λ=0 时,非零向量被送到零向量,而零向量没有方向。三种情况都
可以用“直线映入自身”来统一描述。
一般矩阵不必保持每一条直线。标量矩阵会把每一条直线映入自身,包括零
标量的情形:0I 会把所有向量都映到零向量。固定一个标量 λ 时,
非零特征向量加上零向量才形成对应的特征空间;不同特征值的特征向量合
在一起,未必构成子空间。特征向量很有用,因为矩阵对它的作用完全由一
个标量说明。这个标量记录拉长、缩短、反向或压到零的程度。
寻找矩阵作用最简单的方向
当你解 Ax=b 时,你关心的是整个系统。当你研究特征值时,你关心的是
满足下式的特殊向量 v:
Av=λv.
一旦找出这些向量,之后便能更容易处理矩阵的幂、理解对角化,并描述哪
些方向在变换下是不变的。实际计算时,我们会把它改写成齐次系统,再使用
零空间、秩和行列式的工具。
先看 A=diag(2,3)。水平向量 e1 变成 2e1,竖直向量
e2 变成 3e2,各自仍在原来的直线上;但 (1,1)T 变成 (2,3)T,
已经不在同一条过原点的直线上。因此,寻找特征向量,就是寻找矩阵只用一个
缩放因子便能描述其作用的特殊方向。
定义
特征值与特征向量
设 A 是一个实 n×n 方阵,其中 n≥1。设 λ∈R,
而 v∈Rn 是非零列向量。
若
Av=λv,便称 v 是 A 对应于 特征值 λ 的 特征向量。
这里必须要求 v 非零。因为零向量对每个标量都满足 A0=λ0,
若把它允许进来,每个标量都会看起来像特征值。下面定义的特征空间会包含
零向量,但“特征向量”一词始终只指非零向量。
重点在于,矩阵作用与标量乘法在这个向量上完全一致。向量的直线被保留,
但长度、方向,甚至向量本身是否被压到零,都由 λ 决定。
从不变直线到齐次方程组
方程 Av=λv 等价于 (A−λI)v=0。固定一个特征值 λ 后,非零解是特征向量;把零向量也包括进来,就得到特征空间。
由定义立刻得到的结论
定理
同一个特征向量只会对应一个特征值
设 A 是方阵,而 v 是非零向量。若
Av=λv且Av=μv,则 λ=μ。
把两式相减得到
(λ−μ)v=0.
选择一个非零坐标 vr=0。上式的第 r 个坐标为
(λ−μ)vr=0。由于 vr 是非零实数,可以约去它,得到
λ−μ=0,所以 λ=μ。非零条件正是这一步约去的依据。
定理
特征向量的非零倍数仍是特征向量
若 v 是 A 对应于特征值 λ 的特征向量,则每个非零标量倍
cv 也是 A 的特征向量,而且对应同一个特征值 λ。
事实上,若 c=0,则
A(cv)=c(Av)=c(λv)=λ(cv)
而 cv=0,所以它满足定义。经过一个特征向量的直线上所有非零点都
是特征向量,而不是只有一个孤立的向量。
定理
同一特征值下的线性组合
设 u1,…,uk 都是 A 的特征向量,而且都对应同一个特征值
λ。则每个非零线性组合
α1u1+⋯+αkuk仍然是 A 对应于 λ 的特征向量。
令 w=α1u1+⋯+αkuk。利用线性性,
Aw=j=1∑kαjAuj=j=1∑kαjλuj=λw
若 w=0,它就是特征向量;若各个非零系数的分量恰好抵消而 w=0,
它只属于稍后定义的特征空间,并不是特征向量。因此定理必须写“非零线
性组合”,不能只看系数是否全部为零。
定义
特征空间
设 A 是 n×n 矩阵,而 λ 是一个标量。A 对应于
λ 的 特征空间 定义为
EA(λ)=N(A−λIn).当 λ 是特征值时,这个空间包含所有对应的特征向量,以及零向量;
即使 λ 不是特征值,同一个零空间记号仍给出 EA(λ)={0}。
第一批例子
例题
一个 2×2 矩阵有两个不同特征值
令
A=[13−630−14].先看向量
u1=[5−2].计算得
Au1=[13−630−14][5−2]=[5−2]=1⋅u1.所以 u1 是特征向量,而对应特征值为 1。
再看
u2=[2−1].则
Au2=[−42]=−2[2−1]=−2u2.所以 u2 对应的特征值是 −2。
因此不能把同一特征值的线性组合定理用于 u1+u2,因为这两个向量的
特征值不同。事实上,
u1+u2=[7−3],A(u1+u2)=[10].如果这个和是特征向量,第二个坐标会迫使 λ=0,但第一个坐标又会
要求 1=0⋅7,矛盾。
例题
同一个特征值可以对应多个方向
令
B=211121112.对 x=(x1,x2,x3)T,令 s=x1+x2+x3,则
Bx=(s+x1, s+x2, s+x3)T.向量
u1=111满足 Bu1=4u1,所以 4 是其中一个特征值。
再看
u2=1−10,u3=10−1.二者都满足
Bu2=u2,Bu3=u3.对 λ=1,三个坐标方程都化为 s=0,所以
EB(1)={x∈R3:x1+x2+x3=0}=span⎩⎨⎧1−10,10−1⎭⎬⎫.这两个向量线性无关:若
a(1,−1,0)T+b(1,0,−1)T=0,第二个坐标给出 a=0,第三个坐标给出
b=0。它们构成这个二维特征空间的一组基。
对 λ=4,方程 Bx=4x 给出
s=3x1=3x2=3x3,所以 x1=x2=x3,从而
EB(4)=span⎩⎨⎧111⎭⎬⎫.同一个非零向量不能对应两个不同特征值,但同一个特征值可以有一个平面
的特征向量。平面中的零向量属于特征空间,其非零点才是特征向量。
特征值其实是一个零空间问题
把方程 Av=λv 的所有项移到同一边,便变成
Av−λv=0.
由于 λv=λInv,这等价于
(A−λIn)v=0.
定理
用齐次系统描述特征值
设 A 是 n×n 矩阵,λ 是标量,而 v=0。
下列两件事等价:
- v 是 A 对应于特征值 λ 的特征向量。
- v 是齐次系统
(A−λIn)x=0的一个非平凡解。
正向是把 Av=λv 移项;反向则是把齐次系统的非平凡解移项回去,
得到 Av=λv。因此特征值问题便被改写成普通的线性系统问题。
向量 v 必须落在 A−λIn 的零空间内。
因为它本身是一个零空间,所以 EA(λ) 自动是一个子空间。若
x,y∈EA(λ),则对任意标量 a,b,
(A−λIn)(ax+by)=a(A−λIn)x+b(A−λIn)y=0
所以它对线性组合封闭。特征空间的维数大于零,当且仅当 λ 是
特征值;否则它的维数为零。
dimEA(λ) 称为 λ 的几何重数,表示这个特征值有多
少个线性无关的特征向量方向。在上面的矩阵 B 中,特征值 1 的几何
重数是 2,而特征值 4 的几何重数是 1。这个维数要在找出候选标量
之后,通过零空间方程求得。
概念视角结构
固定特征值的空间不是所有特征向量的总集合
固定一个特征值 λ,方程 (A−λIn)x=0 的解组成一个零空间。
这个解集本来就包含零向量,并且对加法与数乘封闭;它的非零成员才是特征向量。
不同特征值的特征空间取并集,未必对加法封闭:前面的 u1+u2 就不是特征向量。
因此必须区分固定 λ 的特征空间与所有特征值的特征向量集合。
特征值的等价刻画
一旦有了零空间版本,前面学过的可逆性词典便可以直接搬过来。
定理
判断 λ 是否为特征值的等价条件
设 A 是 n×n 矩阵,而 λ 是标量。下列叙述等价:
- λ 是 A 的特征值。
- (A−λIn)x=0 有非平凡解。
- N(A−λIn)={0}。
- A−λIn 不可逆。
- det(A−λIn)=0。
这五个条件来自齐次系统和可逆性词典。非平凡解正是非零零空间向量,所
以 (1)、(2)、(3) 等价;方阵有非零零空间,当且仅当它不可逆,于是
(3) 与 (4) 等价;行列式判准再给出 (4) 与 (5) 等价。特别地,特征值
会迫使 rank(A−λIn)≤n−1。行列式方程用来
找出候选标量,但仍须求零空间,才能找到相应的特征向量。
从方程读出特征空间
找到候选 λ 后,先构造 A−λIn,再把它行化简。所得齐次
系统中的自由变量会参数化 EA(λ) 的全部向量。自由变量的个数
就是这个特征空间的维数,也就是可取得的线性无关特征向量方向数。张成
空间中的每个非零向量都是特征向量;零向量之所以保留,是因为子空间必
须包含零。因此,只有特征值清单而没有特征空间是不完整的。若特征值重
复,行列式记录代数重数,而零空间计算记录几何重数;比较二者是下一节
研究对角化时的关键问题。
例题
用行化简找特征值与特征空间
令
C=[332−2].先解
det(C−λI)=0.得到
det[3−λ32−2−λ]=(3−λ)(−2−λ)−6=λ2−λ−12.因此特征值为
λ=4,λ=−3.对 λ=4,
C−4I=[−132−6]∼[10−20].所以 x1=2x2,可取一组基为
[21].对 λ=−3,
C+3I=[6321]∼[3010].所以 3x1+x2=0,可取一组基为
[1−3].因此
EC(4)=span{[21]},EC(−3)=span{[1−3]}.行列式计算负责找出两个标量;之后的两次行化简分别解决两个齐次系
统。最后代回可核对 C(2,1)T=4(2,1)T 与
C(1,−3)T=−3(1,−3)T,这一步不能由特征多项式单独替代。
对角矩阵与三角矩阵
若 D=diag(d1,…,dn) 是对角矩阵,则特征向量方程
可按坐标分别读取:
(D−λIn)x=0⟺(dj−λ)xj=0对每个 j.
因此
ED(λ)=span{ej:dj=λ},
若没有对角项等于 λ,这个张成空间就是 {0}。重复的对角项
会增加特征空间的维数。对零矩阵,每个非零向量都是特征值 0 的特征
向量;对单位矩阵,每个非零向量都是特征值 1 的特征向量。
若 T 是上三角或下三角矩阵,则 T−λIn 仍是三角矩阵,并且
det(T−λIn)=j=1∏n(tjj−λ).
所以特征值正好是对角项(按代数重数计算)。但对角项不能单独决定特征
向量。考虑上三角矩阵
T=[1042].
有 Te1=e1,所以 ET(1)=span{e1}。当
λ=2 时,(T−2I)x=0 给出 −x1+4x2=0,于是
ET(2)=span{[41]}.
改变右上角的 4 会保留对角项和特征值,却会改变这个特征空间;因此
三角矩阵的对角数据不足以确定全部特征向量。
重要性质
定理
0 是否为特征值,正好等同于可逆与否
对方阵 A,下列两件事等价:
- 0 是 A 的特征值。
- A 不可逆。
等价地,A 可逆当且仅当 0 不是它的特征值。
这其实只是把前面的等价条件代入 λ=0。
具体来说,若 0 是特征值,则存在非零 v 使 Av=0,所以 A 的零空
间非平凡,A 不可逆。反过来,若 A 不可逆,方阵的可逆性词典保证
N(A) 含有非零向量;这个向量满足 Av=0⋅v,所以 0 是特征值。
两个方向都使用同一个齐次系统判准。
定理
简单矩阵运算下的特征值变化
若 λ 是 A 的特征值,则:
- kλ 是 kA 的特征值;
- 对每个正整数 m,λm 是 Am 的特征值;
- λ 也是 AT 的特征值;
- 若 A 可逆,则 λ−1 是 A−1 的特征值。
若 Av=λv,则
(kA)v=kλv,
这证明了 (1),也包括 k=0:零矩阵的每个非零向量都是特征值 0 的特
征向量。对于正整数 m,用归纳法得到 Amv=λmv;归纳步骤
把 A 乘到等式两边,并使用 Av=λv。m=0 是独立情形:
A0=I,同一个非零向量对单位矩阵的特征值是 1,这里不把它写成
00。
对于 (4),若 A 可逆,则 λ=0。在 Av=λv 两边乘
A−1,得到 v=λA−1v,再除以 λ 即得
A−1v=λ−1v。对于 (3),使用行列式恒等式
det(AT−λIn)=det((A−λIn)T)=det(A−λIn)=0.
这通过特征值判准证明了转置结论;不能把 Av=λv 直接转置为
ATv=λv。对矩阵 T=[1042],
Te1=e1,但
TTe1=[1402][10]=[14]=e1.
可选复数扩展:数域很重要
上面的基本定义是在 R 上进行的。考虑实矩阵
R=[11−11]
若 λ 是实数,则
det(R−λI2)=(1−λ)2+1=(λ−1)2+1>0.
因此 R 没有实特征值,也没有实特征向量。若把标量域明确改为
C,根变成 1+i 与 1−i。取
u+=[i1],u−=[−i1],
直接计算可得
Ru+=[i−1i+1]=(1+i)u+,Ru−=[−1−i1−i]=(1−i)u−.
解 (R−(1+i)I)x=0 的第一行得到 x=iy;对 1−i 则得到 x=−iy。
所以在 C 上,两个特征空间分别是
spanC{u+} 与
spanC{u−}。这个例子标出数域的边界:非实
根不是实特征值,但在明确更换标量域后便成为有效的特征值。
常见错误
常见错误
零向量永远不是特征向量
学生常会留意到 A0=λ0 对所有标量都成立,于是误以为零向量是所
有特征值的特征向量。正因如此,定义才必须排除零向量。特征向量代表的
是某个真正的方向,而零向量没有方向可言。
快速检查
思考检查
若 v 是 A 对应于特征值 λ 的特征向量,则 3v 也是吗?
解答 · 答案
是。因为 A(3v)=3Av=3λv=λ(3v),所以 3v 仍是对应于
λ 的特征向量。
思考检查
为什么 0 是特征值会迫使 A 不可逆?
利用系统 (A−0I)x=0 去理解。
解答 · 答案
若 0 是特征值,则 Ax=0 有非平凡解。这表示 A 的零空间非平
凡,因此 A 不可逆。
思考检查
用零空间语言写出 EA(λ)。
解答 · 答案
EA(λ)=N(A−λIn).它就是 A−λIn 的零空间。
练习
思考检查
证明 [11] 是 [2112] 的特征向量,并找出其特征值。
解答 · 引导解答
[2112][11]=[33]=3[11].所以这个向量是特征向量,而对应特征值是 3。
思考检查
求 A=[1004] 对应于特征值 4 的特征空间。
解 (A−4I)x=0。
解答 · 引导解答
A−4I=[−3000].因此 x1=0,x2 自由,所以
EA(4)=span{[01]}.
思考检查
若 A 可逆,则 0 会否成为 AT 的特征值?
解答 · 引导解答
不会。若 0 是 AT 的特征值,则 AT 不可逆。但 A 可逆当且仅
当 AT 可逆,所以这是不可能的。
相关笔记
请把
7.2 行变换、乘积与可逆性
放在旁边,因为测试 det(A−λI)=0 直接依赖行列式与可逆性的理
解。
接着可读
8.2 对角化与相似,
看看一整组特征向量如何改变矩阵的整体形状。
本节的子空间观点也依赖
6.5 基底与维数。