求解线性方程组的一种迭代算法.pdf

( ) 第 25 卷第 3 期 暨南大学学报 自然科学版 　　 Vol . 25 No . 3 　　　 ( ) 　　　 2004 年 6 月 　　　 Journal of Jinan University Natural Science 　　 　Jun . 2004 求解线性方程组的一种迭代算法 常双领 , 　张传林 (暨南大学数学系 ,广东 广州 510632) [摘 　要] 　对系数为对称正定矩阵的线性方程组 ,利用系数矩阵主对角线上元素的和构造一种新 的收敛迭代格式. [ 关键词] 　线性方程组 ; 　迭代算法 ; 　对称正定矩阵 ( ) [ 中图分类号] 　O2416 　　[文献标识码] 　A 　　[文章编号] 　1000 - 9965 2004 03 - 0256 - 04 　　科学研究与生产实践中许多问题都可归结为线性方程组的求解. 随着计算机科学技术的 飞速发展 ,可用直接法求解的线性方程组的规模越来越大 ,但直接法大多数均需对系数矩阵 A 进行分解 ,一般不能保持稀疏性 ,而在实际应用中 ,经常遇到大规模稀疏线性方程组的求解问 题. 对于大规模线性方程组的求解问题 ,特别是大规模稀疏线性方程组 ,迭代法是求解线性方 程组的一种有效方法. 与直接法不同 ,迭代法是从一个初始向量出发 ,按照一定的规则经过迭 代得到一个向量序列 ,逐步逼近线性方程组的解. 迭代方法具有存储量小 、程序结构简单的特 点. 线性方程组 ( ) Ax = b 1 的迭代格式一般按如下方式构造, 首先将矩阵 A 分裂 , 即 A = A 1 + A 2 , 其中A 1 是非奇异的且 ( ) - 1 - 1 ( ) 易于求逆, 则方程组 1 可等价地写成 x = - A 1 A 2 x + A 1 b , 于是方程组 1 的迭代格式可构 造为 : x = A - 1A x + A - 1 b , 一些经典的求解线性方程组的方法如J acobi , Gauss - Seidel 都是 k +1 1 2 k 定常迭代法. 所谓定常迭代法就是迭代格式 x k +1 = Mx k + c (2) 其中, M 为不含参数的常数矩阵, 称为迭代矩阵, 这种方法称为定常迭代法. 不同的分裂方法, 可以构造不同的迭代算法, 文献[ 1 ～3] 中介绍了许多已有结果. 本文利用系数矩阵 A 的对角 线上元素的和给出了系数为对称正定矩阵的线性方程组 Ax = b 的一种新的定常迭代格式, 如 果系数矩阵 A 为可逆的非正定的矩阵, 可以