数值方法第三章节1 线性方程组.pptVIP

基本思想：每次消元之前在系数矩阵中按一定的范围选取绝对值最大的元素作为主元素，以便减少舍入误差的影响。 交换原则：通过方程或变量次序的交换，使在对角线位置上获得绝对值尽可能大的系数作为akk(k)，称这样的akk(k) 为主元素，并称使用主元素的消元法为主元素法 根据主元素选取范围分为：列主元素法、行主元素法、全主元素法 列主元素法：在待消元的所在列中选择主元，经方程的变换，置主元素于对角线位置后进行消元的方法。 全主元素：在全体待选系数中选取主元，则得全主元素法。 记笔记 3.2.4 高斯主元素消去法 主元素法的意义 例3.2 用高斯消去法求下列方程组的解 解： 确定乘数 ，再计算系数 假设计算在4位浮点十进值的计算机上求解,则有 这时方程组的实际形式是 由此回代解出 ,但这个解不满足原方程组,解是错误的。这是因为所用的除数太小使得上式在消元过程中“吃掉”了下式，解决这个问题的方法之一就是采用列选主元高斯消元法。即按列选绝对值大的系数作为主元素，则将方程组中的两个方程相交换，原方程组变为 得到消元后的方程组 这时 因而方程组的实际形式是 由此回代解出 ,这个结果是正确的 可见用高斯消去法解方程组时,小主元可能导致计算失败,因为用绝对值很小的数作除数,乘数很大,引起约化中间结果数量级严重增长,再舍入就使得计算结果不可靠了,故避免采用绝对值很小的主元素。以便减少计算过程中舍入误差对计算解的影响。 每一步选绝对值最大的元素为主元素，保证 。 Step k: ① 选取 ② If ik ? k then 交换第 k 行与第 ik 行; If jk ? k then 交换第 k 列与第 jk 列; ③ 消元 全主元素法不是按列选主元素，而是在全体待选系数中选取，则得全主元素法。 例3.3 用全主元素法解下列线组 10x1 - 19x2 - 2x3=3 (1) -20x1 +40x2 + x3 =4 (2) x1 + 4x2 + 5x3=5 (3) 解：选择所有系数中绝对值最大的40作为主元素，交换第一、二行和交换第一、二列使该主元素位于对角线的第一个位置上，得 40x2 - 20x1 + x3 =4 (4) -19x2+10x1 - 2x3=3 (5) 4x2+ x1 +5x3=5 (6) 记笔记 计算m21=-19/40=0.475，m31=4/40=0.1 (5)- m21(4), (6)- m31(4) 消去x2 得 0.5x1 –1.525 x3 =4.9 (7) 3x1 + 4.9 x3 =4.6 (8) 选4.9为主元素 4.9 x3 + 3x1=4.6 (9) 1.525 x3 +0.5x1=4.9 (10) 计算m32=-1.525/4.9=-0.31122， (10)- m32(9)消去x2得 1.43366x1=6. 33161 (11) 记笔记 保留有主元素的方程 40x2 - 20x1 + x3 =4 (4) 4.9x3 + 3x1=4.6 (9) 1.43366x1=6. 33161 (11) 进行回代 x1=4.41634 x3 =-1.76511 x2=2.35230 3.2.4.1 列主元素法 列主元素法就是在待消元的所在列中选取主元，经方程的行交换，置主元素于对角线位置后进行消元的方法。 即： 在高斯消元第k步之前，做如下的事情： 若 交换 k 行和 j 行 行的交换，不改变方程组的解，同时又有效地克服了高斯消元的缺陷。 例3.4 用列主元素法解下列线性方程组 10x1 - 19x2 - 2x3=3 (1) -20x1 +40x2 + x3 =4 (2) x1 + 4x2 + 5x3=5 (3) 解：选择-20作为该列的主元素， -20x1 +40x2 + x3 =3 (4) 10x1 - 19x2 - 2x3=4 (5) x1 + 4x2 + 5x3=5 (6) 计算m21 =10/-20=-0.5 m31=1/-20=-0.05 (5)- m21(4), (6)- m31(4)得 x2 – 1.5x3=