数值分析实验四.docVIP

数值分析实验一 一、实验名称 解线性方程组的列主元素高斯消去法和LU分解法 二、实验目的 通过数值实验，从中体会解线性方程组选主元的必要性和LU分解法的优点，以及方程组系数矩阵和右端向量的微小变化对解向量的影响。 三、实验内容 解下列两个线性方程组 （1） （2） 四、实验要求 （1）用你熟悉的算法语言编写程序用列主元高斯消去法和LU分解求上述两个方程组，输出 中矩阵A及向量b，A=LU分解的L及U，detA及解向量x （2）将方程组（1）中系数3.01改为3.00，0.987改为0.990，用列主元高斯消去法求解变换后的方程组，输出列主元行交换次序，解向量X及detA，并与（1）中结果比较。 （3）将方程组（2）中的2.099999改为2.1，5.900001改为5.9，用列主元高斯消去法求解变换后的方程组，输出解向量x及detA，并与（1）中的结果比较。 （4）用MATLAB的内部函数inv求出系数矩阵的逆矩阵，再输入命令x=inv（A）*b，即可求出上述各个方程组的解，并与列主元高斯消去法和LU分解法求出的解进行比较，体会选主元的方法具有良好的数值稳定性。用MATLAB的内部函数inv求出系数行列式的值，并与（1）、（2）、（3）中的输出系数行列式的值进行比较。 五、实验结果 （1）第一个方程组 ： A1 = 3.0100 6.0300 1.9900 1.2700 4.1600 -1.2300 0.9870 -4.8100 9.3400 b1 = 1 1 1 L1 = 1.0000 0 0 0.4219 1.0000 0 0.3279 -4.2006 1.0000 U1 = 3.0100 6.0300 1.9900 0 1.6158 -2.0696 0 0 -0.0063 利用列主元素高斯消去法所求x x1 = 1.0e+003 * 1.5926 -0.6319 -0.4936 det1 = -0.0305 第二个方程组 ： A2 = 10.0000 -7.0000 0 1.0000 -3.0000 2.1000 6.0000 2.0000 5.0000 -1.0000 5.0000 -1.0000 2.0000 1.0000 0 2.0000 b2 = 8.0000 5.9000 5.0000 1.0000 L2 = 1.0e+006 * 0.0000 0 0 0 -0.0000 0.0000 0 0 0.0000 -2.5000 0.0000 0 0.0000 -2.4000 0.0000 0.0000 U2 = 1.0e+007 * 0.0000 -0.0000 0 0.0000 0 -0.0000 0.0000 0.0000 0 0 1.5000 0.5750 0 0 0 0.0000 利用列主元素高斯消去法所求x x2 = 0.0000 -1.0000 1.0000 1.0000 det2 = -762.0001 （2）x3 = 119.5273 -47.1426 -36.8403 det3 = -0.4070 （3）x4 = 0.0000 -1.0000 1.0000 1.0000 det4 = -762 （4）x11 = 1.0e+003 * 1.5926 -0.6319 -0.4936 det11 = -0.0305 x22 = 0 -1.0000 1.0000 1.0000 det22 = -7