蒙特卡洛实验报告总结计划.docx

专业：核工程与核技术 实验一 蒙特卡罗方法 一、实验目的 1、了解蒙特卡罗方法方法的基本思想； 2、掌握蒙特卡罗方法计算面积、体积的方法； 3、掌握由已知分布的随机抽样方法。 二、实验原理 Monte Carlo 方法，又称统计模拟方法或计算机随机模拟方法，是一种基于 “随机数”进行数值模拟的方法，一种采用统计抽样理论近似求解物理或数学 问题的方法。 倘若待求量可以表述成某些特征量的期望值、某些事件出现的概率或两者的函数形式，那么可采用蒙特卡罗方法求解。在求解某些特征量的期望值或某 些事件出现的概率时，必须构建合符实际的数学模型。例如采用蒙特卡罗方法计算某函数所围面积时，构建的数学模型是构造一已知面积的可均匀抽样区域， 在该区域投点， 由伯努利定理大数定理可知， 进入待求区域投点的频率依概率 1 收敛于该事件出现的概率（面积之比） 。 由已知分布的随机抽样方法指的是由已知分布的总体中抽取简单子样。具 体方法很多，详见课本第三章。 三、实验内容 1、安装所需计算工具（ MATLAB等）； 以下内容采用工具软件中自带伪随机数发生器进行计算。 2、求解以下区域的面积、体积： 、给定曲线 y =2 – x2 和曲线 y3 = x2，曲线的交点为： P1( – 1，1 ) 、P2( 1， 1 ) 。曲线围成平面有限区域，用蒙特卡罗方法计算区域面积； z x2 y 2 所围体积 、计算 z 1 1 x2 y2 其中{( x, y, z) | 1 x 1, 1 y 1,0 z 2} 。 3、对以下已知分布进行随机抽样： 三、实验报告编写 1、给出各题的抽样程序并解释语句的含义； 2、给出和抽样结果误差随抽样次数的关系图，并解释原因； 表 1 实验记录表 3、给出 3 题的抽样框图、试验累积频率与理论累积频率关系图，并给出抽 样次数（ 106 ）与抽样时间。 程序代码编写如下 ： N=10^6;% 总投点个数 S=0; %记录投点在所围图形中的个数 SS=0; for i=1:N x=2*rand-1;% 产生的随机变量 x,y y=2*rand; ;% 产生 x 和 y 的坐标 if((y=2-x^2)(y^3=x^2))% 判定是否落入所围图像中 S=S+1; %进入则加 1 SS=SS+1^2; end end Area=4*S/N %计算面积 Dev=SS/N-(S/N)^2% 计算方差 A=sqrt(Dev/N)% 计算标准差 toc 实验数据如下： 请输入总投点个数 : 150000 2 1.8 1.6 1.4 1.2 1 0.8 0.6 0.4 0.2 0 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 -1 实 验代码如下： clear; clc; M=0; N= 5*10^4 ； tic; for i=1:N x=2*rand()-1; y=2*rand()-1; z=2*rand(); t=x^2+y^2; s=z^2; if s=t if t=-s+2*z M=M+1; end end end toc MIANJI=M/N*8 clear M N i x y; 计算结果： N=50000 时面积为，计算时间约。 实验数据如下： 程序代码编写如下 ： clear; clc; M = input( 输入所需产生随机变量的个数 :\n); x = zeros(M,1); tic; for i=1:M if(rand()= x(i) = max(rand(),rand()) x(i) = max(x(i),rand()); x(i) = max(x(i),rand()); else x(i) = min(rand(),rand()); x(i) = min(x(i),rand()); end end plot(x,g.) toc clear M; 六、实验心得 通过本次实验后，让我发现这门课非常有趣，并没有 想象的那么枯燥无味，是一门很有实用价值的一门学科。 同时让我学习到 MATLAB的基本操作和用法。