matlab ode45函数

MATLAB的ode45函数是用于求解一阶或二阶常微分方程组的数值解。它使用的是常规的Runge-Kutta方法，适用于解决大多数常微分方程问题。

ode45函数的语法如下：

[t,y] = ode45(odefun,tspan,y0)

其中，t是一个列向量，包含求解的时间点，y是一个矩阵，每一行对应一个时间点的解，odefun是一个函数句柄，用于计算微分方程的右手边，tspan是一个包含求解时间区间的向量，y0是一个列向量，包含初始条件。

例如，下面的代码用ode45函数求解一阶常微分方程 y' = -2y，y(0) = 1：

% 定义微分方程右手边的函数
function dydt = myodefun(t,y)
    dydt = -2*y;
end

% 求解微分方程
[t,y] = ode45(@myodefun,[0 5],1);

% 绘制结果
plot(t,y);
xlabel('t');
ylabel('y');

运行结果如下图所示：