本人matlab小白,但苦于毕业论文需要,实在要用,matlab进行仿真.....借鉴大佬的程序改了下参数,做了个三方演化博弈的仿真,想用其他论文的仿真结果试试这个程序行不行。但得出的图和论文的仿真结果不一样。还想请各位大佬们指教 帮助 谢谢各位大佬了。
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')
论文的结果是:
我的结果是: