Question

下面是有限差分法的一个函数，它是一种非常标准的计算导数的方法，它给出了一些函数f（x），一个网格（np.linspace）以及每一片之间的均匀距离。格（h）。

遇到的问题是当我尝试从0,10穿过网格的已知函数（例如x ** 2）时；我收到特定错误。在代码之后，我将发布遇到的错误。

def finitedifference(f,x,h,n):
"""f : function you are attempting to differentiate.
   x : grid/domain with with you will differentiate.
   h : distance between uniform mesh.
   n : required for loop?"""
    df = np.zeros_like(x)
    for i in range(1,n):
        df[i] = (f[i+1]-f[i-1])/(2*h)
        #end_points
        df[0] = (f[1]-f[0])/h
        df[-1] = (f[-1]-f[-2])/h
    return print(df)

我用什么

f = x**3
x = np.linspace(0,10,11)
h = x[1] - x[0]
finitedifference(f,x,h,11)

我收到错误：

"IndexError: index 11 is out of bounds for axis 0 with size 11"

不幸的是，我不确定这意味着什么，所以也许可以澄清一下错误/补救措施？谢谢！

Answer 1

大小为11的数组（或列表）的索引为0, 1, 2, ..., 10。

如果您在for上进行range(1, 11)循环，则会在1, 2, 3, ..., 10上进行迭代。如果然后尝试访问索引i+1，则当i到达10时，您将最终超出数组。

另一方面，您省略了索引0。因此，您可能希望使用for i in range(n-1)来确保从头开始并保持在限制范围内。

有限差分函数索引错误：

1 个答案: