5、线相交,当N充分大时,可以认为显然,实验次数N越多,p的近似程度越好。需要指出的是,上述由投针试验求得π的近似值的方法,是进行真正的试验,并统计试验结果,要使获得的频率值与概率值偏差小,就要进行大量的试验,这在实际中,往往难以做到。所以,在现代计算机技术出现之前,用频率近似概率的方法—抑或称为雏形时代的蒙特卡罗方法—并没有得到实质上的应用。2.2蒙特卡罗方法与蒲丰问题21/21随着计算机和计算机技术的迅速发展,可以非常方便地利用计算机模拟随机实验。用数值模拟方法代替上述真正的投针实验,是利用均匀分布于(0,1)之间的随机数序列,并构造出随机投针的数学模型,然后进行大量的随机统计并求得到π的近
6、似值。如图1建立坐标系,平面上一根针的位置可以用针中心Ml的坐标x和针与平行线的夹角θ来决定,在y方向上的位置不影响相交性质,任意投针,意味着x与θ都是任意取的.但θ的范围可限于[0,π],x的范围可限于[0,a],在这种情况下,针与平行线相交的数学条件是其次,怎样模拟投针呢?亦即如何产生任意的[x,θ],x在[0,a]任意取值,意味着x在[0,a]上取哪一点的概率都一样,即x的概率密度函数为21/21由此,产生任意(x,θ)的过程就变为由f1(x)抽样x,由f2(θ)抽样θ的过程,容易得到式中,ξ1,ξ2均为(0,1)上均匀分布的随机数,只要随机数的均匀性和独立性良好,如此构造的数值模型就
7、很好地模拟了实际试验中的一次投针,并用下式判断是否相交且记录统计结果:是相交几率p的估计值,这样就实现了用数值方法模拟真正投针实验。3 蒙特卡罗方法的基本思路和特点用蒙特卡罗方法求解问题时,应建立一个概率模型,使待解问题与此概率模型相联系,然后通过随机试验求得某些统计特征值作为待解问题的近似解,与此相似,在一些物理问题,如核裂变、直流气体放电等过程中,粒子的输运过程及粒子输运总效应,也是可以与某些概率过程联系