matlab做演化博弈分析

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

哎，做毕设要用到MATLAB做演化博弈分析，程序已经编好了，只是不知道哪里不对总是显示Missing variable or function.
求诸位大神帮忙看看吧,急求啊！！！！！
jixiao.m
function dy=jixiao(t,y,R,M,A,E,S)
dy=zeros(2,1);
dy(1)=y(1)*(1-y(1))*(y(2)*(S-A-M)-(R-M));
dy(2)=y(2)*(1-y(2))*(y(1)*(A+M)-E);
end


主程序.m
p0=0.2,q0=0.5,R=0.25,M=0.2,A=0.32,E=0.48,S=0.8;
[t,y]=ode45('jixiao',[0:40],[p0,q0],[],R,M,A,E,S);
y1=y(:,1);
plot(t,y1,'b*');
xlabel('t');
ylabel('p');
hold on;

p0=0.4,q0=0.5;
[t,y]=ode45('jixiao',[0:40],[p0,q0],[],R,M,A,E,S);
y1=y(:,1);
plot(t,y1,'b0');
xlabel('t');
ylabel('p');
hold on;

title('演化示意图')

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏28 回帖

关键词：MATLAB matla atlab 演化博弈 博弈分析 matlab

你把

[t,y]=ode45('jixiao',[0:40],[p0,q0],[],R,M,A,E,S);

换成

[t,y]=ode45(@(t,y) jixiao(t,y,R,M,A,E,S),[0:40],[p0,q0]);

再试试看，应该没问题吧

cmwei333 发表于 2016-5-31 14:50
你把

[t,y]=ode45('jixiao',[0:40],,[],R,M,A,E,S);

太谢谢啦！！！简直解救我于水火之中啊，我在这上面已经卡了好长时间了，终于出图了

cmwei333 发表于 2016-5-31 14:50
你把

[t,y]=ode45('jixiao',[0:40],,[],R,M,A,E,S);

p0=0.2,q0=0.5;
R=0.25,M=0.2,A=0.32,E=0.48,S=0.8;
[t,y]=ode45(@(t,y)jixiao(t,y,R,M,A,E,S),[0:50],[p0,q0]);
y1=y(:,1);
plot(t,y1,'b*');
xlabel('t');
ylabel('p');
hold on;

p0=0.4,q0=0.5;
R=0.25,M=0.2,A=0.32,E=0.48,S=0.8;
[t,y]=ode45(@(t,y)jixiao(t,y,R,M,A,E,S),[0:50],[p0,q0]);
y1=y(:,1);
plot(t,y1,'b0');
xlabel('t');
ylabel('p');
hold on;
按照您说的我改了程序，可是，我是想画出两条线对比，可是画出来的图只出一条线，而且纵坐标的最高点是0.25而不是1，这是怎么回事呢？求教呀

wangguanwen 发表于 2016-6-1 09:20
p0=0.2,q0=0.5;
R=0.25,M=0.2,A=0.32,E=0.48,S=0.8;
[t,y]=ode45(@(t,y)jixiao(t,y,R,M,A,E,S),[0:50], ...

你需要把

plot(t,y1,'b0');

换成

plot(t,y1,'bo');


出来的图就应该是这个样的



MATLAB只承认 'o'，而不是 '0'

cmwei333 发表于 2016-6-1 10:59
你需要把

plot(t,y1,'b0');

对的对的，把‘b0’换成‘bo’就出图啦，太谢谢啦

wangguanwen 发表于 2016-6-6 10:09
对的对的，把‘b0’换成‘bo’就出图啦，太谢谢啦

楼主解决了吗，最近在写毕业论文，也是写的这个，不会呀，能教教我吗，谢谢！在线等，急

1392638875 发表于 2017-3-8 15:20
楼主解决了吗，最近在写毕业论文，也是写的这个，不会呀，能教教我吗，谢谢！在线等，急

最近没上线，所以看到的比较慢，主程序不变
p0=0.2,q0=0.5;
R=0.25,M=0.2,A=0.3,E=0.48,S=0.8;
[t,y]=ode45(@(t,y)jixiao(t,y,R,M,A,E,S),[0:50],[p0,q0]);
y1=y(:,1);
plot(t,y1,'b*');
xlabel('t');
ylabel('p');
hold on;

p0=0.4,q0=0.5;
R=0.25,M=0.2,A=0.3,E=0.48,S=0.8;
[t,y]=ode45(@(t,y)jixiao(t,y,R,M,A,E,S),[0:50],[p0,q0]);
y1=y(:,1);
plot(t,y1,'b-');
xlabel('t');
ylabel('p');
hold on;

p0=0.6,q0=0.5;
R=0.25,M=0.2,A=0.3,E=0.48,S=0.8;
[t,y]=ode45(@(t,y)jixiao(t,y,R,M,A,E,S),[0:50],[p0,q0]);
y1=y(:,1);
plot(t,y1,'b^');
xlabel('t');
ylabel('p');
hold on;

p0=0.8,q0=0.5;
R=0.25,M=0.2,A=0.3,E=0.48,S=0.8;
[t,y]=ode45(@(t,y)jixiao(t,y,R,M,A,E,S),[0:50],[p0,q0]);
y1=y(:,1);
plot(t,y1);
xlabel('t');
ylabel('p');
hold on;

title('演化示意图');
legend('p0=0.2','p0=0.4','p0=0.6','p0=0.8')

楼主 能否请教你一下 啊 我最近也是在做演化博弈 不知道怎样用这个matlab啊

楼主您好，能否请教您一下，我用您上面的代码为什么运行不出来呢，总说出错