附件是一个一阶求导或偏导数的代码,有需要的可以参考学习一下。
'一阶中心差商求函数的导数或偏导数
'背景:当工程计算时,有时候由于原函数形式过于复杂从而很难推导出其偏导数的解析形式,从而需要采用插值解法对其进行求导
'然而,利用中心差商时有时候,由于Δx过小,从而无法识别出f(x-Δx)和f(x+Δx)的差异也非常的小,导致求解精度不高甚至无法求导
'以下给出中心差商的迭代形式,效果还算不错
'以EXP(x)函数为例,其在x=1处的导数的精确解为2.71828182845905
'以下方法在Epsilon=1E-06时解为 2.71828200543069,已具备很高的精度,若期望更精确的解只需提高Epsilon的精度即可
附 件:
点击下载此附件学习源代码