matlab 画图，微分方程组

来源：百度知道编辑：UC知道时间：2024/06/25 08:51:19

帮忙画个图，matlab装太麻烦了！
下面的参数可以随便赋初值，请给出源码和截图，谢谢！

先编一个M函数文件

function x=rr(t,u)

%下面的参数值自己设置合适的，否则不一定有解，画不出图来。

b1=2;

d1=-3;

q=-5;

beta=2;

a=-3;

gamma=-5;

x=[b1*u(2)-u(1)*(d1+u(1)+u(2))-q*u(1)*u(3);

beta*u(2)*(1-u(1)-u(2))-a*u(2)*u(3);

gamma*u(3)*(1-u(3)/(u(1)+u(2)))];

再使用命令

t_final=10;%终止时间，可调

u0=[-1;0.7;1]; %0时刻初始值，可调

[t,u]=ode45('rr',[0,t_final],u0);

plot(t,u(:,1),'r',t,u(:,2),'g',t,u(:,3),'k');