第3章 解线性方程组的迭代方法.pptVIP

第1页，共40页，星期日，2025年，2月5日定义:设｛xk｝是Rn上的向量序列，又设x*＝（x1*，x2*,…，xn*)Ｔ是Rn上的向量.则称向量x*是向量序列｛xk｝的极限，若一个向量序列有极限,称这个向量序列是收敛的.向量序列的极限如果向量序列｛xk｝收敛于向量x*的充分必要定理1（i=1,2,…,n）条件是第2页，共40页，星期日，2025年，2月5日矩阵序列的极限定义:设｛Ak｝是上的矩阵序列.若存在矩阵则称矩阵A是矩阵序列｛Ak｝的极限，记为若一个矩阵序列有极限,称这个矩阵序列是收敛的.使得矩阵序列｛Ak｝收敛于矩阵A的充分必要定理2（i,j=1,2,…,n）条件是这里第3页，共40页，星期日，2025年，2月5日证：依次取x为，其中则所以定理3的充要条件是对任何x∈Rn，有设矩阵定理4，则的充要条件是ρ(A)1第4页，共40页，星期日，2025年，2月5日证：矩阵A相似于其Jordan标准型,即存在可逆矩阵P,使得J为对角分块矩阵(Ji称为Jordan块):其中:ni为特征值λi的重数,且n1+n2+…+nr=n由于第5页，共40页，星期日，2025年，2月5日所以而第6页，共40页，星期日，2025年，2月5日一、简单迭代思想设矩阵A可逆，把矩阵A分裂为则迭代过程B称为迭代矩阵。给定初值就得到向量序列第7页，共40页，星期日，2025年，2月5日定义：若称简单迭代法收敛，否则，称逐次逼近法不收敛或发散。问题：是否是方程组x=Bx+f的解？结论1：任意给定初始向量，若由迭代公式(1)产生的迭代序列收敛到，则是方程组x=Bx+f的解。证：又如何判定所给迭代格式(1)是否收敛哪？第8页，共40页，星期日，2025年，2月5日迭代法收敛的条件定理1：对任意初始向量，由（1）得到的迭代序列收敛的充要条件是迭代矩阵的谱半径证：因此结论2：设矩阵，则注：要检验一个矩阵的谱半径小于1比较困难，所以我们希望用别的办法判断迭代格式是否收敛。第9页，共40页，星期日，2025年，2月5日定理2：若迭代法的迭代矩阵满足（矩阵的某一种算子范数），则迭代格式产生的序列收敛于x=Bx+f的精确解x*，且有误差估计式：证：由定理1、结论1和知迭代格式产生的序列收敛于x=Bx+f的精确解x*。且第10页，共40页，星期日，2025年，2月5日整理即得估计式。Remark:因为矩阵范数，都可以直接用矩阵的元素计算，因此，用定理2，容易判别迭代法的收敛性。定理2的条件只是充分的，而不是必要的，也就是说：如果，则迭代法收敛；但若，我们并不能断定迭代法就一定发散，此时需要用定理1来判定迭代法的敛散性。第11页，共40页，星期日，2025年，2月5日迭代格式的收敛速度与初始值x(0)有关，同时也与||B||和?(B)有关，一般来说，||B||和?(B)越小，收敛速度越快。Def：称为迭代法的渐近收敛速度。第12页，共40页，星期日，2025年，2月5日二、Jacobi迭代法例1：用迭代法解方程组解:将方程组化为等价形式:构造迭代格式:取初始值代入计算，得第13页，共40页，星期日，2025年，2月5日注：如何判断迭代过程终止？利用定理2的误差估计式可以判断迭代过程是否可以终止，但这种方法比较麻烦，通常采用的方法是通过前后两次迭代近似值的差来判断，即利用：终止迭代过程。上述这种求解方程组的方法称为Jacobi迭代法。第14页，共40页，星期日，2025年，2月5日Jacobi迭代法的步骤：3、判断迭代格式的收敛性。取初值x(0)带入计算。1、写出等价方程组—即将第i个方程的xi解出。2、写出相应的迭代格式