9讲：ch4-1消元法.pptVIP

回代过程算法 利用MATH验证下列矩阵作用于同阶矩阵后的变换效果 验证效果 第2节 主元素法 主元素法 解:因为a11= 0 ,故此题不能用Gauss消元法求解,但交换方程组的顺序后,就可用Gauss消元法求解了. 问题1:考虑如下线性方程组的Gauss消元法 2x2=1 2x1+3x2=2 问题2:讨论下面方程组的解法 0.0001x1+x2=1 x1+x2=2 假设求解是用四位小数计算机上进行 问 题 0.1000?10-3 x1 + 0.1000 ?101 x2 = 0.1000 ?101 0.1000 ?101 x1 +0.1000 ?101 x2 = 0.2000 ?101 解:本题的计算机机内形式为 因为a11 =0.0001?0,故可用Gauss消元法求解,进行第一次消元时有 a22(1)= 0.1000 ?101 - 104 ? 0.1000 ?101 = 0.00001 ?105 - 0.1000 ?105 (对阶计算) = 0.0000 - 0.1000 ?105 = - 0.1000 ?105 , (m21=a21/a11=1/0.0001= 104) 得三角方程组 问 题 0.1000?10-3 x1 + 0.1000 ?101 x2 = 0.1000 ?101 -0.1000 ?105 x2 = -0.1000 ?105 回代解得x2=1, x1=0 严重失真! (因为本题的准确解为 x1=10000/9999, x2=9998/9999 失真的原因：除数的绝对值远远小于被除数的绝对值 问 题 用高斯消去法求解线性方程组时,应避免小的主 元.在实际计算中,进行第k步消去前,应该在第k列元 素 中找出绝大值最大者,例如 再把第p个方程与第k个方程组进行交换,使 成为主元.我们称这个过程为选主元.由于只在第k列元素中选主元,通常也称为按列选主元(或称列主元素法). 一、列主元的基本思想 如果在第k步消去前,在第k个方程到第n个方程所有的 xk到xn的系数 中,找出 绝对值最大者,例如 再交换第k,p两个方程和第k,q两个未知量的次序, 使 成为主元. 称这个过程为完全选主元或全主元素法. 不论是哪种方式选出主元,而后再按上面介绍的计算步骤进行消去的计算,一般都称为选主元的高斯消去法. 在实际计算中,常用按列选主元的高斯消去法. 完全选主元或全主元素法 用列主元消去法 解方程组 解 第一次消元对 因列主元素为a31(1),故先作行变换r1 r3,然后进行消元计算可得 例2.2 举 例 由此回代,得x=(1.9272,-0.69841,0.90038)T 与精确解x=(1.9273,-0.69850,0.90042)T相比较是 比较准确的. 第二次消元对[A(2) |b(2) ],因列主元素为a32(2) ,故先作行变换r2__r3,然后进行消元计算可得 举 例 用全主元素法求解线性方程组，计算过程保留三位小数 解： 例2.3 二、全主元素法举例 由回代过程得解 注意解的次序 二、全主元素法举例 例2.3的计算结果表明，全主元素法的精度略优于列主元素法，这是由于全主元素法是在全体系数中选主元，故它对控制舍入误差比较有效。 但全主元素法在计算过程中，需同时作行与列的互换，因而程序比较复杂，计算时间较长。列主元素法的精度虽稍低于全主元素法，但其计算简单，工作量大为减少，且计算经验与理论分析均表明，它与全主元素法同样具有良好的数值稳定性，故列主元素法是求解中小型稠密线性方程组的最好方法之一。 二、全主元素法举例 1. 用高斯消元法解方程组（取3位小数） 2. 用列主元素法解方程组： 作 业 作 业 AX=b 直接法 迭代法 是 列主元 消去法 Gauss 消去法 全主元 消去法 否 是 LU分解法 追赶法 A对称 且正定 平方根法