第五十一篇 有限差分法
在有限差分方法中,出现在偏微分方程中的导数被“有限差分”的值近似取代。由于偏微分方程包含了至少两个自变量,有限差分表达式包括了在要求导数的周围的“网格点”上的未知变量的各种组合。
在下图中考虑一个二维笛卡尔解域,它被分割成一个规则的矩形网格,x方向上间距为h, y方向上间距为k。在这种情况下,单个因变量u是两个自变量x和y的函数,需要求导的点是下标i j,其中i在x方向上计数,j在y方向上计数。
正向差分导数公式:
中心差分导数公式:
逆向差分导数公式:
任何偏导数都可以用有限差分表达式来近似。例如,
是一阶导数的中心差分公式。这些公式有时很容易地表示为计算“分子”,例如
结合这两个一阶导数,让h = k,得到
二阶导数也是一样
在h = k的特殊情况下,会得到拉普拉斯式的“分子”
三、四阶导数的中心差分“分子”可以写成这样的形式
“混合”导数(h = k)为
双调和算子(h = k)等于
所有这些例子都使用了最低阶中心差分形式。显然,如果需要,可以设计更高阶的版本,包括正向和逆向的差分。在可能的情况下,最好使用中心差分公式,因为它的准确度较高一些。
当使用这些中心差分公式时,第(i, j)个网格点总是位于“分子”的中间。这是一个近似导数的点,“分