-有限差分法基础.ppt

《材料成形模拟技术》讲义 第2章 有限差分法基础 华中科技大学 周建新 Tel: 027Email: zhoujianxin1975@163.com 第三节 收敛性与稳定性/收敛性（3/6） 相容性不一定能保证收敛性，那么对于一定的差分格式，其解能否收敛到相应微分问题的解？答案是差分格式的解 收敛于微分问题的解是可能的。至于某给定格式是否收敛，则要按具体问题予以证明。下面以一个差分格式为例， 讨论其收敛性： 微分问题 的FTBS格式为 在某结点（xi , tn）微分问题的解为 ，差分格式的解为 ，则离散化误差为 第三节 收敛性与稳定性/收敛性（4/6） 按照截断误差的分析知道 以FTBS格式中的第一个方程减去上式得 或写成 若条件 和 成立，即 ，则 式中 表示在第n层所有结点上 的最大值。 第三节 收敛性与稳定性/收敛性（5/6） 由上式知，对一切i有 故有 于是 综合得 第三节 收敛性与稳定性/收敛性（6/6） 由于初始条件给定函数 的初值，初始离散化误差 。并且 是一有限量，因而 可见本问题FTBS格式的离散化误差与截断误差具有相同的量级。最后得到 这样就证明了，当 时，本问题的RTBS格式收敛。这种离散化误差的最大绝对值趋于零的 收敛情况称为一致收敛。 。 此例介绍了一种证明差分格式收敛的方法，同时表明了相容性与收敛性的关系：相容性是收敛性的必要条件， 但不一定是充分条件，还可能要求其他条件，如本例就是要求 。 第三节 收敛性与稳定性/稳定性（1/2） 首先介绍一下差分格式的依赖区间、决定区域和影响区域。还是以初值问题 为例。先看FTCS格式，如图（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点为例的，事实上对任意层的任一结点，都在初始层上有一对应的依赖区间，而初始层的任一区间都有一对应的决定区域。 (a) FTCS (b) FTFS (c) FTBS 差分格式的依赖区间 第三节 收敛性与稳定性/稳定性（2/2） 随着时间的推移，一点函数值将影响以后某些结点的解。如图，设p为第n层的某结点，当用FTCS格式 计算第n+1层上的结点值时，a、b、c这3点的解必须用到p点的函数值，在第n+2层上则有更多点的解受p点 函数值的影响。所有受p点函数值影响的结点总和为p点的影响区域，如图中阴影所示区域。 FTCS格式 (b) FTFS格式 (c) FTBS格式 差分格式的影响区域 第四节 Lax等价定理（1/2） 前面讨论了差分问题的相容性、收敛性和稳定性。已经知道，相容性是收敛性的必要条件；还发现，稳定性与收敛性有一定的联系。Lax等价定理就是阐述相容性、收敛性和稳定性三者之间关系的。 Lax等价定理：对一个适定的线性微分问题及一个与其相容的差分格式，如果该格式稳定则必收敛，不稳定必不收敛。换言之，若线性微分问题适定，差分格式相容，则稳定性是收敛性的必要和充分条件。这也可表示为 第四节 Lax等价定理（2/2） 根据此定理，在线性适定和格式相容的条件下，只要证明了格式是稳定的，则一定收敛；若不稳定，则不收敛。由于收敛性的证明往往比稳定性更难，故人们就可以把注意力集中在稳定性的研究上。 * * 主要内容 1、差分原理及逼近误差 2、差分方程，截断误差和相容性 3、收敛性与稳定性 4、Lax等价定理 第一节 差分原理及逼近误差/差分原理（1/8） 1．差分原理 设有x的解析函数y=f(x)，从微分学知道函数y对x的导数为 是函数对自变量的导数，