基础知识
偏微分方程的定解问题
各种物理性质的定常（即不随时间变化）过程，都可用椭圆型方程来描述。其最典型、最简单的形式是泊松(Poisson)方程
?2u?2u?u?2?2?f(x,y) （1）
?x?y特别地，当 f ( x, y) ≡ 0 时，即为拉普拉斯(Laplace)方程，又称为调和方程
?2u?2u?u?2?2?0 （2）
?x?y带有稳定热源或内部无热源的稳定温度场的温度分布，不可压缩流体的稳定无旋流动及静电场的电势等均满足这类方程。
Poisson 方程的第一边值问题为
??2u?2u??u??x2??y2?f(x,y)(x,y)?? （3） ??u(x,y)??(x,y)????(x,y)???其 中 Ω 为 以 Γ 为 边 界 的 有 界区 域 ， Γ 为 分 段 光 滑 曲 线， Ω U Γ 称 为 定 解区 域 ，f (x, y),?(x, y) 分别为 Ω,Γ 上的已知连续函数。
第二类和第三类边界条件可统一表示成
??u??0(a?0) （4） ???u???n?(x,y)??其中 n 为边界 Γ 的外法线方向。当α = 0 时为第二类边界条件，α ≠ 0 时为第三类边界条件。
在研究热传导过程，气体扩散现象及电磁场的传播等随时间变化的非定常物理问题时，常常会遇到抛物型方程。其最简单的形式为一维热传导方程
?u?2u?a2?0(a?0) （5） ?t?x方程（5）可以有两种不同类型的定解问题： 初值问题（也称为 Cauchy 问题）
??u?2u??a2?0t?0,???x??? （6） ?x??t?u(x,0)??(x)???x????初边值问题
??u?2u??t?a?x2?00?t?T,0?x?l?? （7） ?u(x,0)??(x)?u(0,t)?g(t),u(l,t)?g(t),0?x?l12???其中?
?(x),g1(x),g2(x)为已知函数，且满足连接条件 ?(0)?g1(0),?(l)?g2(0)问题（7）中的边界条件u(0,t)?g1(t),u(l,t)?g2(t)称为第一类界条件。第二类和第三类边界条件为
??u???(t)u?g1(t),0?t?T1???x??x?0??u???(t)u?g2(t),0?t?T2????x?x?l
（8）
其中?1?0,?2?0。当?1??2?0时，为第二类边界条件，否则称为第三类边界条件。 双曲型方程的最简单形式为一阶双曲型方程
?u?u?a?0?t?x （9）
物理中常见的一维振动与波动问题可用二阶波动方程
?2u?2u?a2 （10） ?t2?x描述，它是双曲型方程的典型形式。方程（10）的初值问题为
??2u?2u?2?a2?t?x???u(x,0)??(x)??u???(x)??t?t?00?t,???x??????x??????x??? （11）
边界条件一般也有三类，最简单的初边值问题为
??2u?2u0?t,0?x?l??t2?a?x2??u?u(x,0)??(x)??(x)0?x?l ??tt?0??u(0,t)?g1(t),u(l,t)?g2(t)0?t?T??如果偏微分方程定解问题的解存在，唯一且连续依赖于定解数据（即出现在方程和定解条件中的已知函数），则此定解问题是适定的。可以证明，上面所举各种定解问题都是适定的。 §2 偏微分方程的差分解法
差分方法又称为有限差分方法或网格法，是求偏微分方程定解问题的数值解中应用最广泛的方法之一。它的基本思想是：先对求解区域作网格剖分，将自变量的连续变化区域用有限离散点（网格点）集代替；将问题中出现的连续变量的函数用定义在网格点上离散变量的函数代替；通过用网格点上函数的差商代替导数，将含连续变量的偏微分方程定解问题化成只含有限个未知数的代数方程组（称为差分格式）。如果差分格式有解，且当网格无限变小时其解收敛于原微分方程定解问题的解，则差分格式的解就作为原问题的近似解（数值解）。因此，用差分方法求偏微分方程定解问题一般需要解决以下问题：
（i）选取网格；
（ii）对微分方程及定解条件选择差分近似，列出差分格式； （iii）求解差分格式；
（iv）讨论差分格式解对于微分方程解的收敛性及误差估计。 下面我们只对偏微分方程的差分解法作一简要的介绍。 2.1 椭圆型方程第一边值问题的差分解法
以 Poisson 方程（1）为基本模型讨论第一边值问题的差分方法。 考虑 Poisson 方程的第一边值问题（3）
??2u?2u?2?2?f(x,y)(x,y)???y ??x?u(x,y)??(x,y)????(x,y)???取 h,τ 分别为 x 方向和 y 方向的步长，以两族平行线 x?xk?kh,y?yj?j?
(k,j?0,?1,?2,?)将 定 解 区 域 剖 分 成 矩 形 网 格 。 节 点 的 全 体 记 为
R?{(xk,yk)|xk?kh,yj?j?, i,j为整数} 。定解区域内部的节点称为内点，记内点集
R??为?h? 。边界Γ与网格线的交点称为边界点，边界点全体记为 Γ
hτ
。与节点
(xk,yj)沿 x 方向或 y 方向只差一个步长的点(xk?1,yj)和 (xk,yj?1) 称为节点 (xk,yj)的相邻节点。如果一个内点的四个相邻节点均属于ΩU Γ，称为正则内点，正则内点的全体记为Ω(1)，至少有一个相邻节点不属于ΩU Γ 的内点称为非正则内点，非正则内点的全体记为Ω(2)。我们的问题是要求出问题（3）在全体内点上的数值解。
为简便记，记(k,j)?(xkyj),u(k,j)?u(xk,yj),fk,j?f(xk,yj)。对正则内点
(k,j)??(1)，由二阶中心差商公式
?2uu(k?1,j)?2u(k,j)?u(k?1,j)??O(h2) 22?x(k,j)h?2uu(k,j?1)?2u(k,j)?u(k,j?1)??O(?2) 22?y(k,j)?Poisson 方程（1）在点(k,j)处可表示为
uk?1,j?2uk,j?uk?1,jh2?fk,j?O(h2??2)?uk,j?1?2uk,j?uk,j?1?2 （12）
在式（12）中略去O(h2??2)，即得与方程（1）相近似的差分方程
uk?1,j?2uk,j?uk?1,jh2?uk,j?1?2uk,j?uk,j?1?2?fk,j （13）
式（13）中方程的个数等于正则内点的个数，而未知数 uk,j, 则除了包含正则内点处解u的近似值，还包含一些非正则内点处u的近似值，因而方程个数少于未知数个数。在非正则内点处 Poisson 方程的差分近似不能按式（13）给出，需要利用边界条件得到。
边界条件的处理可以有各种方案，下面介绍较简单的两种。 (i) 直接转移 (ii) 线性插值
由式（13）所给出的差分格式称为五点菱形格式，实际计算时经常取h =τ ，此时 五点菱形格式可化为
1?uk?1,j?uk?1,j?uk,j?1?uk,j?1?4uk,j??fk,j （14） 2h简记为
1?uk,j?fk,j （15） h2其中 ?uk,j?uk?1,j?uk?1,j?uk,j?1?uk,j?1?4uk,j。
求解差分方程组最常用的方法是同步迭代法，同步迭代法是最简单的迭代方式。除 边界节点外，区域内节点的初始值是任意取定的。
例 1 用五点菱形格式求解 Laplace 方程第一边值问题
??2u?2u?2?2?0?y??x22?u(x,y)?lg(1?x)?y(x,y)???(x,y)????
????1。 3其中??{(x,y)|0?x,y?1}。取 h???当h??时，利用点(k, j),(k ±1, j .1),(k ±1, j +1)构造的差分格式
1?uk?1,j?1?uk?1,j?1?uk?1,j?1?uk?1,j?1?4uk,j??fk,j （16） 2h2称为五点矩形格式，简记为
1?uk,j?fk,j （17） 2h2其中?uk,j?uk?1,j?1?uk?1,j?1?uk?1,j?1?uk?1,j?1?4uk,j。
2.2 抛物型方程的差分解法 以一维热传导方程（5）
?u?2u?a2?0(a?0) ?t?x为基本模型讨论适用于抛物型方程定解问题的几种差分格式。
首先对xt 平面进行网格剖分。分别取h,τ 为x方向与t方向的步长，用两族平行直 线x?xk?kh (k = 0,±1,±2,…) ，t?tj?j?k ( j = 0,1,2, …)，将 xt 平面剖分成矩形网格，节点为(xk,yj) (k = 0,±1,±2, …, j = 0,1,2, …)。为简便起见，记(k,j)?(xk,yj),
u(k,j)?u(xk,yj),?k??(xk),g1j?g1(tj),g2j?g2(tj),?1j??1(tj),
?2j??2(tj)。
2.2.1 微分方程的差分近似
?2u?u在网格内点(k, j)处，对 分别采用向前、向后及中心差商公式，对2采用二
?t?x阶中心差商公式，一维热传导方程（5）可分别表示为
u(k,j?1)?u(k,j)u(k?1,j)?2u(k,j)?u(k?1,j)2?O(??h)2?hu(k,j)?u(k,j?1)u(k?1,j)?2u(k,j)?u(k?1,j)?a?O(??h2) 2?hu(k,j?1)?u(k,j?1)u(k?1,j)?2u(k,j)?u(k?1,j)?a?O(??h2)22?h?a由此得到一维热传导方程的不同的差分近似
uk,j?1?uk,j?uk,j?uk,j?1?auk?1,j?2uk,j?uk?1,jh2?0 （18）
??auk?1,j?2uk,j?uk?1,jh2?0 （19）
uk,j?1?uk,j?12??auk?1,j?2uk,j?uk?1,jh2?0 （20）
2.2.2 初、边值条件的处理
为用差分方程求解定解问题（6），（7）等，还需对定解条件进行离散化。
偏微分方程—matlab(DOC)
