第六章数值分析.ppt

第六章

解线性代数方程组的迭代法§1向量和矩阵序列的极限§2迭代法的基本理论§3几种常用的迭代法三、逐次超松弛迭代法（SOR法）1、迭代格式①SOR迭代格式的分量形式将Gauss-Seidel迭代格式的分量形式改写为：一般一、Jacobi迭代法1、迭代格式将方程组变形为：①Jacobi迭代格式的分量形式该迭代格式称为Jacobi迭代格式。Jacobi迭代格式的分量形式也可写为：②Jacobi迭代格式的矩阵形式设A=D+L+U则因，故D可逆，于是Jacobi迭代格式的矩阵形式为令则得称为Jacobi迭代格式的迭代矩阵。2、收敛条件①充要条件Jacobi迭代关于任意初始向量收敛的充要条件是②充分条件（I）若,则Jacobi迭代关于任意初始向量都收敛。（II）设系数矩阵严格对角占优，即（按行）或（按列）则Jacobi迭代关于任意初始向量都收敛。证明（1）当A为按行严格对角占优，由得即故Jacobi迭代关于任意初始向量都收敛。（2）当A为按列严格对角占优，则可证得（略）故Jacobi迭代关于任意初始向量都收敛。定义设，若A不能经过行置换与相应的列置换化为其中是方阵。则称A为不可约矩阵。或若可找到置换矩阵P，使A呈则称A为可约矩阵。以n阶单位矩阵的n个列向量为列构成的n阶矩阵称为置换矩阵。这里是1,2,……,n的一个排列。如P为四阶置换矩阵定义设，若且其中至少有一个严格不等式成立，则称A为按行弱对角占优。（III）设A为不可约弱对角占优矩阵，则Jacobi迭代关于任意初始向量都收敛。（证明略）二、与Jacobi迭代相应的Gauss-Seidel迭代1、迭代格式①Gauss-Seidel迭代格式的分量形式由Jacobi迭代格式：将Jacobi迭代格式中得到的分量用最新计算的数值取代，即该迭代格式称为Gauss--Seidel迭代格式。Gauss-Seidel迭代格式的分量形式也可写为：②Gauss-Seidel迭代格式的矩阵形式设A=D+L+U则因，故D+L可逆，于是Gauss-Seidel迭代格式的矩阵形式为令则得 称为Gauss-Seidel迭代格式的迭代矩阵。2、收敛条件①充要条件Gauss-Seidel迭代关于任意初始向量都收敛的充要条件是②充分条件(I)若，则Gauss-Seidel迭代关于任意初始向量都收敛。（II）设系数矩阵严格对角占优，则Gauss-Seidel迭代关于任意初始向量都收敛。证明（1）当当A为按行严格对角占优，注意即得故Gauss-Seidel迭代关于任意初始向量都收敛。（2）当A为按列严格对角占优，则可证得（略）故Gauss-Seidel迭代关于任意初始向量都收敛。（III）设A为不可约弱对角占优矩阵，则Gauss-Seidel