如题所述
用ode45()函数求解微分方程组的关键是建立起自定义函数,该函数的主要内容为
function dy=odefun(t,x); %保存为 odefun.m 文件
mu1=1500;mu2=150;
dy(1)=x(1)*(1-x(2)/mu2);
dy(2)=-x(2)*(1-x(1)/mu1);
dy=dy(:);
end
然后用ode45()函数命令,调用odefun()函数。编程后可以得到如下结果。