欢迎访问 生活随笔!

生活随笔

当前位置: 首页 > 编程资源 > 编程问答 >内容正文

编程问答

第二十二讲 对角化分解和幂公式

发布时间:2025/4/16 编程问答 64 豆豆
生活随笔 收集整理的这篇文章主要介绍了 第二十二讲 对角化分解和幂公式 小编觉得挺不错的,现在分享给大家,帮大家做个参考.

一,对角化分解

A=SΛS−1A=S\Lambda S^{-1}A=SΛS1

S−1AS=ΛS^{-1}AS=\LambdaS1AS=Λ
A表示具有n个线性无关的x(特征向量)的矩阵
S表示由x组成的可逆方阵,称作特征向量矩阵
Λ\LambdaΛ表示由A的λ(特征值)作为对角元素的对角矩阵
比较:A=LU(消元化分解),A=QR(正交化分解)

二,矩阵的幂公式

Ak=SΛkS−1A^{k}=S\Lambda ^{k}S^{-1}Ak=SΛkS1

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ΛkS1u0=SΛkS1SΛ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+FkF100F_{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 &amp; 1\\ 1 &amp; 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 &amp; 1\\ 1 &amp; 0\end{bmatrix}A=[1110]

第四步,根据上一节的方法求特解u99=SΛ99Λuu_{99}=S\Lambda ^{99}\Lambda _{u}u99=SΛ99Λu

F100F_{100}F100u99u_{99}u99第一行的值
此例中,∣λ1∣&gt;1\left | \lambda _{1} \right |&gt; 1λ1>1∣λ2∣&lt;1\left | \lambda _{2} \right |&lt; 1λ2<1,因此当k→∞时,λ1k→∞\lambda _{1}^{k}→\inftyλ1k对应的解为稳态解,λ2k→0\lambda _{2}^{k} →0λ2k0对应的解为暂态解

总结

以上是生活随笔为你收集整理的第二十二讲 对角化分解和幂公式的全部内容,希望文章能够帮你解决所遇到的问题。

如果觉得生活随笔网站内容还不错,欢迎将生活随笔推荐给好友。