有限差分近似求导

有限差分近似求导

有限差分法以变量离散取值后对应的函数值来近似微分方程中独立变量的连续取值。在有限差分方法中,我们放弃了微分方程中独立变量可以取连续值的特征,而关注独立变量离散取值后对应的函数值。但是从原则上说,这种方法仍然可以达到任意满意的计算精度。因为方程的连续数值解可以通过减小独立变量离散取值的间格,或者通过离散点上的函数值插值计算来近似得到。这种方法是随着计算机的诞生和应用而发展起来的。其计算格式和程序的设计都比较直观和简单,因而,它在计算数学中使用广泛。

有限差分法的具体操作分为两个部分:

\1. 用差分代替微分方程中的微分,将连续变化的变量离散化,从而得到差分方程组的数学形式;

\2. 求解差分方程组

首先我们直接给出结果,然后在后面进行公式推导:

有限差分近似求导

推导过程:

我们知道任一f(x)均可将其展开成泰勒公式:

一阶泰勒公式展开式为

有限差分近似求导
我们将f(xi+h)和f(xi-h)均用泰勒公式展开:

有限差分近似求导

将(1)和(2)相减,整理可得:

有限差分近似求导

有限差分近似求导

有限差分近似求导