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

数值试验报告分析 一、实验名称：解线性方程组的列主元素高斯消去法和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； 四、实验内容： 解下列两个线性方程组 = 1 \* GB2 （1） = 2 \* GB2 （2） 用你熟悉的算法语言编写程序用列主元高斯消去法和LU分解求解上述两个方程组，输出Ax=b中矩阵A及向量b, A=LU分解的L及U，detA及解向量x. 将方程组 = 1 \* GB2 （1）中系数3.01改为3.00，0.987改为0.990，用列主元高斯消去法求解变换后的方程组，输出列主元行交换次序，解向量x及detA，并与 = 1 \* GB2 （1）中结果比较。 将方程组 = 2 \* GB2 （2）中的2.099999改为2.1，5.900001改为5.9，用列主元高斯消去法求解变换后的方程组，输出解向量x及detA，并与 = 1 \* GB2 （1）中的结果比较。 d、用MATLAB的内部函数inv求出系数矩阵的逆矩阵，再输入命令x=inv(A)*b,即可求出上述各个方程组的解，并与列主元高斯消去法和LU分解法求出的解进行比较，体会选主元的方法具有良好的数值稳定性。用MATLAB的内部函数det求出系数行列式的值，并与（1）、（2）、（3）中输出的系数行列式的值进行比较。 实验结果： a、用你熟悉的算法语言编写程序用列主元高斯消去法和LU分解求解上述两个方程组，输出Ax=b中矩阵A及向量b, A=LU分解的L及U，detA及解向量x. Gauss(高斯）列主元法程序及结果 （１） .A=[3.01 6.03 1.99;1.27 4.16 -1.23;0.987 -4.81 9.34];b=[1;1;1];[RA,RB,n,X]=liezhuY(A,b) D = -0.0305 RA =3 RB =3 n = 3 X = 1.5926 -0.6319 -0.4936 （２）． A=[10 -7 0 1;-3 2.099999 6 2;5 -1 5 -1;2 1 0 2];b=[8;5.900001;5;1];[RA,RB,n,X]=liezhuY(A,b) D =-762.0001 RA = 4 RB =4 n =4 X = 0.0000 -1.0000 1.0000 1.0000 ＬＵ分解法程序及结果 （１） A=[3.01 6.03 1.99;1.27 4.16 -1.23;0.987 -4.81 9.34];h1=zhjLU(A) D = -0.0305 RA =3 U = 3.0100 6.0300 1.9900 0 4.1600 -2.0696 0 0 5.3107 L = 1.0000 0 0 0.4219 1.0000 0 0.3279 -1.6316 1.0000 h1 = 3.0100 4.8635 -0.0305 （２） A=[10 -7 0 1;-3 2.099999 6 2;5 -1 5 -1;2 1 0 2];h1=zhjLU(A) D =