改进Gauss消去法求解线性方程组.docx

文章编号 1000 - 5269 (2004) 02 - 0127 - 05Ξ改进 Ga uss 消去法求解线性方程组胡尧 ,罗文俊( 贵州大学 数学系 ,贵州 贵阳 550025)摘要研究了线性方程组解的误差起源 ,在 Gauss 消元过程中避开除法 ,切断由于消元过程中系数相除所产生的舍入误差 ,用改进的 Gauss 消去法求解线性方程组 ,大大提高了 线性方程组解的精确值 。关键词线性方程组 ; Gauss 消去法 ;精确值 ;除法中图分类号O241文献标识码A引言在科学和工程计算中 ,线性方程组0a ( 1)( 1)( 1)( 1)x 1 + a12 x 2 ++ a1 n x n==b111a ( 1)( 1)( 1)( 1)x 1 + a22 x 2 ++ a2 n x nb221(1)a ( 1)( 1)( 1)( 1)x 1 + a n2 x 2 ++ a n n x n=b nn1数值解的重要性是显而易见 ,而求解线性方程组 ( 1) 的计算方法一般分为两类 : 直接法和迭代法 。不管是哪一类算法都只能在预定的计算步数内或给定的精度内得到近似解 ,很少得到精确解 ,也就是说所有 的算法都有误差问题 。误差的来源是什么 ? 就是在求解线性方程组的过程中 ,系数相除所产生的舍入误差累积带入了未知 量的直接求解式 ,导致了线性方程组解的误差 ,即求解线性方程组 ———误差来源于除法 。如果我们在求解 过程中不使用或尽可能少使用除法 ,或对于除法采取分数代入 (因为计算机的字长总是有限的) ,误差就可 以完全消除 ,或达到误差最小 。本文在不考虑计算量的情况下 ,将求最大公因式中的辗转相除法融入到 Gauss 消去法中 ,在归一消元化为等价同解的上三角形方程组a ( 1)( 1)( 1)( 1)x 1 + a12 x 2 ++ a1 n x n = b111ωa ( n)( n)n n x n= b n的过程中 ,将系数相除取整 ,避开了除法运算所产生的舍入误差 ,消除了消元过程中除法造成的误差累积 ,大大提高了线性方程组解的精确值 。1 传统的 Ga uss 消去法求解线性方程组文献1 Gauss消去法求解线性方程组的步骤如下 :收稿日期 :2003 - 12 - 29作者简介 :胡 尧( 1972 - ) ,男 ,贵州大学数学系讲师。研究方向 :数值算法 ;罗文俊( 1966 - ) ,男 ,贵州大学数学系副教授。研究方向 :计算机算法、编码理论、密码学等。Ξ贵州大学学报 (自然科学版)第 21 卷·128 ·akj ]a , ( j = k + 1 ,, n)kjakk①对于 k 从 1 到 n - 1归一化bk ]bkakkaij - ai k ·akj ] aij , ( i , j = k + 1 ,, n )消去bi - aik ·bk ] bi , ( i = k + 1 ,, n)bn ]②回代xna n nnbi - ∑aij ·x j ] x i , ( i =n - 1 ,, 2 , 1)j = i +1Gauss 消去法在归一过程中每执行一次就要作 n - k + 1 次系数相除 , 相应地产生了n - k + 1 个舍入误差 ,产生的舍入误差又带入随后的乘法运算 , 累积所产生的误差个数就为 ( n -k + 1) ( n - k) 。因此 ,仅归一nn1消元过程就产生了 ∑( n -k + 1) ( n - k) 个累积误差 ,再加上回代过程中所作的 ∑( n - k) =n ( n - 1)2k = 1k = 13次乘法所产生的累积误差 ,整个 Gauss 消去法中总共产生的误差个数(加减法不会产生误差) 为 n + n2 -n 。33文献2 尽管每次乘除所产生的误差符号不会都相同 ,从而使误差在累积时发生相互抵消现象 ,但累积所产生的误差是相当大的 。现仅以 2 阶线性方程组为例说明 :7 x 1 + 2 x 2 = 11例 1方程组消元Gauss 消去法求解 :3 x 1 + 4 x 2 = r 2 - 0143 r 1 721172116127 a21 3( 因为=≈0143)=704 - 0132 ×211 - 0143 ×1103114a11×11996x 2 = 6127≈11996 , x 1 = 11 - 2回代得≈1 .方程组的精确解为 x 1 = 1 , x 2 = 2 , 仅 2 阶就产生了这么大的误差原因何在 ? 就在于系数 a21 与a11不能整除或在计算机允许的范围内不能除尽 , 将所产生的舍入误差带入了