蒙特卡洛实验(一)报告

蒙特卡洛实验(一)报告

ID:35236150

大小:342.00 KB

页数:7页

时间:2019-03-22

蒙特卡洛实验(一)报告_第1页
蒙特卡洛实验(一)报告_第2页
蒙特卡洛实验(一)报告_第3页
蒙特卡洛实验(一)报告_第4页
蒙特卡洛实验(一)报告_第5页
资源描述:

《蒙特卡洛实验(一)报告》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库

1、姓名:王蕴楠学号:1020330210专业:核工程与核技术实验一蒙特卡罗方法一、实验目的1、了解蒙特卡罗方法方法的基本思想;2、掌握蒙特卡罗方法计算面积、体积的方法;3、掌握由已知分布的随机抽样方法。二、实验原理MonteCarlo方法,又称统计模拟方法或计算机随机模拟方法,是一种基于“随机数”进行数值模拟的方法,一种采用统计抽样理论近似求解物理或数学问题的方法。倘若待求量可以表述成某些特征量的期望值、某些事件出现的概率或两者的函数形式,那么可采用蒙特卡罗方法求解。在求解某些特征量的期望值或某些事件出现的概率时,必须构建合符实际的数学模型。例如采用蒙特卡罗方法计算某函

2、数所围面积时,构建的数学模型是构造一已知面积的可均匀抽样区域,在该区域投点,由伯努利定理大数定理可知,进入待求区域投点的频率依概率1收敛于该事件出现的概率(面积之比)。由已知分布的随机抽样方法指的是由已知分布的总体中抽取简单子样。具体方法很多,详见课本第三章。三、实验内容1、安装所需计算工具(MATLAB等);以下内容采用工具软件中自带伪随机数发生器进行计算。2、求解以下区域的面积、体积:2.1、给定曲线y=2–x2和曲线y3=x2,曲线的交点为:P1(–1,1)、P2(1,1)。曲线围成平面有限区域,用蒙特卡罗方法计算区域面积;2.2、计算所围体积其中。3、对以下已

3、知分布进行随机抽样:三、实验报告编写1、给出各题的抽样程序并解释语句的含义;2、给出2.1和2.2抽样结果误差随抽样次数的关系图,并解释原因;表1实验记录表序号1234567试验次数1031×1045×1041.2×1051.5×1051.8×1062.0×107试验时间计算结果实验误差3、给出3题的抽样框图、试验累积频率与理论累积频率关系图,并给出抽样次数(>106)与抽样时间。2.1程序代码编写如下:N=10^6;%总投点个数S=0;%记录投点在所围图形中的个数SS=0;fori=1:Nx=2*rand-1;%产生的随机变量x,yy=2*rand;;%产生x和y的

4、坐标if((y<=2-x^2)&(y^3>=x^2))%判定是否落入所围图像中S=S+1;%进入则加1SS=SS+1^2;endendArea=4*S/N%计算面积Dev=SS/N-(S/N)^2%计算方差A=sqrt(Dev/N)%计算标准差toc实验数据如下:序号1234567试验次数1031×1045×1041.2×1051.5×1051.8×1062.0×107试验时间0.000000s0.031000s1.182000s17.609000s29.453000s310.23515s4.4532523s计算结果2.13202.13522.14442.12922.

5、13942.12072.1258实验误差0.01580.00500.00220.00140.00134.198e-0045.568e-004请输入总投点个数:1500002.2实验代码如下:clear;clc;M=0;N=5*10^4;tic;fori=1:Nx=2*rand()-1;y=2*rand()-1;z=2*rand();t=x^2+y^2;s=z^2;ifs>=tift<=-s+2*zM=M+1;endendendtocMIANJI=M/N*8clearMNixy;计算结果:N=50000时面积为3.1350,计算时间约0.282s。实验数据如下:序号12

6、34567试验次数1035×1045×1041.2×1051.5×1051.8×1062.0×107试验时间0.000000s0.016000s0.000000s0.016000s0.204000s0.265000s2.609000s计算结果2.80803.15943.15363.13953.14013.14713.1419实验误差0.3336-0.0178-0.01200.00210.0015-0.0055-3.0735e-0042.3程序代码编写如下:clear;clc;M=input('输入所需产生随机变量的个数:');x=zeros(M,1);tic;fo

7、ri=1:Mif(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.')tocclearM;六、实验心得通过本次实验后,让我发现这门课非常有趣,并没有想象的那么枯燥无味,是一门很有实用价值的一门学科。同时让我学习到MATLAB的基本操作和用法。

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。