初值的线性微分方程组求解 | 例题

陌雨’

于 2022-05-05 16:00:47 发布

阅读量557

点赞数

分类专栏：常微分方程文章标签：微分方程

抓住一只白嫖怪(｀・ω・´)

本文链接：https://blog.csdn.net/Infinity_07/article/details/124581775

版权

线性微分方程组特征值特征向量解法矩阵指数函数

关键词由CSDN通过智能技术生成

常微分方程专栏收录该内容

6 篇文章 4 订阅

订阅专栏

问题重述

本博客旨在求解带初值的线性微分方程组问题。

$A(t)x+f(t),\quad x(t_0)=x_0$

理论证明这里就不想写了，只写一下具体怎么求解。

求解公式

$x(t)=\exp(A(t-t_0))x_0+\exp(At)\int_{t_0}^t\exp(-As)f(s)\,{\rm d}s$

示例

$x'=\begin{pmatrix}1& 0& 0\\ 2 & 1 & -2\\ 3 & 2 & 1\end{pmatrix}x+\begin{pmatrix}0\\0\\e^t\cos(2t)\end{pmatrix},\quad x(0)=\begin{pmatrix}0\\1\\1\end{pmatrix}$

首先求解 $\Phi (t)$

$A$ 的特征值和特征向量为 $1, 1 - 2 i, 1 + 2 i$ ，对应的特征向量为 $[1,-\frac{3}{2},1]^T$ , $0,-i,1]^T$ , $0,i,1]^T$ 。于是

$x_1(t)=\begin{bmatrix}1\\-\frac{3}{2}\\1\end{bmatrix}e^t$

$\begin{aligned}\bar{x}_{2,3}&=\begin{bmatrix}0\\-i\\\ 1\end{bmatrix}e^{(1-2i)t}\\ &=e^t\begin{bmatrix}0\\-ie^{-2it}\\e^{-2it}\end{bmatrix}\\ &=e^t\begin{bmatrix}0\\ -i(\cos2t-i\sin2t)\\\cos2t-i\sin2t\end{bmatrix}\\ &=e^t\begin{bmatrix}0\\-\sin2t\\\cos2t\end{bmatrix}+ie^t\begin{bmatrix}0\\-\cos2t\\-\sin2t\end{bmatrix}\end{aligned}$

于是就能得到，

$x_2(t)=e^t\begin{bmatrix}0\\-\sin2t\\\cos2t\end{bmatrix},\quad,x_3(t)=e^t\begin{bmatrix}0\\-\cos2t\\-\sin2t\end{bmatrix}$

于是 $\Psi(t)=e^t\begin{bmatrix}1 & 0 & 0\\ -\frac{3}{2} &-\sin2t&-\cos2t\\ 1 & \cos2t &-\sin2t\end{bmatrix}$

$\begin{aligned} \exp At&=\Psi(t)*\Psi^{-1}(0)\\ &=e^t\begin{bmatrix}1 & 0 & 0\\ -\frac{3}{2} &-\sin2t&-\cos2t\\ 1 & \cos2t &-\sin2t\end{bmatrix}*\begin{bmatrix} 1 & 0 & 0\\ -1 & 0 & 1\\ -\frac{2}{3} & -1 & 0 \end{bmatrix}\\ &= e^t\begin{bmatrix} 1 & 0 & 0\\ \frac{{\left(2\,\cos \left(2\,t\right)+3\,\sin \left(2\,t\right)-2\right)}}{3} & \cos2t & -\sin2t \\ \frac{{\left(2\,\sin \left(2\,t\right)-3\,\cos \left(2\,t\right)+3\right)}}{3} &\sin2t & \cos2t \end{bmatrix}\\ \end{aligned}$

于是

$\begin{aligned} x(t)&=\exp(A(t-t_0))x_0+\exp(At)\int_{t_0}^t\exp(-As)f(s)\,{\rm d}s\\ &=\exp(At)\begin{bmatrix}0\\1\\1\end{bmatrix}+\exp(At)\int_0^t\exp-st\begin{bmatrix}0\\0\\e^s\cos 2s\\ \end{bmatrix}{\rm d}s\\ &=\exp At\begin{bmatrix}0\\1\\1\end{bmatrix}+\exp At\int_0^t\begin{bmatrix}0\\\cos2s\sin2s\\\cos^22s\end{bmatrix}{\rm d}s\\ &=\exp At*\left(\begin{bmatrix}0\\1\\1\end{bmatrix}+\begin{bmatrix} 0\\ \frac{{\sin \left(2\,t\right)}^2 }{4}\\ \frac{t}{2}+\frac{\sin \left(4\,t\right)}{8} \end{bmatrix} \right)\\ &=e^t\begin{bmatrix} 0\\ \cos \left(2\,t\right)-\sin \left(2\,t\right)-\frac{t\,\sin \left(2\,t\right)}{2}\\ \cos \left(2\,t\right)+\frac{5\,\sin \left(2\,t\right)}{4}+\frac{t\,\cos \left(2\,t\right)}{2} \end{bmatrix}\end{aligned}$