解线性方程组的列主元素高斯消去法和LU分解法.pdfVIP

. 数值试验报告分析 一、实验名称：解线性程组的列主元素高斯消去法和 LU 分解法 二、实验目的及要求： 通过数值实验，从中体会解线性程组选主元的必要性和 LU 分解法的优点，以及程组系 数矩阵和右端向量的微小变化对解向量的影响。 三、算法描述： 本次试验采用的是高斯列主元消去法和 LU 分解法求解线性程组的解。 其中， 高斯消去法的基本思想是避免接近于零的数作分母； 能进行到底的条件 :当 A 可逆 时，列主元 Gauss( 高斯）消去法一定能进行到底。 优点 :具有很好的数值稳定性；具有与顺序 Gauss 消去法相同的计算量。列主元 Gauss( 高 斯）消去法的精度显著高于顺序 Gauss( 高斯）消去法。 注意 :省去换列的步骤，每次仅选一列中最大的元。 矩阵的三角分解法是 A=LU,L 是下三角阵， U 是上三角阵 ,Doolittle 分解 :L 是单位下三角 阵， U 是上三角阵； Crout 分解 :L 是下三角阵， U 是单位上三角阵。矩阵三角分解的条件 是矩阵 A 有唯一的 Doolittle 分解的充要条件是 A 的前 n-1 顺序主子式非零；矩阵 A 有唯 一的 Crout 分解的充要条件是 A 的前 n-1 顺序主子式非零。三角分解的实现是通过 (1)Doolittle 分解的实现； (2)Doolittle 分解的缺点：条件苛刻，且不具有数值稳定性。 （3 ）用 Doolittle 分解求解程组 : AX=b LUX=b LY=b A=LU UX=Y ； 四、实验容： 解下列两个线性程组 3.01 6.03 1.99 x1 1 （1 ） 1.27 4 .16 1.23 x 2 1 0.987 4 .81 9.34 x 3 1 10 7 0 1 x1 8 3 2.099999 6 2 x2 5.900001 （2 ） 5 1 5 1 x3 5 2 1 0 2 x4 1 a 、用你熟悉的算法语言编写程序用列主元高斯消去法和 LU 分解求解上述两个程组，输出 ..页脚 . . Ax=b 中矩阵 A 及向量 b, A=LU 分解的 L 及 U，detA 及解向量 x. b 、将程组（ 1）中系数 3.01 改为 3.00 ，0.987 改为 0.990 ，用列主元高斯消去法求解变换 后的程组，输出列主元行交换次序，解向量 x 及 detA ，并与（ 1 ）中结果比较。 c 、将程组（ 2 ）中的 2.099999 改为 2.1 ， 5.900001 改为 5.9 ，用列主元高斯消去法求解变 换后的程组，输出解向量 x 及 detA ，并与（ 1）中的结果比较。 d 、用 MATLAB 的部函数 inv 求出系数矩阵的逆矩阵，再输入命令 x=inv(A)*b, 即可求出上 述各个程组的解，并与列主元高斯消去法和 LU 分解法求出的解进行比较，体会选主元的法 具有良好的数值稳定性。用