matlab做演化博弈分析

9049

收藏 2016-05-31

哎，做毕设要用到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('演化示意图')

扫码加我拉你入群

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

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

全部回复

cmwei333

2016-5-31 14:50:08

你把

[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]);

再试试看，应该没问题吧

扫码加我拉你入群

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

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

wangguanwen

2016-6-1 08:43:00

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

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

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

扫码加我拉你入群

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

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

wangguanwen

2016-6-1 09:20:11

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，这是怎么回事呢？求教呀

扫码加我拉你入群

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

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

cmwei333

2016-6-1 10:59:16

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'

扫码加我拉你入群

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

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

wangguanwen

2016-6-6 10:09:50

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

plot(t,y1,'b0');

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

扫码加我拉你入群

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

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

点击查看更多内容…

1392638875

2017-3-8 15:20:13

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

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

扫码加我拉你入群

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

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

wangguanwen

2017-4-21 09:44:12

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')

扫码加我拉你入群

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

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

nanianhuakai@

2017-5-28 10:37:09

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

扫码加我拉你入群

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

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

金坛映黼帷6

2017-8-28 09:46:45

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

扫码加我拉你入群

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

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

深潭微澜

2017-9-14 10:06:51

楼主你好，我运行您的程序的时候怎们总显示未定义jixiao变量

扫码加我拉你入群

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

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

深潭微澜

2017-9-14 10:10:37

金坛映黼帷6 发表于 2017-8-28 09:46
楼主您好，能否请教您一下，我用您上面的代码为什么运行不出来呢，总说出错

你好，请问楼主的代码您运行出来了吗，因为最近我也要用到演化博弈这方面的仿真在学习，不会啊

。如果您运行出来了，能不能教教我怎么运行出来的。

扫码加我拉你入群

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

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

小小数的欢乐

2019-3-2 16:08:30

请问ode5中yo初始参数如何设置呢

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群