数值分析几种常用的迭代法.ppt

6.3几种常用的迭代法记，A非奇异，且对角元，可以把A分解为其中雅可比迭代法华长生制作*方程组Ax=b等价于由此构造迭代公式：其中迭代距阵和向量为称之为Jacobi迭代法（简称J法），称为雅可比迭代矩阵。华长生制作*雅可比法的分量形式为由前面的定理知雅可比迭代关于任意初始向量收敛的充要条件为，充分条件为利用这些判别J法的收敛性，有时不太方便，对于大型方程组，要求出迭代矩阵谱半径是不容易的。下面给出一些容易验证收敛性的充分条件，先讨论对角占优矩阵的性质。华长生制作*定义1若满足则称A为严格对角占优矩阵。若满足且其中至少有一个严格不等式成立，则称A为弱对角占优矩阵。华长生制作*定义2设，若A不能经过行置换与相应的列置换化为其中和均为方阵，则称A为不可约的，否则称A为可约的。定理若A为严格对角占优矩阵，或不可约的弱对角占优矩阵，则解方程组的J法关于任意初始向量收敛。设，这里只给出A为严格对角占优阵时的证明。对J法，迭代矩阵，易得。由A的严格对角占优性，得到，所以J法收敛。证华长生制作*与雅可比法相应的高斯-赛德尔迭代法在J法中，计算时，分量已经算出，所以可考虑在J法中的求和分成两部分，从而得到与雅可比迭代法相应的高斯-赛德尔迭代法为这就是Gauss-Seidel迭代法，简称GS法。华长生制作*将上式写成距阵形式整理为简单迭代的形式其中迭代矩阵和向量为Jacobi迭代法和Gauss-Seidel迭代法的分量形式供计算编程用，它们的矩阵形式供研究迭代序列是否收敛等理论分析用。华长生制作*例用Ｊ法和ＧＳ法分别求解方程组。解用J法计有华长生制作*用GS法计算有取，J法迭代4次的计算结果是GS法迭代4次的计算结果是精确解为（1，1，1），从计算结果看，本例用GS法显然比用J法收敛快，但并不是任何时候GS法都比J法快，甚至有J法收敛而GS法不收敛的例子。华长生制作*显然，高斯-赛德尔法关于任意初始向量收敛的充要条件是另外与雅可比法相仿有如下结论：定理若A为严格对角占优矩阵，或不可约的弱对角占优矩阵，则解方程组Ax=b的GS法关于任意初始向量收敛。华长生制作*例.判别下列方程组用J法和G-S法求解是否收敛解:(1)求Jacobi法的迭代矩阵华长生制作*因此不能用范数判断所以即Jaobi迭代法收敛(2)求Gauss-Seidel法的迭代矩阵华长生制作*所以Gauss-Seidel迭代法发散华长生制作*无论是解线性方程组的Jacobi迭代法和G—S迭代法都涉及到收敛速度问题如何加快迭代法的速度呢？如何改善迭代法的适用范围呢？逐次超松弛（SOR）迭代法华长生制作*考虑解线性方程组的Gauss-Seidel迭代法------(1)华长生制作*令因此------(2)华长生制作*上式称为逐次超松弛法(SOR迭代法),逐次超松弛法(SOR迭代法)的矩