前几天学习了嵩天老师的一节挺有意思的课,讲的是通过蒙特卡罗方法来计算圆周率。蒙特卡罗方法在这里就不多描述了,有兴趣的朋友可以看一下。蒙特卡罗方法
主要思想是:在一个正方形中按照r/2为半径画圆,那么圆的面积跟单位正方形面积之比就是圆周率的相关数据。
代码如下:
# -*- coding: utf-8 -*-
import random
import time
num = 1000*1000 #撒点总数为num
count = 0
start = time.perf_counter() #通过 time.perf_counter()函数来计算程序运行时间
for i in range(1,num+1):
x,y = random.random(),random.random() #随机生成两个个0-1之间的数
length = pow(x*x+y*y,0.5) #得出某个点距离圆心(坐标轴原点)的距离
if length <= 1:
count += 1
pi = 4 * (count/num)
print("圆周率={}".format(pi))
print("运行时间共计{:.5f}s".format(time.perf_counter()-start))