数值分析习题.docVIP

定义2 若某方法在节点值上有大小为的摄动，而其后各节点上的误差均不超过，则称此方法是稳定的。 例5 对方程考察Euler公式和后退的Euler公式的稳定性。 解 对应的Euler公式为。 若在上有摄动值，而它使产生的摄动值为，则 。 显然，即时，Euler公式稳定，称之为条件稳定。 后退的Euler公式为，，而，即后退的Euler公式无条件稳定。 定理1 若(1)对任意，有(压缩映射)；(2)存在正数，使对任意，有(Lipschitz条件)，则对任意，迭代序列收敛于的唯一根，且有误差估计式 。 证 令，则由，得，且。 若，则即为的根，否则由零点定理，在内至少有一根，故在上至少有一根。（存在性） 若有两个根，即，从而有 ， 矛盾，故。（唯一性） 。 由于，，故 即。（收敛性） 。 ，得。 令，则有。 定义2 设迭代公式收敛于的根，如迭代误差 满足，则称为p阶收敛。 ，称之为牛顿公式。 牛顿法对附近的初值收敛，即局部收敛。再由定理3知，牛顿法为二阶收敛的。 例5 应用牛顿法于方程和，分别导出求的迭代公式，讨论迭代公式的收敛性，并求。 解 牛顿公式为。 因，即有下界。 又，即单调下降，故收敛。 。 三对角方程追赶法 计算机解法：(1) A为一般方阵时，用LU分解(Gauss消去法)；(2) A对称正定时，用Cholesky分解。 定理6 A的谱半径不超过A的范数，即。 证 设为A的任一特征值，为相应的特征向量，即，则 (相容性)，故，即。 定理7 若，则可逆，且。 证 若不可逆，即，有非零解，亦即有，使，，，矛盾，故可逆。 ， ，得。 §5. 误差分析和矩阵的条件数 1. 右端向量误差对解的影响 设的右端向量有误差，相应的解为，则 ，。 又，故。 2. 系数矩阵误差对解的影响 设的系数矩阵有误差，相应的解为，则 。 由定理7，若，则可逆，， ， 且，故 迭代法收敛的充要条件为。 (1) 若A为严格对角占优矩阵，即，则Jacobi迭代和G-S迭代收敛； (2) 若A对称正定，则G-S迭代收敛。 计算表明，当阶数较高时，G-S迭代收敛速度仍然较慢，可在G-S迭代基础上对其加速收敛—超松弛法。 从G-S迭代公式求得的结果不作为第次近似解，而仅作为中间结果，然后再将与进行加权平均后作为第次近似解，即。 时，即为G-S迭代，时称为超（低）松弛法。 《数值分析》 第五章 常微分方程数值解法 1