Question

我创建了一个简单的python程序，通过计算适合圆圈的1个像素点的数量来计算pi。起初我使用以下值运行程序：

from tkinter import *

tk = Tk()
canvas = Canvas(tk, width=400, height=400)
canvas.pack()

canvas.create_oval(100,100,400,400)

global pi_count
pi_count = 0

for x in range(100,400):
    for y in range(100,400):
        point = canvas.create_oval(x,y,(x+1),(y+1))
        c = canvas.coords(point)
        sc = list(c)
        xc = sc[1]
        yc = sc[2]
        fxc = int(xc)
        fyc = int(yc)
        dist = ((fxc - 250)**2 + (fyc - 250)**2)**0.5
        if dist > 150:
            print("outside of the circle")
        else:
            pi_count += 1
        print(pi_count)

pi = (pi_count/150**2)
print(pi)

这给了我输出3.1412888889。

对此结果比较满意，我将椭圆尺寸更改为 (100,100,500,500) 并相应地更改了所有其他值（半径，循环等）。

然而，这个较大的圆圈和可能更精确的区域产生了更加不准确的估计值3.140675。

为什么会这样，以及如何修复计算器以提供更准确的估算？

修改

以下是经过编辑和改进的代码，更易于测试：

from tkinter import *

MIN_POS = 100
MAX_POS = 300
center = (MAX_POS + MIN_POS)/2
radius = (MAX_POS - MIN_POS)/2

tk = Tk()
canvas = Canvas(tk, width=MAX_POS, height=MAX_POS)
canvas.pack()

canvas.create_oval(MIN_POS,MIN_POS,MAX_POS,MAX_POS)

global pi_count
pi_count = 0

for x in range(MIN_POS,MAX_POS):
    for y in range(MIN_POS,MAX_POS):
        point = canvas.create_oval(x,y,(x+1),(y+1))
        c = canvas.coords(point)
        sc = list(c)
        xc = sc[1]
        yc = sc[2]
        fxc = int(xc)
        fyc = int(yc)
        dist = ((fxc - center)**2 + (fyc - center)**2)**0.5
        if dist > radius:
            print("outside of the circle")
        else:
            pi_count += 1
        print(pi_count)

pi = (pi_count/radius**2)
print(pi)

如果拥有快速计算机的任何人都可以运行此程序并且以{100}的步长提升Max_pos，我将不胜感激。

Answer 1

此方法非常不准确。测试一个半径的结果可能比1个或更少的结果更准确或更不准确：

100 -> 3.141700 ~ 1.000034
125 -> 3.140928 ~ 0.999788
150 -> 3.141378 ~ 0.999932
175 -> 3.141518 ~ 0.999976
200 -> 3.140725 ~ 0.999724
225 -> 3.141393 ~ 0.999936
250 -> 3.141136 ~ 0.999855
275 -> 3.140826 ~ 0.999756
300 -> 3.141078 ~ 0.999836
325 -> 3.141311 ~ 0.999910
350 -> 3.140939 ~ 0.999792
375 -> 3.141582 ~ 0.999997
400 -> 3.141406 ~ 0.999941

原因是您只测试积分值，因此这些值往往会向下舍入。虽然增加半径确实提高了准确性，但是在你超越Python的浮点精度本身的限制之前，应该不可能得到足够小数的π - 我不知道半径需要多大。

对于您的价值150，测试149和151以及显示两者都“更好”; 两个都差得多！

150 -> 3.141378 ~ 0.999932
 -1 -> 3.140804 ~ 0.999749
 +1 -> 3.140783 ~ 0.999742

即使对于像10,000这样的非常大的半径，你仍然得到3.141591~0.999999（需要一段时间来计算）。

使用以下代码进行测试（不使用图形显示）。请注意，代码中存在一个小错误，这也会影响结果！您运行(100,400)范围，您可以将其写为(0,300)。减去半径，你得(-150,150) - 但Python的range函数从start to less than stop运行。这意味着您正在计算-150..-1，然后计算0，然后计算1..149，并且左侧有一个小偏差。

（但是，考虑到这一点并不能明显提高准确性。）

import math

def calc_pi(radius):
    pi_count = 0
    rad2 = radius*radius
    for x in range(-radius,radius+1):
        x2 = x*x
        for y in range(-radius,radius+1):
            dist = x2 + y*y
            if dist <= rad2:
                pi_count += 1
    return float(pi_count)/radius**2

for i in range(100,401,25):
    result = calc_pi(i)
    print("%d -> %f ~ %f" % (i, result, result/math.pi))

Pi计算器输出问题

1 个答案: