微分方程如下:
function dxdt=green(t,x)
dxdt=[x(1)*(1-x(1))*(500*x(2)-250);x(2)*(1-x(2))*(400*x(1)-300)];
end
主函数如下:
for i=0.2:0.2:0.8
for j=0.2
[T,Y]=ode45('green',[0 0.08],[i j]);
figure(2)
grid on
plot(T,Y(:,1));
xlabel('t')
ylabel('x')
hold on
end
end
for i=0.2:0.2:0.8
for j=0.8
[T,Y]=ode45('green',[0 0.08],[i j]);
figure(3)
grid on
plot(T,Y(:,1));
xlabel('t')
ylabel('x')
hold on
end
图2是符合预期的
图3 i取0.2的时候应该是趋向于1,结果是趋向于0,请问各位是什么原因?