楼主: 115861
82913 91

[问答] 【求教】运用Matlab做演化博弈仿真分析   [推广有奖]

  • 4关注
  • 35粉丝

副教授

37%

还不是VIP/贵宾

-

威望
0
论坛币
125362 个
通用积分
402.3751
学术水平
15 点
热心指数
36 点
信用等级
13 点
经验
28084 点
帖子
515
精华
0
在线时间
1070 小时
注册时间
2010-3-31
最后登录
2024-4-14

相似文件 换一批

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
最近写演化博弈方面的论文,最后一部分需要用Matlab软件做一个数值分析,请大神指导!下面是我从百度中搜索到的别人分享的程序,但是没学过Matlab,看不懂,运行之后总不能出现结果

%函数

function dxdt=differential(t,x)

dxdt=[x(1)*(1-x(1))*(5-6*x(2));x(2)*(1-x(2))*(1-6*x(1))];

end

%主程序

clear

%y-x

for i=0:0.1:1

    for j=0:0.1:1

        [T,Y]=ode45('differential',[0 5],[i j]);

        figure(1)

        grid on

        plot(Y(:,1),Y(:,2));

        hold on

    end

end

%dx/dt-t

for i=0:0.1:1

    for j=0:0.1:1

        [T,Y]=ode45('differential',[0 5],[i j]);

        figure(2)

        grid on

        plot(T,Y(:,1));

        hold on

    end

end

%dy/dt-t

for i=0:0.1:1

    for j=0:0.1:1

        [T,Y]=ode45('differential',[0 5],[i j]);

        figure(3)

        grid on

        plot(T,Y(:,2));

        hold on

    end

end

主要想求教两个问题:1,如何把这个程序在Matlab运行处结果?2,这个程序每一步什么意思?如果我需要把自己数据带入,需要修改哪些内容?

这个程序的结果

Y-X图像

QQ图片20141229164412.jpg

QQ图片20141229164342.png

dx/dt-t图像

QQ图片20141229164537.jpg

QQ图片20141229164515.png

dy/dt-t图像


QQ图片20141229164729.jpg

求大神指教,不胜感激!


二维码

扫码加我 拉你入群

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

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

关键词:MATLAB matla atlab 仿真分析 演化博弈 演化博弈 数值分析 仿真分析 Matlab博弈

QQ图片20141229164718.jpg (70.06 KB)

QQ图片20141229164718.jpg

回帖推荐

walqy222 发表于4楼  查看完整内容

把前四行和后面的分别保存成文件就可以运行。前一个是函数,后一个才是主程序。 接着是ode45,这是数值求解常微分方程的。('differential',[0 5],),第一部分是函数,[0 5]是变量区间,是初值。至于函数为啥长那个样子,你肯定比谁都明白,我只能帮你到这了。
已有 1 人评分经验 收起 理由
我的素质低 + 20 鼓励积极发帖讨论

总评分: 经验 + 20   查看全部评分

沙发
115861 发表于 2014-12-29 16:53:21 |只看作者 |坛友微信交流群
搜索好久也没找到一个具体的操作方案,现在主要两问题,一个就是这个程序不知道怎么在Matlab中实现;另一个就是不太明白程序的意思,不知道如果如果换成自己的数据,需要修改哪些方面?

使用道具

藤椅
115861 发表于 2014-12-29 16:54:39 |只看作者 |坛友微信交流群
自己也下载Matlab尝试了,但是运行不出结果,不知道是程序问题,还是自己操作问题

使用道具

板凳
walqy222 发表于 2014-12-29 17:16:47 |只看作者 |坛友微信交流群
把前四行和后面的分别保存成文件就可以运行。前一个是函数,后一个才是主程序。

接着是ode45,这是数值求解常微分方程的。('differential',[0 5],[i j]),第一部分是函数,[0 5]是变量区间,[i j]是初值。至于函数为啥长那个样子,你肯定比谁都明白,我只能帮你到这了。
已有 5 人评分论坛币 学术水平 热心指数 信用等级 收起 理由
sohji + 1 + 1 + 1 精彩帖子
gobaw + 5 + 1 + 1 + 1 精彩帖子
115861 + 5 + 1 + 1 热心帮助其他会员
crystal8832 + 36 + 1 + 1 热心帮助其他会员
statax + 10 + 1 精彩帖子

总评分: 论坛币 + 56  学术水平 + 4  热心指数 + 5  信用等级 + 2   查看全部评分

使用道具

报纸
115861 发表于 2014-12-29 17:26:58 |只看作者 |坛友微信交流群
walqy222 发表于 2014-12-29 17:16
把前四行和后面的分别保存成文件就可以运行。前一个是函数,后一个才是主程序。

接着是ode45,这是数值求 ...
谢谢啦~~我运行中还遇到一点问题,我把前四行定义为函数了,后面的以什么格式保存?我直接把后面的粘在命令窗口,但是运行之后,不明白为什么提示:未定义函数或变量 'differential'。

使用道具

地板
115861 发表于 2014-12-29 17:34:02 |只看作者 |坛友微信交流群
walqy222 发表于 2014-12-29 17:16
把前四行和后面的分别保存成文件就可以运行。前一个是函数,后一个才是主程序。

接着是ode45,这是数值求 ...
QQ图片20141229173240.jpg

使用道具

7
walqy222 发表于 2014-12-29 17:40:40 |只看作者 |坛友微信交流群
1.把前四行保存成文件'differential.m'
2.把后面的保存成文件'XXX.m',名称任意
3.运行'XXX.m'文件
已有 4 人评分论坛币 学术水平 热心指数 信用等级 收起 理由
妄语9 + 1 谢谢楼主
蓝天明 + 1 + 1 + 1 热心帮助其他会员
Xaero + 40 + 1 + 1 精彩帖子
statax + 10 + 1 精彩帖子

总评分: 论坛币 + 50  学术水平 + 2  热心指数 + 4  信用等级 + 1   查看全部评分

使用道具

8
115861 发表于 2014-12-29 17:52:31 |只看作者 |坛友微信交流群
walqy222 发表于 2014-12-29 17:40
1.把前四行保存成文件'differential.m'
2.把后面的保存成文件'XXX.m',名称任意
3.运行'XXX.m'文件
大神厉害!终于成功运行出来结果了,非常感谢

使用道具

9
115861 发表于 2014-12-29 18:12:32 |只看作者 |坛友微信交流群
经过大神指点,终于运行出想要的结果了在这里我把自己的操作总结一下,一方面表达对@walqy222的感谢,另一方面也和大家分享一下(1)这段程序分为两个部分:第一部分前面四行;第二部分从第五行到最后
(2)在matlab界面点击新建——函数,然后把前面四行粘进去,关闭编辑器页面,软件提示是否保存?将文件保存为differential.m
(3)在matlab界面点击新建——函数,然后把后面内容粘进去,关闭编辑器页面,软件提示是否保存?将文件保存为XXX.m
(4)注意将differential.m和XXX.m放在同一个文件夹,例如都放在”Matlab学习“
(5)将你打算调用的M文件所在的目录选为当前工作目录。(我这里Untitled4.m为主程序文件,即XXX.m)
QQ图片20141229180905.jpg
(6)右击文件XXX.m,点击”运行“,结果就出来了
最后再次感谢walqy222的帮助!
已有 2 人评分论坛币 学术水平 热心指数 信用等级 收起 理由
mermaid1258 + 1 + 1 + 1 + 1 精彩帖子
statax + 10 + 1 精彩帖子

总评分: 论坛币 + 11  学术水平 + 1  热心指数 + 2  信用等级 + 1   查看全部评分

使用道具

10
Xaero 发表于 2015-1-1 20:48:59 |只看作者 |坛友微信交流群
“但是没学过Matlab,看不懂,运行之后总不能出现结果”
如果是这种情况,建议还是先学习一些基础知识吧。不然大家不知怎么帮忙啊。

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-4-19 22:43