sboyer的个人博客分享 http://blog.sciencenet.cn/u/sboyer

博文

蒙特卡洛方法——求圆周率

已有 9545 次阅读 2010-11-30 22:10 |个人分类:科研笔记|系统分类:科研笔记| 蒙特卡洛

半径为1的圆(假定圆心为坐标原点)面积为(以下用pi表示),作圆外接正方形(各边与坐标轴平行),则正方形面积为4。设(X,Y)服从该正方形中的二维均匀分布,则随机分布点位于圆内的概率P=pi/4,于是pi=4*P,对于概率P可用蒙特卡洛方法模拟计算,下面给出一个简例。

n=input('请输入随机数个数');
x=2*(rand(1,n)-0.5);
y=2*(rand(1,n)-0.5);
r=sqrt(x.^2+y.^2);
r(r<=1)=1;
r(r>1)=0;
pai=4*sum(r)/n;



https://blog.sciencenet.cn/blog-114642-389035.html

上一篇:博士生3.0
下一篇:机器人相关期刊
收藏 IP: .*| 热度|

0

发表评论 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-12-31 01:27

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部