(Monte—Carlo)方法求圆周率π.docVIP

蒙特卡洛（Monte—Carlo）方法求圆周率π

蒙特卡罗(MonteCarlo)方法，或称计算机随机模拟方法，是一种基于“随机数”的计算方法。设表示针的中点到最近的一条平行线的距离，表示针与平行线的夹角（如图1所示），如果，或时，针与一条直线相交。

图1蒙特卡洛法求π示意图

由于向桌面投针是随机的，所以用来确定针在桌面上位置的是二维随机向量。并且在上服从均匀分布，在上服从均匀分布,与相互独立。由此可以写出的联合概率密度函数：

于是，所求概率为：

由于最后的结果与π有关，因此有些人想利用它来计算π的值。其方法是向桌面投针n次，若针与直线相交k次，则针与直线相交的频率为，以频率代替概率，则有，所以。

namespace是指标识符的各种可见范围

srand函数是随机数发生器的初始化函数。

原型：voidsrand(unsignedseed);

用法：它需要提供一个种子，这个种子会对应一个随机数，如果使用相同的种子后面的rand()函数会出现一样的随机数。如：srand(1);直接使用1来初始化种子。不过为了防止随机数每次重复常常使用系统时间来初始化，即使用time函数来获得系统时间，它的返回值为从00:00:00GMT,January