首页 > 编程语言 > 详细

【Python】圆周率的计算

时间:2020-01-31 16:23:21      阅读:96      评论:0      收藏:0      [点我收藏+]

1.公式法

技术分享图片

 

 代码:

#CalPiV1.py
pi=0
N=100
for k in range(N):
    pi+=1/pow(16,k)*(        4/(8*k+1)-2/(8*k+4)-        1/(8*k+5)-1/(8*k+6))
print("圆周率的值是:{}".format(pi))    

技术分享图片

 

 2.蒙特卡洛方法:

撒点方法

取一个正圆和一个正方形的1/4形成一个单位方形

单位方形中有1/4个圆

圆的面积和单位方形之比即为圆周率的相关数据

然后进行撒点

技术分享图片

 

 

#CalPiV2.py
from random import random#调用random库中的random
from time import perf_counter#调用time库中的perf_counter,用来进行程序计时
DARTS=1000*1000#在当前区域中抛洒点的数量(数量越大,精度越高,就看你电脑的性能怎么样了!):100W
hits=0.0#目前在圆内部点的数量
start=perf_counter()#启动计时,初始值为当前系统的时间
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#撒点后圆内部的点的数量
pi=4*(hits/DARTS)#用hits的值比上全部撒点的值为1/4π
print("圆周率值是:{}".format(pi))
print("运行时间是:{:.5f}s".format(perf_counter()-start))

    

 技术分享图片

 

 技术分享图片

 

 

技术分享图片

 

【Python】圆周率的计算

原文:https://www.cnblogs.com/HGNET/p/12245701.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!