第二十二讲 对角化分解和幂公式
一,對角化分解
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?+Fk?Fk+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+2?Fk+1??]=[11?10?][Fk+1?Fk??]
第三步,化為標準型: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+2?Fk+1??],uk=[Fk+1Fk]u_{k}=\begin{bmatrix}F_{k+1}\\ F_{k}\end{bmatrix}uk?=[Fk+1?Fk??],A=[1110]A=\begin{bmatrix}1 & 1\\ 1 & 0\end{bmatrix}A=[11?10?]
第四步,根據上一節的方法求特解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對應的解為暫態解
總結
以上是生活随笔為你收集整理的第二十二讲 对角化分解和幂公式的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 第二十一讲 特征值和特征向量
- 下一篇: 第二十六讲 有特殊特征值的微分方程组