线性方程组的四种数值解法.pdfVIP

线性方程组的四种数值解法 ( 电子科技大学物理 电子学院，四川 成都 610054) 摘要： 本文介 绍 了四种求解 线性方程 组 的数 值 解法 : 雅克比迭代法、高斯 赛德 尔迭代法、高斯消去法 和改 进 的平方根法的基本原理和算法流程，通 过求解具体方程， 对 四种求解方法 进 行了 对 比。 对 于雅克比 迭代法和高斯 赛德 尔迭代法，研究了两种算法 对求解同一方程 组 的迭代效率差异， 结果表明高斯 赛德 尔迭 代法达到同 样 精度所需迭代次数 较少。 对 于高斯消去法，通 过选择 列主元的方法提高算法的准确度， 计算 结果表明高斯消去法 计 算精确，且运算复 杂度也不是很高。 对于改 进 的平方根法，其运算复 杂度低，但 对 于给 定的方程 组有着 严 苛的要求。 关键词： 雅克比迭代法；高斯 赛 德尔 迭代法；高斯消去法；改 进 的平方根法； 线性方程 组 引言 线性方程组的求解在日常生活和科研中有着极其重要的应用， 但在实际运算中， 当矩阵 的维数较高时， 用初等方法求解的计算复杂度随维数的增长非常快， 因此， 用数值方法求解 线性方程组的重要性便显现出来。 经典的求解线性方程组的方法一般分为两类： 直接法和迭 代法。前者例如高斯消去法 , 改进的平方根法等，后者的例子包括雅克比迭代法，高斯赛德 尔迭代法等。这些方法的计算复杂度在可以接受的范围内，因此被广泛采用。一般来说，直 接法对于阶数比较低的方程组比较有效； 而后者对于比较大的方程组更有效。 在实际计算中， 几十万甚至几百万个未知数的方程组并不少见。在这些情况下，迭代法有无可比拟的优势。 另外， 使用迭代法可以根据不同的精度要求选择终止时间， 因此比较灵活。 在问题特别大的 时候，计算机内存可能无法容纳被操作的矩阵， 这给直接法带来很大的挑战。 而对于迭代法， 则可以将矩阵的某一部分读入内存进行操作， 然后再操作另外部分。 本文使用上述四种算法 求解对应的方程组，验证各种算法的精确度和计算速度。 1 算法介绍 1.1 雅克比迭代法 1.1.1 算法理论 设线性方程组 Ax b (1) a11 ,a22 ,...,ann 的系数矩阵 A 可逆且主对角元素 均不为零 , 令 D diag a11 ,a22 ,...,ann 并将 A 分解成 A A D D (2) 从而 (1) 可写成 Dx D A x b 令 x B x f 1 1 其中 1 1 B I D A, f D b 1 1 . (3) 以 1 为迭代矩阵的迭代法