偏微分方程数值解解析.ppt

于是有 矛盾． （负的极小），则一定在边界上取． 推论4.2 若网函数 注　本定理说明在定理中所给条件下，若有正的极大 推论4.1 差分方程有唯一解． 满足 则 定理4.2 设 是两个网函数，满足 则 证明　由条件可知 由 Cor4.2 知，结论成立． 推论4.3 的解 证明　设 由定理4.2知： 满足 是如下问题的解 若 则 若 由定理4.1知， 的最大值只能在 达到，因此 从而结论成立． 4.2 五点格式的敛速估计 设 是Poisson方程第一边值问题的解， 差分格式的解，即 令 则 是五点 截断误差 令 无论 则设 设 R 是以 (0,0) 为心包含在内部的最小圆域 的半径 . 令 是正则内点还是非正则内点， 于是 由定理4.2知： 作业：P 105 1 ，2 从而 若 则 一致收敛到 u，且有如上敛速估计． 第三章 椭圆形方程的有限差分法 有限差分法和有限元方法是解偏微分方程的两种 用数值微商或数值积分导出相应的线性代数方程组. 主要数值方法. 两种方法的主要差别：离散化的第二步. 有限差分法：从定解问题的微分或积分形式出发， 但基函数要按特定方式选取. 有限元方法：从定解问题的変分形式出发，用 Ritz-Galerkin 方法导出相应的线性代数方程组， 差分法求解的主要步骤： (3) 差分方程的解法. (1) 对求解区域做网格剖分. 一维：将区间分成等距或不等距的小区间单元. 二维：将区域分割成一些均匀或不均匀的矩形，其边与 坐标轴平行，或分割成一些三角形或一些凸四边形等. (2) 构造逼近微分方程定解问题的差分格式：三种方法 直接差分化法、积分插值法、有限体积法(或广义差分法). 差分解的存在唯一性、差分性及稳定性的研究. §1 差分逼近的基本概念 考虑二阶微分方程边值问题 ① 将其分成等分，分点为 ?将方程 (1.1) 在节点 其中 q，f 为 [ a , b ] 上的连续函数， 为给定常数. 处再离散化. 由 Taylor 展开得 其中 于是在 其中 舍去 将方程 (1.1) 写成 点取值. 表示方括号内的函数在 得逼近方程 (1.1) 的差分方程为： 称 形成关于 注 此方程组尽管是高阶方程组，但每个方程未知数 ④ 对方程组 (1.4)～(1.5) 的解分析需要考虑以下几个问题： 最多有3个易于求解. 为差分方程 (1.3) 的截断误差. 的线性代数方程组 (a) 解是否惟一？ 以 定义在 对 (b) 当网格无限加密时，即 时，差分解 是否收敛到真解 (c) 在何种度量下收敛？ (d) 收敛速度如何？ 为了解决如上问题，需要给出如下说明： 内点和界点 表示网格内点 的集合， 表示网格 集合. 上的网函数. 或 上的函数 称为 上的网函数引进如下范数： 其中 定义1.1 设 U 是某一充分光滑的函数类， 称差分算子 断误差定义的网格函数. 为相容条件. 是由截 恒有 若对任何 逼近微分算子 L ，并称 (1.6) 注 当用 定义1.2 当 h 充分小时，若 (1.4) ～ (1.5) 的解 称 如果存在与网格 阶也就不同. 且按某一范数 定义1.3 对于差分方程 数 M 和 逼近 L 时，选择网函数的范数不同，逼近的 存在， 有 收敛到边值问题的解 u . 无关的常数 及右端 使 称差分方程关于右端稳定. 相同的收敛阶. 定理1.1 若边值问题的解 u 充分光滑，差分方程按 满足相容条件，且关于右端稳定，则差分解 收敛到边值问题的解 u ，且有与 按 §2 两点边值问题的差分格式 考虑两点边值问题 其中 是给定的常数. 2.1 直接差分法 (1) 取 N+1 个节点将 I =[a, b] 分成 N 个小区间: 于是，得到 I 的一个网格剖分. (2) 对 I = [a, b] 进行对偶剖分 取 称为半整数点，则 的中点 构成 I 的一个对偶剖分. (3) 将方程 (2.1) 在内点 (2.3) (2.4) 处离散化. (2.5) 将 (2.5) 减 (2.4) 再除以 得 于是得逼近方程 (2.1)~(2.2) 的差分方程： 误差为 2.2 积分插值法 在 [a, b] 内任一小区间 考虑守恒型微分方程 上积分得： 其中 令 取 则 或 于是 又 沿 积分，得 得守恒型差分方程： 其中 (2.6) 2.3 边值条件的处理 考虑第二、第三边值条件 (2.7) (2.8) 由于 不妨设 (2.9) (2.10) 取 得 而 故 于是逼近边值条件 (2.9)~(2.10) 的差分方程为： (2.11) 可以证明：