线性方程组解法综述.pdfVIP

线性方程组解法的研究综述 摘要： 这篇论文在说明了线性方程组的应用目的的基础上，提出了线性方程组 求解的研究现状，并列举了常用的求解方法，同时说明了它们的应用条件，剖析 了各种方法的不足之处。 关键词 高斯消元 迭代 病态方程组 一、问题提出 在自然科学和工程实际应用中，有许多问题的求解最终都转化为线性方程组 的求解问题。例如, 电学中的网络问题，曲线拟合中常用的最小二乘法、样条函 数插值、解非线性方程组、求解偏微分方程的差分法、有限元法和边界元法以及 目前工程实践中普遍存在的反演问题等。特别是在图像恢复、模型参数估计、解 卷积、带限信号外推、地震勘探等众多领域，都需要求解线性方程组。 由于线性方程组问题在理论上的重要性和在工程实际应用中的大量存在，多 年来人们在这方面做了广泛深入的研究和探讨，并取得了许多有价值的成果. 由 于模型误差、测量误差、计算误差等各种误差的存在，常常使得线性方程组中的 系数矩阵和非齐次项信息具有某种程度的近似性(即扰动性) ，这种近似性显然会 使得线性方程组的求解不容易得到真实的理论解。此时，不同的求解方法由于运 算机理不一样，求解过程中误差积累程度就不一样，因此必然会使得不同的求解 方法得到的解具有不同的逼近真解的误差程度，尤其对具有病态性的方程组而言 由于病态线性方程组的条件数很大，数据误差以及计算过程中引入的舍入误差往 往会使线性方程组的解不稳定，即不管原始数据的误差多么小，都可能造成解的 很大变化，使线性方程组的解严重失真。因此，许多现有的方法都是无效的，病 态线性方程组的求解变得相当困难。求解线性方程组的最常用的方法主要有直接 法和迭代法两大类，其中直接法中最常用的方法是高斯消元法。但是，该方法求 解病态线性方程组时不能得到合理的解，误差很大。 二、研究现状 目前关于线性方程组的数值解法一般有两大类。一类是直接方法，另一类是 迭代方法。直接方法最基本的是高斯消元法及其变形，这类方法是解低阶稠密矩 阵方程组的有效方法，近十几年来直接法在求解具有较大型稀疏矩阵方程组方面 取得了较大进展。迭代法就是用某种迭代过程去逐步逼近线性方程组的精确解 迭代法具有需要计算机的存储单元较少，程序设计简单，原始系数矩阵在计算过 程中始终不变等优点，但存在收敛性及收敛速度问题。迭代法是解大型稀疏矩阵 方程组的重要方法。当前对迭代算法的研究已经较为成熟，但如何使之适合新体 系模型，以获得更好的性能加速一直是应用和体系设计者关心的问题。 三、常用方法比较 1.直接方法 直接方法是指假设计算过程中不产生舍入误差，经过有限次运算可求得方程组的 精确解的方法。事实上，由于舍入误差的存在，用直接法一般也只能求得方程组 的近似解。直接方法中主要有三种方法：克拉默法则、高斯消元法、LU 分解 法。 （1）克拉默法则 设有线性方程组（ n 个未知数 n 个方程） a x  a x  ... a x  b  11 1 12 2 1n n 1  a x  a x  ... a x  b  2 1 1 22 2 2n n 2     a x  a x  ... a x  b  n1 1 n 2 2 nn n n 其矩阵形式为 Ax = b 其中 a11 a12  a1n  x1  b1        a a  a x b  2 1 22 2n   2   2  A  , x  ,b                a a  a x b  n1 n 2 nn   n   n  如果线性方程组的系数行列式不为零，即 det( A) ≠ 0 ,则该方程组有唯一解