一,对角化分解
A=SΛS−1
S−1AS=Λ
A表示具备n个线性无关的x(特征向量)的矩阵
S表示由x组成的可逆方阵,称做特征向量矩阵
Λ表示由A的λ(特征值)做为对角元素的对角矩阵
比较:A=LU(消元化分解),A=QR(正交化分解)html
二,矩阵的幂公式
Ak=SΛkS−1
Akx=λkx
A和
Ak具备相同的特征向量
若是全部λ均知足|λ|<1,那么当k→∞时,
Ak→0web
三,从特征值判断矩阵是否可对角化
性质1:若是A的λ各不相同,那么A的x都线性无关,可对角化
性质2:若是A有相同的λ,那么没法肯定A的x是否线性无关
若是A是单位矩阵,λ都=1,而x都线性无关
若是A是退化矩阵(上一讲),λ相同,但x只有一个app
四,应用:求解差分方程
uk+1=Auk
前提:A可对角化
根据方程规律可导出:
uk+1=Auk=Ak+1u0
所以,方程的通解:
uk=Aku0
将初始条件
u0带入,便可得特解:
第一步:根据上一讲的方法,求出A的Λ和S
第二步:将
u0分解成
u0=SΛu,
Λu表示
u0特征值对角矩阵
第三步:根据已知条件S,求出
Λu
第四步:根据幂公式,
Aku0=SΛkS−1u0=SΛkS−1SΛu=SΛkΛu
第五步:得特解,
uk=SΛkΛusvg
五,应用:求解斐波那契数列通项公式
Fk+2=Fk+1+Fk,
F100的值
第一步,创建方程组:
{Fk+2=Fk+1+FkFk+1=Fk+1
第二步,化为矩阵:
[Fk+2Fk+1]=[1110][Fk+1Fk]
第三步,化为标准型:
uk+1=Auk
uk+1=[Fk+2Fk+1],
uk=[Fk+1Fk],
A=[1110]spa
第四步,根据上一节的方法求特解
u99=SΛ99Λu
F100为
u99第一行的值
此例中,
∣λ1∣>1,
∣λ2∣<1,所以当k→∞时,
λ1k→∞对应的解为稳态解,
λ2k→0对应的解为暂态解orm