最新数值分析大作业35843584.docVIP

数 值 分 析 大 作 业 姓名：黄 晨 晨 学号： 学院：储运与建筑工程学院学院 班级：储建研17-2 实验3.1 Gauss消去法的数值稳定性实验 实验目的： 理解高斯消元过程中出现小主元即很小时引起方程组解数值不定性 实验内容： 求解方程组Ax=b，其中 （1）A1=0.3×10-1559.14315.291-6.130-1211b1=59.1746.7812 ; （2）A2=10-701-32.099999999999625-15-10102，b2=85.90000000000151 ; 实验要求： (1)计算矩阵的条件数，判断系数矩阵是良态的还是病态的 (2)用Gauss列主元消去法求得L和U及解向量x1,x2∈R4 (3)用不选主元的高斯消去法求得L和U及解向量x1,x2∈R4 (4)观察小主元并分析对计算结果的影响 （1）计算矩阵的条件数，判断系数矩阵是良态的还是病态的 代码： format longe format compact A1=[0.3*10^-15,59.14,3,1;5.291,-6.130,-1,2;11.2,9,5,2;1,2,1,1] b1=[59.17;46.78;1;2] n=4 C1=cond(A1,1) %C1为A1矩阵1范数下的条件数 C2=cond(A1,2) %C2为A1矩阵2范数下的条件数 C3=cond(A1,inf) %C3为1矩阵谱范数下的条件数 结果： C1 = 1.362944708720448e+02 C2 = 6.842955771253409e+01 C3 = 8.431146020518001e+01 显然A1矩阵为病态矩阵 将矩阵A2，b2输入上述代码中求得A2矩阵的条件数为： C1 = 1.928316831682894e+01 C2 = 8.993938090170119e+00 C3 = 1.835643564356072e+01 显然A2矩阵也为病态矩阵 (2)用Gauss列主元消去法求得L和U及解向量x1,x2∈R4 代码： for k=1:n-1 a=max(abs(A1(k:n,k))) if a==0 return end for i=k:n if abs(A1(i,k))==a y=A1(i,:) A1(i,:)=A1(k,:) A1(k,:)=y x=b1(i,:) b1(i,:)=b1(k,:) b1(k,:)=x break end end if A1(k,k)~=0 A1(k+1:n,k)=A1(k+1:n,k)/A1(k,k) A1(k+1:n,k+1:n)=A1(k+1:n,k+1:n)-A1(k+1:n,k)*A1(k,k+1:n) else break end end L=tril(A1,0); for i=1:n L(i,i)=1; end L U=triu(A1,0) y1=L\b1 x1=U\y1 得到如下结果： x1 = 3.845714853511634e+00 1.609517394778522e+00 -1.547605454206655e+01 1.041130489899787e+01 将A2=[10,-7,0,1;-3,2.099999999999,6,2;5,-1,5,-1;0,1,0,2] b2=[8;5.900000000001;5;1] 代入上述代码求得结果如下： X2 = 4.440892098500626e-16 -9.999999999999993e-01 9.999999999999997e-01 1.000000000000000e+00 (3)用不选主元的高斯消去法求得L和U及解向量x1,x2∈R4 代码： format longe format compact A1=[0.3*10^-15,59.14,3,1;5.291,-6.130,-1,2;11.2,9,5,2;1,2,1,1] b1=[59.17;46.78;1;2] [L,U]=lu(A1) y1=L\b1 x1=