5.1Guass消去法.pptVIP

华长生制作 * 第五章 解线性代数方程组的直接法 5.1 高斯消去法 　　在科学与工程计算中，大量的问题归结为求解线性代数方程组Ax=b，其中 分别称为方程组的系数矩阵，右端向量和解向量，若A可逆，则方程组存在唯一解。 然而，如果使用克莱姆法则求解n阶方程组所需的乘除运算量为(n+1)!(n-1)+n次，当n=20时，用每秒上亿次运算速度的计算机也要30万年之多！所以必须研究好的算法来解决这个问题，下面介绍高斯消去法。 我们把方程组Ax=b写成 设方程组的系数矩阵A非起奇异,记 这样,方程组又可写成 。消元过程就是要按确定的计算过程对方程组进行初等行变换,将方程组化为上三角方程组. 第一步消元:假设 　　,作初等行变换运算 将 　变换为 它对应的方程组 　与原方程组等价,而在第2至第n个方程中,含 的项已经消去. 　　第k步消元:设消去法已进行k-1步,得到方程组 ,此时对应的增广矩阵是 假设 令 将 变换为 　， 　中第1至第k行的元素与中的对应元素相同，它对应的方程组 　与原方程组等价，而在k+1至第n个方程中，含 　的项已经消去。 上述过程可进行n-1步，得到方程组 其中系数矩阵是一个上三角阵，原方程组等价转换为上三角方程组 这就完成了消元过程。 　　因为A非奇异，所以可求解上三角方程组，通过逐次代入计算可得方程组的解，其计算公式为 求解上式的过程称为回代过程。 　　以上由消去过程和回代过程合起来求解方程组的过程就称为Gauss消去法，或称为顺序Gauss消去法。 由前面可知,消元过程的第k步需出发运算n-k次，乘法和减法运算各需（n-k）（n+1-k）次，所以消元过程共需乘除法运算的次数为 需加减法运算的次数 回代过程需乘除法运算n（n+1）/2次，加减法运算n（n-1）/2次。所 以，Gauss消去总共需乘除法预示内的次数为 该运算量较克莱姆法则法要少很多！ 例 用Gauss消去法解方程组 解 第一步消元，令 　　得增广矩阵 第二步消元，令 　　得增广矩阵 利用回代公式依次得到 　　在这个例子中我们写出的是分数运算的结果。如果在计算机上进行计算，系数矩阵和中间结果都用经过舍入的机器数表示，中间结果和方程组的解都会有误差。 例1. 用Gauss消去法解线性方程组(用3位十进制浮 点数计算） 解: 本方程组的精度较高的解为 用Gauss消去法求解(用3位十进制浮点数计算) Gauss列主元消去法的引入 9999 回代后得到 与精确解相比,该结果相当糟糕 究其原因,在求行乘数时用了很小的数0.0001作除数 主元 如果在求解时将1,2行交换,即 0.9999 回代后得到 这是一个相当不错的结果 例2. 解线性方程组(用8位十进制尾数的浮点数计算) 解: 这个方程组和例1一样,若用Gauss消去法计算会有 小数作除数的现象,若采用换行的技巧,则可避免 绝对值最大 不需换行 经过回代后可得 事实上,方程组的准确解为 例2所用的方法是在Gauss消去法的基础上,利用换行 避免小主元作除数,该方法称为Gauss列主元消去法 行变换相 当于左乘 初等矩阵 　　列主元素消去法也称按列部分主元的消去法。一般地，在完成 了第k-1步消元运算后，在 的第k 列元素 之下的所有 元素中选一个绝对值最大的元素作为主元素，即若 则以 为主元素，这里 ,且 由于非奇异， 有 　.这样， 有达到控制舍入误差的作用。 　　选出主元素后，若则进行顺序Gauss消去法的第k步若 ，则将 　的第 行与第k行交换，然后进行消元运算。 　　完成了n-1步主元，换行与消元运算后，得到