蒙特卡诺法计算不规则图形的面积

结课大作业题目自拟:

利用蒙特卡诺法计算不规则图形的面积,由于题目自拟,所以选择用椭圆来表示待测面积的不规则图形。

代码实现:

Matlab 复制代码
%%先绘制一个椭圆
a=5;%长轴
b=2;%短轴
f=@(x,y) x.^2/a^2 + y.^2/b^2-1;%椭圆方程
fimplicit(f,'LineWidth',2)%绘制椭圆 
设置窗口
set(gcf,'units','normalized','position',[0.2 0.2 0.6 0.6]); %设置figure窗口位置和尺寸
grid on
axis equal
axis([-(a+0.5),(a+0.5),-(b+0.5),(b+0.5)])
pause(2)
hold on
绘制椭圆的包络矩形
rectangle('Position',[-a,-b,2*a,2*b],'EdgeColor','r','LineWidth',2)
A=4*a*b;%矩形面积
pause(2)
投放N个均匀分布的随机点
N=500;
xk=-a+(a+a)*rand(1,N);   %x坐标
yk=-b+(b+b)*rand(1,N);   %y坐标
scatter(xk,yk,'g.')
统计落在椭圆里面的随机点个数
r = xk.^2/a^2+yk.^2/b^2;
m=find(r<=1);
n=length(m);
S=(n/N)*A

效果演示:

注:由于计算过程比较复杂,所以我是选择在matlab中的实时脚本里写的。

注:为了图像效果较好,结果图中只生成了500个随机点,结果不是很准确,增加至1000以上时,随N增大,面积值变化不再明显。

网上也有其他类似的代码,可以多方借鉴。

相关推荐
hmbbcsm8 分钟前
python学习之路(二)
学习
冷崖8 分钟前
定时器的学习(二)
linux·c++·学习
好奇龙猫1 小时前
【学习AI-相关路程-mnist手写数字分类-一段学习的结束:自我学习AI-复盘-代码-了解原理-综述(5) 】
人工智能·学习·分类
小苏兮1 小时前
【C++】priority_queue和deque的使用与实现
开发语言·c++·学习
Yurko132 小时前
【C语言】环境安装(图文)与介绍
c语言·开发语言·学习
十安_数学好题速析2 小时前
数论探秘:如何用模4思想破解平方数谜题
笔记·学习·高考
在繁华处3 小时前
C语言初步学习:数组的增删查改
c语言·数据结构·学习
呵呵哒( ̄▽ ̄)"3 小时前
专项智能练习(科尔伯格道德发展阶段理论)
学习
武文斌773 小时前
项目学习总结:CAN总线、摄像头、STM32概述
linux·arm开发·stm32·单片机·嵌入式硬件·学习·c#
SPFFC189380330533 小时前
AI玩具排线专业生产与全球营销策略
人工智能·学习·智能手机·显示器·智能手表·平板·游戏机