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

结课大作业题目自拟:

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

代码实现:

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

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

相关推荐
Mephisto.java30 分钟前
【大数据学习 | kafka高级部分】kafka中的选举机制
大数据·学习·kafka
南宫生1 小时前
贪心算法习题其三【力扣】【算法学习day.20】
java·数据结构·学习·算法·leetcode·贪心算法
武子康2 小时前
大数据-212 数据挖掘 机器学习理论 - 无监督学习算法 KMeans 基本原理 簇内误差平方和
大数据·人工智能·学习·算法·机器学习·数据挖掘
使者大牙2 小时前
【大语言模型学习笔记】第一篇:LLM大规模语言模型介绍
笔记·学习·语言模型
As977_3 小时前
前端学习Day12 CSS盒子的定位(相对定位篇“附练习”)
前端·css·学习
ajsbxi3 小时前
苍穹外卖学习记录
java·笔记·后端·学习·nginx·spring·servlet
Rattenking3 小时前
React 源码学习01 ---- React.Children.map 的实现与应用
javascript·学习·react.js
顶呱呱程序3 小时前
2-143 基于matlab-GUI的脉冲响应不变法实现音频滤波功能
算法·matlab·音视频·matlab-gui·音频滤波·脉冲响应不变法
dsywws3 小时前
Linux学习笔记之时间日期和查找和解压缩指令
linux·笔记·学习
道法自然04023 小时前
Ethernet 系列(8)-- 基础学习::ARP
网络·学习·智能路由器