Python之蒙特卡罗法求圆周率

个人感觉把蒙特卡罗法理解了,对于理解机器学习的有关知识很有帮助

蒙特卡罗法--4*(圆内点数/总点数)=pi

from random import random
import time
DARTS=1000*1000 #在正方形区域中总共撒下的点数
hits=0.0  #定义一个原点
start=time.perf_counter()#计时函数
for i in range(1,DARTS+1): 
    x, y=random(),random()#生成两个随机数,模拟向方形区域撒点。x,y是该点的坐标
    distance=pow((x**2+y**2),0.5)#求出该点到原点的距离
    if(distance<=1.0):#距离<=1说明这些点落在了圆形区域内部
        hits=hits+1
pi=4*(hits/DARTS)#利用公式求出近似值pi
print("the value of pi is:{}".format(pi))
print("run time is:{:.5}s".format(time.perf_counter()-start))

 

你可能感兴趣的:(python)