实验4求解线性方程组的高斯消去法 参考答案 L.带选主元步骤的高斯消去法:functionX=uptrbk(A,B)(见代码文件),其中需要调用求 解上三角形方程组的后代法:functionX=backsub(A,B)(见代码文件)。 2.令,为配制这种油漆所需要的一加仑不受欢迎的油漆j的百分比。那么,新油漆中色素 A、B、C和D的总量各自为 80x. +30x.+10x.=40 80x2+10x,+10x=27 16x+20x,+60x+72x,=31 +8x4=2 写成矩阵相乘的形式为AX=B: [80030107[x1「40 0801010x2 27 16206023 4008Lx4」L2」 利用第1题的程序求解得到:x=0.40.30.250.05 3.不带选主元步骤的追赶法:function[X,Y,L,U=chase((A,B)(见代码文件) 利用追赶法程序求解得到方程组的解为:x=[132-2] 需要注意的是,追赶法木质上是没有洗主元的高斯消去法。如果采用选主元(不管是列 选主元、按比例列选主元还是全选主元),则破坏了方程组的三对角性质,就不能用追赶法 求解,而只能按 般的方程组求解。由于追赶法中没有采用选主元,因此,对于 一般的三对 角线线性代数方程组来说,追赶法的计算过程是不稳定的。但是,当三对角矩阵满足对角占 优条件时,即 laul>la21>0 lal>a+a>0i=2.3.....n-1 a>a> 追赶法的计算过程中不会出现中间结果数量级的巨大增长和舍入误差的严重积累。而在工程 实际中的许多三对角线线性代数方程组是满足上述条件的,因此,追赶法是有实用价值的。 实验 4 求解线性方程组的高斯消去法 参考答案 1. 带选主元步骤的高斯消去法:function X = uptrbk(A,B)(见代码文件),其中需要调用求 解上三角形方程组的后代法:function X = backsub(A,B)(见代码文件)。 2. 令 xj 为配制这种油漆所需要的一加仑不受欢迎的油漆 j 的百分比。那么,新油漆中色素 A、B、C 和 D 的总量各自为 1 3 4 234 1 2 3 4 1 4 80 30 10 40 80 10 10 27 16 20 60 72 31 4 8 2 x x x x x x x x x x x x + + = + + = + + + = + = 写成矩阵相乘的形式为 AX=B: 1 2 3 4 80 0 30 10 40 0 80 10 10 27 16 20 60 72 31 4 0 0 8 2 x x x x = 利用第 1 题的程序求解得到:x = [0.4 0.3 0.25 0.05]’. 3. 不带选主元步骤的追赶法:function [X,Y,L,U]=chase(A,B) (见代码文件) 利用追赶法程序求解得到方程组的解为:x = [1 3 2 -2]’. 需要注意的是,追赶法本质上是没有选主元的高斯消去法。如果采用选主元(不管是列 选主元、按比例列选主元还是全选主元),则破坏了方程组的三对角性质,就不能用追赶法 求解,而只能按一般的方程组求解。由于追赶法中没有采用选主元,因此,对于一般的三对 角线线性代数方程组来说,追赶法的计算过程是不稳定的。但是,当三对角矩阵满足对角占 优条件时,即 11 12 , 1 , 1 , 1 0 0, 2,3, , 1 0 ii i i i i nn n n a a a a a i n a a − + − + = − 追赶法的计算过程中不会出现中间结果数量级的巨大增长和舍入误差的严重积累。而在工程 实际中的许多三对角线线性代数方程组是满足上述条件的,因此,追赶法是有实用价值的