对角化
相似矩阵
相似:A、B是n阶方阵,存在n阶可逆矩阵P,如果:P−1AP=B,则:A、B相似(A ~ B)
- 反身性,A ~ A,E−1AE=B
- 对称性,A ~ B -> B ~ A,同时左乘P以及右乘P−1,即:PP−1APP−1=PBP−1,得PBP−1=A,转换下:(P−1)−1B(P−1)=A
- 传递性,A ~ B 并且 B ~ C -> A ~ C,P−1AP=B,Q−1BP=C,代入B:Q−1P−1APQ=C,(PQ)−1A(PQ)=C
性质1:如果A ~ B,那么A、B有相同的特征值并且|A|=|B|(对角线相乘)、tr(A)=tr(B)
证明:∣λE−B∣=∣λE−P−1AP∣=∣λP−1EP−P−1AP∣=∣p−1∣×∣λ−A∣×∣P∣=∣λ−A∣
性质2:如果A ~ B,那么A可逆的充要条件是B可逆,并且如果A、B都可逆,那么A−1 ~ B−1
证明:根据性质1可知,如果A ~ B,那么|A|=|B|,①|A|、|B|均不等于0,则A、B可逆,②|A|、|B|均等于0,则A、B不可逆,即:A、B同时可逆或同时不可逆
B−1=(P−1AP)−1=P−1A−1P,所以A−1 ~ B−1
性质3:如果A ~ B,那么Am ~ Bm
证明:Bm=(P−1AP)m=P−1AP×P−1AP×P−1AP...=P−1AmP
A ~ B 性质总结:
- tr(A) = tr(B)
- |A| = |B|
- 均可逆(或均不可逆)
- A−1=B−1
- Am=Bm
- 特征值相同
- r(A) = r(B),矩阵左乘可逆矩阵或右乘可逆矩阵或左右都乘可逆矩阵,矩阵的秩不变
计算题(相似求未知量)
A=⎝⎛20001201a⎠⎞B=⎝⎛10003−102b⎠⎞ A ~ B,求:a、b的值
用tr(A) = tr(B) 和 |A| = |B|,解方程即可
对角形
P−1AP = ∧(对角形)
与对角形相似的条件
定理:A ~ ∧ <=> A有n个线性无关的特征向量 =>
证明必要性:P−1AP = ∧ = ⎝⎛λ1⋮00⋱00⋮λn⎠⎞P=(α1,...,αn)
P−1AP=∧,同时左乘P,得:AP=P∧,代入向量组,A(α1,...,αn)=(α1,...,αn)⎝⎛λ1⋮00⋱00⋮λn⎠⎞ = Aα1,...,Aαn=λ1α1,...,λnαn(根据分块矩阵定义)
所以,Aα1=λ1α1,...,Aαn=λnαn
因为P可逆,所以α1,...,αn!=0
所以,满足特征向量的性质:Aα=λα且α != 0
因为,r(P)=n,所以r(α1,...,αn)=n(矩阵的秩等于矩阵的列秩),所以α1,...,αn线性无关
证明充分性:把必要性的步骤反过来证明即可
推论:如果A有n个互异的特征根(充分条件),那么A ~ ∧
证明:λ1,...,λn是互异的话,那么每一个特征值找一个特征向量出来(α1,...,αn),它们必定是线性无关的(条件强于上面的定理)
- 如果矩阵所有的根都是单根(必定线性无关),一定相似于对角形
- 如果矩阵有重根,需要找到n(n等于重根的个数)个线性无关的特征向量,否则不能不相似于对角形
定理:A ~ ∧ <=> 每个ri重的特征根,其基础解系有ri个解
计算题(对角形相似)
A=⎝⎛3−232−26−12−1⎠⎞,是否与对角形相似?如果相似,P=?,∧=?
∣λE−A∣=∣∣∣∣∣∣λ−32−3−2λ+2−61−2λ+1∣∣∣∣∣∣,第3列加到第1列,∣∣∣∣∣∣λ−20λ−2−2λ+2−61−2λ+1∣∣∣∣∣∣=(λ−2)2(λ+4)=0
λ1=−4λ2=λ3=2
当λ1=−4时,对(λ1E−A)x=0,⎝⎛−72−3−2−2−61−2−3⎠⎞−>⎝⎛100010−31320⎠⎞
{x1=31x3x2=−32x3
令x3=1,得:⎝⎛31−321⎠⎞=α1
当λ2=λ3=4时,得:⎝⎛−210⎠⎞=α2⎝⎛101⎠⎞=α3
P=(α1,α2,α3)=⎝⎛31−321−210101⎠⎞
∧ = ⎝⎛−400020002⎠⎞(把特征值依次带入即可)
注意:特征值和特征向量的顺序必须是对应的
计算题(相似反推矩阵)
3阶方阵A,λ1=1,λ2=2,λ3=3,α1=(1,1,1)T,α2=(1,2,3)T,α3=(1,3,6)T,求:A、AT的特征值和特征向量
P=(α1,α2,α3)=⎝⎛111123136⎠⎞P−1AP=⎝⎛100020003⎠⎞=∧
同时左乘P,右乘P−1,得:A=P∧P−1,得:A=⎝⎛0031−1−9027⎠⎞
AT的特征值不变,所以还是λ1=1,λ2=2,λ3=3
P−1AP=∧,(P−1AP)T=PTAT(p−1)T=∧(对角形的转置等于它本身)
逆矩阵的转置等于转置的逆矩阵,所以PTAT(PT)−1=∧
整理为定理:[(PT)−1]−1AT(PT)−1=∧,求得:(PT)−1=⎝⎛3−31−35−21−21⎠⎞
(PT)−13个列向量就是AT的特征向量
计算题(利用对角形求高次幂)
A=⎝⎛0−2−1101−120⎠⎞,求:A100
先求出特征值,λ1=−1,λ2=0,λ3=1
再求出特征向量并构成矩阵P,P=⎝⎛101111143⎠⎞
P−1AP=∧,同时左乘P和右乘P−1,得:A=P∧P−1
P∧P−1P∧P−1...P∧P−1 = P∧的100次方P−1
计算出结果...
实对称矩阵的对角化(预热)
所有的实对称矩阵都能对角化(也就是说:所有的实对称矩阵都能相似于一个对角形)
内积:两个向量对应的分量相乘后累加,得到的值就是内积(内积是个数)
α=⎝⎛a1a2a3⎠⎞β=⎝⎛b1b2b3⎠⎞
(α,β) = a1b1+a2b2+a3b3=num
我们习惯上用"行"乘"列"的方式来表示内积。如:α、β都是列向量,那么可以表示为:αTβ
- (α,α) >= 0,因为:a12+a22+a32肯定大于等于0,所以:(α,α) = 0 <=> α = 0
- 对称性:(α,β) = (β,α)
- 齐次性:(kα,β)=k(α,β),(α,kβ)=k(α,β),(kα,kβ)=k2(α,β)
- (α+β,γ)=(α,γ)+(β,γ), (α,β+γ)=(α,β)+(α,γ)
向量的长度(或者说:范数、模),∣∣α∣∣=√(α,α),两边同时平方:(α,α)=∣∣α∣∣2
α=(−1,−1,5),那么∣∣α∣∣=√(−1)2+(−1)2+(5)2=√27,点到原点的距离
如果||α||=1,则α是单位向量
单位化(标准化),∣∣α∣∣1α=β,乘以倒数后的向量的长度为1
证明:∣∣β∣∣=√(∣∣α∣∣1α,∣∣α∣∣1α)=√∣∣α∣∣21(α,α)=∣∣α∣∣1√(α,α)=∣∣α∣∣1∣∣α∣∣=1
性质1:∣∣α∣∣>=0,||α||=0 <=> α=0
=>
性质2:||kα||=|k| * ||α||(kα的长度=k的绝对值乘以α的长度)
证明:∣∣kα∣∣=√(kα,kα)=√k2(α,α)=∣k∣×∣∣α∣∣
柯西施瓦茨不等式,|(α,β)| <= ||α|| * ||β||
∣a1b1+a2b2+...+anbn∣<=√a12+a22+...+an2×√b12+b22+...+bn2
三角不等式,||α+β|| <= ||α|| + ||β||,两边之和大于第三边,在向量中如果同方向就可以等于
正交
正交(垂直),(α,β) = 0,那么 α ⊥ β,0向量和任何向量都正交,和自身正交的向量只能是0向量
正交向量组:给一组向量(不含0向量),两两都正交,α1,α2,...,αs中随便取2个向量都是正交的
标准正交向量组,每一个向量都是单位向量,{(αi,αi)=1(αi,αj)=0
定理:如果α1,α2,...,αs是正交向量组,那么α1,α2,...,αs一定线性无关(线性相关须在同一直线上)
证明: 设:k1α1+k2α2+...+ksαs=0
(α1,k1α1+k2α2+...+ksαs)=(α1,0)=0
k1(α1,α1)+k2(α1,α2)+...+ks(α1,αs)=0
根据正交的性质,后面的都等于0,只剩:k1(α1,α1)=0,因为正交向量组不含0向量,所以(α1,α1)!=0,所以:k1=0
用同样的方式,求出k2,...,ks=0,所以线性无关
施密特正交化,给一组线性无关的α1,α2,...,αs,求一组与之等价的并且是正交的β1,β2,...,βs
β1=α1
β2=α2−(β1,β1)(α2,β1)×β1
β3=α3−(β1,β1)(α3,β1)×β1−(β2,β2)(α3,β2)×β2
...
计算题(施密特正交化)
线性无关的向量,α1=(1,1,1,1)Tα2=(1,−2,−3,−4)Tα3=(1,2,2,3)T
一般流程:α1,α2,...αs 正交化-> β1,β2,...βs 单位化-> η1,η2,...ηs
β1=α1
β2=α2−(β1,β1)(α2,β1)β1
β3=α3−(β1,β1)(α3,β1)β1−(β2,β2)(α3,β2)β2
η1=∣∣β1∣∣1β1η2=∣∣β2∣∣1β2η3=∣∣β3∣∣1β3
正交矩阵
正交矩阵:A是n阶方阵,ATA=E
性质1:如果A是正交矩阵,那么|A| = 1或-1
证明:∣AT∣×∣A∣=1,则:∣A∣2=1
性质2:如果A是正交矩阵,那么A−1=AT且A−1和AT均为正交矩阵
证明:因为ATA=E,所以A−1=AT
(A−1)TA−1=(A−1)TAT=(AA−1)T=E
(AT)TAT=AAT=AA−1=E
性质3:A、B正交,那么AB也正交
证明:(AB)TAB=BTATAB=B−1A−1AB=E
性质4:如果A是正交矩阵,α、β是n维列向量,那么(Aα,Aβ)=(α,β)
证明:(Aα)TAβ=αTATAβ=αTβ=(α,β)
定理:A是正交矩阵 <=> A的列(行)向量组是标准正交向量组 =>
证明:设A的列向量为(α1,α2,...,αn)
向量可以看作是一个分块矩阵,分块矩阵转置:先把子块元素求转置,再对每个子块求转置
ATA=(α1,α2,...,αn)T(α1,α2,...,αn)
=⎝⎜⎜⎛α1Tα2T⋮αnT⎠⎟⎟⎞(α1,α2,...,αn)
=⎝⎜⎜⎛α1Tα1α2Tα1⋮αnTα1α1Tα2α2Tα2⋮αnTα2⋯⋯⋱⋯α1Tαnα2Tαn⋮αnTαn⎠⎟⎟⎞
=⎝⎜⎜⎛(α1α1)(α2α1)⋮(αnα1)(α1α2)(α2α2)⋮(αnα2)⋯⋯⋱⋯(α1αn)(α2αn)⋮(αnαn)⎠⎟⎟⎞
因为是正交矩阵,所以最终结果等于单位阵E,即:⎝⎜⎜⎛10⋮001⋮0⋯⋯⋱⋯00⋮1⎠⎟⎟⎞
(αi,αi)=∣∣α∣∣2=1,(ai,aj)=0
所以,满足必要性(标准正交向量组)
充分性的证明逻辑跟必要性类似
计算题(正交矩阵)
三阶非0矩阵A,aij=Aij,求:①|A|=1,②A是正交矩阵
证明:A∗=⎝⎛A11A12A13A21A22A23A31A32A33⎠⎞=⎝⎛a11a12a13a21a22a23a31a32a33⎠⎞=AT
ATA=A∗A=∣A∣E
两边同时取行列式
∣AT∣×∣A∣=∣∣A∣∣
- ∣AT∣=∣A∣
- 行列式提取k要提n次
所以,∣A∣2=∣A∣3
以上式子可以写成:∣A∣2(∣A∣−1)=0,有2个解:1和0
因为A是三阶非0矩阵,所以必定有元素!=0,假设:a11!=0
那么,按第一行展开后,得:∣A∣=a11A11+a12A12+a13A13=a112+a122+a132>0(aij=Aij)
所以,|A| != 0,只有一个解:|A| = 1
实对称矩阵的对角化
实数并且是对称矩阵的对角化
n阶实对称矩阵的n个特征值都是实数并且特征向量都是实向量
定理1:实对称矩阵A的不同特征值对应的特征向量一定正交
证明:AT=A
λ1,λ2(λ1!=λ2)
Aα1=λ1α1Aα2=λ2α2
只要证出内积(α1,α2)=0,就表示正交
(Aα1,α2)=(λ1α1,α2)=λ1(α1,α2)
(Aα1,α2)=(Aα1)Tα2=α1TATα2=α1TAα2=λ2α1Tα2=λ2(α1,α2)
所以,λ1(α1,α2)=λ2(α1,α2)
所以,(λ1−λ2)(α1,α2)=0
因为,λ1!=λ2
所以,(α1,α2)=0,所以特征向量一定正交
正交相似:如果A、B是同阶方阵,如果存在正交矩阵P,使得 P−1AP=B
定理2:假设A是实对称矩阵,一定存在正交矩阵Q,使得:Q−1AQ=∧=⎝⎜⎜⎛λ10⋮00λ2⋮0⋯⋯⋱⋯00⋮λn⎠⎟⎟⎞
这个定理说明:n阶实对称矩阵A一定有n个线性无关的特征向量
推论:假设A是实对称矩阵,λi是它的一个ri重特征根,那么A对应于ri的对应λi的特征向量一定有ri个(实对称矩阵一定能对角化)
计算步骤总结
3阶实对称矩阵A,求正交矩阵Q,使得Q−1AQ=∧
证明:①求特征值 ②求特征向量 ③特征向量正交化、单位化 ④特征向量做成列构成Q、特征值构成∧
注意:对角形的λ顺序要跟特征向量Q的λ一致
情况1,设:λ1=a,λ2=b,λ3=c(三个根都是单根)
根据上面的定理1:α1,α2,α3一定正交,直接单位化即可:η1=∣∣α1∣∣1α1,η2=∣∣α2∣∣1α2,η3=∣∣α3∣∣1α3
Q=(η1,η2,η3)
Q−1AQ=⎝⎛λ1000λ2000λ3⎠⎞
情况2,设:λ1=a,λ2=λ3=b(一个单根、一个二重根)
α2,α3做施密特正交化,α1不用做施密特正交化,因为一定是正交的,但是为了符号统一,可以写成β1=α1
再单位化后,得到η1,η2,η3
情况3,设:λ1=λ2=λ3=a(一个三重根)
α1,α2,α3都得做施密特正交化,得到:β1,β2,β3
再单位化后,得到η1,η2,η3