《专业：核工程与核技术》由会员分享，可在线阅读，更多相关《专业：核工程与核技术（11页珍藏版）》请在金锄头文库上搜索。
1、专业：核工程与核技术实验一蒙特卡罗方法一、实验目的1、了解蒙特卡罗方法方法的基本思想;2、掌握蒙特卡罗方法计算面积、体积的方法;3、掌握由已知分布的随机抽样方法。二、实验原理MonteCarlo方法，又称统计模拟方法或计算机随机模拟方法，是一种基于“随机数”进行数值模拟的方法，一种采用统计抽样理论近似求解物理或数学问题的方法。倘若待求量可以表述成某些特征量的期望值、某些事件出现 的概率或两者的函数形式，那么可采用蒙特卡罗方法求解。在求 解某些特征量的期望值或某些事件出现的概率时，必须构建合符 实际的数学模型。例如采用蒙特卡罗方法计算某函数所围面积时,构建的数学模型是构造一已知面积的可均匀抽样区
2、域，在该区域 投点，由伯努利定理大数定理可知，进入待求区域投点的频率依概率1收敛于该事件出现的概率（面积之比）。由已知分布的随机抽样方法指的是由已知分布的总体中抽取简单子样。具体方法很多，详见课本第三章。三、实验内容1、安装所需计算工具（MATLAB等）;以下内容采用工具软件中自带伪随机数发生器进行计算。2、求解以下区域的面积、体积：2.1、给定曲线y=2 - X2和曲线y3=x2，曲线的交点为：P, _ 1,1）、P（1，1 ）。曲线围成平面有限区域，用蒙特卡罗方法计算区2域面积；所围体积2.2、计算；z x2+y2z 1 + 抻1 - X 2 y 2其中 Q = （x, y, z） I 1
3、 x 1,1 y 1,0 z 23、对以下已知分布进行随机抽样：三、实验报告编写1、给出各题的抽样程序并解释语句的含义;2、给出2.1和2.2抽样结果误差随抽样次数的关系图，并解释原因；表1实验记录表3、给出3题的抽样框图、试验累积频率与理论累积频率关系 图，并给出抽样次数（106）与抽样时间。2.1程序代码编写如下:N=106%总投点个数S=0； %记录投点在所围图形中的个数SS=0；fori=l：Nx=2*rand-1；%产生的随机变量x,y-来源网络，仅供个人学习参考y=2*rand； ；%产生X和y的坐标f (y=x八2)%判定是否落入所围图像中S=S+1；%进入则加1SS二SS+2；
4、endendArea=4*S/N%计算面积Dev二SS/N-(S/N厂2%计算方差A=sqr t(Dev/N)%计算标准差toe实验数据如下:请输入总投点个数：15000021.81.61.41.210.80.60.40.20-1-0.8-0.6-0.4-0.200.20.40.60.812.2实验代码如下：clear；clc；M=0；N=5*104；t ic；fori=1：Nx=2*rand()T；y=2*rand()T；z=2*rand()；t=x2+y2；s=z2；ifs=tift=-s+2*zM=M+1；endendendtocMIANJI=M/N*8clearMNixy；计算结果：N
5、=50000时面积为3. 1350，计算时间约0282s。实验数据如下：2.3程序代码编写如下：clear；clc；M=inp ut(输入所需产生随机变量的个数:n)；x二zeros(M,1)；t ic；fori=l：M辻(rand()=0.5)x(i)=max(rand(),rand()x(i)=max(x(i),rand();x(i)=max(x(i),rand();elsex(i)二min(rand(),rand()；x(i)二min(x(i),rand()；endendplot (x,g )toeclearM；六、实验心得通过本次实验后，让我发现这门课非常有趣，并没有想象的那么枯燥无味，是一门很有实用价值的一门学科。同时让我学习到MATLAB的基本操作和用法。
