Matlab计算常微分方程

在 MATLAB 中，可以使用 `ode45` 函数来求解常微分方程。

例如，考虑求解一阶常微分方程 $y' = f(t,y)$，其中 $f(t,y)$ 是一个已知的函数。可以使用以下代码来求解：

% 定义函数 f(t,y)
f = @(t,y) t - y;

% 定义初始条件
t0 = 0;
y0 = 1;

% 定义时间范围
tspan = [0 5];

% 使用 ode45 函数求解
[t,y] = ode45(f, tspan, y0);

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

这段代码中，首先定义了函数 `f(t,y)`，然后定义了初始条件 `t0` 和 `y0`，以及时间范围 `tspan`。接着使用 `ode45` 函数求解，得到时间数组 `t` 和对应的解 `y`。最后使用 `plot` 函数绘制结果。

需要注意的是，`ode45` 函数的第一个参数是一个函数句柄，即一个指向函数的指针。在本例中，使用了匿名函数的方式定义了函数 `f(t,y)`，也可以将其定义为一个独立的函数文件。