§8.3迭代法的收敛性(缩减).ppt

* ? 2009, Henan Polytechnic University * §3 迭代法的收敛性 第八章 解线性方程组的迭代法 * * ? 2009, Henan Polytechnic University * §3 迭代法的收敛性 第八章 解线性方程组的迭代法 计算方法 第三节 迭代法的收敛性 　迭代矩阵B满足什么条件时，由迭代法产生的向量序列{x(k)}收敛到x*. 引进误差向量 可得误差向量的递推公式 则研究迭代法收敛性问题就是要研究迭代矩阵B满足什么条件时，有. 8.3.1 迭代法的收敛性 定义8.2 设有矩阵序列 Ak=(aij(k))∈Rn×n 及A=(aij)∈Rn×n，如果n2个数列极限存在且有 则{Ak}称收敛于A，记为lim (k→∞). 　　例 设有矩阵序列{Ak}, 其中Ak=Bk,而 且设|λ|1，考查矩阵序列极限. 　解 显然, 当|λ|1时, 则有 　　矩阵序列极限概念可以用矩阵算子范数来描述. 定理8.1 　　　　　　　　　　　　　其中||·||为矩阵的任意一种算子范数. 定理8.2 设B=(bij)∈Rn×n，则 limBk=0 (k→∞)(零矩阵)的充分必要条件是矩阵B的谱半径?(B)1. 定理8.3 (迭代法基本定理) 设有方程组 x=Bx+f. 及一阶定常迭代法 x(k+1)=Bx(k)+f. 对任意选择初始向量x(0)，迭代法收敛的充要条件是矩阵B的谱半径?(B)1. 定理 8.2 和定理 8.3 的结论和起来即为 (1) 迭代法x(k+1)=Bx(k)+f 收敛?limBk=O； (2) 迭代法x(k+1)=Bx(k)+f 收敛??(B)1. 推论 设Ax=b，其中A=D-L-U为非奇异矩阵且D非奇异矩阵，则有 (1) Jacobi迭代法收敛??(J)1, 其中J=D-1(L+U). (2) G-S迭代法收敛??(G)1, 其中G=(D-L)-1U. (3) SOR迭代法收敛??(Lω)1, 其中 Lω=(D-ωL)-1[(1-ω)D+ωU]. 例： 考察用Jacobi方法解系数矩阵为A的方程组的收敛性. 解 迭代矩阵J为 得迭代矩阵J的特征方程为 解得 即?(J)1. 所以用Jacobi方法解方程组是收敛的. 例：当方程组的系数矩阵为 可以验证：Jacobi迭代收敛，而Gauss-Seidel迭代不收敛。 例：当方程组的系数矩阵为 可以验证：Jacobi迭代不收敛，而Gauss-Seidel迭代收敛。 例 考察用迭代法解方程组的收敛性. 其中 解 方程组的迭代矩阵B的特征方程为 矩阵B的特征值为 　　　　即?(B)1. 这说明用迭代法解此方程组不收敛. 迭代法的基本定理在理论上是重要的，根据谱半径的性质?(B)≤||B||，下面利用矩阵B的范数建立判别迭代法收敛的充分条件. 8.3.2 关于解某些特殊方程组迭代法的收敛性 在科学及工程计算中，要求解方程组Ax=b，其矩阵A常常具有某些特性. 例如，A具有对角占优性质或A为不可约阵，或A是对称正定阵，下面讨论用基本迭代法解这些方程组的收敛性. 定义8.4 (对角占优阵) 设A=(aij)n×n . (1) 如果A的元素满足 称A为严格(按行)对角占优阵. (2) 如果A的元素满足 且上式至少有一个不等式成立，称A为弱(按行)对角占优阵. 定理 8.6 (对角占优定理) 如果A=(aij)n×n为严格对角占优矩阵,则A为非奇异矩阵. 定理 8.7 设方程组Ax=b，如果 A为严格对角占优阵，则解Ax=b的Jacobi迭代法, Gauss-Seidel 迭代法均收敛. 定义8.3 迭代法的渐近收敛速度定义为 计算方法 * * ? 2009, Henan Polytechnic University * §3 迭代法的收敛性 第八章 解线性方程组的迭代法 * * ? 2009, Henan Polytechnic University * §3 迭代法的收敛性 第八章 解线性方程组的迭代法 计算方法