matlab求解微分方程组

在Matlab中，可以使用ode45函数求解微分方程组。以下是一个求解常微分方程组的例子：

假设有一个常微分方程组：

dy1/dt = 2*y1 - 3*y2
dy2/dt = -y1 + 4*y2

其初始条件为y1(0) = 1，y2(0) = 0。

代码如下：

% 定义微分方程组
f = @(t,y) [2*y(1) - 3*y(2); -y(1) + 4*y(2)];

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

% 定义时间区间
tspan = [0 10];

% 求解微分方程组
[t, y] = ode45(f, tspan, y0);

% 绘图
plot(t, y(:,1), t, y(:,2));
legend('y1', 'y2');

解释一下：

- 第一行定义了微分方程组，用匿名函数表示；
- 第三行定义了初始条件；
- 第五行定义了时间区间；
- 第七行用ode45函数求解微分方程组，返回时间和解；
- 最后一行将结果绘制出来。