原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
博弈论赌博仿真
%PROJECT1%%function project1(N,count,mode,pb)N=100;count=500;mode=2;pb=1;% 随即过程project1仿真程序% N 人数% count 次数% mode 1-零和 2-囚徒% pb 问题序号% 26-10-2008if nargin~=4;end%if nargin2|isempty(p), p=1/length(x)*ones(size(x)); end;%根据N产生概率空间P=rand(1,N);Y=zeros(N,count);%单词赌博尽收益BOSS=zeros(N/2,count);%庄家在每次赌博中的收益%**********求得所有玩家的策略均值,以便后续策略调整*******%eP=0;for i=1:N eP=eP+P(i);endeP=eP/N;%*****************赌博过程仿真******************%for j=1:count suffixn=zeros(1,N); suffixn(1)=ceil(N.*rand(1)); %产生随机下标,模拟抽签for i=2:N suffixn(i)=ceil(N.*rand(1)); %产生随机下标,模拟抽签 k=1; while(ki) %检测重复,若果重复,则重新产生 if (suffixn(i)==suffixn(k)) suffixn(i)=ceil(N.*rand(1)); k=1; else k=k+1; end end %对当前i抽签结束 end %让下标近邻的两个人赌博 1-2 3-4 ... if(pb==1) for i=1:N/2 %抽签决定座位 R:-1 L:1, 小于0.5,坐右边R pos=round(rand(1)); if(pos==0) pos=-1; else pos=1; end %产生一个随机数,如果小于p,则出红桃,否则,出黒桃 pp=rand(1); if(ppP(suffixn(i))) %出红桃 x1=-1; else x1=1; end pp=rand(1); if(ppP(suffixn(i+N/2))) %出红桃 x2=-1; else x2=1; end if(mode==1) Y(suffixn(i),j)=pos*x1*x2; Y(suffixn(i+N/2),j)=pos*(-1)*x1*x2; end if(mode==2) if(x1==-1 x2==-1) Y(suffixn(i),j)=1; Y(suffixn(i+N/2),j)=1; BOSS(i,j)=-2; elseif(x1==-1 x2==1) Y(suffixn(i),j)=pos*(-1)*4; Y(suffixn(i+N/2),j)=pos*4; BOSS(i,j)=0; elseif(x1==1 x2==-1) Y(suffixn(i),j)=pos*4; Y(suffixn(i+N/2),j)
文档评论（0）