MIT 线性代数 22 对角化和A的幂，差分方程的线性代数解法

2022-05-19  本文已影响0人  光能蜗牛

一.对角化和矩阵的幂次计算

假设矩阵A有n个线性无关的特征向量X_1,X_2,X_3,......,X_n，这些特征向量按列组成特征向量矩阵S=\begin{bmatrix}X_1&X_2&X_3&....&X_n\end{bmatrix}
则AS=A\begin{bmatrix}X_1&X_2&X_3&....&X_n\end{bmatrix}=\begin{bmatrix}\lambda_1X_1&\lambda_2X_2&\lambda_3X_3&....&\lambda_nX_n\end{bmatrix}

=\begin{bmatrix}X_1&X_2&X_3&....&X_n\end{bmatrix} \begin{bmatrix} \lambda_1& & &\cdots&0 \\0&\lambda_2& &\cdots&0 \\0& &\lambda_3&\cdots&0 \\\vdots& & &\ddots&\vdots \\0& & &\cdots& \lambda_n \end{bmatrix}
=S \Lambda （\Lambda称为对角特征值矩阵）

强调一遍，并不是所有的矩阵都会有n个线性无关的特征向量，关于这一点，我们在上一节就发现了一些无法求解出多个不同的线性无关特征向量的矩阵

也就是AS=S\Lambda
我们可以写成如下形式
S^{-1}AS=\Lambda 这就是矩阵A的对角化方法

或
A=S\Lambda S^{-1}

其中我们把A=S\Lambda S^{-1} 这是一个新的矩阵分解方法，我们之前消元的时候有LU分解，以及施密特正交化的QR分解

如果我们尝试对A的幂次方做计算，比如A的平方
A^2=S\Lambda S^{-1}S\Lambda S^{-1}=S\Lambda^2 S^{-1}
于是同理可得有
A^k=S\Lambda^k S^{-1}
也就是A的这种特征值对角化分解方式对于求A的幂次计算非常方便

二.差分方程线性代数解法

给定差分方程u_{K+1}=Au_{k}
u_{1}=Au_{0}
u_{2}=Au_{1}
u_{k}=A^ku_{0}
这样的问题，当A可对角化时，即存在n个 不同的特征向量以及对应的特征值的情况下，即A是满秩的，也是可逆的
那么很显然可以把u_0 用A的特征向量进行线性组合表示，

于是
u_0=c_1X_1+c_2X_2+...+c_nX_n 其中X_1,X_2.....X_n是矩阵A的特征向量
Au_0=c_1\lambda_1 X_1+c_2\lambda_2 X_2+...+c_n\lambda_n X_n,其中\lambda_1,\lambda_2......\lambda_n表示矩阵A的各个特征向量对应的特征值
于是我们得到
A^ku_0=c_1\lambda_1^k X_1+c_2\lambda_2^k X_2+...+c_n\lambda_n^k X_n
=S\Lambda^kC

这里举了一个斐波那契数列的例子
众所周知，斐波那契数列排列方式是0,1,1,2,3,5,8,13,21.........，

即有F_{k+2}=F_{k+1}+F_k-------(k>=0)

问题是如何求解F_{100}

参考前面的例子我们可以这么写
追加一个等式得到方程组
F_{k+2}=F_{k+1}+F_k
F_{k+1}=F_{k+1}
如果我们令向量u_k=\begin{bmatrix}F_{k+1}\\F_k\end{bmatrix}
于是就会有
u_{k+1}=\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}

也就是u_{k+1}=\begin{bmatrix}1&1\\1&0\end{bmatrix}u_k

其中特征矩阵c

令|A-\lambda I|=0
马上可以求出特征值是\lambda_1 =\frac{1+\sqrt{5}}{2} =1.618.....,\lambda_2=\frac{1-\sqrt{5}}{2}=-0.618.....
其特征向量即A-\lambda I的零空间 有
X_1=\begin{bmatrix}\lambda_1 \\1\end{bmatrix}
X_2=\begin{bmatrix}\lambda_2 \\1\end{bmatrix}

由前面的推论我们知道,特征向量X_1和X_2,进行线性组合能得到u_0
所以u_{0}=c_1 X_1+c_2X_2
而A又是可对角化的矩阵，所以
u_{1}=Au_0=c_1 AX_1+c_2AX_2=c_1\lambda_1 X_1+c_2\lambda_2 X_2
u_{100}=A^{100}u_0=c_1\lambda_1^{100} X_1+c_2\lambda_2^{100} X_2

我们注意到第二项\lambda_2为负数，意味着平方后这一项c_2\lambda_2^{100} X_2整体收敛到无穷小，这意味着u_k数组系列的取值基本上由\lambda_1决定

因为u_{0}=c_1 X_1+c_2X_2=\begin{bmatrix}F_{1}\\F_0\end{bmatrix}=\begin{bmatrix}1\\0\end{bmatrix}

于是可以解出

c_1=\frac{1}{\sqrt{5}} ,,c_2=\frac{-1}{\sqrt{5}}

于是我们甚至求出了斐波那契数列的通项公式

u_k=\begin{bmatrix}F_{k+1}\\F_k\end{bmatrix}=c_1\lambda_1^{k} X_1+c_2\lambda_2^{k} X_2
=\begin{bmatrix}X_1 &X_2\end{bmatrix}\begin{bmatrix}\lambda_1^k &0\\0&\lambda_2^k \end{bmatrix}\begin{bmatrix}c_1 \\c_2 \end{bmatrix}
=\begin{bmatrix}\lambda_1 &\lambda_2 \\1&1\end{bmatrix}\begin{bmatrix}\lambda_1^k &0\\0&\lambda_2^k \end{bmatrix}\begin{bmatrix}c_1 \\c_2 \end{bmatrix}

我们只看下面的F_k的展开

F_k=c_1\lambda_1^k+c_2\lambda_2^k=\frac{1}{\sqrt{5}}[(\frac{1+\sqrt{5}}{2})^k-(\frac{1-\sqrt{5}}{2})^k]

这就是斐波那契数列的通项公式

上一篇下一篇

猜你喜欢

热点阅读