首先看看上面这个问
题
快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题
。这个问题是我在一个MATLAB交流群里碰到的提问,计算阴影部分面积。
什么是蒙特卡罗在这里我就不多做介绍了,感兴趣的朋友可以自己去查阅相关资料,相信可以得到全面的解释,在这里我只介绍如果用蒙特卡罗
方法
快递客服问题件处理详细方法山木方法pdf计算方法pdf华与华方法下载八字理论方法下载
来计算上图中阴影部分的面积,注意这只是蒙特卡罗方法的一个应用而已。
记得第一次接触到蒙特卡罗是在一次数学建模
培训
焊锡培训资料ppt免费下载焊接培训教程 ppt 下载特设培训下载班长管理培训下载培训时间表下载
中。当时我们老师给我们讲了一个故事,故事的全部我已经记不清了,大概
内容
财务内部控制制度的内容财务内部控制制度的内容人员招聘与配置的内容项目成本控制的内容消防安全演练内容
是:
一个古人要求一个图形的面积,他把图形画在一块方形布上,然后找来一袋豆子,然后将所有豆子洒在布上,落在图形内豆子的重量比上那块布上所有豆子的重量再乘以布的面积就是他所要求的图形的面积。
这确实是一个求面积的好方法,这是我听到这个故事后的第一反应。从此我就记住了这个方法,记得很深刻。所以当群里有人问如何求上面这个图形的面积的时候我马上就回想起用蒙特卡罗方法来计算。
仔细思考后,以我的知识面我能找到两种编程思路来计算这个面积:
方法一:将整个坐标轴看成一个边长为12的正方形,然后均匀的这个正方形分成N(N的大小取决于划分的步长)个点,然后找出N个点中有多少个点是属于阴影部分中,假设这个值为k,则阴影部分的面积为:k/N*12^2
方法二:将整个坐标轴看成一个边长为12的正方形,然后在(-6,6)中随机出N(N越大越好,至少超过1000)个点,然后找出这N个点中有多少个点在阴影区域内,假设这个值为k,则阴影部分的面积为:k/N*12^2。然后重复这个过程100次,求出100次面积计算结果的均值,这个均值为阴影部分面积。
对比
分析
定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析
:以上两个方法都是利用蒙特卡罗方法计算阴影部分面积,只是在处理的细节有一点区别。前者是把豆子均匀分布在布上;后者则是随机把豆子仍在布上。就计算结果的精度而言,前者取决点的分割是否够密,即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