python 图中的大舍入误差

Question

我尝试 plot 下面的简单序列

a_n=\frac{3^n+1}{7^n+8}

它应该趋于 0，但是 plot 对接近 20 的 $n$ 值显示出奇怪的效果....

我用代码

import numpy as np
import matplotlib.pyplot as plt

def f(n):
    return (3**n+1)/(7**n+8)

n=np.arange(0,25, 1)

plt.plot(n,f(n),'bo-')

另一方面，对上述序列进行数值计算并没有找到这么大的值

for i in range(0,25):
    print([i,f(i)])

[0, 0.2222222222222222]
[1, 0.26666666666666666]
[2, 0.17543859649122806]
[3, 0.07977207977207977]
[4, 0.034039020340390205]
[5, 0.014510853404698185]
[6, 0.0062044757218015075]
[7, 0.0026567874970706124]
[8, 0.0011382857610720493]
[9, 0.00048778777316480816]
[10, 0.00020904485804220367]
[11, 8.958964415487241e-05]
[12, 3.8395417418579486e-05]
[13, 1.6455158259653074e-05]
[14, 7.05220773432529e-06]
[15, 3.022374322043928e-06]
[16, 1.295303220696569e-06]
[17, 5.551299431298911e-07]
[18, 2.3791283154177113e-07]
[19, 1.0196264191387531e-07]
[20, 4.3698275080881505e-08]
[21, 1.872783217393992e-08]
[22, 8.026213788319863e-09]
[23, 3.439805909206865e-09]
[24, 1.4742025325067883e-09]

为什么会这样？

Answer 1

问题不matplotlib ，而在于arange生成的数字的数据类型。 您没有指定dtype ，因为在arange的文档中，它声明是从输入中推断出来的。 你的输入是整数，所以它必须假设它们是 32 位整数，因为dtype是未修改的，所以当我检查类型时：

print(type(n[0]))
<class 'numpy.int32'>

如果我将 dtype 更改为单精度浮点数，我们会得到您期望的行为：

n = np.arange(0,25,1, dtype=np.float32)

print(type(n[0]))
<class 'numpy.float32'>

plt.plot(n,f(n),'bo-')

或者，您可以在1 -> 1.后面放一个句点来暗示您需要双精度浮点数（即使结果数组包含整数式数字 [0., 1., 2., ...]）