- 1、本文档共6页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
概率论与数理统计论文
基于matlab的蒙特卡罗
定积分的实现
班级:
姓名:
基于matlab的蒙特卡罗定积分的实现
摘 要:在对蒙特卡罗方法概念学习的基础上,利用计算机产生了随机数序列.在此基础上,研究了蒙特卡罗积分,并应用matlab加以实现.
关键词:蒙特卡罗;MATLAB;定积分
0 引言
随着电子计算的出现和发展,近年来用概率模型来作近似计算的方法得到了很大的发展,即蒙特卡罗(Monte—Garlo)方法.它是一种采用统计抽样理论近似的求解数学与物理问题的方法,它既可以用来研究概率问题,也可以用来解决非概率问题.蒙特卡罗法已被广泛地运用到各个领域中,如高维数学问题求解、医学技术中的诊断识别、大型系统的可靠性分析等.一般蒙特卡罗方法在数学中最常见的应用就是蒙特·卡罗积分.对于那些由于计算过于复杂而难以得到解析解或者根本没有解析解的问题,蒙特卡罗方法能够提供一种有效而可行的解决方法.
1 蒙特卡罗方法的基本原理
用蒙特卡罗方法解决数学分析问题时,基本思想是:首先建立与描述该问题有相似性的概率模型,利用这些相似把使该模型的一些数字特征(如概率或均值)与数学分析的解答(如积分值或微分方程的解)联系起来;然后对该模型进行随机模拟或统计抽样;最后利用所得的结果求出这些特征的统计估计值作为原来问题的近似解.蒙特卡罗方法计算的结果收敛的理论依据来自于大数定律,且结果渐进地服从正态分布的理论依据是中心极限定理.
本文将对蒙特卡罗方法计算定积分,并用matlab加以实现.
2 蒙特卡罗方法求定积分的一个实例
2.1 采用均匀随机数的蒙特卡罗方法计算和 matlab 实现
为了使用蒙特卡罗方法求该积分,首先选定一个矩形区域 D =[0,1;O,1].然后抽取n个随即点的坐标(i=1,2,?,n),他们服从区域D上的二维均匀分布.计算落在区域D内的随机点数设为m,即满足0≤≤.则由强大数定律,要求的定积分近似的等于两个区域随机落点的比值 与矩形面积的乘积.
matlab程序及其运行结果:
2.2 matlab生成随机掷点效果图
x1=rand(1,2000);
y1=rand(1,length(x1));
x=0:0.01:1;
y2=0;y3=x.^3;
subplot(2,1,1),plot(x,y2,k,x,y3,k,x1,y1,.)
xlabel(x轴),ylabel(y轴)
m=1;
for j=10:10:30000;
a(m)=j;
x0=rand(1,j);
y0=rand(1,length(x0));
s=0;n=length(x0);
for i=1:n y4(i)=0;y5(i)=x0(i).^3;
if y0(i)y5(i)y0(i)y4(i) s=s+1;
end
end
p(m)=s/n;m=m+1;
end
t=mean(p);
t
xx=0;length(x0);yy=1/4;
subplot(2,1,2),plot(a,p,b,xx,yy,b)
t ≈ 0.2499
2.3 matlab均匀随机数法计算定积分
由以上运行结果可知,均匀随机数法求的定积分
≈ 0.2499
由定积分的几何意义可知,它是三次曲线与直线y = 0所围的平面图形的面积,由微积分基本公式可知,.模拟值与精确值的误差率为0.4%.
3 小结
Monte Carlo方法的实质是通过大量随机试验,利用概率论解决问题的一种数值方法,基本思想是基于概率和体积间的相似性,因此程序十分简单;蒙特卡罗方法在解决确定性问题时,具有较高的精确度;误差主要取决于样本的容量而与维数无关,随着抽取点增加,近似面积也将逼近真实面积;对问题求解的过程取决于所建立的概率模型.
参 考 文 献:
[1] 概率论与数理统计教程/茆诗松,程依明,濮晓龙编著. —2版. —北京:高等教育出版社,2011.2
[2] MATLAB教程及实训/曹弋主编. —北京:机械工业出版社,2008.4
您可能关注的文档
- 听障儿童听觉、言语培建目标.doc
- 周其凤校长在北京大学2010年本科毕业典礼上的讲话.doc
- 周杰伦所有专辑名称及内容歌名.doc
- 呼啸山庄男女主角的爱情.doc
- 和田玉相关文案集.doc
- 咏花卉草木诗词之三.doc
- 品一壶诗歌的香茗.doc
- 品德教养,让人生更精彩.doc
- 品牌推广是指企业塑造自身及产品品牌形象使广大消费者甚至精准的目标消费者广泛的认同的系列动作.doc
- 品读永州道州篇庄村遗韵.doc
- 10《那一年,面包飘香》教案.docx
- 13 花钟 教学设计-2023-2024学年三年级下册语文统编版.docx
- 2024-2025学年中职学校心理健康教育与霸凌预防的设计.docx
- 2024-2025学年中职生反思与行动的反霸凌教学设计.docx
- 2023-2024学年人教版小学数学一年级上册5.docx
- 4.1.1 线段、射线、直线 教学设计 2024-2025学年北师大版七年级数学上册.docx
- 川教版(2024)三年级上册 2.2在线导航选路线 教案.docx
- Unit 8 Dolls (教学设计)-2024-2025学年译林版(三起)英语四年级上册.docx
- 高一上学期体育与健康人教版 “贪吃蛇”耐久跑 教案.docx
- 第1课时 亿以内数的认识(教学设计)-2024-2025学年四年级上册数学人教版.docx
最近下载
- 7.2 类比推理及其方法-高中政治课件 (统编版选择性必修3).pptx VIP
- 《数学物理方程-福州大学-江飞》作业chapter1.pdf VIP
- 重庆渝北中交·中央公园 C96, C98-1 地块山地新中式商业街项目 GOA.pdf
- 2024年江苏省高考物理真题试卷含答案.pdf VIP
- 《数学物理方程-福州大学-江飞》数学物理方程A.doc VIP
- 《数学物理方程-福州大学-江飞》作业chapter2.ppt VIP
- 《数学物理方程-福州大学-江飞》第四章.doc VIP
- 《数学物理方程-福州大学-江飞》数学物理方程A答案.doc VIP
- 2023年辽宁省检察系统招聘聘用制书记员考试真题及答案.docx VIP
- 2024年高考真题——物理(河北卷)含答案.pdf VIP
文档评论(0)