浅析蒙特卡洛方法原理及应用.docx

浅析蒙特卡洛方法原理及应用于希明（英才学院1236103班测控技术与仪器专业6120110304）摘要：本文概述了蒙特卡洛方法产生的历史及基本原理，介绍了蒙特卡洛方法的最初应用一 一蒲丰投针问题求圆周率，并介绍了蒙特卡洛方法在数学及生活中的一些简单应用，最后总 结了蒙特卡洛方法的特点关键词：蒙特卡洛方法蒲丰投针生活应用蒙特卡洛方法（Monte Carlo method），也称统计模拟方法，是二十世纪四 十年代中期由于科学技术的发展和电子计算机的发明，而被提出的一种以概率统 计理论为指导的一类非常重要的数值计算方法它是以概率统计理论为基础，依 据大数定律（样本均值代替总体均值），利用电子计算机数字模拟技术，解决 一些很难直接用数学运算求解或用其他方法不能解决的复杂问题的一种近似计 算法蒙特卡洛方法在金融工程学，宏观经济学，计算物理学（如粒子输运计算、 量子热力学计算、空气动力学计算）等领域应用广泛一蒙特卡洛方法的产生及原理、蒙特卡洛方法于20世纪40年代美国在第二次世界大战中研制原子弹的“曼 哈顿计划”计划的成员S.M.乌拉姆和J.冯•诺伊曼首先提出数学家冯•诺伊 曼用驰名世界的赌城一摩纳哥的Monte Carlo一来命名这种方法，为它蒙上了一 层神秘色彩。

在这之前，蒙特卡洛方法就已经存在1777年，法国数学家蒲丰 （Georges Louis Leclere de Buffon，1707一1788）提出用投针实验的方法求 圆周率n这被认为是蒙特卡洛方法的起源其基本原理如下：由概率定义知，某事件的概率可以用大量试验中该事件发 生的频率来估算，当样本容量足够大时，可以认为该事件的发生频率即为其概率 因此，可以先对影响其可靠度的随机变量进行大量的随机抽样，然后把这些抽样 值一组一组地代入功能函数式，确定结构是否失效，最后从中求得结构的失效概 率蒙特卡洛法正是基于此思路进行分析的设有统计独立的随机变量Xi（i=1，2，3，…，k），其对应的概率密度函数 分别为fx1，fx2，…，fxk，功能函数式为Z=g（x1，x2，…，xk）首先根据各 随机变量的相应分布，产生N组随机数x1，x2,…，xk值，计算功能函数值 Zi=g（x1，x2，…，xk）（i=1，2，…，N），若其中有L组随机数对应的功能函数 值ZiW0，则当N—8时，根据伯努利大数定理及正态随机变量的特性有：结构 失效概率，可靠指标蒲丰投针问题 •-、作为蒙特卡洛方法的最初应用，是解决蒲丰投针问题。

1777年，法国数学 家蒲丰提出利用投针实验求解圆周率的问题设平面上等距离（如为2a）画有一 些平行线，将一根长度为2l（l

我们定 义一个随机变量n为吁1 = r⑴.则显然有n的期望值等于积分值I只要抽取足够多的随机点，即取随机点数足够大时，咱勺平均值f (E )就是积分I的一个无偏估计值n的方差显然V{n }依赖于被积函数f(x)在积分域上的方差当f(x)在乂的定义域内变化 平坦，即和1的差处处都较小时，方差也小；反之，则方差较大欧)从这里可以看出：尽量减小被积函数在积分域上的方差，可以减小积分估计值的 方差，加速收敛推而广之来说，就是要减少模拟量在模拟范围内的方差根据这样的原则，当被积函数f(x)在积分域内的方差较大时，可以采用各 种抽样技巧如采用重要抽样法，将f(x )的方差吸收到g(x)中去，这样模拟量一 记录函数f*(x)=f(x)/g(x)在定义域内相当平坦，则我们将积分式的计算变为pl fU _J = £溥加"J,翥号⑴成=I f仙3）虫若选取n'为服从分布密度函数g(x)的函数f*(x)的抽样值这里g(x)称为偏倚 分布密度函数我们得到1=聊｝ •因此它的平均值给出了i的一个无偏估计值这时的方差为:E=fl/我-妙血*=£膏］血5；在实际计算中，方差通过下式得到计算结果：式中角型括号表示对括号内所有可能的［0, 1］区间，按g(x)分布的随机坐标数序 列｛xj对应的数值求平均。

方程右边第一项对｛&(%)｝求平均 f )，第二项表2示求｛f*(x”平均值的平方(f忝)上式可以经推导得到:由此我们看出其误差平方与f*在［0, 1］区间的方差成正比，并且X 1/打这 与中心极限定理所得到的结果一致四、蒙特卡洛方法在学校分班问题中的应用学生分班是学校管理中的一项经常而重要的工作.每个新学期伊始,学校要 对新招取的新生进行分班.操作时需要综合考虑男女生比例、各科成绩、生源分 布、是否住校、学生干部安排等要素.假设分班的结果要求满足以下条件：A. 班级规模尽量相当，人数相差不超过1人；B. 各班学生总分的均分差值在2分内；C. 各班的各科平均分差值在3分内；D. 各班的住校人数相差最多3人；E. 入学前任过学生干部的人数相差不超2人.现要求条件A〜B必须满足,C〜E尽量满足(有时还需设定男女生人数之差的上 限等).采用蒙特卡罗算法来实现的过程为：a. 随机产生一个分班的方案使符合条件B ;b. 检验是否符合条件A ,若符合执行步骤c ,否则返回执行步骤a ;c. 计算条件C〜E的参数，作为目标函数的参考值.在此可以通过构造加权系数 来确定条件C ,D ,E的优先级，作为最终参考目标；d. 通过大量模拟得到多个目标值,选择其中最佳目标值作为最优方案.五、总结 蒙特卡洛方法的解题过程可以归结为三个主要步骤：构造或描述概率过程；实现 从已知概率分布抽样；建立各种估计量。

蒙特卡洛方法解题过程的三个主要步骤:（1） 构造或描述概率过程对于本身就具有随机性质的问题，如粒子输运问题，主要是正确描述和模拟这个 概率过程，对于本来不是随机性质的确定性问题，比如计算定积分，就必须事 先构造一个人为的概率过程，它的某些参量正好是所要求问题的解即要将不具 有随机性质的问题转化为随机性质的问题2） 实现从已知概率分布抽样构造了概率模型以后，由于各种概率模型都可以看作是由各种各样的概率分布构 成的，因此产生已知概率分布的随机变量（或随机向量），就成为实现蒙特卡洛 方法模拟实验的基本手段，这也是蒙特卡洛方法被称为随机抽样的原因最简单、 最基本、最重要的一个概率分布是（0，1）上的均匀分布（或称矩形分布）随 机数就是具有这种均匀分布的随机变量随机数序列就是具有这种分布的总体的 一个简单子样，也就是一个具有这种分布的相互独立的随机变数序列产生随机 数的问题，就是从这个分布的抽样问题在计算机上，可以用物理方法产生随机 数，但价格昂贵，不能重复，使用不便另一种方法是用数学递推公式产生这 样产生的序列，与真正的随机数序列不同，所以称为伪随机数，或伪随机数序列 不过，经过多种统计检验表明，它与真正的随机数，或随机数序列具有相近的性 质，因此可把它作为真正的随机数来使用。

由已知分布随机抽样有各种方法，与 从（0,1）上均匀分布抽样不同，这些方法都是借助于随机序列来实现的，也就是 说，都是以产生随机数为前提的由此可见，随机数是我们实现蒙特卡洛模拟的 基本工具3） 建立各种估计量一般说来，构造了概率模型并能从中抽样后，即实现模拟实验后，我们就要确定 一个随机变量，作为所要求的问题的解，我们称它为无偏估计建立各种估计量， 相当于对模拟实验的结果进行考察和登记，从中得到问题的解作为一种解决物理数学问题和系统性质分析的近似计算法，蒙特卡洛方法 和传统方法相比，具有思路新颖，直观性强，简便易行的优点特别是借助电子 计算机可以在很大程度上模拟许多大型的、难以实现的复杂实验或社会行为过程, 使蒙特卡洛方法逐渐成为重要的计算方法尤其在复杂系统的性能评价上给出量 化指标，利用蒙特卡洛方法对系统进行模拟几乎是必不可少的参考文献[1] 吴海霞 刘潞锋 蒙特卡罗方法在实际问题中的应用 太原师范学院学报（自然科学版）2009.3 Vol.8 No.3[2] 杨莉军赵贤淑蒙特卡洛方法及在二维随机游动问题中的应用初探北京印刷学院学报2001.9 Vol.9 No.3[3] 柳海东蒙特卡洛方法在概率计算中的应用苏州职业大学学报2004.8 Vol.15 No.3[4] 百度百科。