圆周率使用 Python

Question

我是 Python 的初学者。我尝试用 Python 语言评估 pi。我在这里使用Archimedes method。

我的代码是：

import math
from decimal import *
getcontext().prec = 50

pi = 0
PolySides = 6
SideLen = 1
Perim = 6
Dia = 2

while PolySides < 10000000000:
    SideLen2 = Decimal(SideLen) / Decimal(2)
    radius_a = Decimal(math.sqrt(1 - SideLen2**2))
    radius_b = Decimal(1) - Decimal(radius_a)
    SideLenNew = Decimal(math.sqrt(SideLen2**2 + radius_b**2))
    PolyCircum = Decimal(PolySides * SideLen)
    pi = PolyCircum / Dia
    print('Polygon Sides:', PolySides, 'Pi = ', pi)
    SideLen = SideLenNew
    PolySides = 2 * PolySides

但出于某种原因，此代码将 pi 打印到小数点后 14 位。我不明白我在哪里失去了精度。谢谢你。

Answer 1

您在 float 和 decimal 之间穿梭。 math.sqrt 需要一个 float 和 returns 一个 float，所以会有转换，你的精度被截断为 float。

这些更改就足够了：

while PolySides < 10000000000000000000:
    SideLen2 = Decimal(SideLen) / 2
    radius_a = (1 - SideLen2**2).sqrt()
    radius_b = 1 - radius_a
    SideLenNew = (SideLen2**2 + radius_b**2).sqrt()
    PolyCircum = Decimal(PolySides * SideLen)
    pi = PolyCircum / Dia
    print('Polygon Sides:', PolySides, 'Pi = ', pi)
    SideLen = SideLenNew
    PolySides = 2 * PolySides

Answer 2

您必须使用 .sqrt() 小数法：

from decimal import *
getcontext().prec = 50

pi = 0
PolySides = 6
SideLen = 1
Perim = 6
Dia = 2

while PolySides < 1000000000000:
    SideLen2 = Decimal(SideLen) / Decimal(2)
    radius_a = (1 - SideLen2**2).sqrt()
    radius_b = Decimal(1) - Decimal(radius_a)
    SideLenNew = (SideLen2**2 + radius_b**2).sqrt()
    PolyCircum = Decimal(PolySides * SideLen)
    pi = PolyCircum / Dia
    print('Polygon Sides:', PolySides, 'Pi = ', pi)
    SideLen = SideLenNew
    PolySides = 2 * PolySides

输出：

Polygon Sides: 824633720832 Pi =  3.1415926535897932384626357839340399511415560167917

有 22 位小数。

Answer 3

确保变量 pi 本身也是一个 Decimal，并避免使用 @Amadan 所写的 math.sqrt()。使用 Decimal.sqrt()

import math
from decimal import *
getcontext().prec = 50

pi = Decimal(0)
PolySides = 6
SideLen = 1
Perim = 6
Dia = 2

while PolySides < 10000000000:
    SideLen2 = Decimal(SideLen) / Decimal(2)
    radius_a = Decimal((Decimal(1) - SideLen2**2).sqrt())
    radius_b = Decimal(1) - Decimal(radius_a)
    SideLenNew = Decimal((SideLen2**2 + radius_b**2).sqrt())
    PolyCircum = Decimal(PolySides * SideLen)
    pi = PolyCircum / Dia
    print('Polygon Sides:', PolySides, 'Pi = ', pi)
    SideLen = SideLenNew
    PolySides = 2 * PolySides

给出：

Polygon Sides: 6 Pi =  3
Polygon Sides: 12 Pi =  3.1058285412302491481867860514885799401888268158392
Polygon Sides: 24 Pi =  3.1326286132812381971617494694917362446497769154815
Polygon Sides: 48 Pi =  3.1393502030468672071351468212084211891503505893626
Polygon Sides: 96 Pi =  3.1410319508905096381113529264596601070364122161628
Polygon Sides: 192 Pi =  3.1414524722854620754506093089612256452476623045496
Polygon Sides: 384 Pi =  3.1415576079118576455164633451298595415043764795884
Polygon Sides: 768 Pi =  3.1415838921483184086689696037211533505200449157810
Polygon Sides: 1536 Pi =  3.1415904632280500957384585059309517235542823086758
Polygon Sides: 3072 Pi =  3.1415921059992715505447766406101173531274972549662
Polygon Sides: 6144 Pi =  3.1415925166921574475928740847688319059677188923680
Polygon Sides: 12288 Pi =  3.1415926193653839551895493120653190422221927255946
Polygon Sides: 24576 Pi =  3.1415926450336908966721415089192384127226230112419
Polygon Sides: 49152 Pi =  3.1415926514507676517042536404922190204484274723852
Polygon Sides: 98304 Pi =  3.1415926530550368416911231804154742022572768058866
Polygon Sides: 196608 Pi =  3.1415926534561041392646431596150783313543414874292
Polygon Sides: 393216 Pi =  3.1415926535563709636628233165541133642749135041839
Polygon Sides: 786432 Pi =  3.1415926535814376697626683659225751788703495073714
Polygon Sides: 1572864 Pi =  3.1415926535877043462876483788980471857502408672475
Polygon Sides: 3145728 Pi =  3.1415926535892710154188945540564999738004063788314
Polygon Sides: 6291456 Pi =  3.1415926535896626827017061710907747199859793791592
Polygon Sides: 12582912 Pi =  3.1415926535897605995224090799271347533561151459628
Polygon Sides: 25165824 Pi =  3.1415926535897850787275848074223367208751396830886
Polygon Sides: 50331648 Pi =  3.1415926535897911985288787393140192102034265840860
Polygon Sides: 100663296 Pi =  3.1415926535897927284792022222880574573760314838880
Polygon Sides: 201326592 Pi =  3.1415926535897931109667830930316368707217160322736
Polygon Sides: 402653184 Pi =  3.1415926535897932065886783107175360897801705020852
Polygon Sides: 805306368 Pi =  3.1415926535897932304941521151390111674024112028327
Polygon Sides: 1610612736 Pi =  3.1415926535897932364705205662443799538615730707255
Polygon Sides: 3221225472 Pi =  3.1415926535897932379646126790207221515422136434928
Polygon Sides: 6442450944 Pi =  3.1415926535897932383381357072148077010289894182968

圆周率使用 Python

Pi using Python

python

pi

python-3.x