大型线性代数方程组解法的探讨-河海期刊网.pdfVIP

第 卷第 期 河 海 大 学 学 报 !% ! /:? ’!% C: ’! 年 月 !$ # J4KEC6L 4A M4M6N KCN/OE)NPQ *9 ’!$ 大型线性代数方程组解法的探讨 $ ! 徐小明 ，汪德 （ 河海大学数学物理系，江苏 南京 ； 河海大学水资源水文及环境学院，江苏 南京 ） $’ !$%( ! ’ !$%( 摘要：针对)*+,-./0,*,- 方程组离散后所形成的线性代数方程组的求解问题，将1*233 列主元消去法 与压缩存贮技术相结合，提出了存贮单元少、舍入误差小且数值计算稳定的计算方法，使得河网非 恒定流的数值计算更加高效，并且计算的精度可得到充分的保证’ 关键词：河网非恒定流；列主元消去法；压缩存贮 中图分类号： 文献标识码： 文章编号： （ ） 4!5$’ 6 $ $%( !$ ! #( 7 ! ! ! 利用)*+,-./0,*,- 方程组求解河网非恒定流时，通常有两种解法：一种是由荷兰水力学家89:,;093 提出 ［，］ ［］ ［，］ 的后经许多学者发展的分级解法 $ ! ；另一种是 +=*0? 6@0+, # 提出的由A90*B 5 7 等人发展的用 C0D-:,. 方法直接求解非线性方程组的方法 无论是哪一种方法，对具有 个断面的单一河道，最终都要求 E*F3:, ’ ! 解一系列的线性代数方程组 （） ! G # $ 其中 为每行最多只有 个非零元的 阶五对角阵 对方程组（）的求解，目前常用的方法有两种：其一 ! 5 ! ! $ ［］ ［］ 5 ! 为消元法 ；其二为追赶法 但这两种方法都有个共同的缺点，这就是在消元或追赶过程中，当出现分母为 零或机器零的数时，求解就不能顺利进行下去 以往补救的办法是选择一个适当的正数 ，当某个作为分母 ! 的数 的绝对值小于 时，就令 （）（其中 为符号函数），以使方程组的求解能进行下去 但当 # # $ 3,H # 3H, ! !