有限差分方法基础-研究生课程.pptVIP

有限差分方法应用-研究生课程讲义 有限差分法基础 材料学院 周建新 Tel: 027Email: zhoujianxin1975@163.com 第三节 收敛性与稳定性/收敛性（3/6） 相容性不一定能保证收敛性，那么对于一定的差分格式，其解能否收敛到相应微分问题的解？答案是差分格式的解 收敛于微分问题的解是可能的。至于某给定格式是否收敛，则要按具体问题予以证明。下面以一个差分格式为例， 讨论其收敛性： 微分问题 的FTBS格式为 在某结点（xi , tn）微分问题的解为 ，差分格式的解为 ，则离散化误差为 （3-6） （3-5） （3-4） 第三节 收敛性与稳定性/收敛性（4/6） 按照截断误差的分析知道 以FTBS格式中的第一个方程减去上式得 或写成 若条件 和 成立，即 ，则 式中 表示在第n层所有结点上 的最大值。 （3-7） （3-8） （3-9） （3-10） 第三节 收敛性与稳定性/收敛性（5/6） 由上式知，对一切i有 故有 于是 综合得 （3-11） （3-13） （3-12） （3-14） 第三节 收敛性与稳定性/收敛性（6/6） 由于初始条件给定函数 的初值，初始离散化误差 。并且 是一有限量，因而 可见本问题FTBS格式的离散化误差与截断误差具有相同的量级。最后得到 这样就证明了，当 时，本问题的RTBS格式收敛。这种离散化误差的最大绝对值趋于零的 收敛情况称为一致收敛。 。 （3-15） （3-16） 此例介绍了一种证明差分格式收敛的方法，同时表明了相容性与收敛性的关系：相容性是收敛性的必要条件， 但不一定是充分条件，还可能要求其他条件，如本例就是要求 。 第三节 收敛性与稳定性/稳定性（1/8） 首先介绍一下差分格式的依赖区间、决定区域和影响区域。还是以初值问题 为例。先看FTCS格式，如图3-1（a）欲计算第二层p点的函数值，必先知道第一层上a、b、c这3点的函数值，故说p点的解依赖于a、b、c这3点的解。而a点的解又依赖于第0层（初值线）上A、d、e的初值，b点的解依赖于d、e、f的初值，c点的解依赖于e、f、B的初值。因此p点的解依赖于初值线AB段上所有结点的初值，故称AB段上所有结点为p点的依赖区间。又，三角形pAB区域内任一结点的依赖区间都包含在AB之内，即该区域内任一结点上的解都由AB段上某些结点的初值所决定，而与AB以外结点的初值无关 ，故称此三角形区域为AB区间所决定的区域。这里为方便起见，是以第二层的p点为例的，事实上对任意层的任一结点，都在初始层上有一对应的依赖区间，而初始层的任一区间都有一对应的决定区域。 （3-17） (a) FTCS (b) FTFS (c) FTBS 图3-1 差分格式的依赖区间 第三节 收敛性与稳定性/稳定性（2/8） 随着时间的推移，一点函数值将影响以后某些结点的解。如图3-2，设p为第n层的某结点，当用FTCS格式 计算第n+1层上的结点值时，a、b、c这3点的解必须用到p点的函数值，在第n+2层上则有更多点的解受p点 函数值的影响。所有受p点函数值影响的结点总和为p点的影响区域，如图3-2中阴影所示区域。 FTCS格式 (b) FTFS格式 (c) FTBS格式 图3-2 差分格式的影响区域 第三节 收敛性与稳定性/稳定性（3/8） 例如微分问题 其解为零，即 。若用FTBS格式计算，且计算中不产生任何误差，则结果也是零，即 由上可知，同一微分问题，当采用不同差分格式时，其依赖区间、决定区域和影响区域可以是不一样的。 依赖区间、决定区域和影响区域是由差分格式本身的构造所决定的，并与步长比 有关。 （3-18） （3-19） 假设在第k层上的第j点，由于计算误差得到 。不妨设k=0, j=0, ，即相当于FTBS格式写成 （3-20） 第三节 收敛性与稳定性/稳定性（4/8） 现分别取 和2，列表计算如下： （1） 4 3 2 1 0 -1 -2 -3 -4 i 0 0 0 0 1 0 0 0 0 0 0 0 0 1 2 1 2 0 0