MonteCarlo方法及其简单应用(图文)

1、MonteCarlo方法及其简单应用(图文)论文导读：本文介绍了Monte Carlo方法的思想，主要从在定积分计算方面介绍了随机投点法和平均值法，并将其推广到二重积分、三重积分和多重积分情形，最后以棋手分奖金问题介绍了 Monte Carlo方法在古典概率问题中的应用.分析了误差，介绍了减少误差的方法. 给出这些方法的实例及其Mathematica实现程序.关键词：MonteCarlo方法，积分计算，古典概率，模拟1 引言Monte Carlo方法，源于第二次世界大战美国关于研制原子弹的曼哈顿方案;.该方案的主持人之一、数学家冯-诺伊曼用著名世界的赌城摩纳哥的Monte Carlo来命名这种

2、方法,为它蒙上了一层神秘色彩.Monte Carlo方法研究的问题大致可分为两种类型:一种是问题本身就是随机的,另一种本身属于确定性问题,但可以建立它的解与特定随机变量或随机过程的数字特征或分布函数之间的联系,因而也可用随机模拟方法解决.文【1】-【7】 介绍了Monte Carlo方法的思想，但没有给出具体的实例及实现过程。发表论文。本文介绍了MonteCarlo方法的思想，从计算定积分和古典概率两方面的应用进行研究，给出了实例及其Mathematica实现程序.2 Monte Carlo方法2.1 Monte Carlo方法思想概述Monte Carlo方法，有时也称随机模拟(Random

3、Simulation)方法或统计试验(Statistical Testing)方法.它的根本思想是:首先建立一个概率模型或随机过程,使它的参数等于问题的解；然后通过对模型或过程的观察、抽样来计算所求参数的统计特征；最后给出所求解的近似值,而解的精度可用估计值的标准误差来表示.假设所求的量是随机变量的数学期望,那么近似确定的方法是对进行重复抽样,产生相互独立的值的序列并计算其算术平均值:根据大数定理，当充分大时,以概率1成立，即可用作为的估计值.Monte Carlo方法以概率统计理论为根底,以随机抽样(随机变量的抽样)为手段，在很多方面有重要的应用.它的优点表现在三个方面：方法和程序的结构简单

4、，易分析、易理解；收敛的概率性和收敛速度与问题的维数无关，很好的防止了维数问题；受问题条件限制的影响较小，很好的提高可行性.使用Monte Carlo方法的步骤如下：(l)构造或描述概率过程(2)实现从概率分布中抽样(3)建立各种估计量2.2 Monte Carlo方法的可行性从Monte Carlo方法的根本思想可以得到它通常的做法，利用数学或物理方法产生中均匀分布的随机数，在变换得到任意分布的随机数.随机数个数很大时，可以由大数定理，求出事件的概率值.这种做法的可行性主要依据下面的事实：即服从上的均匀分布.因此是服从分布函数的随机变量.所以我们只要能够产生中均匀分布的随机变量的子样，那么通

5、过2式我们就可以得到任意分布函数的随机变量的子样.再结合大数定理、就可以运用Monte Carlo方法进行随机模拟，解决一些实际的问题.3 Monte Carlo方法在定积分中的应用3.1随机投点法对于定积分.为使计算机模拟简单起见，设，有限，令，并设是在上均匀分布的二维随机变量，其联合密度函数为.那么是中曲线下方的面积如图2.图2假设我们向中进行随机投点.假设点落在下方即称为中的，否那么称不中.那么点中的概率为，假设我们进行次投点，其中次中的.那么可以得到的一个估计该方法的具体计算步骤为： 独立地产生2个随机数,i=1,n； 计算， ，和； 统计的个数； 用由矩法，假设有个那么可给出的一个矩

6、估计，这便是样本平均值法的根本原理.假设，有限，可取.设是那么的一个估计为该方法的具体计算步骤为： 独立地产生个随机数； 计算和 ，；用其中为S维单位立方体，在上有：.很明显.此时积分对于这种较为一般形式的多重积分计算问题，采用的还是随机投点.具体步骤如下：首先产生个随机数i=1,2,，及，构造维随机向量 ，然后检验是否落后在V中，同理可以推论.检验是否成立，如果在构成的个随机向量中，有个随机向量落于V中，那么取作为积分的近似值，即，如果积分区域及被积函数不满足上述条件，那么可以通过变换便可到达所希望的条件.方法二：其中积分区域包含在维多面体中，此多面体决定于个不等式.设函数在内连续且满足条件

7、：，是在维多面体中均匀分布的随机质点的个数，是在个随机点之中落入以维区域V为底以为顶之曲顶柱体内的随机点的个数.这里表示由不等式和决定的维多面体.那么重积分的Monte Carlo近似计算公式为：=例2 在三维空间中，由三个圆柱面：，围成一个立体，利用Monte Carlo方法求它的体积.分析：据题意，所求体积,其中，且，.记，,考虑在空间内随机的产生个点，落在空间内有个，那么.在Mathematica中模拟程序见附录2.5 在古典概率问题中的应用下面的例子说明了Monte Carlo方法在古典概率中的应用.例3 甲乙两位棋手棋艺相当，现他们在一项奖金为1000元的比赛中相遇，比赛为五局三胜制

8、，已经进行了三局的比赛，结果为甲三胜一负，现因故要停止比赛，问应该如何分配这1000元比赛奖金才算公平？分析：平均分对甲欠公平，全归甲那么对乙欠公平.合理的分法是按一定的比例分配.现在我们用计算机模拟两位棋手后面的比赛，是否就可以知道奖金分配方案.由于两位棋手的棋艺相当，可以假定他们在以下每局的比赛胜负的时机各半.Mathematica中函数产生随机数0或1，0与1出现的时机各占一半，可以用随机数1表示甲棋手胜，而随机数0表示乙胜.最终以甲分到；乙分到.即甲750元，乙250元.6 误差分析6.1 收敛性蒙特卡罗方法是由随机变量的简单子样的算术平均值：作为所求解的近似值.由大数定律可知，如独立

9、同分布，且具有有限期望值6.2 误差蒙特卡罗方法的近似值与真值的误差问题，概率论的中心极限定理给出了答案.该定理指出，如果随机变量序列，独立同分布，且具有有限非零的方差，即是的分布密度函数.那么当N充分大时，有如下的近似式其中称为置信度，1称为置信水平.这说明，不等式近似地以概率1成立，且误差收敛速度的阶为.通常，Monte Carlo方法的误差定义为上式中与置信度是一一对应的，根据问题的要求确定出置信水平后，查标准正态分布表，就可以确定出.关于蒙特卡罗方法的误差需说明两点：第一，蒙特卡罗方法的误差为概率误差，这与其他数值计算方法是有区别的.第二，误差中的均方差是未知的，必须使用其估计值来代替

10、，在计算所求量的同时，可计算出.例4 求用平均值法估计圆周率，并考虑置信度为5%，精度要求为0.01的情况下所需的试验次数.解：易知，故考虑令，令，其期望值为，因此=，其中是区间上均匀分布的随机数.此时，所以次.6.3 减小方差的各种技巧显然，当给定置信度后，误差由和N决定.要减小，或者是增大N，或者是减小方差.在固定的情况下，要把精度提高一个数量级，试验次数N需增加两个数量级.因此，单纯增大N不是一个有效的方法.另一方面，如能减小估计的均方差，比方降低一半，那误差就减小一半，这相当于N增大四倍的效果.因此降低方差的各种技巧，引起了人们的普遍注意.总的来说，增大样本的值对计算机要求较高；减小方

11、差的技巧都只具有指导思想上的意义.对于实际的计算问题，往往要求对涉及的随机变量有先验的了解，或者对发生的物理过程的性态有一定的认识.通过利用这些预知的信息采取相应的手段减小误差，提高精度.附录1.1 n=1000;p=Do;y=Random;If,k,1,n;AppendTo,t,1,10;Print;Sum,t,1,10/102 n=10000;p=Do;y=Random;If,k,1,n;AppendTo,t,1,10;Print;Sum,t,1,10/103 n=100000;p=Do;y=Random;If,k,1,n;AppendTo,t,1,10;Print;Sum,t,1,10/102. n=1000;p=Do;y=Random;z=Random;If,k,1,n;AppendTo,t,1,10;Print;Sum,t,1,10/103. n=1000;p=Do+2;y=Random+1;If,k,1,n;AppendTo,t,1,20Print;Sum,t,1,20/20,1000-Sum,t1,20/20参考文献【1】 徐钟济蒙特卡罗方法上海：上海科学技术出版社，1985：171-188.【2】 茆诗松，王静龙，濮晓龙高等数理统计.北京：高等教育出版社，2006：415454【3】 周铁，徐树方