用MATLAB的ode45求解三元二阶微分方程代码

以下是使用MATLAB的ode45求解三元二阶微分方程的代码示例：

function dydt = ode_fun(t, y)
% y(1): x
% y(2): y
% y(3): z
% y(4): dx/dt
% y(5): dy/dt
% y(6): dz/dt

% Define parameters
a = 0.7;
b = 0.8;
c = 3.5;

% Define equations
dydt = [y(4); y(5); y(6);
    -a*y(4) + b*y(5)*y(6);
    -b*y(4)*y(6) + c*y(5);
    -y(4)*y(5)];

% Set initial conditions
y0 = [0; 0; 0; 1; 1; 1];

% Define time interval
tspan = [0 10];

% Solve ODE using ode45
[t, y] = ode45(@ode_fun, tspan, y0);

% Plot results
plot3(y(:,1), y(:,2), y(:,3));
xlabel('x');
ylabel('y');
zlabel('z');

在上面的代码中，首先定义了一个包含三个未知数和它们的导数的向量`y`。然后定义了微分方程的右侧，即`dydt`。对于这个特定的微分方程，`dydt`包含了三个未知数的导数，以及它们的表达式。接下来，设置了初始条件`y0`和时间间隔`tspan`。最后，使用MATLAB的`ode45`函数解决微分方程，并用`plot3`函数绘制结果。