matlab ode45输出数值解,ode45求微分方程组数值解的问题

%函数

function dxdt=differential(t,x)

dxdt=[x(1)*(1-x(1))*(195.75-223.23*x(2));x(2)*(1-x(2))*(51.75-38.525*x(1))];

end复制代码

clear

%y-x

for i=0:0.1:1

for j=0:0.1:1

[T,Y]=ode45('differential',[0 5],[i j]);

figure(1)

grid on

plot(Y(:,1),Y(:,2));

xlabel('α');

ylabel('β');

hold on

title('演化示意图')

end

end复制代码这是做演化博弈时用到的一段程序，i,j表示的是概率不会超过1，运行发现初值为i=0.7,j=0.2和i=0.9,j=0.2时，出的图是错的。

i=0.9,j=0.2.jpg (23.99 KB, 下载次数: 4)

i=0.9,j=0.2

2020-2-29 15:10 上传