蒙特卡罗方法及应用实验讲义2016资料.pdf

蒙特卡罗方法及应用 实验讲义 东华理工大学核工系 2016.8 1 实验一 蒙特卡罗方法基本思想 一、实验目的 1、了解蒙特卡罗方法方法的基本思想； 2、掌握蒙特卡罗方法计算面积、体积的方法； 3、掌握由已知分布的随机抽样方法。 二、实验原理 Monte Carlo 方法，又称统计模拟方法或计算机随机模拟方法，是一种基于 “随机数”进行数值模拟的方法，一种采用统计抽样理论近似求解物理或数学问 题的方法。 如待求量可以表述成某些特征量的期望值、某些事件出现的概率或两者的函 数形式，那么可采用蒙特卡罗方法求解。在求解某些特征量的期望值或某些事件 出现的概率时，必须构建合符实际的数学模型。例如采用蒙特卡罗方法计算某函 数所围面积时，构建的数学模型是构造一已知面积的可均匀抽样区域，在该区域 投点，由伯努利定理大数定理可知，进入待求区域投点的频率依概率 1 收敛于该 事件出现的概率（面积之比）。 由已知分布的随机抽样方法指的是由已知分布的总体中抽取简单子样。具体 方法很多，详见教材第三章。 三、实验内容 1、安装所需计算工具（MATLAB、fortran、C++等）； 2、学习使用rand(m,n)、unifrnd(a,b,m,n)函数 3、求解下列问题： 3.0、蒲丰氏投针求圆周率。 2 3 2 3.1、给定曲线y =2 – x 和曲线y = x ，曲线的交点为：P ( – 1 ，1 )、P ( 1，1 )。曲线 1 2 围成平面有限区域，用蒙特卡罗方法计算区域面积；  2 2 z x y 3.2、计算 所围体积 2 2 z1 1x y 其中{(x, y, z) | 1x1, 1y1,0 z2} 。 4 、对以下已知分布进行随机抽样： 2 3   3  2   4.1、f x 2x  1x ,x 0,1 ； 2   1     4.2、f x   f x ,x1,2E1 E k E 1    E1x2 1 1 1 f x     1  Ex  x x x2 3 其中   。     2 E1   1 4 1 k E  1 2 ln 2E1     2   E  2 E 2 2E1 四、实验报告编写 1、给出各题的抽样程序并解释语句的含义； 2、给出3.1和3.2抽样结果误差随抽样次数的关系图，并解释原因； 表1 实验记录表 序号 1 2 3 4 5 6 7