第四章迭代求解线性方程组概述.pptVIP

单步线性定常迭代算法 雅可比(Jacobi)迭代法 雅可比迭代法的分量形式 雅可比迭代程序设计 高斯-赛德尔(Gauss-Seidel)迭代法 G-S迭代法的分量形式 G-S迭代程序设计 超松驰(Successive Overrelaxation)迭代法 SOR迭代法的分量形式 G-S迭代程序设计 * 河南师大数学与信息科学学院 线性方程组迭代解法 本章讨论解线性方程组 古典迭代解法。 主要内容包括： 雅可比(Jacobi)迭代法 高斯-赛德尔(Gauss-Seidel)迭代法 超松驰(SOR)迭代法 收敛性分析 编制计算程序时,应注意以下几个问题: 算法的终止 程序设计 点击对象展开程序 算法的实现 编制计算程序时,应注意以下几个问题: 算法的终止 程序设计 点击对象展开程序 算法的实现 编制计算程序时,应注意以下几个问题: 算法的终止 程序设计 点击对象展开程序 算法的实现 注意: 迭代阵M不唯一,影响收敛性。 单步线性定常迭代法的收敛性 定义：矩阵M的谱半径为： 定义：如果存在矩阵G使：G(I – M)=A，Gg=b.则称迭代法与方程组Ax=b是相容的。 引理1 迭代法 收敛当且仅当 引理4.3.2 设 ，则有 1. 对 上的任意矩阵范数 ，有 2. 对任给的 ,存在 上的算子范数 使得 注: 引理的详细证明过程,请阅读第四章,第3节. 定理1 当且仅当 定理2 单步线性定常迭代法收敛当且仅当 例 1 当方程组的系数矩阵为 可以验证：Jacobi迭代收敛，而Gauss-Seidel迭代不收敛。 例 2 当方程组的系数矩阵为 可以验证：Jacobi迭代不收敛，而Gauss-Seidel迭代收敛。 例3 对于方程组 Jacobi迭代公式为 ，迭代矩阵为 Gauss-Seidel迭代公式为 ，迭代矩阵为 显然，两种迭代都不收敛。