小知识,大挑战!本文正在参与“程序员必备小知识”创作活动。
前言
Hello!小伙伴!
非常感谢您阅读海轰的文章,倘若文中有错误的地方,欢迎您指出~
自我介绍 ଘ(੭ˊᵕˋ)੭
昵称:海轰
标签:程序猿|C++选手|学生
简介:因C语言结识编程,随后转入计算机专业,有幸拿过一些国奖、省奖...已保研。目前正在学习C++/Linux/Python
学习经验:扎实基础 + 多做笔记 + 多敲代码 + 多思考 + 学好英语!
机器学习小白阶段
文章仅作为自己的学习笔记 用于知识体系建立以及复习
知其然 知其所以然!
5.2 方阵的特征值与特征向量
定义6
设A是n阶矩阵,如果数λ和n维非零列向量x使关系式
Ax=λx(1)
成立,那么数λ称为矩阵A的特征值,非零向量x称为A的对应特征值λ的特征向量
(1)式也可以写为
(A−λE)x=0(2)
(2)式有非零解说明∣A−λE∣=0
方程Ax=0有非零解时,说明R(A)<n,即∣A∣=0
即
⎣⎡a11−λa21..an1a12a22−λ..an2.........a1na2n..ann−λ⎦⎤=0(3)
(3)式是以λ为未知数的一元n次方程,称为矩阵A的特征方程;∣A−λE∣是λ的n次多项式,记作f(λ),称为矩阵A的特征多项式
- A的特征值就是特征方程的解
- 特征方程在复数范围内一定有解,其解的个数为方程的次数
- n阶矩阵A在复数范围内有n个特征值
性质1:设n阶矩阵A=(aij)的特征值为λ1,λ2,...,λn,则有
- λ1+λ2+...+λn=a11+a22+...+ann
- λ1λ2....λn=∣A∣
性质2:设λ=λi为矩阵A的一个特征值,则由方程
(A−λiE)x=0
求得非零解x=pi,那么pi就是A的对应于特征值λi的特征向量
显然pi是矩阵A对应于特征值λi的特征向量
那么kpi也是矩阵A对应于特征值λi的特征向量
定理2
设λ1,λ2,...,λm是方阵的m个特征值,p1,p2,...,pm依次是与之对应的特征向量,如果λ1,λ2,...,λm各不相等,则p1,p2,...,pm线性无关
证明:使用数学归纳法
当m=1时,特征向量为p1,则只含有一个向量p1的向量组一定是线性无关
p1是特征向量,那么p1肯定是非空向量(由定义而来)
当m>1时
假设m=k−1时,结论成立,即
λ1,λ2,...,λk−1是方阵的(k−1)个特征值,那么对应特征向量p1,p2,...,pk−1线性无关
那么只需要证明m=k时结论同样成立即可
令
x1p1+x2p2+....+xk−1pk−1+xkpk=0(1)
xi(i∈[1,k])表示一个常数
用A左乘(1)式,得
x1Ap1+x2Ap2+....+xk−1Apk−1+xkApk=0(2)
因为Api=λipi,得
x1λ1p1+x2λ2p2+....+xk−1λk−1pk−1+xkλkpk=0(3)
式(3)−λk(2),得
(x1λ1p1+x2λ2p2+....+xk−1λk−1pk−1+xkλkpk)−λk(x1p1+x2p2+....+xk−1pk−1+xkpk)=0
化简,得
x1(λ1−λk)p1+x2(λ2−λk)p2+....+xk−1(λk−1−λk)pk−1=0
根据假设p1,p2,...,pk−1线性无关
所以
x1(λ1−λk)=x2(λ2−λk)=...=xk−1(λk−1−λk)=0
即
xi(λi−λk)=0(i∈[1,k−1])
因为λ1,λ2,...,λk各不相等,得
(λi−λk)=0
所以
xi=0
再将xi=0(i∈[1,k−1])代入(1)式,得到
xkpk=0
因为pk=0,故
xk=0
从而
x1=x2=....=xk−1=xk=0
即p1,p2,...,pk线性无关
从而当m=k时,结论同样成立
证明完成!
举例
例5
求矩阵A=[3−1−13]的特征值和特征向量
解答:
A的特征多项式为
∣A−λE∣=A=[3−λ−1−13−λ]=(3−λ)2−1=8−6λ+λ2=(4−λ)(2−λ)
解得
λ1=2,λ2=4
当λ1=2时,对应的特征向量满足
(A−λ1E)x=0(2)
即
A=[3−2−1−13−2]x=[1−1−11][x1x2]=0
解得
x1=x2
所以λ1=2对应的特征向量可以取p1=[11]
同理
当λ2=4时,解得
x1=−x2
所以λ2=4对应的特征向量可以取p2=[−11]
综上,pi是对应于λi的特征向量,那么kpi(k=0)同样也是对应于λi的特征向量
例6
求矩阵A=⎣⎡−1−41130002⎦⎤的特征值和特征向量
解答
A的特征多项式为
∣A−λE∣=∣∣−1−λ−4113−λ0002−λ∣∣=(2−λ)(1−λ2)
(2−λ)(1−λ2)=0,解得
λ1=2,λ2=λ3=1
当λ1=2时,解方程(A−2E)x=0
A−2E=⎣⎡−3−41110000⎦⎤∼r⎣⎡100010000⎦⎤
得到基础解系
p1=⎣⎡001⎦⎤
kp1(k=0)是对应于λ1=2的全部特征向量
当λ2=λ3=1,解方程(A−E)x=0
A−E=⎣⎡−2−41120001⎦⎤∼r⎣⎡100010120⎦⎤
得基础解系
p2=⎣⎡−1−21⎦⎤
kp2(k=0)是对应于λ2=λ3=1的全部特征向量
例8
设λ是方阵A的特征值,证明
(1)λ2是A2的特征值
(2)当A可逆时,λ1是A−1的特征值
证明(1):
因为λ是方阵A的特征值
则存在非零向量x,使得
Ax=λx
A2x=A(Ax)=A(λx)=λAx=λ(λx)=λ2x
即
A2x=λ2x
所以λ2是A2的特征值
证明(2):
因为λ是方阵A的特征值
则存在非零向量x,使得
Ax=λx
当A可逆时,则A−1存在
等式两边同时左乘A−1
A−1Ax=A−1λx=λA−1x
即
x=λA−1x
因为x=0,所以λ=0
等式两边再同时除以λ,得
λ1x=A−1x
综上:λ1是A−1的特征值
例9
设3阶矩阵A的特征值为1,−1,2,求A∗+3A−2E的特征值
解答:
A的特征值全不为0,得A可逆
若至少有一个特征值为0
那么∣A∣=λ1λ2...λn=0
说明R(A)<n
推出A不可逆
相反,特征值全不为0
得∣A∣=λ1λ2...λn>0、R(A)=n
知A可逆
所以
A∗=∣A∣A−1
又因为
∣A∣=λ1λ2λ3=−2
得到
A∗+3A−2E=−2A−1+3A−2E
令
F(A)=A∗+3A−2E=−2A−1+3A−2E
假设A的特征值是λ
那么F(A)的特征值就是F(λ)=−2λ1+3λ−2
结合上述,得F(A)的特征值就是
- F(1)=−1
- F(−1)=−3
- F(2)=3
如果A的特征值是λ1,λ2,...,λm
那么F(A)的特征值是F(λ1),F(λ2),...,F(λm)
例10
设λ1和λ2是矩阵A的两个不同的特征值,对应的特征向量依次为p1和p2,证明p1+p2不是A的特征向量
证明:
因为λ1和λ2是矩阵A的两个不同的特征值
所以
Ap1=λ1p1
Ap2=λ2p2
使用反证法 假设p1+p2是是A的特征向量
那么就存在一个λ,使得
A(p1+p2)=λ(p1+p2)
又
A(p1+p2)=Ap1+Ap2=λ1p1+λ2p2
有
λ1p1+λ2p2=λ(p1+p2)
移项得
(λ1−λ)p1+(λ2−λ)p2=0
由定理2可知,p1,p2线性无关
则有
λ1−λ=0
λ2−λ=0
即
λ1=λ=λ2
但与题意λ1=λ2相矛盾
故假设不成立
即p1+p2不是A的特征向量
结语
说明:
- 参考于 课本《线性代数》第五版 同济大学数学系编
- 配合书中概念讲解 结合了自己的一些理解及思考
文章仅作为学习笔记,记录从0到1的一个过程
希望对您有所帮助,如有错误欢迎小伙伴指正~
我是 海轰ଘ(੭ˊᵕˋ)੭
如果您觉得写得可以的话,请点个赞吧
谢谢支持 ❤️