Question

求助MATLAB高手，给出一个三维的微分方程组，条件齐全，求在一个图内显示三条原方程轨迹及方程组的相轨线

初始条件及各常数设置如图所示，大神们救命~
要做两个图，图1是三条方程轨迹，图2是相轨线

Answer 1

 

clear all
clc
r=1;a=0.1;m=0.5;b=0.02;c=0.1;n=0.6;d=0.05;x0=100;y0=30;z0=5;
f=@(t,y)([r*y(1)-a*y(1)*y(2);
         -m*y(2)+b*y(1)*y(2)-c*y(2)*y(3);
          -n*y(3)+d*y(2)*y(3)]);
[t y]=ode45(f,0:0.5:50,[x0 y0 z0]);
plot(t,y(:,1),t,y(:,2),t,y(:,3))
legend('t-x','t-y','t-z')
xlabel('t')

Answer 2

采用龙格库塔数值算法求解，然后作图。