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

结课大作业题目自拟:

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

代码实现:

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增大,面积值变化不再明显。

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

相关推荐
CV学术叫叫兽2 分钟前
一站式学习:害虫识别与分类图像分割
学习·分类·数据挖掘
我们的五年13 分钟前
【Linux课程学习】:进程程序替换,execl,execv,execlp,execvp,execve,execle,execvpe函数
linux·c++·学习
一棵开花的树,枝芽无限靠近你37 分钟前
【PPTist】添加PPT模版
前端·学习·编辑器·html
VertexGeek1 小时前
Rust学习(八):异常处理和宏编程:
学习·算法·rust
二进制_博客2 小时前
Flink学习连载文章4-flink中的各种转换操作
大数据·学习·flink
火山口车神丶2 小时前
某车企ASW面试笔试题
c++·matlab
codebolt2 小时前
ADS学习记录
学习
今天吃饺子2 小时前
2024年SCI一区最新改进优化算法——四参数自适应生长优化器,MATLAB代码免费获取...
开发语言·算法·matlab
Komorebi.py3 小时前
【Linux】-学习笔记05
linux·笔记·学习
朝九晚五ฺ10 小时前
【Linux探索学习】第十四弹——进程优先级:深入理解操作系统中的进程优先级
linux·运维·学习