我正在使用以下代码绘制一个函数。
from __future__ import division
from scipy.misc import comb
import matplotlib.pyplot as plt
import numpy as np
N = 100
X = np.arange(2,N)
def k_loop(w,n):
K = np.arange(0, w+1)
return (comb(w,K)*(comb(w,K)/2**w)**(3*float(n)/np.log(n))).sum()
def w_loop(n):
v = [comb(n,w)*k_loop(w,n) for w in range(1,n+1)]
return sum(v)
Y = [w_loop(n) for n in X]
plt.plot(X,Y)
plt.show()
但是我想确定没有数字错误。是否有可能让 python 准确地进行计算并且只在绘图时的最后一刻转换为 float ?日志可能不准确,但我希望尽可能保持其他一切准确。
最佳答案
如果您希望使用以 10 为基数的数字进行尽可能精确的计算,请尝试使用 decimal模块。
不过,这将需要重写您的代码。
decimal
有一个ln
和log10
函数,以及一个power
函数。- 但您可能需要根据
decimal
重写comb
,或者使用exact=True
并将结果转换为十进制
。
关于python - 精确的计算和绘图,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/20593177/