数值分析-第六章-解线性方程组的迭代法.pptVIP

*/68 郑州大学研究生2011-2012学年课程 数值分析 Numerical Analysis §6.5 迭代法收敛判定的应用举例 证: 因为B=I-A, 故?(B)= ?(I)- ?(A)=1 - ?(A) ?(A) + ?(B) = 1 由于已知?(A) 和 ?(B)全为正数，故 0?(B)1 ,从而? (B) 1 所以该迭代格式收敛。 */68 郑州大学研究生2011-2012学年课程 数值分析 Numerical Analysis §6.5 迭代法收敛判定的应用举例 例6.5.3 设 方程组 写出解方程组的Jacobi迭代公式和迭代矩阵 并讨论迭代收敛的条件。 写出解方程组的Gauss-Seidel迭代矩阵,并讨 论迭代收敛的条件。 解 ① Jacobi迭代公式和Jacobi矩阵分别为 */68 郑州大学研究生2011-2012学年课程 数值分析 Numerical Analysis §6.5 迭代法收敛判定的应用举例 当时 时,Jacobi矩阵??GJ??∞1,对初值x(0)均收敛 */68 郑州大学研究生2011-2012学年课程 数值分析 Numerical Analysis 写出解方程组的Gauss-Seidel迭代矩阵，并讨论 迭代收敛的条件。 解 ② Gauss-Seidel矩阵为 当时 时,Gauss-Seidel矩阵 ??Gs??∞1, 所以对任意初值x(0)均收敛。 §6.5 迭代法收敛判定的应用举例 */68 郑州大学研究生2011-2012学年课程 数值分析 Numerical Analysis 解： 先计算迭代矩阵 例6.5.4 讨论用雅可比迭代法和高斯-塞德尔迭代 法解线性方程组Ax=b的收敛性。 */68 郑州大学研究生2011-2012学年课程 数值分析 Numerical Analysis 求特征值 雅可比矩阵 ? ( B ) = 1 ∴用雅可比迭代法求解时，迭代过程不收敛 ?1 = - 1， ?2,3 = 1/2 */68 郑州大学研究生2011-2012学年课程 数值分析 Numerical Analysis 求特征值 高斯-塞德尔迭代矩阵 ? (G1) = 0.3536 1 ∴用高斯-塞德尔迭代法求解时，迭代过程收敛 ?1=0, */68 郑州大学研究生2011-2012学年课程 数值分析 Numerical Analysis §6.5 迭代法收敛判定的应用举例 求解AX=b,当?取何值时迭代收敛？ 解:所给迭代公式的迭代矩阵为 例6.5.5 给定线性方程组 AX= b 用迭代公式 X(K+1)=X(K)+?(b-AX(K)) (k=0,1,…） */68 郑州大学研究生2011-2012学年课程 数值分析 Numerical Analysis 即 ?2-(2-5 ?)?+1- 5 ?+4 ?2=0 ?2-(2-5 ?)?+(1- ? )(1-4?)=0 [?-(1-?)][?- (1-4?)]=0 ?1=1-? ?2=1-4? ?(B)=max{|1- ?|, |1-4?|}1 取0 ?1/2迭代收敛 */68 郑州大学研究生2011-2012学年课程 数值分析 Numerical Analysis 内容小结 本章介绍了解线性方程组 迭代法的一些基本理论和具体方法。迭代法是一种逐次逼近的方法，即对任意给定的初始近似解向量，按照某种方法逐步生成近似解序列，使解序列的极限为方程组的解。注意到在使用迭代法 解方程组时，其迭代矩阵B和迭代向量f在计算过程中始终不变,迭代法具有循环的计算公式,方法简单，程序实现方便，它的优点是能充分利用系数的稀疏性,适宜解大型稀疏系数矩阵的方程组。 */68 郑州大学研究生2011-2012学年课程 数值分析 Numerical Analysis 迭代法不存在误差累积问题。使用迭代法的关键问题是其收敛性与收敛速度，收敛性与迭代初值的选取无关，这是比一般非线性方程求根的优越之处。在实际计算中，判断一种迭代格式收敛性较麻烦，由于求迭代的谱半径时需要求特征值，当矩阵的阶数较大时，特征值不易求出，通常采用矩阵的任一种范数都小于1或