当前位置:首页 >> 三年级数学 >>

蒙特卡罗算法举例


首先看看上面这个问题。这个问题是我在一个 MATLAB 交流群里碰到的提问,计算阴影部 分面积。 什么是蒙特卡罗在这里我就不多做介绍了,感兴趣的朋友可以自己去查阅相关资料,相信可 以得到全面的解释, 在这里我只介绍如果用蒙特卡罗方法来计算上图中阴影部分的面积, 注意这 只是蒙特卡罗方法的一个应用而已。 记得第一次接触到蒙特卡罗是在一次数学建模培训中。当时我们老师给我们讲了一个故事, 故事的全部我已经记不清了,大概内容是: 一个古人要求一个图形的面积,他把图形画在一块方形布上,然后找来一袋豆子,然后将所 有豆子洒在布上, 落在图形内豆子的重量比上那块布上所有豆子的重量再乘以布的面积就是他所 要求的图形的面积。 这确实是一个求面积的好方法,这是我听到这个故事后的第一反应。从此我就记住了这个方 法, 记得很深刻。 所以当群里有人问如何求上面这个图形的面积的时候我马上就回想起用蒙特卡 罗方法来计算。 仔细思考后,以我的知识面我能找到两种编程思路来计算这个面积: 方法一:将整个坐标轴看成一个边长为 12 的正方形,然后均匀的这个正方形分成 N(N 的大小 取决于划分的步长)个点,然后找出 N 个点中有多少个点是属于阴影部分中,假设这个值为 k, 则阴影部分的面积为:k/N*12^2 方法二:将整个坐标轴看成一个边长为 12 的正方形,然后在(-6,6)中随机出 N(N 越大越 好,至少超过 1000)个点,然后找出这 N 个点中有多少个点在阴影区域内,假设这个值为 k, 则阴影部分的面积为:k/N*12^2。然后重复这个过程 100 次,求出 100 次面积计算结果的均 值,这个均值为阴影部分面积。

对比分析: 以上两个方法都是利用蒙特卡罗方法计算阴影部分面积, 只是在处理的细节有一点区 别。前者是把豆子均匀分布在布上;后者则是随机把豆子仍在布上。就计算结果的精度而言,前 者取决点的分割是否够密,即 N 是否够大;后者不仅仅通过 N 来控制精度,因为随机的因素会 造成单次计算结果偏高和偏小,所以进行反复多次计算最后以均值来衡量阴影部分面积。 附上 MATLAB 程序: 方法一: clear x=-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 a1<1 if a2<1 if a3<9 k=k+1; end end end end end mj=(12^2)*k/zs; 运行结果: mj = 7.2150 方法二: clear N=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 a1<1 if a2<1 if a3<9 k=k+1; end end end end m(j)=(12^2)*k/N; end mj=mean(m); 运行结果: mj = 7.2500


相关文章:
蒙特卡罗 算法
蒙特卡罗 算法 - 1、蒙特卡罗定位 足球机器人中自定位方法是由 Fox 提出的蒙特卡罗定位。这是一种概率方法,把足球机器 人当前位置看成许多粒子的密度模型。 每个...
蒙特卡罗算法
蒙特卡罗算法 - 第一阶段:算法 ——蒙特卡洛算法 蒙特卡罗方法(MC) (Monte Carlo) : 蒙特卡罗(Monte Carlo)方法,或称计算机随机模拟方法,是一种基于“随 机数”...
蒙特卡罗算法与matlab(精品教程)
蒙特卡罗算法与matlab(精品教程)_理学_高等教育_教育专区。蒙特卡罗算法与matlab...的一些基础知识,另外包括一个最简单的用 Monte Carlo 方法计算数值积分的例子。...
蒙特卡罗算法【数学建模】
蒙特卡罗算法【数学建模】 - 首先看看上面这个问题。这个问题是我在一个 MATLAB 交流群里碰到的提问,计算阴影部 分面积。 一个古人要求一个图形的面积,他把图形画...
蒙特卡洛算法详讲
蒙特卡洛算法详讲 - Monte Carlo 法§8.1 概述 Monte Carlo 法不同于前面几章所介绍的确定性数值方法,它是用来解决数 学和物理问题的非确定性的(概率统计的或...
蒙特卡罗方法模拟计算圆周率
蒙特卡罗方法模拟计算圆周率 - 蒙特卡罗方法模拟计算圆周率 摘要:蒙特卡罗方法是一种新型计算方法,它需要真 实的随机数,在统计学方面有强有力的应用,随着高性能计 算...
蒙特卡洛方法模拟小例子
蒙特卡洛方法模拟小例子_理学_高等教育_教育专区。蒙特卡洛方法模拟小例子 ...可以说,计算高维积分是蒙特卡洛方法最成功和典型 的应用。 基本的蒙特卡洛法具有...
蒙特卡罗算法计算Pi
蒙特卡罗算法计算Pi - 蒙特卡罗算法 以概率和统计理论方法为基础的一种计算方法。将所求解的问题同一定的概率模型相联系,用电子计算机实现统计模拟或抽样,以获得问题...
蒙特卡洛算法
蒙特卡洛算法 - 根据我的理解简单的说就是以部分估计整体, 利用随机数来解决问题的方法称为 蒙特卡罗算法,记得课本上讲了个例题: 在数值积分法中,我们利用求单位...
蒙特卡罗算法举例
蒙特卡罗算法举例 - 首先看看上面这个问题。这个问题是我在一个 MATLAB 交流群里碰到的提问,计算阴影部 分面积。 什么是蒙特卡罗在这里我就不多做介绍了,感兴趣的...
更多相关标签: