p0、q0、z0对应论文里的x、y、z
function dy=jixiao(t,y,g,r,c,k,v,s,n,o,p)
dy=zeros(3,1);
dy(1)=y(1)*(1-y(1))*(y(2)*g*-1+g+r-c);
dy(2)=y(2)*(1-y(2))*(y(3)*(y(1)*k*v+s-p)-v-n-y(1)*g);
dy(3)=y(3)*(1-y(3))*(y(1)*k*v*1-y(2)-o);
end
p0=0.2,q0=0.4,z0=0.8;
s=5,v=10,o=2,n=3,r=5,k=0.5,g=8,c=3,p=2;
[t,y]=ode45(@(t,y)jixiao(t,y,g,r,c,k,v,s,n,o,p),[0:10],[p0,q0,z0]);
y1=y(:,1);
plot(t,y1,'b*');
xlabel('t');
ylabel('p');
hold on;
title('演化示意图');
legend('p0=0.2')
论文的结果是:
我的结果是: