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

结课大作业题目自拟:

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

代码实现:

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

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

相关推荐
Hello_Embed2 小时前
STM32HAL 快速入门(二十):UART 中断改进 —— 环形缓冲区解决数据丢失
笔记·stm32·单片机·学习·嵌入式软件
咸甜适中2 小时前
rust语言 (1.88) 学习笔记:客户端和服务器端同在一个项目中
笔记·学习·rust
Magnetic_h3 小时前
【iOS】设计模式复习
笔记·学习·ios·设计模式·objective-c·cocoa
研梦非凡4 小时前
ICCV 2025|从粗到细:用于高效3D高斯溅射的可学习离散小波变换
人工智能·深度学习·学习·3d
limengshi1383925 小时前
机器学习面试:请介绍几种常用的学习率衰减方式
人工智能·学习·机器学习
知识分享小能手5 小时前
React学习教程,从入门到精通,React 组件核心语法知识点详解(类组件体系)(19)
前端·javascript·vue.js·学习·react.js·react·anti-design-vue
周周记笔记6 小时前
学习笔记:第一个Python程序
笔记·学习
优雅鹅7 小时前
ARM、AArch64、amd64、x86_64、x86有什么区别?
arm开发·学习
..过云雨7 小时前
05.【Linux系统编程】进程(冯诺依曼体系结构、进程概念、进程状态(注意僵尸和孤儿)、进程优先级、进程切换和调度)
linux·笔记·学习
咸甜适中7 小时前
rust语言 (1.88) egui (0.32.2) 学习笔记(逐行注释)(二十八)使用图片控件显示图片
笔记·学习·rust·egui