数值分析（河海大学）ch6解线性方程组的迭代法.pptVIP

第6章 解线性方程组的迭代法6．1 引言 线性方程组的迭代解法的优点： 能有效利用矩阵的稀疏特点； 可用于提高已知近似解的精度。 迭代法解线性方程组时，需解决三个问题： (1)迭代格式的构造； (2)迭代格式的判敛； (3)迭代格式的收敛速度估计。 6．2 基本迭代法 与之等价的方程组为 它们是从原方程组中的第 个方程解出 。 方程的右端可用以构造迭代函数。 定义 取定初值 ，若用公式 逐步代入求得近似解的方法称为迭代 法（或一阶线性定常迭代法）。 称为迭代矩 阵。若 ，则称迭代法收敛， 为方 程组的解，否则称迭代法发散。 J-迭代的分量形式 G-S迭代的分量形式 称为逐次超松驰（SOR） 迭代格式， 称为松驰因子， ，称为超 松驰法， ，称为低松驰法。 ，是G-S迭代格式。 SOR迭代的矩阵形式 6．3 迭代法的收敛性1.一阶线性定常迭代法的基本定理 设一阶线性定常迭代格式为 基本定理 迭代法收敛的充要条件为 。 定理 设方程组 ，一阶线性定常迭 代格式 ，且对某一 范数有 ： ，则任意初始向量 有 迭代收敛，且方程组的解为 ； ； ； 。 2.关于解某些特殊方程组迭代法的收敛性 定义（对角占优）若 的元素满足 则称 为行严格对角占优阵，若 的元素满足 则称 为列严格对角占优阵。若以上严格不等式 改为≥,则称 为行(列)弱对角占优阵。 定义（可约）若存在置换阵 使 则称 为可约矩阵，否则称 为不可约矩阵。 定理 设 为严格对角占优或弱对角占优且不可约矩阵，则 有唯一解，且J-迭代和G-S迭代均收敛。 证明（仅证 为行严格对角占优阵） 定理（SOR迭代收敛的必要条件）SOR迭代 收敛 。 定理 若 （1） 是对称正定矩阵， （2） ， 则 SOR迭代收敛。 定理 若 （1） 是严格对角占优或弱对角占优且不可约矩阵， （2） ， 则SOR迭代收敛。 6．4 分块迭代法 方程组 ， ， 注意 ， G-S迭代收敛。 证毕 证明 因为迭代矩阵 令 ， 若收敛 ， ， ，即： 。 证毕 证明 设 为 的特征值和特征向量， （ 不一定是实的，可以是复的。） ， 左乘 记 取共轭再转置，得 ， ， 所以，SOR迭代收敛。 证毕 （证略） 例 试对方程组 建立收敛的G-S迭代格式，其中 ， 。 解 即 变换后的方程组其系数阵是弱对角占优且不可约矩阵， 所以，G-S迭代格式收敛。 河海大学理学院《数值分析》 数 值 分 析 Computational Method Chapter 6 解线性方程组的迭代法 设线性方程组： 分量形式 从第 个方程解出 的等价方程的分量形式 例 解线性方程组 （精确到第五位小数） 解 等价方程组为 要求精确到第五位小数即允许的迭代误差限为 （雅可比迭代法） 取 ，计算结果见下表： （高斯-塞德尔迭代法） 取 ，计算结果见下表： 解毕 一般, ，改写成： 设 ， ， ， 等价方程的矩阵形式 : J-迭代的矩阵形式 迭代矩阵 G-S迭代的矩阵形式 迭代矩阵 若记 作为修正量， SOR迭代的分量形式 SOR迭代矩阵 上例的SOR迭代格式 证明 迭代法收敛即 ， 令 ，则 ， 为固定向量。所以， 迭代法收敛 ，