我使用以下代码计算一阶导数:
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
变得很大。
这是为什么?这是怎么回事
本文收集自互联网,转载请注明来源。
如有侵权,请联系 [email protected] 删除。
我来说两句