本文已参与「新人创作礼」活动,一起开启掘金创作之路。
特征值、特征向量
定义:设A是n阶矩阵,α是n维非0列向量,且
Aα=λα
则称λ是矩阵A的特征值,α是矩阵A对应于特征值λ的特征向量
由Aα=λα,α=0⇒(λE−A)α=0,(λE−A)x=0⇒α是齐次方程组(λE−A)x=0的非0解
-
由∣λE−A∣=0求特征值λi,共n个(含重根)
-
对(λiE−A)x=0,求基础解系,即特征值λi的线性无关的特征向量,写通解得λi所有的特征向量
定理:如果α1,α2都是矩阵A对应于特征值λ的特征向量,则当k1α1+k2α2=0时,k1α1+k2α2仍是矩阵A关于特征值λ的特征向量
定理:如果λ1与λ2是A不同的特征值,对应的特征向量分别是α1和α2,则α1,α2必定线性无关
定理:设A是n阶矩阵,特征值是λ1,λ2,⋯,λn,则有
例:求A=⎝⎛17−2−2−214−4−2−414⎠⎞特征值,特征向量
三阶行列式主对角线元素都含有未知数,直接展开要解三次方程,一般考虑通过行列加减,找出某一行或某一列所有的元素都含有含未知数的公因式或0
由A的特征多项式
∣λE−A∣=∣∣λ−17222λ−14424λ−14∣∣=∣∣λ−17202λ−1418−λ24λ−18∣∣=∣∣λ−17204λ−10024λ−18∣∣=(λ−18)∣∣λ−1724λ−10∣∣=(λ−18)(λ2−27λ+162)=(λ−18)2(λ−9)
因此λ1=λ2=18,λ3=9
当λ=18时,(18E−A)x=0
⎝⎛122244244⎠⎞→⎝⎛100200200⎠⎞
得基础解系:α1=(−2,1,0)T,α2=(−2,0,1)T,因此特征向量k1α1+k2α2,k1,k2不同时为0
当λ=9时,(9E−A)x=0
⎝⎛−8222−5424−5⎠⎞→⎝⎛200010−1−10⎠⎞
得基础解系α3=(1,2,2)T,因此特征向量k3α3,k3=0
对于一个矩阵若,秩等于1即存在一行能表示其他所有行,秩等于2即存在两行能表示其他所有行,之后同理,类似最大线性无关组
⎝⎛−8222−5424−5⎠⎞求基础解系,本题题解省略的化简步骤,其实化简步骤并不容易,但此处可以用其他思路。
由于已知∣λE−A∣=0,即该矩阵的秩一定小于3;随便选两行,这里选二三行,发现二者线性无关(不成比例),因此该矩阵秩大于1,可得该矩阵秩为2。
由于这两行能线性表示另一行,因此可以构造新的矩阵
⎝⎛220−5404−50⎠⎞
可以理解为由于这两行能线性表示另一行,因此另一行一定能被全消为0
此时只需对新的矩阵行化简即可
设Aα=λα,α=0,则有
例:A为3阶矩阵,特征值是−1,0,4,如果A+B=2E,则B的特征值为()
设Aα=λα,α=0
由A+B=2E,有B=2E−A,则
Bα=(2E−A)α=2α−Aα=(2−λ)α
因此B的特征值为3,2,−2
例:A是3阶矩阵,A2+2A−3E=0,证明矩阵A的特征值只能是1或−3
设λ是A的任一特征值,对应的特征向量是α,即Aα=λα,α=0,那么
A2α=λ2α
由A2+2A−3E=0
有
A2α+2Aα−3α(λ2+2λ−3)αλ2+2λ−3=0=0,α=0=0
因此λ=1或λ=−3
从题目条件只能推出矩阵特征值只能是1或−3
例如:⎝⎛1 1 1⎠⎞,⎝⎛3 3 3⎠⎞,⎝⎛1000−1−20−2−1⎠⎞
例:已知α=(1,1,−1)T是A=⎝⎛25−1−1ab23−2⎠⎞的一个特征向量,则a=(),b=()
设Aα=λα,α=0
⎝⎛25−1−1ab23−2⎠⎞⎝⎛11−1⎠⎞=λ⎝⎛11−1⎠⎞
有
⎩⎨⎧2−1−2=λ5+a−3=λ−1+b+2=−λ
可得λ=−1,a=−3,b=0
相似矩阵
设A,B都是n阶矩阵,如果存在可逆矩阵P使
P−1AP=B
就称矩阵A相似与矩阵B,B是A的相似矩阵,记作A∼B
相似的基本性质
-
A∼A
-
如果A∼B,则B∼A
-
如果A∼B,B∼C,则A∼C
如果A∼B
-
∣λE−A∣=∣λE−B∣⇒λA=λB
证明:
∣λE−B∣=∣λE−P−1AP∣=∣P−1(λE−A)P∣=∣P−1∣∣λE−A∣∣P∣=∣λE−A∣
-
r(A)=r(B)
证明:
如果A可逆,有r(AB)=r(B),r(BA)=r(B),则
r(B)=r(P−1AP)=r(AP)=r(A)
-
∣A∣=∣B∣
∣B∣=∣P−1AP∣=∣P−1∣∣A∣∣P∣=∣A∣
-
∑aii=∑bii
-
An∼Bn
A+kE∼B+kE
A−1∼B−1
相似对角化
如果A∼Λ,则称矩阵A可相似对角化
定理:A∼Λ⇔A有n个线性无关的特征向量
证明:
如果Aα1=λ1α1,Aα2=λ2α2,Aα3=λ3α3,且α1,α2,α3线性无关,则
A(α1,α2,α3)=(Aα1,Aα2,Aα3)=(λ1α1,λ2α2,λ3α3)=(α1,α2,α3)⎝⎛λ1 λ2 λ3⎠⎞
令P=(α1,α2,α3),由于α1,α2,α3线性无关,则P可逆,有
APP−1AP=PΛ=Λ=⎝⎛λ1 λ2 λ3⎠⎞
必要性得证
如果P−1AP=Λ,有
APA(α1,α2,α3)(Aα1,Aα2,Aα3)=PΛ=(α1,α2,α3)⎝⎛a1 a2 a3⎠⎞=(a1α1,a2α2,a3α3)=(a1α1,a2α2,a3α3)
因此有
Aα1=a1α1,Aα2=a2α2,Aα3=a3α3
又因为P=(α1,α2,α3)可逆,则α1,α2,α3线性无关,充分性得证、
推论:如果A有n个不同的特征值,则A∼Λ
定理:A∼Λ⇔λ是A的k重特征值,则λ有k个线性无关的特征向量
相似对角化解题步骤:
-
求特征值λ1,λ2,λ3
-
求特征向量α1,α2,α3
-
构造可逆矩阵P=(α1,α2,α3),则P−1AP=⎝⎛λ1 λ2 λ3⎠⎞
例:已知A=⎝⎛0 01010100⎠⎞,求可逆矩阵P,使P−1AP=Λ
由特征多项式
∣λE−A∣=∣∣λ0−10λ−10−10λ∣∣=(λ−1)2(λ+1)
则A的特征值1,1,−1
当λ=1时,由(E−A)x=0
⎝⎛10−1000−101⎠⎞→⎝⎛100000−100⎠⎞
特征向量α1=(0,1,0)T,α2=(1,0,1)T
当λ=−1时,由(−E−A)x=0
⎝⎛−10−10−20−10−1⎠⎞→⎝⎛100010100⎠⎞
特征向量α3=(−1,0,1)T
令
P=(α1,α2,α3)=⎝⎛010101−101⎠⎞
有
P−1AP=⎝⎛1 1 −1⎠⎞
例:已知A=⎝⎛20000101x⎠⎞,B=⎝⎛2000y000−1⎠⎞相似
由于A∼B,有
∑aii=∑bii∣A∣=∣B∣∣λE−A∣=∣λE−B∣⇒2+0+x=2+y+(−1)(1)⇒−2=−2y(2)⇒∣∣λ−2000λ−10−1λ−x∣∣=∣∣λ−2000λ−y000λ+1∣∣⇒(λ−2)(λ2−xλ−1)=(λ−2)[λ2+(1−y)λ−y]⇒{−x=1−y−1=−y(3)
可以选(1),(2)或单独选(3),一般选(1),(2),计算较简单
或
B的特征值2,y,−1,则λ=−1是A的特征值,有
∣−E−A∣=∣∣−3000−1−10−1−1−x∣∣=−3∣∣−1−1−1−1−x∣∣=−3x=0
若用λ=2,有
∣2E−A∣=∣∣00002−10−12−x∣∣=0
若用λ=y,有
∣yE−A∣=∣∣y−2000y−10−1y−x∣∣=(y−2)(y2−xy−1)=0
该式需要结合上面的(1),(2),(3)使用,或者选择两个不相关的使用,例如λ=0,λ=y搭配
该方法当出现本题λ=2时∣2E−A∣本身有一行就为0或两行成比例,则需要换λ或用(1),(2),(3)方法
解得x=0,y=1,代入题设有
A=⎝⎛200001010⎠⎞∼⎝⎛2 1 −1⎠⎞=B
此处省略步骤,可得
P=⎝⎛10001101−1⎠⎞
实对称矩阵
向量的内积
定义:设α=(a1,a2,⋯,an)T,β=(b1,b2,⋯,bn)T,有
(α,β)=a1b1+a2b2+⋯+anbn=αTβ=βTα
如果(α,β)=0,称α与β正交
(α,α)=a12+a22+⋯+an2,称a12+a22+⋯+an2为向量α的长度,记作∣∣α∣∣
内积的性质
-
(α,β)=(β,α)
-
(kα,β)=(α,kβ)=k(λ,β)
-
(α+β,γ)=(α,γ)+(β,γ)
-
(α,α)≥0,等号当且仅当α=0时成立
例:求与α1=(−1,−1,1)T,α2=(1,2,1)T都正交的向量
设α=(x1,x2,x3)T与α1,α2都正交,则αTα1=0,αTα2=0
{−x1−x2+x3=0x1−2x2−x3=0
有
(−11−1−21−1)→(1001−10)
基础解系(1,0,1)T,因此α=k(1,0,1)T
正交矩阵
定义:A是n阶矩阵,若AAT=ATA=E,则称A是正交矩阵
如果A是正交矩阵
⇔AT=A−1
⇔A的列向量都是单位向量且两两正交
如果A是正交矩阵⇒∣A∣为1或−1
证明:
AAT=E⇒∣AAT∣∣A∣⋅∣AT∣∣A∣2=∣E∣=1=1
因此∣A∣为1或−1
实对称矩阵
定理:实对称矩阵必可相似对角化
定理:实对称矩阵不同特征值所对应的特征向量相互正交
定理:实对称矩阵必存在正交矩阵Q,使Q−1AQ=QTAQ=Λ
例:A为3阶实对称矩阵,满足A2=A,如果r(A−E)=2,则A的特征值()
设Aα=λα,α=0,则
A2α=λ2α
由A2=A,有
λ2α(λ2−λ)α=λα=0⇒λ2−λ=0
因此λ是1或0,又由于A是实对称矩阵,r(A−E)=2,有
A∼Λ⇒A−E∼Λ−E⇒r(Λ−E)=2
Λ与A相似,即特征值相同,又有r(Λ−E)=2,则可构造
Λ=⎝⎛1 0 0⎠⎞
因此A的特征值为1,0,0
求Q−1AQ=QTAQ=Λ
-
求出A的特征值λ1,λ2,λ3
-
求出对应的特征向量α1,α2,α3
-
该特征向量为γ1,γ2,γ3
1. 如果特征值不同,只需单位化(实对称矩阵不同特征值所对应的特征向量相互正交)
2. 若特征值有重根
1. 如果特征向量已正交,只需单位化
2. 如果特征向量不正交,需施密特(Schmidt)正交化
- 构造正交矩阵Q=(γ1,γ2,γ3),Q−1AQ=Λ=⎝⎛λ1 λ2 λ3⎠⎞
例:已知A=⎝⎛1a−1a31−111⎠⎞,若r(A)=2,求a,求正交矩阵Q和对角矩阵Λ,使Q−1AQ=Λ
求a可以用行变换,这里用秩的定义,即行列式不为0。由于A中有二阶子式
∣∣3111∣∣=0
于是r(A)=0⇔∣A∣=0
∣∣1a−1a31−111∣∣=∣∣0a+10a31−111∣∣=−(a+1)2
因此a=−1,有
A=⎝⎛1−1−1−131−111⎠⎞
由A的特征多项式
∣λE−A∣=∣∣λ−1111λ−3−11−1λ−1∣∣找除主对角线外对应位置两数字和为0=∣∣λ110λ−3−1λ−1λ−1∣∣=λ(λ−1)(λ−4)
因此A的特征值为0,1,4
对λ=1
(E−A)=⎝⎛0111−2−11−10⎠⎞→⎝⎛100010110⎠⎞⇒α1=(−1,−1,1)T
对λ=4
(4E−λ)=⎝⎛31111−11−13⎠⎞→⎝⎛1000101−20⎠⎞⇒α2=(−1,2,1)T
对λ=0
(0E−A)=⎝⎛−1111−3−11−1−1⎠⎞→⎝⎛100010−100⎠⎞⇒α3=(1,0,1)T
实对称矩阵,特征值不同,特征向量已经正交,只需单位化
γ1=31⎝⎛−1−11⎠⎞,γ2=61⎝⎛−121⎠⎞,γ3=21⎝⎛101⎠⎞
令
Q=(γ1,γ2,γ3)=⎝⎛ −31−3131−61626121021⎠⎞
则
Q−1AQ=Λ=⎝⎛1 4 0⎠⎞