蒙特卡罗法求PI

蒙特·卡罗方法(Monte Carlo method),也称统计模拟方法,是二十世纪四十年代中期由于科学技术的发展和电子计算机的发明,而被提出的一种以概率统计理论为指导的一类非常重要的数值计算方法。是指使用随机数(或更常见的伪随机数)来解决很多计算问题的方法。与它对应的是确定性算法。为象征性地表明这一方法的概率统计特征,故借用赌城蒙特卡罗命名。

对于求圆周率PI的问题,可以模拟在一个长宽均为1的正方形里面均匀的随机落入点,然后根据落在半径为1的圆中的点的个数来求得Pi的大小,点的数目越多,结果越准确。


#include <iostream>
#include <time.h>
using namespace std;

#define  MAX_NUM  5000000

int g_nInCircle = 0;

int main()
{
	double nX;
	double nY;
	double dDis = 0;
	double Pi = 0;
	srand(time(NULL));

	for (long int i=0; i< MAX_NUM; i++)
	{
		nX = (double)rand()/RAND_MAX;  //可能生成都最大随机数
		nY = (double)rand()/RAND_MAX;

		dDis = pow(nX, 2) + pow(nY, 2);
		if (dDis < 1.0)
		{
			g_nInCircle++;

		}
	}
	Pi =(double) g_nInCircle / MAX_NUM * 4;
	cout<<"Pi = "<<Pi<<endl;

	system("pause");

}


评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值