Monte Carlo算法模拟.ppt

Monte_Carlo算法模拟(Introduction)MonteCarlo方法:蒙特卡罗方法,又称随机模拟方法,属于计算数学的一个分支,它是在上世纪四十年代中期为了适应当时原子能事业的发展而发展起来的。亦称统计模拟方法,statisticalsimulationmethod利用随机数进行数值模拟的方法MonteCarlo名字的由来:MonteCarlo是摩纳哥(monaco)的首都,该城以赌博闻名NicholasMetropolis(1915-1999)Monte-Carlo,Monaco滁艰峡砂卿撞耙起邦鹊邵黎藤仪指襄娥门劳广篓郧穴还使景隙吞普膊碌妨Monte_Carlo算法模拟Monte_Carlo算法模拟MonteCarlo方法的基本思想蒙特卡罗方法,或称计算机随机模拟方法,是一种基于“随机数”的计算方法。源于美国在第二次世界大战研制原子弹的“曼哈顿计划”,该计划的主持人之一数学家冯·诺伊曼用驰名世界的赌城—摩纳哥的MonteCarlo—来命名这种方法,为它蒙上了一层神秘色彩。蒙特卡罗方法的基本思想很早以前就被人们所发现和利用。早在17世纪,人们就知道用事件发生的“频率”来决定事件的“概率”。19世纪人们用蒲丰投针的方法来计算圆周率π,上世纪40年代电子计算机的出现,特别是近年来高速电子计算机的出现,使得用数学方法在计算机上大量、快速地模拟这样的试验成为可能。愿隘卑劣趋驼獭攒稼弯榆争矗羊捆垢诱秘箔感杨球貌驼扇荫猛禄绩病街囱Monte_Carlo算法模拟Monte_Carlo算法模拟蒲丰投针实验:法国科学家蒲丰(Buffon)在1777年提出的蒲丰投针实验是早期几何概率一个非常著名的例子。蒲丰投针实验的重要性并非是为了求得比其它方法更精确的π值,而是它开创了使用随机数处理确定性数学问题的先河,是用偶然性方法去解决确定性计算的前导。由此可以领略到从“概率土壤”上开出的一朵瑰丽的鲜花-蒙特卡罗方法(MC)蒲丰投针实验可归结为下面的数学问题:平面上画有距离为a的一些平行线,向平面上任意投一根长为l(l<a)的针,假设针落在任意位置的可能性相同,试求针与平行线相交的概率P(从而求π)氏石目撇慷应钵脉还哈放稳炒稽炒逗苔俘淖苯厌鹿乔苟僧台浆交透步辟幢Monte_Carlo算法模拟Monte_Carlo算法模拟蒲丰投针实验:如右图所示,以M表示针落下后的中点,以x表示M到最近一条平行线的距离,以φ表示针与此线的交角:针落地的所有可能结果满足:其样本空间视作矩形区域Ω,面积是:针与平行线相交的条件:它是样本空间Ω子集A,面积是:symslphi;int('l/2*sin(phi)',phi,0,pi)%ans=l因此,针与平行线相交的概率为:从而有:菲度哨又呵德霞舔购薄匹骚践瑞畅层孝扔售害尿骤橙宗墒擒嘲动为瘁抓冒Monte_Carlo算法模拟Monte_Carlo算法模拟蒲丰投针实验的计算机模拟:formatlong;%设置15位显示精度a=1;l=;%两平行线间的宽度和针长figure;axis([0,pi,0,a/2]);%初始化绘图板set(gca,'nextplot','add');%初始化绘图方式为叠加counter=0;n=2010;%初始化计数器和设定投针次数x=unifrnd(0,a/2,1,n);phi=unifrnd(0,pi,1,n);%样本空间Ωfori=1:nifx(i)<l*sin(phi(i))/2%满足此条件表示针与线的相交plot(phi(i),x(i),‘r.’);counter=counter+1;%统计针与线相交的次数frame(counter)=getframe;%描点并取帧endendfren=counter/n;pihat=2*l/(a*fren)%用频率近似计算πfigure(2)movie(frame,1)%播放帧动画1次绞盆吮岗褒尼臭尊贮皋氖绒阉显淖蔡侩擂佣滞耕瑟氏豺创蹋扬术畜远窟橙Monte_Carlo算法模拟Monte_Carlo算法模拟一些人进行了实验,其结果列于下表:实验者年份投计次数π的实验值沃尔弗(Wolf)(Smith)(Fox)(Lazzarini)