python 圆周率,用python计算圆周率π

用python计算圆周率π

一、要求

1.要起义能计算到圆周率后面越多位越好。

2.用进度条显示计算的进度，能给出越多种进度条越好。

3.要求给出圆周率π的具体公式或者算法说明。

二、算法

1、马青公式：马青公式每计算一项可以得到1.4位的十进制精度。

2、蒙特卡罗法(我们使用这种方法)

一个正方形内部相切一个圆，圆和正方形的面积之比是π/4。

在这个正方形内部，随机产生n个点(这些点服从均匀分布)，计算它们与中心点的距离是否大于圆的半径，以此判断是否落在圆的内部。

统计圆内的点数，与n的比值乘以4，就是π的值。理论上，n越大，计算的π值越准

三、进度条

1、time库：

(1)perf_counter( )函数，用于测量时间，返回一个CPU级别的精确时间计数值，单位为妙，由于这个计数值起点不确定，连续调用使用差值才有意义。

(2)sleep(s)函数，产生 s妙的休眠时间，可以是浮点数，如time.sleep(3.5)。

(3)time.clock()函数，统计cpu时间 的工具，这在统计某一程序或函数的执行速度最为合适。两次调用time.clock()函数的插值即为程序运行的cpu时间。

2、tqdm库：Tqdm 是一个快速，可扩展的Python进度条，可以在 Python 长循环中添加一个进度提示信息，用户只需要封装任意的迭代器 tqdm(iterator)。

步骤如下：

(1)安装tqdm

pip install tqd

(2)导入tqdm库

from tqdm import tqdm

(3)tqdm示例

for i in tqdm(range(1000000)):

pass

四、算法实现

1、time库

from random import random

from time import clock

DARTS=100000

hits = 0.0

a=1

clock()

for i in range(1,DARTS+1):

x,y = random(),random()

dist = pow(x ** 2 + y ** 2,0.5)

if dist<=1.0:

hits = hits + 1

if i== DARTS*0.01*a :

print("\r%{} [{}->{}]".format(a,'*'*a,'-'*(10-a)),end="")

a+=3

pi = 4* (hits/DARTS)

print("\nPi=={:.7f}".format(pi))

print("运行时间为:{:.4f}s".format(clock()))

效果如下：

2、tqdm库

from math import *

from tqdm import tqdm

from time import *

total,s,n,t=0.0,1,1.0,1.0

clock()

while(fabs(t)>=1e-6):

total+=t

n+=2

s=-s

t=s/n

k=total*4

print("π值是{:.10f} 运行时间为{:.4f}秒".format(k,clock()))

for i in tqdm(range(101)):

print("\r{:3}%".format(i),end="")

sleep((clock())/100) #用执行程序的总时间来算出进度条间隔的时间