r₁、r₂、N₁、N₂、σ₁、σ₂值必须已知,这样才能用ode45(或其他数值方法)来求解,其
微分方程组。求解方法,
function eq=fun(t,x)
eq = [ r1* x(1) * (1-x(1)/N1+sigma1*x(2)/N2); r2* x(2) * (1-x(2)/N2+sigma2*x(1)/N1)];
end
x0 = [xi;xj]; %初值
[t,x]=ode45(@fun,[0 10],x0);
plot(t,x(:,1),'--rs',t,x(:,2),'k-');
温馨提示:答案为网友推荐,仅供参考