2.4　有限差分法求解偏微分方程

自然科学与工程技术中种种运动发展过程与平衡现象各自遵守一定的规律。这些规律的定量表述一般地呈现为关于含有未知函数及其导数的方程。只含有未知多元函数及其偏导数的方程，称为偏微分方程。初始条件和边界条件称为定解条件，未附加定解条件的偏微分方程称为泛定方程。对于一个具体的问题，定解条件与泛定方程总是同时提出。定解条件与泛定方程作为一个整体，称为定解问题。

2.4.1　椭圆方程

由于大多数工程问题都是二维问题，所以得到的微分方程一般都是偏微分方程，对于一维问题得到的是常微分方程，解法与偏微分方程类似，为了不是一般性，这里只讨论偏微分方程。由于工程中高阶偏微分较少出现，所以本节仅仅给出二阶偏微分方程的一般形式，对于高阶的偏微分，可进行类似地推广。二阶偏微分方程的一般形式如下

A Φ_x，_x+BΦ_xy+CΦ_yy=f（x，y，Φ，Φ_x，Φ_y）

其中，Φ表示一个连续函数。当A，B，C都是常数时，上式成为准线性，有三种准线性方程形式：

（1）如果Δ=B²-4AC＜0，则称为椭圆型方程。

（2）如果Δ=B²-4AC=0，则称为抛物型方程。

（3）如果Δ=B²-4AC＞0，则称为双曲型方程。

椭圆方程是工程技术应用中所涉及的偏微分方程里最为普遍的一种形式。根据椭圆方程的具体形式又可以将其分为以下三种形式：

（1）拉普拉斯（Laplace）方程：=0；

（2）泊松（Poisson）方程：=g（x，y）；

（3）亥姆霍兹（Helmholtz）方程：+λu=0。

其中，u是关于x和y的二元函数。

2.4.2　有限差分法

差分方法又称为有限差分方法或网格法，是求偏微分方程定解问题的数值解中应用最广泛的方法之一。它的基本思想是：先对求解区域作网格剖分，将自变量的连续变化区域用有限离散点（网格点）集代替；将问题中出现的连续变量的函数用定义在网格点上离散变量的函数代替；通过用网格点上函数的差商代替导数，将含连续变量的偏微分方程定解问题化成只含有限个未知数的代数方程组（称为差分格式）。如果差分格式有解，且当网格无限变小时其解收敛于原微分方程定解问题的解，则差分格式的解就作为原问题的近似解（数值解）。因此，用差分方法求偏微分方程定解问题一般需要解决以下问题：

（1）选取网格；

（2）对微分方程及定解条件选择差分近似，列出差分格式；

（3）求解差分格式；

（4）讨论差分格式解对于微分方程解的收敛性及误差估计。

下面就以拉普拉斯方程的数值解法为例来演示一下有限差分法的基本思路。首先写出完整的拉普拉斯方程如下

现在的问题其实是要求在一个给定的二维区域中求解满足方程的每一点（x，y）。一些区域中的点将被用来给出边界条件（hold boundary conditions）。

于是将整个二维区域离散化成若干个点，其中的5个相邻点如图2-3所示。

根据偏导数的定义则有

同理可得

将上述两个结果代入拉普拉斯方程可得

2.4.3　方程组求解

回想雅可比迭代法，假设有一个由n个线性方程组成的系统（也就是线性方程组）

Ax=b

那么雅可比迭代可以描述为

其中，k表示第k轮迭代。

注意在2.4.2节最后得出的拉普拉斯方程离散化形式给出了（离散化后）区域上众多点中的一个点的求解方程，所有点的求解方程合在一起就构成了一个大的方程组。把求解某点（x，y）的方程重写成雅可比迭代的形式，则有

重复应用上述迭代式，最后方程就会收敛到解的附近。

本来连续的一个区域经过离散化处理之后就变成了一个网格结构，假设网格的大小是n²，标签为x₁，x₂，…，如图2-4所示。

上面这种自然排列的点序可以得出不超过n²个五元线性方程

x_i-n+x_i-1-4x_i+x_i+1+x_i+n=0

注意，一般不对处于边界上的点（如x₁，x₂，x₃，x₄，x₅，x₈，…）应用上述方程。最后将得到一个大型的稀疏线性方程组

除了使用雅可比迭代之外，还可以采用高斯迭代加速方程组解的收敛速度。高斯迭代的一般形式为

此时，拉普拉斯方程需用下式进行求解，即

偏微分方程在图像处理中有非常重要的应用。本节主要是以拉普拉斯方程的数值解为例来讨论的，而前面也提到过椭圆方程中除了拉普拉斯方程之外，还有一类叫做泊松方程。图像处理中基于泊松方程的算法构成了一大类的具有广泛应用的算法，可以用于图像融合、图像去雾、图像拼接等。例如，图2-5就是基于解泊松方程的方法实现的图像泊松编辑的效果图。本书后面还会详细介绍这种算法的原理。