Python数值微分和h的最小值

我使用以下代码计算一阶导数：

def f(x):
   f = np.exp(x)
   return f

def dfdx(x):
   Df = (f(x+h)-f(x-h)) / (2*h)
   return Df

例如,对于x == 10,这很好.但当我将h设置为10E-14或更低时,Df开始
获得真正远离预期值f(10)的值,并且期望值和Df之间的相对误差变得很大.

这是为什么？这里发生了什么？

解决方法:

您可能遇到了一些数值不稳定性,对于x = 10和h = ~1E-13,np.exp的参数非常接近10,无论h是加还是减,因此np.exp值的小近似误差通过非常小的2 * h的划分显着缩放.

标签：python,numerical-methods
来源： https://codeday.me/bug/20190607/1193579.html