大型线性代数方程组解法的探讨.pdfVIP

第 29 卷第 2 期 河 海 大 学 学 报 Vol. 29 No. 2 2001 年 3 月 JOURNAL OF HOHAI UNIVERSITY Mar. 2001 大型线性代数方程组解法的探讨 1 2 徐小明 ,汪德　 ( 1. 河海大学数学物理系 ,江苏 南京 　210098 ; 2. 河海大学水资源水文及环境学院 ,江苏 南京 　210098) 摘要 :针对 SaintVenant 方程组离散后所形成的线性代数方程组的求解问题 ,将 Gauss 列主元消去法 与压缩存贮技术相结合 ,提出了存贮单元少、舍入误差小且数值计算稳定的计算方法 ,使得河网非 恒定流的数值计算更加高效 ,并且计算的精度可得到充分的保证. 关键词 :河网非恒定流 ;列主元消去法 ;压缩存贮 ( ) 中图分类号:O241. 6 　　　文献标识码 :A 　　　文章编号 2001 02 003805 利用 SaintVenant 方程组求解河网非恒定流时 ,通常有两种解法 :一种是由荷兰水力学家 Dronkers 提出 的后经许多学者发展的分级解法[1 ,2 ] ; 另一种是 Michael Amein[3 ] 提出的由 Fread [4 ,5 ] 等人发展的用 Newton Raphson 方法直接求解非线性方程组的方法. 无论是哪一种方法 ,对具有 N 个断面的单一河道 ,最终都要求 解一系列的线性代数方程组 Ax = b ( 1) ( ) 其中 A 为每行最多只有 4 个非零元的 2 N 阶五对角阵. 对方程组 1 的求解, 目前常用的方法有两种:其一 为消元法[4] ; 其二为追赶法[2] . 但这两种方法都有个共同的缺点, 这就是在消元或追赶过程中, 当出现分母为 ε 零或机器零的数时, 求解就不能顺利进行下去. 以往补救的办法是选择一个适当的正数 , 当某个作为分母 ε ε ( ) ( ) 的数 d 的绝对值小于 时, 就令 d = sng d 其中 sgn 为符号函数 , 以使方程组的求解能进行下去. 但当 A 的条件数cond ( A) = ‖A ‖‖A - 1 ‖较大时, 即对应的线性代数方程组为病态方程组时, 经这样一改动所得 的方程组的解与原方程组的解将会相差甚远, 毫无实际意义. 在河网非恒定流的数值模拟中, 对每一计算时 ( ) ( ) 刻的每一条河道都要求解形如 1 的方程组. 因此整个计算过程将要求解大量的形如 1 的方程组, 由于河网 水流及边界条件的复杂性, 难免出现 A 的条件数较大及运算过程中出现分母为零或接近于零的情形, 从而 产生误差较大的解甚至出现计算中断. 这种情况在实际问题的数值模拟中曾出现过, 它给数值计算提出了一 个不可回避的问题. 如何处理好这一问题, 将直接关系到河网水力数值模拟的成败. 针对这一问题, 本文充分 考虑到所求解的线性方程组系数矩阵的特点