微分方程


简述

这篇博文主要将一些微分方程的解法记录在这里,会持续更新。

一

v(t)=exp(At)v(0)v(t) = exp(At)v(0)是微分方程ddtv(t)=Av(t)\frac{d}{dt}v(t) = Av(t)的解。

其中AA是矩阵,当然也可以是一个常数。v(0)v(0)是vv的初始值。

证明

先对v(t)=exp(At)v(0)v(t) = exp(At)v(0)进行泰勒展开,然后对v(t)v(t)进行微分。

etAv(0)=t0A0v(0)+t1A1v(0)+t22A2v(0)+t36A3v(0)+⋯+tnn!Anv(0)+⋯ddtv(t)=ddtetAv(0)=t0A1v(0)+t1A2v(0)+t22A3v(0)+t36A4v(0)+⋯+t(n−1)(n−1)!Anv(0)+⋯=A(t0A0v(0)+t1A1v(0)+t22A2v(0)+t36A3v(0)+⋯+t(n−1)(n−1)!An−1v(0)+⋯ )=AetAv(0)\begin{align} & e^{tA}v(0) = t^0A^0v(0) + t^1A^1v(0) + \frac{t^2}{2}A^2v(0) + \frac{t^3}{6}A^3v(0) + \cdots + \frac{t^n}{n!}A^nv(0) + \cdots \notag \\ & \frac{d}{dt}v(t) \notag \\ &= \frac{d}{dt}e^{tA}v(0) \notag \\ &= t^0A^1v(0) + t^1A^2v(0) + \frac{t^2}{2}A^3v(0) + \frac{t^3}{6}A^4v(0) + \cdots + \frac{t^{(n-1)}}{(n-1)!}A^{n}v(0) + \cdots \notag \\ &= A(t^0A^0v(0) + t^1A^1v(0) + \frac{t^2}{2}A^2v(0) + \frac{t^3}{6}A^3v(0) + \cdots + \frac{t^{(n-1)}}{(n-1)!}A^{n - 1}v(0) + \cdots) \notag \\ &= Ae^{tA}v(0)\notag \end{align}

所以可得,

ddtv(t)=AetAv(0)=Av(t)\begin{align} \frac{d}{dt}v(t) = Ae^{tA}v(0) = Av(t) \notag \end{align}

即当v(t)=eAtv(0)v(t)=e^{At}v(0)时,v(t)v(t)的微分是Av(t)Av(t)

所以倒退,v(t)=exp(At)v(0)v(t) = exp(At)v(0)是微分方程ddtv(t)=Av(t)\frac{d}{dt}v(t) = Av(t)的解。


Author: CAI Ming
Reprint policy: All articles in this blog are used except for special statements CC BY 4.0 reprint policy. If reproduced, please indicate source CAI Ming !
评论
  TOC