2归纳.4差分方程的相容性、收敛性和稳定性归纳.pptVIP

差分方程收敛性有两种证明方法，直接证明法和数值试验法。 一、直接证明法 对流方程 的FTBS差分格式为： (a) 设求解区域内任意一点 ， ，它的微分方程精确解为u，差分方程解为 ，则离散化误差为 ，把差分方程和微分方程相减可得离散化误差方程： (b) 由（b）式可以看出离散化误差方程在形式上和差分方程是完全相同的，由此可以得到： 设a≥0， ≤1，则0≤ ≤1，于是有： (c) 式中 表示在n层的所有节点上离散化误差 绝对值最大值，对于所有节点j有： 于是有： … 由此可得到： (d) 在t=0时，差分方程的初始条件应该是完全准确的，即： 即： 即差分方程离散化误差和截断误差是相同数量级，因此，若 →0，则： (f) 由此可知，FTBS格式在a＞0， 时，是收敛的。 (e) 二、数值试验法 数值试验法基本思想是用差分方程求出FTBS数值解，然后和微分方程精确解进行比较，确定差分方程是否收敛。 直接证明法比较简单，但是只有很少几个差分方程可以采用直接证明法来证明其收敛性，而数值试验法又非常麻烦，一般来说，很难用数值试验结果