蒙特卡罗算法举例.doc_第1页
蒙特卡罗算法举例.doc_第2页
蒙特卡罗算法举例.doc_第3页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

首先看看上面这个问题。这个问题是我在一个MATLAB交流群里碰到的提问，计算阴影部分面积。 什么是蒙特卡罗在这里我就不多做介绍了，感兴趣的朋友可以自己去查阅相关资料，相信可以得到全面的解释，在这里我只介绍如果用蒙特卡罗方法来计算上图中阴影部分的面积，注意这只是蒙特卡罗方法的一个应用而已。 记得第一次接触到蒙特卡罗是在一次数学建模培训中。当时我们老师给我们讲了一个故事，故事的全部我已经记不清了，大概内容是： 一个古人要求一个图形的面积，他把图形画在一块方形布上，然后找来一袋豆子，然后将所有豆子洒在布上，落在图形内豆子的重量比上那块布上所有豆子的重量再乘以布的面积就是他所要求的图形的面积。 这确实是一个求面积的好方法，这是我听到这个故事后的第一反应。从此我就记住了这个方法，记得很深刻。所以当群里有人问如何求上面这个图形的面积的时候我马上就回想起用蒙特卡罗方法来计算。 仔细思考后，以我的知识面我能找到两种编程思路来计算这个面积：方法一：将整个坐标轴看成一个边长为12的正方形，然后均匀的这个正方形分成N（N的大小取决于划分的步长）个点，然后找出N个点中有多少个点是属于阴影部分中，假设这个值为k，则阴影部分的面积为：k/N*122方法二：将整个坐标轴看成一个边长为12的正方形，然后在（-6，6）中随机出N（N越大越好，至少超过1000）个点，然后找出这N个点中有多少个点在阴影区域内，假设这个值为k，则阴影部分的面积为：k/N*122。然后重复这个过程100次，求出100次面积计算结果的均值，这个均值为阴影部分面积。对比分析：以上两个方法都是利用蒙特卡罗方法计算阴影部分面积，只是在处理的细节有一点区别。前者是把豆子均匀分布在布上；后者则是随机把豆子仍在布上。就计算结果的精度而言，前者取决点的分割是否够密，即N是否够大；后者不仅仅通过N来控制精度，因为随机的因素会造成单次计算结果偏高和偏小，所以进行反复多次计算最后以均值来衡量阴影部分面积。附上MATLAB程序：方法一：clearx=-6:0.01:6;y=x;s=size(x);zs=s(1,2)2;k=0;for i=1:s(1,2) for j=1:s(1,2) a1=(x(i)2)/9+(y(j)2)/36; a2=(x(i)2)/36+y(j)2; a3=(x(i)-2)2+(y(j)+1)2; if a11 if a21 if a39 k=k+1; end end end endendmj=(122)*k/zs;运行结果：mj = 7.2150方法二：clearN=10000;n=100;for j=1:n k=0;for i=1:N a=12*rand(1,2)-6; x(i)=a(1,1); y(i)=a(1,2); a1=(x(i)2)/9+(y(i)2)/36; a2=(x(i)2)/36+y(i)2; a3=(x(i)-2)2+(y(i)+1)2; if a11 if a21 if a39 k
人人文库> 全部分类> 行业资料 > 管理策划

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论