第3讲蒙特卡洛方法基本思想.pptVIP

functionliti1(n,p,mm)pro=zeros(1,mm);randnum=binornd(n,p,1,mm)a=0;fori=1:mma=a+randnum(1,i);pro(i)=a/i;endpro=pronum=1:mm;plot(num,pro)在Matlab中编辑.m文件输入以下命令：第30页，共61页，星期日，2025年，2月5日在Matlab命令行中输入以下命令：liti1(1,0.5,1000)第31页，共61页，星期日，2025年，2月5日在Matlab命令行中输入以下命令：liti1(1,0.5,10000)第32页，共61页，星期日，2025年，2月5日练习频率的稳定性1事件的频率R=binornd(N,P,mm,nn)在一组不变的条件下，重复作n次试验，记m是n次试验中事件A发生的次数。频率f=m/n2.频率的稳定性练习掷一枚不均匀硬币，正面出现概率为0.3，记录前1000次掷硬币试验中正面频率的波动情况，并画图。第33页，共61页，星期日，2025年，2月5日在Matlab命令行中输入以下命令：liti1(1,0.3,1000)第34页，共61页，星期日，2025年，2月5日例2掷两枚不均匀硬币，每枚正面出现概率为0.4，记录前1000次掷硬币试验中两枚都为正面频率的波动情况，并画图。在Matlab中编辑.m文件输入以下命令：functionliti2(n,p,mm)pro=zeros(1,mm);randnum=binornd(n,p,2,mm)；a=0;fori=1:mma=a+randnum(1,i)*randnum(2,i);pro(i)=a/i;endpro=pro，num=1:mm;plot(num,pro)第35页，共61页，星期日，2025年，2月5日熊宇乐y=zeros(1,1000);a=binornd(1,0.4,1,1000);b=binornd(1,0.4,1,1000);c=0;d=0;fori=1:1000c=c+a(1,i).*b(1,i);y(i)=c/i;endy=y;num=1:1000;plot(num,y)第36页，共61页，星期日，2025年，2月5日孟亚functionbino(n,p,m)x=binornd(n,p,1,m);y=binornd(n,p,1,m);fori=1:mifx(i)==1y(i)==1s(i)=1;elses(i)=0;endendfori=1:my(i)=sum(s(1,1:i))/i;endplot(y)第37页，共61页，星期日，2025年，2月5日liti2(1,0.4,100)第38页，共61页，星期日，2025年，2月5日liti2(1,0.4,10000)第39页，共61页，星期日，2025年，2月5日在一袋中有10个相同的球，分别标有号码1,2,…,10。每次任取一个球，记录其号码后不放回袋中，再任取下一个。这种取法叫做“不放回抽取”。今不放回抽取3个球，求这3个球的号码均为偶数的概率。（用频率估计概率）例3:解：令A={不放回抽取3个球，求这3个球的号码均为偶数}={(2,4,6)，(2,4,8)，….，(6,8,10)}第40页，共61页，星期日，2025年，2月5日functionproguji=liti3(nn,num,mm)%nn是每盒中的火柴数%num是剩余的火柴数%mm是随机实验次数frq=0;randnum=binornd(1,0.5,mm,2*nn);proguji=0;fori=1:mma1=0;a2=0;j=1;while(a120)(a220)ifrandnum(i,j)==1a1=a1+1;elsea2=a2+1;endj=j+1;endifabs(a1-a2)=5frq=frq+1;end%a1=a1,a2=a2,frq%pauseendproguji=frq/mm第41页，共61页，星期日，2025年，2月5日functionproguji=liti3(nn,num,m