CH02抛物方程差分法CH2.5-2.8 微分方程数值解ppt培训资料.ppt

我们先看一个数值例子，考虑初边值问题;利用显式差分格式(2.29)，即;的精确解。黑点是用差分格式在 时算出的相应各层上的近似值。二者符合得很好，由于对称性我们只给出一半图形。 图2.10是当 时差分方程解和微分方程精确解的图示，黑点仍表示差分方程解，其中 分别为在 时的计算结果。从图中看出，随着 的增大，差分方程的解越来越远离微分方程的解。 由此可见， 值的不同，得出的结果有很大的差别，如 的结果是可用的，但是 时的结果就完全没有用。 当然上面各种情况所得的差分方程解是由计算机得到的，不可能是差分方程理论上的准确解 ; ，而是差分方程的近似解，我们用 表示。显然 与 之间存在着差别，差分方程的准确解 与微分方程的解 之间，如前所述,也是有差别的。因而从计算机上解得的差分方程近似解 与微分方程解 之间的差别实质上包括两方面的差别，即;第一层的差分方程解 时，用到第0层上的 值，也就是初始值。由于计算机存储数据为二进制数位的限制, 不可能完全精确地存储在机器中，也就是计算 用到的是带有误差的初始值 。一般来说，在计算 时又出现了误差，因此 中包括了由于 参加运算而出现的误差，即初始误差的传递，以及本身计算过程中出现的误差。这样，在第 时间层计算 时得到的 是由于前面的误差传递和本身计算中出现的误差引起的。 下面我们给出研究差分格式稳定性的最直接的方法，就是在第0层的一个结点上给出一个误差 ，然后研究这个误差的发展情况，即 图方法。 ;表2.3; 由表2.3可知，用显式差分格式(2.29)( )计算时，由初始数据的误差，在以后各层所引起的误差是逐层减小的，这说明差分格式(2.29)当 时是稳定的。; 可见，用显式差分格式(2.29)( )计算时，由初始数据的误差所引起的误差在以后各层的计算中逐层迅速增大，以致不能控制，因此差分格式(2.29)在 时是不稳定的。 用 图方法讨论格式的稳定性能直观地看到差分格式是稳定性，缺点是必先固定 。 ;2.5.3 稳定性定义、稳定性分析的矩阵方法 ;其中; 设 是初始值引进的误差向量???而在边值以及其他各层计算中未引入其它任何误差。由于 的引入，差分方程的解为 。; 设向量 ，则常用的向量范数有：; (1) ，其中 ，为 的共轭 转置矩阵， 为 的最大特征值； ;对所有的范数都有 ;因此，可推得; 设 为 的特征值，用 表示 的最大值，即 的谱半径，则有 定理 2.1 差分格式(2.60)稳定的必要条件是，存在与 无关的常数 ，使矩阵 的谱半径满足;定理 2.3 若在差分格式(2.60)中， 为正规矩阵，即其满足 ，则条件(2.62)是差分格式(2.60)按欧几里德范数稳定的充分条件。 ;有关矩阵特征值计算的几个结论：;; 利用古典显式差分格式解定解问题(2.26)，(2.3)，(2.4)的相应差分方程组为;写成矩阵形式为;显然，矩阵的特征值为;(2.66);写成矩阵形式;令;方程(2.66)写为; 定解问题(2.26),(2.3),(2.4)的加权六点格式相应的差分问题是;写成矩阵形式;可以写为;其特征值为 ;即;所以当 时，加权六点格式稳定的条件为 ; 由上面的稳定性分析可知：;2.5.4 Gerschgorn定理及其在分析差分格式稳定性中的应用;定理 2.6 (Gerschgorn圆盘定理或Brauer定理) 设 为任意 复矩阵，则 的特征值都在复平面上的 个圆;考虑Crank-Nicolson隐式差分格式(2.66)。如前，这时稳定性分析归结为计算矩阵 的特征值，由 ;如果 的每一特征值的模不超过1，则格式将是稳定的。这相当于要求 ;圆盘定理，有 ;其中 是常数， 。;差分格式写成矩阵形式为;2.5.5 稳定性分析的 级数法( 方法) ; 假定微分方程是常系数的，因此差分方程也是常系数的，考虑两层常系数差分方程;在结点 上差分方程的近似解 与其理论解 之差 满足下面的差分方程 ;将 展开成 级数的复数形式，有;(2.76);因此差分方程(2.71)按 范数稳定的充分必要条件是;