|||
半径为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;
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-12-4 22:40
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社