MATLAB: ode45 求解常微分方程

引入

ode45 是 MATLAB 中用于求解非刚性常微分方程（ODE）的数值方法。它基于 Runge-Kutta 方法，并具有自适应步长调整机制，能够在一定误差控制范围内高效地计算 ODE 的数值解。

下面我们通过这个包含详细注释的代码，一起学习一下这个函数的使用：

使用 ode45 求解非线性常微分方程并绘制数值解与精确解对比图

方程模型

我们考虑以下非线性常微分方程：

$y^{\prime \prime }=-y^{\prime }+\cos (t)-3\sin (t)$

##初始条件

选择初始条件：
$y(0.1)=\cos (0.1)+2\sin (0.1)$
y ′