代码之家 › 专栏 › 技术社区 › Dessie

如何提高此python代码中的计算精度

precision python

Dessie · 技术社区 · 2 年前

我正试图重现《小数学快报》第81页上的一个情节(顺便说一句,这是一本很有趣的数学书)。该图是使用一对迭代应用的简单函数创建的(即,每个循环的输出是下一个循环的输入)。在每个循环中绘制输出x,y点应该会产生一个美丽的复杂点云。然而,当我实现下面的代码时,在第五十次迭代后,它会在4个相同的点之间循环。我认为这是由于我的变量定义精度不够而导致的四舍五入。我该如何改进?

import matplotlib.pyplot as plt

x = 0
y = 0
vecx = []
vecy = []

for i in range(5000):
    x = 1.3 + 0.3*x + 0.6*x*y - 0.6*y - y**2
    y = 0.1 - 0.7*x + 0.5*(x**2) - 0.8*x*y + 0.1*y - 0.6*(y**2)
    vecx.append(x)
    vecy.append(y)
    
plt.plot(vecx, vecy, 'b.')

3 回复 | 直到 2 年前

Navkar Jain 2 年前

我尝试了100次迭代,并且工作正常:

x = 0
y = 0
vecx = []
vecy = []

for i in range(100):
    new_x = 1.3 + 0.3*x + 0.6*x*y - 0.6*y - y**2
    new_y = 0.1 - 0.7*x + 0.5*(x**2) - 0.8*x*y + 0.1*y - 0.6*(y**2)
    
    x, y = new_x, new_y  # Update x and y for the next iteration        

    vecx.append(x)
    vecy.append(y)

    print('x', x)
    print('y', y)

chubercik 2 年前

您可以使用 decimal Python模块:

from decimal import Decimal, getcontext

# Set the precision for decimal calculations
getcontext().prec = 50  # You can adjust the precision as needed

x = Decimal('0')
y = Decimal('0')
vecx = []
vecy = []

for i in range(5000):
    x_new = Decimal('1.3') + Decimal('0.3') * x + Decimal('0.6') * x * y - Decimal('0.6') * y - y**2
    y_new = Decimal('0.1') - Decimal('0.7') * x + Decimal('0.5') * (x**2) - Decimal('0.8') * x * y + Decimal('0.1') * y - Decimal('0.6') * (y**2)
    x, y = x_new, y_new
    vecx.append(x)
    vecy.append(y)

-1

dawid.pawlowski 2 年前

您应该尝试使用十进制库进行精确的浮点数字计算。 from decimal import Decimal 并以字符串形式提供数字 y = Decimal('1.005')

推荐文章

BlurKid · R中for循环时结果的奇怪差异

1 年前

SapereAude · 在Python中,为什么numpy数组的预分配无法限制其打印精度?

1 年前

RuralAnemone · 需要帮助了解Google Sheets百分比

2 年前

Dakolas · 对python十进制库感到困惑

2 年前

IIXenon · for循环的printf仅跳过0[已关闭]

2 年前

despinxz · MIPS:如何在使用双打时保持精度或将其四舍五入到特定的小数位数?

2 年前

pmor · 如何计算转换规范%.*f的精度以保持浮点值的精度?

2 年前

Dessie · 如何提高此python代码中的计算精度

2 年前

ramgorur · 字符串流、setw和setprecision的C++困难[重复]

3 年前

Kevin Spaghetti · 使用numpy.sin时出现精度问题

3 年前