第二十二讲 对角化分解和幂公式
一,对角化分解
A=SΛS−1A=S\Lambda S^{-1}A=SΛS−1
S−1AS=ΛS^{-1}AS=\LambdaS−1AS=Λ
A表示具有n个线性无关的x(特征向量)的矩阵
S表示由x组成的可逆方阵,称作特征向量矩阵
Λ\LambdaΛ表示由A的λ(特征值)作为对角元素的对角矩阵
比较:A=LU(消元化分解),A=QR(正交化分解)
二,矩阵的幂公式
Ak=SΛkS−1A^{k}=S\Lambda ^{k}S^{-1}Ak=SΛkS−1
Akx=λkxA^{k}x=\lambda ^{k}xAkx=λkx
A和AkA^{k}Ak具有相同的特征向量
如果所有λ均满足|λ|<1,那么当k→∞时,AkA^{k}Ak→0
三,从特征值判断矩阵是否可对角化
性质1:如果A的λ各不相同,那么A的x都线性无关,可对角化
性质2:如果A有相同的λ,那么无法确定A的x是否线性无关
如果A是单位矩阵,λ都=1,而x都线性无关
如果A是退化矩阵(上一讲),λ相同,但x只有一个
四,应用:求解差分方程uk+1=Auku_{k+1}=Au_{k}uk+1=Auk
前提:A可对角化
根据方程规律可导出:uk+1=Auk=Ak+1u0u_{k+1}=Au_{k}=A^{k+1}u_{0}uk+1=Auk=Ak+1u0
因此,方程的通解:uk=Aku0u_{k}=A^{k}u_{0}uk=Aku0
将初始条件u0u_{0}u0带入,即可得特解:
第一步:根据上一讲的方法,求出A的Λ和S
第二步:将u0u_{0}u0分解成u0=SΛuu_{0}=S\Lambda _{u}u0=SΛu,Λu\Lambda _{u}Λu表示u0u_{0}u0特征值对角矩阵
第三步:根据已知条件S,求出Λu\Lambda _{u}Λu
第四步:根据幂公式,Aku0=SΛkS−1u0=SΛkS−1SΛu=SΛkΛuA^{k}u_{0}=S\Lambda ^{k}S^{-1}u_{0}=S\Lambda ^{k}S^{-1}S\Lambda _{u}=S\Lambda ^{k}\Lambda _{u}Aku0=SΛkS−1u0=SΛkS−1SΛu=SΛkΛu
第五步:得特解,uk=SΛkΛuu_{k}=S\Lambda ^{k}\Lambda _{u}uk=SΛkΛu
五,应用:求解斐波那契数列通项公式Fk+2=Fk+1+FkF_{k+2}=F_{k+1}+F_{k}Fk+2=Fk+1+Fk,F100F_{100}F100的值
第一步,建立方程组:{Fk+2=Fk+1+FkFk+1=Fk+1\left\{\begin{matrix}F_{k+2}=F_{k+1}+F_{k}\\ F_{k+1}=F_{k+1}\end{matrix}\right.{Fk+2=Fk+1+FkFk+1=Fk+1
第二步,化为矩阵:[Fk+2Fk+1]=[1110][Fk+1Fk]\begin{bmatrix}F_{k+2}\\ F_{k+1}\end{bmatrix}=\begin{bmatrix}1 & 1\\ 1 & 0\end{bmatrix}\begin{bmatrix}F_{k+1}\\ F_{k}\end{bmatrix}[Fk+2Fk+1]=[1110][Fk+1Fk]
第三步,化为标准型:uk+1=Auku_{k+1}=Au_{k}uk+1=Auk
uk+1=[Fk+2Fk+1]u_{k+1}=\begin{bmatrix}F_{k+2}\\ F_{k+1}\end{bmatrix}uk+1=[Fk+2Fk+1],uk=[Fk+1Fk]u_{k}=\begin{bmatrix}F_{k+1}\\ F_{k}\end{bmatrix}uk=[Fk+1Fk],A=[1110]A=\begin{bmatrix}1 & 1\\ 1 & 0\end{bmatrix}A=[1110]
第四步,根据上一节的方法求特解u99=SΛ99Λuu_{99}=S\Lambda ^{99}\Lambda _{u}u99=SΛ99Λu
F100F_{100}F100为u99u_{99}u99第一行的值
此例中,∣λ1∣>1\left | \lambda _{1} \right |> 1∣λ1∣>1,∣λ2∣<1\left | \lambda _{2} \right |< 1∣λ2∣<1,因此当k→∞时,λ1k→∞\lambda _{1}^{k}→\inftyλ1k→∞对应的解为稳态解,λ2k→0\lambda _{2}^{k} →0λ2k→0对应的解为暂态解
总结
以上是生活随笔为你收集整理的第二十二讲 对角化分解和幂公式的全部内容,希望文章能够帮你解决所遇到的问题。
- 上一篇: 第二十一讲 特征值和特征向量
- 下一篇: 第二十六讲 有特殊特征值的微分方程组