python如何计算π

发布网友发布时间：2022-04-06 04:13

我来回答

共3个回答

懂视网时间：2022-04-06 08:35

python利用公式计算π的方法：首先导入数学模块及时间模块；然后计算Pi精确到小数点后几位数，代码为【print(' {:=^70}'.format('计算开始'))】；最后完成计算，代码为【print(' {:=^70}'】。

【相关学习推荐：python教程】

python利用公式计算π的方法：

一、π的简介

π的介绍

　　圆周率用希腊字母 π（读作pài）表示，是一个常数（约等于3.141592654），是代表圆周长和直径的比值。它是一个即无限不循环小数，在日常生活中，通常都用3.14代表圆周率去进行近似计算。

π的求解历程

　　1965年，英国数学家约翰·沃利斯（John Wallis）出版了一本数学专著，其中他推导出一个公式，发现圆周率等于无穷个分数相乘的积。

　　2015年，罗切斯特大学的科学家们在氢原子能级的量子力学计算中发现了圆周率相同的公式。

　　2019年3月14日，谷歌宣布圆周率现已到小数点后31.4万亿位。

　　此处用一个自我感觉‘良好’的公式进行求解，说良好是因为计算结果相对准确，但计算过程用时较长，一起来学习吧~~~

二、π的近似计算

　　1. 计算公式

　　2. 方法讲解

　　所用公式等式右边分子都为1，分母为递增数列，从第一项开始，奇数项符号为正，偶数项符号为负。等式右边的分母越大，越小，圆周率π计算的值越精确；换个角度讲，就是等式右边的项越多，计算的值越精确。

3. 代码实现(python)

 1 from math import fabs  #导入数学模块
 2 from time import perf_counter #导入时间模块
 3 
 4 def Bar(i):  #动态文本条
 5 N = pow(10,level)
 6 a = int((i/N)*50)
 7 b = 50 - a
 8 Y , N = '*' * a , '.' * b
 9 print("
计算中：{:3.0f}% [{}->{}] {:.2f}s"
10  .format(2*a,Y,N,perf_counter()),end='')
11 
12 level = eval(input('计算Pi精确到小数点后几位数：'))
13 print('
{:=^70}'.format('计算开始'))
14 a,b,pi,tmp = 1,1,0,1
15 i = 0
16 '''
17 a 分子 | b 分母 | pi 圆周率
18 tmp 存储a/b的值 | i 执行进度
19 '''
20 perf_counter() #开始计时
21 while (fabs(tmp) >= pow(10,-level)): #计算Pi
22 pi += tmp
23 b += 2
24 a = -a
25 tmp = a/b
26 i += 2
27 Bar(i)  #调用函数，实时显示计算进度
28 
29 print('
{:=^70}'.format('计算完成'))
30 print('
Pi的计算值为：{}'.format(round(pi*4,level))) #输出计算结果

4. 图片示例

　　由上面3张图片可知，精确到小数点后4位只要14.07秒，精确到小数点后6位也需要124.61秒，而精确到小数点后8位就需要 850 / 8% = 10625秒，约为 177 分钟，也就是2.95个小时。这种方法固然好，但计算起来还是需要很长一段时间的。

热心网友时间：2022-04-06 05:43

#coding=utf-8
'''
Created on 2014-11-04

@author: Neo
'''

import sys
import math
from decimal import *

def bbp(n):
    pi=Decimal(0)
    k=0
    while k < n:
        pi+=(Decimal(1)/(16**k))*((Decimal(4)/(8*k+1))-(Decimal(2)/(8*k+4))-(Decimal(1)/(8*k+5))-(Decimal(1)/(8*k+6)))
        k+=1
    return pi

def main(argv):
    if len(argv) !=2:
        exit('Usage: BaileyBorweinPlouffe.py <prec> <n>')

    getcontext().prec=(int(sys.argv[1]))
    my_pi=bbp(int(sys.argv[2]))
    accuracy=100*(Decimal(math.pi)-my_pi)/my_pi

    print "Pi is approximately "+str(my_pi)
    print "Accuracy with math.pi: "+str(accuracy)

if __name__=="__main__":
    main(sys.argv[1:])

result:

d:\workspace\PyDemo>python test.py 10 10

Pi is approximately 3.141592653

Accuracy with math.pi: 1.877369797E-8

d:\workspace\PyDemo>python test.py 25 25

Pi is approximately 3.141592653589793238462644

Accuracy with math.pi: -3.898171852150198570978563E-15

d:\workspace\PyDemo>python test.py 40 40

Pi is approximately 3.141592653589793238462643383279502884195

Accuracy with math.pi: -3.898171832519375446564294863483061824361E-15

d:\workspace\PyDemo>

热心网友时间：2022-04-06 07:01

#马青公式
n = int(input('请键入想要计算到小数点后的位数n:'))
w = n+10
b = 10**w
x1 = b*4//5
x2 = b// -239
he = x1+x2
n *= 2
for i in range(3,n,2):
    x1 //= -25
    x2 //= -57121
    x = (x1+x2) // i
    he += x
    pai = he*4
    pai //= 10**10
print(pai)