边值问题的有限差分法详解.ppt

第九章 边值问题的有限差分法;一、有限差分法的基本概念;1、剖分区域－－建立差分网格;;式中， 于 的 近似， 称为差分方程 （1.4)的截断误差。显然，差分方程 （1.4)当 时成立，加上边值条件 就得到关于 的线性代数方程组 ;为了求解方便，把差分方程（1.5）改写成;;定义1 如果当网格无限加密时，差分格式的 解 存在，且按某一范数 有;为例。记方程;定义2、设差分方程为;; 稳定性概念不仅在理论研究而且在实际应用 中都有重要意义。实际上，由于有实测误差和舍 入误差。右端数据不可能准确给出。如果小的右端 误差会引起解的很大偏离，即差分方程不稳定，则 如此的差分方程便没有实际意义。 ;其中;我们只考虑属于 的节点。如果两个节点沿x ???（或沿y轴）方向只相差一个步长，则称这两个 节点是相邻的；以 表示所有属于G 内部的节点集合，并称如此的节点为内点。以 表示网线 或 与G的交点集合，并称如此 的点为边界点。如果内点的四个相邻点都属于 则称如此的点为正则内点；否则称为非正则内点。;若考虑的微分方程是Laplace方程，则相应的差分 方程为;;可以证明，差分格式（2.8）的截断误差为;其中，;对于;例、用差分法解边值问题;方程的右端项可从边界条件的直接转移法得到：;用Gauss-Seidel方法求解上述线性代数方程组，其结果如下：;上例的解析解为