蒙特卡罗方法及应用实验讲义2016

《蒙特卡罗方法及应用实验讲义2016》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、蒙特卡罗方法及应用实验讲义东华理工大学核工系2016.815实验一蒙特卡罗方法基本思想一、实验目的1、了解蒙特卡罗方法方法的基本思想；2、掌握蒙特卡罗方法计算面积、体积的方法；3、掌握由已知分布的随机抽样方法。二、实验原理MonteCarlo方法，又称统计模拟方法或计算机随机模拟方法，是一种基于“随机数”进行数值模拟的方法，一种采用统计抽样理论近似求解物理或数学问题的方法。如待求量可以表述成某些特征量的期望值、某些事件出现的概率或两者的函数形式，那么可采用蒙特卡罗方法求解。在求解某些特征量的期望值或

2、某些事件出现的概率时，必须构建合符实际的数学模型。例如采用蒙特卡罗方法计算某函数所围面积时，构建的数学模型是构造一已知面积的可均匀抽样区域，在该区域投点，由伯努利定理大数定理可知，进入待求区域投点的频率依概率1收敛于该事件出现的概率（面积之比）。由已知分布的随机抽样方法指的是由已知分布的总体中抽取简单子样。具体方法很多，详见教材第三章。三、实验内容1、安装所需计算工具（MATLAB、fortran、C++等）；2、学习使用rand(m,n)、unifrnd(a,b,m,n)函数3、求解下列问题：3.

3、0、蒲丰氏投针求圆周率。3.1、给定曲线y=2–x2和曲线y3=x2，曲线的交点为：P1(–1，1)、P2(1，1)。曲线围成平面有限区域，用蒙特卡罗方法计算区域面积；3.2、计算所围体积其中。4、对以下已知分布进行随机抽样：154.1、；4.2、其中。四、实验报告编写1、给出各题的抽样程序并解释语句的含义；2、给出3.1和3.2抽样结果误差随抽样次数的关系图，并解释原因；表1实验记录表序号1234567试验次数1032×1045×1041.2×1051.5×1051.8×1062.0×107试验时

4、间计算结果实验误差3、给出4.1和4.2的抽样框图、试验累积频率与理论累积频率关系图，并给出抽样次数（>106）与抽样时间。15实验二由已知分布的随机抽样方法一、实验目的1、掌握由已知分布的随机抽样方法。2、用编程语言实现某具体随机抽样方法。二、实验原理由已知分布的随机抽样方法指的是由已知分布的总体中抽取简单子样。具体方法很多，本实验综合直接抽样方法、挑选抽样方法和替换抽样方法，以散射方位角余弦分布的抽样为例。实验原理详见教材对应章节。1.连续型分布的直接抽样方法对于连续型分布，如果分布函数F(x)

5、的反函数F－1(x)存在，则直接抽样方法是：2.挑选抽样方法为了实现从己知分布密度函数f(x)抽样，选取与f(x)取值范围相同的分布密度函数h(x)，如果则挑选抽样方法为：3.替换法抽样方法为了实现某个复杂的随机变量y的抽样，将其表示成若干个简单的随机变量x1，x2，…，xn的函数得到x1，x2，…，xn的抽样后，即可确定y的抽样，这种方法叫作替换法抽样。15三、实验内容蒲丰氏问题的算法如何产生任意的（x,θ）？x在［0，a］上任意取值，表示x在［0，a］上是均匀分布的，其分布密度函数为：类似地，θ

6、的分布密度函数为：因此，产生任意的（x,θ）的过程就变成了由f1(x)抽样x及由f2(θ)抽样θ的过程了。由此得到：其中ξ1，ξ2均为（0,1）上均匀分布的随机变量。每次投针试验，实际上变成在计算机上从两个均匀分布的随机变量中抽样得到（x,θ），然后定义描述针与平行线相交状况的随机变量s(x,θ)，为如果投针Ｎ次，则是针与平行线相交概率Ｐ的估计值。事实上，于是有15四、实验报告1、给出源程序程序并解释语句的含义；2、作出抽样框图、试验累积频率与理论累积频率关系图，并给出抽样次数（>106）与抽样时间

7、。15实验三MCNP方法在实验核物理中的应用一、实验目的1、了解MCNP程序运行流程；2、掌握MCNP输入文件编写规范；3、理解模拟内容、并能编写输入文件、运行，并获得计算结果；二、实验原理MCNP是一种常见的粒子输运模拟软件，软件的安装、运行和输入文件编写方法详见相关参考资料。MCNP输入文件编写完成后，先确认输入模型是否正确，在DOS环境下进行，打开运行DOS环境，进行以下操作：DOS命令操作命令含义Mcnpi=name.inpo=name.o打开画图框PXvx输出模型在x=vx面上的切面PYv

8、y输出模型在y=vy面上的切面PZvz输出模型在z=vz面上的切面FACTORm将输出图放大1/m倍Extentab切面沿两坐标轴方向分别放大ORIGINXYZ定义画图中心位置（X,Y,Z）三、实验内容1、学习MCNP程序常见各种运行方法；2、编写以下问题的输入文件；2.1对课堂讲解的实例，模拟溴化镧探测器对点源的能谱，实验做一遍。2.2有一HPGe探测器，结构如图1所示。分别给出位于探测器轴心、距离探测器晶体中心25cm处的137Cs源、60Co、131I源对应特征