MATLAB R2016a实现三方演化博弈

%脚本1
functiondydt=cjd(t,y,a,C11,C12,L11,L12,L13,P11,P12,C21,L21,L22,P21,P22,P23,C31,L31,L32,P31,P32)

dydt=zeros(3,1);
dydt(1)= y(1)*(1-y(1))*((y(2)*((1-a)*P11+a*(1-y(3))*(C11+L13)+L12))+(1-a)*y(3)*P12+(a*y(3)-1)*(C11+L13)+y(3)*L12+C12);
dydt(2)= y(2)*(1-y(2))*(y(1)*(P23+L21-L22)+y(3)*P22+P21-C21-L21);
dydt(3)= y(3)*(1-y(3))*(y(1)*(P32-L32)+P31-C31+y(2)*L31);

end

%脚本2
a=0.1,C11=7,C12=1.5,L11=3,L12=5,L13=10,P11=2,P12=2,C21=1,L21=1,L22=1,P21=1.5,P22=1.5,P23=1.5,C31=1.5,L31=1.5,L32=1.5,P31=3,P32=2;

for i=0.1:0.2:0.9
figure(1)
[t,y]=ode45(@(t,y) cjd(t,y,a,C11,C12,L11,L12,L13,P11,P12,C21,L21,L22,P21,P22,P23,C31,L31,L32,P31,P32),[0 10],[i j k]);
plot3(y(:,1),y(:,2),y(:,3),'rh-');
hold on
end

for j=0.1:0.2:0.9
[t,y]=ode45(@(t,y) cjd(t,y,a,C11,C12,L11,L12,L13,P11,P12,C21,L21,L22,P21,P22,P23,C31,L31,L32,P31,P32),[0 10],[i j k]);
plot3(y(:,1),y(:,2), y(:,3),'bo-');
hold on
end

for k=0.1:0.2:0.9
[t,y]=ode45(@(t,y) cjd(t,y,a,C11,C12,L11,L12,L13,P11,P12,C21,L21,L22,P21,P22,P23,C31,L31,L32,P31,P32),[0 10],[i j k]);
plot3(y(:,1),y(:,2),y(:,3),'mx-');
axis([0 1 0 1 0 1])
set(gca,'XTick',[0:0.1:1],'YTick',[0:0.1:1] ,'ZTick',[0:0.1:1])
xlabel('$x$', 'interpreter','latex');
ylabel('$y$','interpreter','latex');
zlabel('$z$','interpreter','latex','rotation',360);
grid on
hold on
end

最后

以上就是执着高跟鞋为你收集整理的MATLAB R2016a实现三方演化博弈的全部内容，希望文章能够帮你解决MATLAB R2016a实现三方演化博弈所遇到的程序开发问题。

如果觉得靠谱客网站的内容还不错，欢迎将靠谱客网站推荐给程序员好友。