猎狗追兔问题.doc

文档介绍：蒲丰投针――MonteCarlo算法背景:蒙特卡罗方法(MonteCarlo),也称统计模拟方法,是在二次世界大战期间随着科学技术的发展和电子计算机的发明,而被提出的一种以概率统计理论为基础的一类非常重要的数值计算方法。蒙特卡罗方法在应用物理、原子能、固体物理、化学、生态学、社会学以及经济行为等领域中得到广泛利用。蒙特卡罗方法的名字来源于世界著名的赌城——摩纳哥的蒙特卡罗。其历史起源可追溯到1777年法国科学家蒲丰提出的一种计算圆周的方法——随机投针法,即著名的蒲丰投针问题。问题:设在平面上有一组平行线,间距为d,把一根长L的针随机投上去,则这根针和平行线相交的概率是多少?(其中L<d)分析:由于L<d,所以这根针至多只能与一条平行线相交。设针的中点与最近的平行线之间的距离为y,针与平行线的夹角为q(0£q£p)。相交情形不相交情形易知针与平行线相交的充要条件是:由于,且它们的取值均满足平均分布。建立直角坐标系,则针与平行线的相交条件在坐标系下就是曲线所围成的曲边梯形区域(见右图)。所以有几何概率可知针与平行线相交的概率是 MonteCarlo方法:随机产生满足平均分布的y和q,其中,判断y是否在曲边梯形内。重复上述试验,并统计y在曲边梯形内的次数m,其与试验次数n的比值即为针与平行线相交的概率的近似值。clear;n=100000;L=1;d=2;m=0;fork=1:ntheta=rand(1)*pi;y=rand(1)*d/2;ify<sin(theta)*L/2m=m+1;endendfprintf('针与平行线相交的概率大约为%f\n',m/n)计算π的近似值利用该方法可以计算p的近似值:下面是一些通过蒲丰投针实验计算出来的p的近似值: **********.,而是在于它是第一个用几何形式表达概率问题的例子。计算π的这一方法,不但因其新颖,奇妙而让人叫绝,而且它开创了使用随机数处理确定性数学问题的先河,是用偶然性方法去解决确定性计算的前导。