第十章偏微分方程数值解法偏微分方程问题，其求解十分困难。除少数特殊情况外，绝大多数情况均难以求出精确解。因此，近似解法就显得更为重要。本章仅介绍求解各类典型偏微分方程定解问题的差分方法。§1差分方法的基本概念1.1几类偏微分方程的定解问题椭圆型方程：其最典型、最简单的形式是泊松（Poisson）方程),(2222yxfyuxuu特别地，当0),(yxf时，即为拉普拉斯（Laplace）方程，又称为调和方程02222yuxuuPoisson方程的第一边值问题为),(),(),(),(),(2222yxyxuyxyxfyuxuyx其中为以为边界的有界区域，为分段光滑曲线，称为定解区域，),(yxf，),(yx分别为，上的已知连续函数。第二类和第三类边界条件可统一表示为),(),(yxuuyxn其中n为边界的外法线方向。当0时为第二类边界条件，0时为第三类边界条件。抛物型方程：其最简单的形式为一维热传导方程220(0)uuaatx方程可以有两种不同类型的定解问题：初值问题xxxuxtxuatu)()0,(,0022初边值问题221200,0(,0)()0(0,)(),(,)()0uuatTxltxuxxxlutgtultgttT其中)(x，)(1tg，)(2tg为已知函数，且满足连接条件)0()(),0()0(21glg边界条件)(),(),(),0(21tgtlutgtu称为第一类边界条件。第二类和第三类边界条件为)()()()(22101tgutxutgutxulxxTt0其中0)(1t，0)(2t。当0)()(21tt时，为第二类边界条件，否则称为第三类边界条件。双曲型方程:最简单形式为一阶双曲型方程0xuatu物理中常见的一维振动与波动问题可用二阶波动方程22222xuatu描述，它是双曲型方程的典型形式。方程的初值问题为xxtuxxuxtxuatut)()()0,(,0022222边界条件一般也有三类，最简单的初边值问题为2222201200,0(,0)(),()0(0,)(),(,)()0tuuatTxltxuuxxxxltutgtultgttT1.2差分方法的基本概念差分方法又称为有限差分方法或网格法，是求偏微分方程定解问题的数值解中应用最广泛的方法之一。它的基本思想是：先对求解区域作网格剖分，将自变量的连续变化区域用有限离散点（网格点）集代替；将问题中出现的连续变量的函数用定义在网格点上离散变量的函数代替；通过用网格点上函数的差商代替导数，将含连续变量的偏微分方程定解问题化成只含有限个未知数的代数方程组（称为差分格式）。如果差分格式有解，且当网格无限变小时其解收敛于原微分方程定解问题的解，则差分格式的解就作为原问题的近似解（数值解）。因此，用差分方法求偏微分方程定解问题一般需要解决以下问题：（1）选取网格；（2）对微分方程及定解条件选择差分近似，列出差分格式；（3）求解差分格式；（4）讨论差分格式解对于微分方程解的收敛性及误差估计。下面，用一个简单的例子来说明用差分方法求解偏微分方程问题的一般过程及差分方法的基本概念。设有一阶双曲型方程初值问题。)()0,(,00xxuxtxuatu（1）选取网格：－2h-h0h2h3h首先对定解区域}0,),{(txtxD作网格剖分，最简单常用一种网格是用两族分别平行于x轴与t轴的等距直线khxxk，(0,1,2,0,1,2,)jttjkj将D分成许多小矩形区域。这些直线称为网格线，其交点称为网格点，也称为节点，h和分别称作x方向和t方向的步长。这种网格称为矩形网格。（2）对微分方程及定解条件选择差分近似，列出差分格式：如果用向前差商表示一阶偏导数，即2211(,)12(,)(,)(,)(,)2(,)(,)(,)2kjkjkjkjkjxxtkjkjkjtxtuxtuxtuhuxhtxhuxtuxtuuxtt其中1,021。方程0uuatx在节点),(jktx处可表示为htxutxuatxutxujkjkjkjk),(),(),(),(11),(2),(21222jkxjktthxuahtxu),2,1,0,,2,1,0(),(jktxRjk其中(,0)()(0,1,2,)kkuxxk。由于当,h足够小时，在式中略去),(jktxR，就得到一个与方程相近似的差分方程1,1,,,0kjkkkjjjuuuuah此处，jku,可看作是问题的解在节点),(jktx处的近似值。同初值条件),2,1,0()(0,kxukk结合，就得到求问题的数值解的差分格式。式)(),(2),(2),(1222hOthxuahtxutxRjkxjktjk称为差分方程的截断误差。如果一个差分方程的截断误差为)(pqhOR，则称差分方程对t是q阶精度，对x是p阶精度的。显然，截断误差的阶数越大，差分方程对微分方程的逼近越好。若网格步长趋于0时，差分方程的截断误差也趋于0，则称差分方程与相应的微分方程是相容的。这是用差分方法求解偏微分方程问题的必要条件。如果当网...
