有限差分方法(2010-07-19).docVIP

椭圆型微分方程的有限差分法 主讲： 谭 林 基本思想（步骤）： 将求解区域（无限个点）限制在有限个离散点上，一般可通过网格剖分获得。 在离散点处，将求微分问题（无限计算问题）近似化为求若干（相邻）离散点上函数值的线性组合问题（有限计算问题），一般利用数值微商（分）（不同有限元法）。形成所谓的差分方程。 差分方程的适定性、收敛性和稳定性分析。 差分方程的解法。 下面以两点边值问题为例介绍有限差分法全过程 一、常见的有限差分方法 直接差分法 模型问题1：椭圆型方程第一边值问题。 其中， 模型问题2： 其中， eq \o\ac(○,1)首先对模型问题1 讨论其有限差分方法的基本步骤 求解区域的离散化 做均匀网格剖分： 其中，分点 剖分步长 在节点处，对微分方程离散化 有 其中 记在节点数值解为 ， 则有 (*1) 比较知 所以 表示用差分算子代替微分算子产生的误差 称之为（局部）截断误差。这里关于的阶为 。 注意 所以 由此知：（局部）截断误差可视为差分格式（方程）(*1)，将数值解换成相应真解值后，左端减右端，再做Taylor展式获得的（可作为计算公式）。 方程 (*1) 的联立形式（中心差分格式） 矩阵形式 (其中 是三对角矩阵) eq \o\ac(○,2)接下来对模型问题2 进行讨论 在节点处，对微分方程离散化 关键，建立 的离散（近似）公式。 做一般网格剖分(非均匀)： 其中，分点 第个剖分单元的剖分步长 引入对偶剖分 其中，分点 建立离散公式的基本思想是一阶中心差商去近似一阶导数 令 Flux Function. 第一次用一阶中心差商近似一阶导数 第二次用一阶中心差商近似一阶导数 最后的格式为 (局部)截断误差的推导和模型问题1类似. 习题.在一般网格剖分下，对微分方程 建立中心差分格式. 积分差分法（有限体积法） 考虑微分方程节点处的离散化 模型问题3： 其中， Step 1. 将微分方程 在节点所属的对偶单元 上积分得：积分方程 Step 2. 对含未知函数的三项 注意：一般要求函数连续，但 允许间断。所以不对直接离散化。 以为例讨论。注意 得 因此 而 类似 而 综上得微分方程在节点处的离散化格式（差分方 程）为 特别利用数值积分，可得： 直接差分方法的结果! 积分差分法（有限体积法）的优点 对系数的光滑性要求低 变于推广到任意网格、边界条件 保守恒性（对收敛性重要） 变分--差分法 类似有限元法，不介绍 二、差分方程的适定性、收敛性和稳定性分析 考虑模型问题3： 其中， 差分方程（格式）为 （*） 等价形式 其中， 矩阵形式 M矩阵，它们满足 （1） （2）对角占优性 极值定理 若（），对， 则不可能在内点取正的极大（负的极小）， 除非常数。 证明. （），对。 反证法，设在某内点取正的极大M, 且常数，则 存在某内点，使 或 这时 矛盾。 推论 1. 差分方程（*）是适定的（存在唯一性）。 证明. 注意差分方程（*）的矩阵形式 差分方程（*）是存在唯一性，只需证明矩阵 非奇异。或奇次方程 （**） 只有零解。 注意，差分方程（*）的奇次情形 （***） 其矩阵形式为（**），所以只需证明差分方程 （***）的解是零解。而由极值定理知只能在边界点上取正的极大和负的极小。所以（***）的解为零解。 推论 2. 若，对，且 ，则 ，对。 即M矩阵的逆矩阵是非负矩阵。 （因为 其中，。 若对任意非负向量, 均为非负向量，则 M矩阵的逆矩阵只能为非负矩阵。这点与其相应的连续算子，即Laplace算子的性质非常类似，亦为正算子。） 比较定理.设序列（网函数），满足 ，且 则 推论 3（关于边界值的稳定性）.差分方程 的解满足估计式 差分方程（*）(或（**）)的收敛性估计（同时，回答了差分方程的解关于右值的稳定性）。 对特殊模型问题 均匀剖分 令 满足差分方程 其中，当时 选取优势函数 直接验算 由比较定理 定理. 设,则 习题： 推导上述未证明的推论或定理。 注.利用嵌入法（对有限元法类似，需用LAX-Milgram Th.）可证明（饱和阶） 设 令 则有 利用上述定理的同样推导可证 从而证之。