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

熊宇乐 y=zeros(1,1000); a=binornd(1,0.4,1,1000);b=binornd(1,0.4,1,1000); c=0;d=0; for i=1:1000 c=c+a(1,i).*b(1,i); y(i)=c/i; end y=y; num=1:1000; plot(num,y) liti2(1,0.4,100) liti2(1,0.4,10000) function proguji=liti5(mm) %mm 是随机实验次数 frq=0; randnum1=unifrnd(0,60,mm,1); randnum2=unifrnd(0,60,mm,1); randnum=randnum1-randnum2; proguji=0; for ii=1:mm if abs(randnum(ii,1))=20 frq=frq+1; end end proguji=frq/mm 练习 频率的稳定性 1 事件的频率 R = binornd(N,P,mm,nn) 在一组不变的条件下，重复作n次试验，记m是n次试验中事件A发生的次数。 频率 f=m/n 2.频率的稳定性 练习掷一枚不均匀硬币，正面出现概率为0.3，记录前1000次掷硬币试验中正面频率的波动情况，并画图。 在Matlab命令行中输入以下命令： liti1(1,0.3,1000) 例2 掷两枚不均匀硬币，每枚正面出现概率为0.4，记录前1000次掷硬币试验中两枚都为正面频率的波动情况，并画图。 在Matlab中编辑.m文件输入以下命令： function liti2(n,p,mm) pro=zeros(1,mm); randnum = binornd(n,p,2,mm)；a=0; for i=1:mm a=a+randnum(1,i)*randnum(2,i); pro(i)=a/i; end pro=pro，num=1:mm;plot(num,pro) 孟亚 function bino (n,p,m) x=binornd(n,p,1,m); y=binornd(n,p,1,m); for i=1:m if x(i)==1y(i)==1 s(i)=1; else s(i)=0; end end for i=1:m y(i)=sum(s(1,1:i))/i; end plot(y) 在一袋中有10 个相同的球，分别标有号码1,2,…,10。每次任取一个球，记录其号码后放回袋中，再任取下一个。这种取法叫做“有放回抽取”。今有放回抽取3个球，求这3个球的号码均为偶数的概率。（用频率估计概率） 例3: 解：令A={有放回抽取3个球，求这3个球的号码均为偶数}={(2,2,2)，(2,2,4)，….，(10,10,10)} function proguji=liti3(n,mm) frq=0; randnum=unidrnd(n,mm,3);proguji=0; for i=1:mm a=(randnum(i,1)+1)*(randnum(i,2)+1)*(randnum(i,3)+1)； if mod(a,2)==1 frq=frq+1 end end; proguji=frq/mm 例4 两盒火柴，每盒20根。每次随机在任一盒中取出一根火柴。问其中一盒中火柴被取完而另一盒中至少还有5根火柴的概率有多大？（用频率估计概率） liti4(20,5,100) proguji = 0.4800 liti4(20,5,1000) proguji = 0.4970 liti4(20,5,10000) proguji = 0.4910 liti4(20,5,100000) proguji = 0.4984 function proguji=liti4(nn,num,mm) %nn 是每盒中的火柴数 %num 是剩余的火柴数 %mm 是随机实验次数 frq=0; randnum=binornd(1,0.5,mm,2*nn);proguji=0; for i=1:mm a1=0;a2=0;j=1; while (a120)(a220) if randnum(i,j)==1 a1=a1+1; else a2=a2+1; end j=j+1; end if abs(a1-a2)=5 f