数学建模-蒙特卡罗方法.ppt

《数学建模-蒙特卡罗方法.ppt》由会员上传分享，免费在线阅读，更多相关内容在PPT专区-天天文库。

1、2017美赛培训：蒙特卡罗方法主讲人：老教练报告人：陈蓉蓉陈雨路杨园园什么叫蒙特卡罗方法？蒙特卡罗方法又称统计模拟法、随机抽样技术，是一种随机模拟方法，以概率和统计理论方法为基础的一种计算方法，是使用随机数（或伪随机数）来解决很多计算问题的方法。将所求解的问题同一定的概率模型相联系，用电子计算机实现统计模拟或抽样，以获得问题的近似解。为象征性地表明这一方法的概率统计特征，故借用赌城蒙特卡罗命名。基本思想当所求问题的解是某个事件的概率，或者是某个随机变量的数学期望，或者是与概率，数学期望有关的量时，通过某种试验的方法，得出该事件发生

2、的概率，或者该随机变量若干个具体观察值的算术平均值，通过它得到问题的解。当随机变量的取值仅为1或0时，它的数学期望就是某个事件的概率。或者说，某种事件的概率也是随机变量（仅取值为1或0）的数学期望。蒙特卡罗方法的特点优点：1、能够比较逼真地描述具有随机性质的事物的特点及物理实验过程2、受几何条件限制小3、收敛速度与问题的维数无关4、具有同时计算多个方案与多个未知量的能力5、误差容易确定6、程序结构简单，易于实现缺点：1收敛速度慢2误差具有概率性3在粒子输运问题中，计算结果与系统大小有关所以在使用蒙特卡罗方法时，要“扬长避短”，只对

3、问题中难以用解析（或数值）方法处理的部分，使用蒙特卡罗方法计算，对那些能用解析（或数值）方法处理的部分，应当尽量使用解析方法主要应用范围粒子输运问题（实验物理，反应堆物理，高能物理）统计物理典型数学问题真空技术激光技术以及医学生物探矿什么是随机数？在连续型随机变量的分布中，最简单而且最基本的分布是单位均匀分布。由该分布抽取的简单子样称为随机数序列，其中每一个体称为随机数符号：两个特点：独立性，均匀性产生随机数随机数表方法物理方法随机数表随机数表是由0,1,2,3,4,5,6,7,8,9十个数字组成，每个数字以0.1的等概率出现，数

4、字之间相互独立，这些数字序列叫作随机数字序列。（如果要得到n位有效数字的随机数，只需将表中每n个相邻的随机数字合并在一起，且在最高位的前边加上小数点即可。例如，某随机数表的第一行数字为763425891...，要想得到三位有效数字的随机数一次为0.763，0.425，0.891...）物理方法利用某些物理现象，在计算机上增加些特殊设备，可以在计算机上直接产生随机数。作为随机数发生器的物理源主要有两种：一种是根据放射性物质的放射性，另一种是利用计算机的固有噪声。一般情况下，任意一个随机数在计算机内总是用二进制的数表示的：其中或者为0

5、，或者为1。因此，利用物理方法在计算机产生随机数，就是要产生只取0或1的随机数字序列，数字之间相互独立，每个数字取0或1的概率均为0.5缺点随机数表需在计算机中占有很大内存，而且也难以满足蒙特卡罗方法对随机数需要量非常大的要求，因此，该方法不适于在计算机上使用。物理方法产生的随机数序列无法重复实现，不能进行程序复算。给验证结果带来很大困难。而且增加随机数发生器和电路联接等附加设备，费用昂贵。因此该方法也不适合在计算机上使用。伪随机数用递推公式产生随机数序列。伪随机数存在的两个问题递推公式和初始值确定后，整个随机数序列便被唯一确定。

6、不满足随机数相互独立的要求。由于随机数序列是由递推公式确定的，而在计算机上所能表示的[0，1]上的数又是有限的，因此，这种方法产生的随机数序列就不可能不出现重复。随机数序列出现周期性的循环现象。解决方案第一个问题：不能从本质上加以改变，但只要递推公式选的比较好，随机数间的相互独立性是可以近似满足的。第二个问题：因为用蒙特卡罗方法解任何具体问题时，所使用的随机数的个数总是有限的，只要所用随机数的个数不超过伪随机数序列出现循环现象时的长度就可以了。应用：蒙特卡罗方法计算积分可以通俗地说，蒙特卡罗方法是用随机试验的方法计算积分，即将所要

7、计算的积分看作服从某种分布密度函数f(r)的随机变量ｇ(r)的数学期望通过某种试验，得到Ｎ观察值r1，r2，…，rN（用概率语言来说，从分布密度函数f(r)中抽取Ｎ个子样r1，r2，…，rN，），将相应的Ｎ个随机变量的值g(r1)，g(r2)，…，g(rN)的算术平均值作为积分的估计值（近似值）。求积分（2.1）蒙特卡罗方法步骤如下：1、在区间【a,b】上利用计算机均匀产生n个随机数x1，x2·····xn，这个可以在MATLAB软件中用unifrnd命令实现。2、计算每一个随机数相应的被积函数值f（x1），f（x2）····f（

8、xn）。3、计算被积函数值的平均值4、所以2.1式的值≈简单定积分例子：用蒙特卡罗方法求首先我们进行严格的数学计算，便于后面与蒙特卡洛计算方法所得结果形成对比：已知的原函数是，那么定积分值就是：我们可以在Matlab中输入以下代码进行精确计算：ex