发帖

【优惠升级】华秋PCB首单最高立减100元，SMT免费贴片！！！

[经验]

基于蒙特卡洛法算圆的面积

2014-9-10 15:16:26

6263

1

这是一个比较简单的入门程序，是通过概率来算面积的比如说一个面积为1的正方形如果里面有n1个点，一个圆里面有n2个点，则圆的面积是n2/n1.大家可以试验一下。
%蒙特卡洛法求面积
sita=0:0.01:2*pi;
x=sin(sita);
y=cos(sita)
m=0
x1=2*rand(1000,1)-1
y1=2*rand(1000,1)-1;
N=1000
for n=1:
p1=x1(1:n);
q1=y1(1:n);
if (x1(n)*x1(n)+y1(n)*y1(n))<1
m=m+1
end
plot(p1,q1,'.',x,y,'-k',[-1 -1 1 1 -1],[-1 1 1 -1 -1],'-k');
axis equal
axis([-2 2 -2 2])
text(-1,-1.2,['試驗總次數 n=',num2str(n)])
text(-1,-1.4,['落入圓中數 m=',num2str(m)]);
text(-1,-1.6,['近似圓面積 S_c=',num2str(m/n*4)]);
set(gcf,'DoubleBuffer','on')
drawnow
end

只看该作者 淘帖0 举报

相关推荐

评论

声明：本文内容及配图由入驻作者撰写或者入驻合作网站授权转载。文章观点仅代表作者本人，不代表电子发烧友网立场。文章及其配图仅供工程师学习之用，如有内容图片侵权或者其他问题，请联系本站作侵删。侵权投诉

精选推荐

请问simulink的s-function模块如何添加多输入输出接口

1349 浏览 2 评论
Matlab中如何将不同的时间的数据进行同时间对比？

1184 浏览 3 评论
使用simulink进行三相短路故障分析时，各参数应该如何设置

1746 浏览 1 评论
想请教一下图中是simulink的什么模块，需要这种三段斜率函数模块但没找到在哪

1823 浏览 1 评论
怎么修改matlab与comsol联合的代码啊

2744 浏览 1 评论

快速回复 返回顶部 返回列表

关注微信公众号

电子发烧友网

电子发烧友论坛

社区合作: 刘勇; 联系电话：15994832713; 邮箱地址：liuyong@huaqiu.com

社区管理: elecfans短短; 微信：elecfans_666; 邮箱：users@huaqiu.com

【优惠升级】华秋PCB首单最高立减100元，SMT免费贴片！！！

返回 Matlab论坛

回复

关闭

站长推荐 /7

快速回复 返回顶部 返回列表

- 威廉希尔官方网站社区: HarmonyOS威廉希尔官方网站社区

RISC-V MCU威廉希尔官方网站社区

FPGA开发者威廉希尔官方网站社区

- OpenHarmony开源社区: OpenHarmony开源社区

- 嵌入式论坛: ARM威廉希尔官方网站论坛

STM32/STM8威廉希尔官方网站论坛

嵌入式威廉希尔官方网站论坛

单片机/MCU论坛

RISC-V威廉希尔官方网站论坛

瑞芯微Rockchip开发者社区

FPGA|CPLD|ASIC论坛

DSP论坛

- 电路图及DIY: 电路设计论坛

DIY及创意

电子元器件论坛

专家问答

- 电源威廉希尔官方网站论坛: 电源威廉希尔官方网站论坛

无线充电威廉希尔官方网站

- 综合威廉希尔官方网站与应用: 机器人论坛

USB论坛

电机控制

模拟威廉希尔官方网站

音视频威廉希尔官方网站

综合威廉希尔官方网站交流

上位机软件（C/Python/Java等）

- 无线通信论坛: WIFI威廉希尔官方网站

蓝牙威廉希尔官方网站

天线|RF射频|微波|雷达威廉希尔官方网站

- EDA设计论坛: PCB设计论坛

DigiPCBA论坛

Protel|AD|DXP论坛

PADS威廉希尔官方网站论坛

Allegro论坛

multisim论坛

proteus论坛|仿真论坛

KiCad EDA 中文论坛

DFM|可制造性设计论坛

- 测试测量论坛: LabVIEW论坛

Matlab论坛

测试测量威廉希尔官方网站

传感威廉希尔官方网站

- 招聘/交友/外包/交易/杂谈: 项目外包

供需及二手交易

工程师杂谈|交友

招聘|求职|工程师职场

- 官方社区: 发烧友官方/活动

华秋商城

华秋电路

time

recommend

hot

post

—
—
—

版
块
导
航