赌徒输光问题的理论分析与模拟.doc

模拟赌徒输光问题的输光概率和博弈持续时间的程序代码：ticp=0.3;a=10;z=5;l=0;m=0;generation=10000;num=0;for i=1:generation flag=1; z=5; while(flag) r = rand(); if r>1-p z=z+1; else z=z-1; end if z==0 l=l+1; flag=0; end if z==a m=m+1; flag=0; end endendwin=m/(m+l)lost=l/(m+l)toc求不同下的输光概率clearp=0.1:0.1:1;a=10;z=5; generation=10000;for i=1:length(p)r=0;for j=1:generation m(j)=0;l(j)=0;flag=1;z=5;while(flag)r=rand(); if r=a flag=0; end m(j)=m(j)+1;endflag=1;u=10;endk(i)=mean(m);endplot(b,k);hold onx1=ones(size(b));x2=(1./b);x=[x1,x2];y=k;[b1,bint,r,rint,stats]=regress(y,x,0.05);plot(b,x*b1,r); legend(模拟的情况,拟合的情况);title(赌注与赌博持续时间的关系);xlabel(赌注的变化);ylabel(赌博的持续时间);。