我使用以下代码计算一阶导数:
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
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。