对病态方程组的处方法研究.pdf

对病态方程组的处理方法研究 蓝醒龙 （广西民族大学数学与计算机科学学院03数本2班 530006） 摘 要： 对病态线性方程组解法研究是数值计算方法的一个重要研究课题。本文分析 了病态方程组的特点，介绍了几种有效的解法。 关键词： 病态线性方程组；条件数；预处理；迭代 StudyingTheAlgorithmForSolvingIll-conditioned SystemOfEquations Abstract： Studying the algorithm for solving ill-conditioned system of equations is an important issue. This paper analyses the equations characteristic, and introduces several effectivealgorithms. Keywords：ill-conditionedsystemofequations conditionnumber pretreatment iteration 1 问题的提出 一个线性方程组 若右端向量 或系数矩阵 的微小变化就会引起方程组的 AX b b A 解发生很大的变化，则称 为病态方程组。方程组的系数矩阵 的条件数 AX b A Cond A  A A1 刻画了方程组的性态，若Cond A 1，则称AX b为 “病态”方     程组；若Cond A 相对较小，则称AX b为 “良态”方程组。良态方程组用GAUSS消   去法和JACOBI等简单的迭代法就可以得到比较好的计算解 而对于病态方程组 一般的直 接法和迭代法会有较大的误差 甚至严重失真。所以 在解方程组时 有必要先对方程组的 性态进行研究 采用相应的算法 才能得到比较精确的计算解。利用方程组的条件数来判断 就是一个很好的办法。 下面的一些直观的现象可作为判别病态矩阵的参考： （1）在主元消去法的过程中出现小主元，则 有可能是病态矩阵，但病态矩阵未必一 A 定有这种小主元； 1 （2）若解方程组时出现很大的解，则 有可能是病态矩阵，但病态矩阵也可能有一个 A 小解； （3）从矩阵本身来看，若元素间数量级相差很大且无一定规律；或矩阵的某些行（列） 近似线性相关，即矩阵的行列式接近于0，这样的矩阵就有可能是病态的。 当然 这些现象只能帮助我们做初步的判断 并且很多病态矩阵也不一定会出现这些现 象。所以，最可靠的判别方法是求出矩阵的条件数。 在文[1]中，引用了Forsythe和Moler 的一种估计Cond A 的简便方法：若残向量        bAx x表示用Gauss消去法求解AX b 的计算解 y 表示用Gauss消去法求解 AY  的计算解，则 ||y ||ˆ  i Cond (A) 10