1 / 11
文档名称:

蒙特卡洛实验(一)报告.doc

格式:doc   大小:149KB   页数:11页
下载后只包含 1 个 DOC 格式的文档,没有任何的图纸或源代码,查看文件列表

如果您已付费下载过本站文档,您可以点这里二次下载

分享

预览

蒙特卡洛实验(一)报告.doc

上传人:guoxiachuanyue005 2021/11/30 文件大小:149 KB

下载得到文件列表

蒙特卡洛实验(一)报告.doc

相关文档

文档介绍

文档介绍:大型企业经典管理资料模板, WORD文档,欢迎下载交流
: .
分享一个苹果,各得一个苹果,分享一种思想,各得两种思想。分享是件快乐的事件,乐于分享的人,事业更容易成功。
经典资料,WORD文档,可编辑修改,欢迎下载交流。
专业:核工程与核技术
实验一蒙特卡罗方法
一、 实验目的
1、 了解蒙特卡罗方法方法的基本思想;
2、 掌握蒙特卡罗方法计算面积、体积的方法;
3、 掌握由已知分布的随机抽样方法。
二、 实验原理
Monte Carlo方法,又称统计模拟方法或计算机随机模拟方法,是一种基于“随机数”进行数值模拟的 方法,一种采用统计抽样理论近似求解物理或数学问题的方法。
(面
倘若待求量可以表述成某些特征量的期望值、某些事件出现的概率或两者的函数形式,那么可采用蒙 特卡罗方法求解。在求解某些特征量的期望值或某些事件出现的概率时,必须构建合符实际的数学模型。 例如采用蒙特卡罗方法计算某函数所围面积时,构建的数学模型是构造一已知面积的可均匀抽样区域,在 该区域投点,由伯努利定理大数定理可知, 进入待求区域投点的频率依概率 1收敛于该事件出现的概率
积之比)。
由已知分布的随机抽样方法指的是由已知分布的总体中抽取简单子样。具体方法很多,详见课本第三
早。
三、 实验内容
1、安装所需计算工具(MATLAB等);
以下内容采用工具软件中自带伪随机数发生器进行计算。
2、求解以下区域的面积、体积:
、给定曲线y =2 -x2和曲线y3 = x2,曲线的交点为:R( -1,1 )、P2( 1,1 )。曲线围成平面有限区域, 用蒙特卡罗方法计算区域面积;
、计算
所围体积
z _、x2 y2
z _1 、1 -x2 - y2
其中门二{(x,y,z)|—1 空 x 乞1,—1^y 汨,0 乞 z 空 2}。
3、对以下已知分布进行随机抽样:
3 2
f x i=2x3 - 1 - x ,x 〔0,11
2
经典资料,WORD文档,可编辑修改,欢迎下载交流。
大型企业经典管理资料模板, WORD文档,欢迎下载交流
、实验报告编写
1给出各题的抽样程序并解释语句的含义;
2、,并解释原因;
表1 实验记录表
3、给出3题的抽样框图、试验累积频率与理论累积频率关系图, 并给出抽样次数(>106)与抽样时间。

N=10八6;%总投点个数
S=0; %记录投点在所围图形中的个数
SS=0;
for i=1:N
分享一个苹果,各得一个苹果,分享一种思想,各得两种思想。分享是件快乐的事件,乐于分享的人,事业更容易成功。
大型企业经典管理资料模板, WORD文档,欢迎下载交流
x=2*rand-1;%产生的随机变量 x,y
y=2*rand; ;%产生x和y的坐标
if((yv=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)A2%+ 算方差
A=sqrt(Dev/N)%计算标准差
toc
实验数据如下:
请输入总投点个数:
150000
2




1




0
i I I I 耳 n
-1 - - - - 0 1

码如下:
clear;
clc;
M=0;
N= 5*10八4 ;
tic;
for i=1:N
x=2*ra nd()-1;
y=2*ra nd()-1;
z=2*ra nd();
t=xA2+yA2;
s=zA2;
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;