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

第三节 迭代法的收敛性 * ? 2009, Henan Polytechnic University * §3 迭代法的收敛性 第八章 解线性方程组的迭代法 * 　迭代矩阵B满足什么条件时，由迭代法产生的向量序列{x(k)}收敛到x*. 引进误差向量 可得误差向量的递推公式 则研究迭代法收敛性问题就是要研究迭代矩阵B满足什么条件时，有. 8.3.1 迭代法的收敛性 定义 设有矩阵序列 Ak=(aij(k))∈Rn×n 及A=(aij)∈Rn×n，如果n2个数列极限存在且有 则{Ak}称收敛于A，记为lim (k→∞). 　　例 设有矩阵序列{Ak}, 其中Ak=Bk,而 且设|λ|1，考查矩阵序列极限. 　解 显然, 当|λ|1时, 则有 　　矩阵序列极限概念可以用矩阵算子范数来描述. 定理1 　　　　　　　　　　　　　其中||·||为矩阵的任意一种算子范数. 　　定理2 定理3 设B=(bij)∈Rn×n，则limBk=0 (k→∞)(零矩阵)的充分必要条件是矩阵B的谱半径?(B)1. 　　证明 由矩阵B的若当标准形,存在非奇异矩阵P使 其中若当(Jordan)块 且　　　　，显然有 其中 显然有, Et,0=I, Et,k=0(当k≥t)，(Et,1)k= Et,k. 由于Ji=λI+Et,1 ，因此 　　下面考查Jik的情况. 引进记号 其中 定理4(迭代法基本定理) 设有方程组 x=Bx+f. 及一阶定常迭代法 x(k+1)=Bx(k)+f. 对任意选择初始向量x(0)，迭代法收敛的充要条件是矩阵B的谱半径?(B)1. 　　证明 充分性. 设?(B)1，则(I-B)可逆，(I-B)x=f有唯一解，记为x*，则 x*=Bx*+f. 误差向量 由设?(B)1，应用定理3，有　 　. 于是对任意x(0)有 ，即 . 其中x(k+1)=Bx(k)+f . 显然，极限x*是方程组的解，且对任意x(0)有 　　必要性. 设对任何x(0)有 由定理2知 再由定理3，即得?(B)1. 　　定理4是一阶定常迭代法的基本理论. 定理3和定理4的结论和起来即为 (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方法解方程组是收敛的. 例 考察用迭代法解方程组的收敛性. 其中 解 方程组的迭代矩阵B的特征方程为 矩阵B的特征值为 　　　　即?(B)1. 这说明用迭代法解此方程组不收敛. 迭代法的基本定理在理论上是重要的，根据谱半径的性质?(B)≤||B||，下面利用矩阵B的范数建立判别迭代法收敛的充分条件. 定理5(迭代法收敛的充分条件) 设有方程组 x=Bx+f， B=(bij)∈Rn×n， 及一阶定常迭代法 x(k+1)=Bx(k)+f. 如果有B的某种算子范数||B||=q1 ，则 (1) 迭代法收敛，即对任取x(0)有 　　证明 (1)由基本定理4结论(1)是显然的. 　　(2) 显然有关系式x*-x(k+1)=B(x*-x(k)) 及 x(k+1) –x(k)=B(x(k) –x(k-1)). 于是有 (a) ||x(k+1) –x(k