高斯消去法课稿.ppt

★ 高斯消元法运算量： 消元过程： 乘： 除： 回代过程： 运算量： ★ 高斯消元法的局限性： 在高斯消元过程中，我们假定了对角元素 由于每次消元时是按未知量的自然顺序进行的，而顺序消元不改变 A 的主子式的值，因此高斯消元法可行的充分必要条件为 A 的各阶主子式不为零。实际上只要detA≠0，方程组就有解。 2. 即使高斯消元法可行，如果 很小，运算中用它作分母会导致其它元素数量计的严重增长和舍入误差的扩散。 例 解方程组 精确解为：x1=1/3，x2=2/3。计算取5位有效数字。 解1 顺序消元： 解2 交换方程的顺序： 经高斯消元： 　　完成了n-1步主元，换行与消元运算后，得到 ，这是与原方程组等价的方程组,　　是一个上三角阵,再代回求解.这就是列主元素消去法的计算过程. 除了列主元素消去法外,还有一种完全主元素消去法.在其过程的第k 步 ，不是按列来选主元，而是在　　　右下角的n-k+1阶子阵中选主元 ,　即 然后将 　的第 行与第k行交换，将第 列与第k列交换,　同时将自变量 与　 的位置交换并记录自变量的排列次序.　直到消去法完成后,再按记录恢复自变量为自然次序.完全主元法比列主元法运算量大得多,由于列主元法的舍入误差一般已较小,所以在实际计算中多用列主元法. 　　例5.3 用列主元素消去法解方程组Ax=b,计算过程中五位有效数字进行运算,其中 解 记 　　. 第一步选列主元为 　，交换第1行与第3行，再消元计算得 第二步选列主元为 　,交换第2行与第3行，再消元计算得 消去过程至此结束。回代计算依次得到解 这个例题的精确解是 而用不选住主元的顺序Gauss消去法，则解得 这个结果误差较大，这是因为消去法的第1步中， 按绝对值比其他元素小很多所引起的。从此例看到去列主元素消法是有效的方法。 解：精确解为（舍入值）： 例5.4：用列主元消去法去解方程组 本例是具有舍入的4位浮点数进行运算，所得的计算解还是比较准确的。 列主元素消去法步骤： 设Ax=b。对于具有行交换的列主元素消去法，消元结果冲掉A，乘数mij冲掉aij，计算解X冲掉常数项b,行列式存在det. 1, 1→det, k=1 （对k=1—n-1 做2—7步。） 3， 若aik,k=0, 则 0→det. 计算终止。 5， 计算乘数mik, mik= aik/ akk→aik.(i=k+1—n) (| mik|≤1) 6, 消元计算： aij - mik aij→ aij (i,j=k+1,…,n) bi - mik bk→ bi (i=k+1,…,n) 优点： 数值稳定。 修正方法： 列主元高斯-约当（Gauss-Jordam）消去法。 缺点： 既消元；又回代。 5.2.4** Gauss-Jordan消元法（可不讲，自学） 高斯消去法有消元和回代两个过程，消去的是对角线下方的元素。当对消元过程稍加改变便可使方程组 化为对角阵 （5.2.9） 这时求解就不需要回代了,这种将主元素化为1,并用主元将其所在列的冗余元素全都消为0,即消去对角线上方与下方的元素,这种方法称为高斯-约当消去法,这时等号右端即为方程组的解 在第k步计算时，考虑对上述矩阵的第k行上、下都进行消元计算（k=1,2…n） 2，换行（当ik≠k）。交换[A,b]第k行与第ik行元素。 3，计算乘数 mik = -aik/akk (i=1—n,i ≠k), mkk =1/ akk (mik可存放在aik的单元中) （| mik |≤1） 上述过程全部执行完后有： 这表明用G—J方法将A约化为单位矩阵，计算解就在常数项位置得到，因此用不着回代求解。用G—J方法接方程组的计算量大约需要 次乘除法，要比Gauss消去法大些。 例5.4 用高斯-约当(Jordan)消去法求方程组的解 解 方程组相应的增广矩阵 列选主元 第五章方程组的直接解法 5.2 Gauss消去法 5.2.4 Gauss-Jordan消元法 5.2.3 主元素消去法