MATLAB常微分方程组画图

在MATLAB中，常微分方程组（ODEs）的可视化通常涉及到使用`ode45`或`odeset`等函数来求解方程组，然后用`plot`或`ezplot`等图形函数来绘制解的路径。以下是一个简单的步骤：

1. 定义方程组：首先，你需要定义你要解决的一组常微分方程。例如，对于二维方程组 `dy/dt = f(t, y)`，你可以定义`f`函数。

function dydt = myODE(t, y)
    dydt = [y(2); -0.1*y(1)]; % 假设这是一个简单的一阶线性方程组
end

2. 求解方程组：使用`ode45`函数，提供初始条件和时间范围。

tspan = [0 10]; % 时间范围
y0 = [1; 0]; % 初始条件
[t, y] = ode45(@myODE, tspan, y0);

3. 画出解：最后，使用`plot`或`hold on`来在同一个图上叠加多个曲线。

figure;
plot(t, y(:, 1), 'r', 'LineWidth', 2); % 第一个方程的解
hold on;
plot(t, y(:, 2), 'b', 'LineWidth', 2); % 第二个方程的解
xlabel('Time');
ylabel('Solution');
legend('Equation 1', 'Equation 2');