蒙特卡罗方法介绍及其建模应用Part_II.pptVIP

Monte-Carlo 方法介绍及其建模应用 朱连华 Tel南京信息工程大学数学与统计学院 E-mail:ahualian@126.com 课程说明 公用邮箱：ahualian2008@126.com key：ahualian2008 参考书目: 黄燕、吴平.SAS统计分析及应用，机械工业出版社. 陈杰. Matlab宝典，电子工业出版社. 张文彤等.SPSS11.0统计分析教程，北京希望电子出版社. 薛益、陈立萍.统计建模与R软件，清华大学出版社. 主要内容 蒙特卡洛方法应用实例 频率的稳定性模拟 频率: 在一组不变的条件下，重复作n次试验，记m是n次试验中事件A发生的次数，频率 f=m/n 频率的稳定性： f----P 例1：掷一枚均匀硬币，记录掷硬币试验中频率P的波动情况 function liti21(p,mm) pro=zeros(1,mm); randnum = binornd(1,p,1,mm) a=0; for i=1:mm a=a+randnum(1,i); pro(i)=a/i; end pro=pro num=1:mm; plot(num,pro) 古典概率模拟 例3: 在一袋中有10 个相同的球，分别标有号码1,2,…,10。每次任取一个球，记录其号码后放回袋中，再任取下一个。这种取法叫做“有放回抽取”。今有放回抽取3个球，求这3个球的号码均为偶数的概率。（用频率估计概率） 解：有放回取3个球, 所有取法有103种; 有放回取3个偶数号码的球, 所有取法有53种. 所以 古典概率模拟 例4: 两盒火柴，每盒20根。每次随机在任一盒中取出一根火柴。问其中一盒中火柴被取完而另一盒中至少还有5根火柴的概率有多大？（用频率估计概率） function proguji=liti24_0(mm) %mm 是随机实验次数 frq=0; randnum=binornd(1,0.5,mm,2*20); 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 frq=frq+1; end end; proguji=frq/mm 古典概率模拟 例4: 两盒火柴，每盒n根。每次随机在任一盒中取出一根火柴。问其中一盒中火柴被取完而另一盒中至少还有k根火柴的概率有多大？（用频率估计概率） function proguji=liti24_1(n,k,mm) %n是每盒中的火柴数 %k 是剩余的火柴数 %mm 是随机实验次数 frq=0;randnum=binornd(1,0.5,mm,2*n);proguji=0; for i=1:mm a1=0;a2=0;j=1; while (a1n)(a2n) if randnum(i,j)==1 a1=a1+1; else a2=a2+1; end j=j+1; end if abs(a1-a2)=k , frq=frq+1; end % a1=a1,a2=a2,frq % pause end; proguji=frq/mm 几何概率模拟 1.定义 向任一可度量区域G内投一点，如果所投的点落在G中任意可度量区域g内的可能性与g的度量成正比，而与g的位置和形状无关，则称这个随机试验为几何型随机试验。或简称为几何概型。 2.概率计算 P(A)=[A的度量]/[S的度量] 例5: 两人约定于12点到1点到某地会面，先到者等20分钟后离去，试求两人能会面的概率？ 解：设x, y分别为甲、乙到达时刻(分钟) 令A={两人能会面}={(x,y)||x-y|≤20，x≤60,y≤60} P(A)=A的面积/S的面积=（602-402）/602=5/9=0.5556 几何概率模拟 function proguji=liti25(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