数值分析实验作业,gauss消去法的数值稳定性分析.docxVIP

实验3.1 Gauss消去法的数值稳定性试验 实验目的： 观察和理解Gauss消元过程中出现小主元 观察和理解Gauss消元过程中出现小主元（即 (k)| akk | 很小）时引起的方程 组解的数值不稳定性。 实验内容: 求解方程组Ax =b，其中 -0.3 0-155.29159.14-6.1303-1A1 =11.2 - 0.3 0-15 5.291 59.14 -6.130 3 -1 A1 = 11.2 9 5 1 2 1 _10 -7 (1) 2 2 1 0 11 b1 = 59.17 46.78 -3 2.099999999999 A2 T 5 -1 -1 b2 = 5.900000000001 (2) 实验要求: （1） 计算矩阵的条件数，判断系数矩阵是良态的还是病态的。 _ 4 （2） 用Gauss列主元消去法求得L和U及解向量人兀，R . （3） 用不选主元的Gauss消去法求得L和U及解向量x1,x2 R （4） 观察小主元并分析其对计算结果的影响. 程序如下:计算矩阵条件数及Gauss列主元消去法： format longeng A1=[0.3e-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; k2=cond(A1) %k2为矩阵的条件数； for k=1: n-1 a=max(abs(A1(k: n,k))); [p,k]=fi nd(A 1==a); B=A1(k,:);c=b1(k); A1(k, :)=A1(p,:);b1(k)=b1(p); A1(p,:)=B;b1(p)=c; 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 L1=tril(A1,0); for i=1: n L1(i,i)=1; end L=L1 U=triu(A1,0) for j=1: n-1 b1(j)=b1(j)/L(j,j); b1(j+1: n)=b1(j+1: n)-b1(j)*L(j+1: n,j); end b1( n)=b1( n)/L(n,n); for j=n:-1:2 b1(j)=b1(j)/U(j,j); b1(1:j-1)=b1(1:j-1)-b1(j)*U(1:j-1,j); end b1(1)=b1(1)/U(1,1); x1=b1 运行结果如下： K2=68.43; - 1 0 0 01 26.79X0 48 1 0 0 L = 0.4724 -0.1755 1 0 _ 0.0893 -0.0202 -0.4929 1 一 11.2 9 5 2 1 0 59.14 3 1 U = 0 0 -2.835 1.231 _ 0 0 0 0.801 一 X1 =[18.9882;3.3378;-34.747;-33.9865] 不选主元的Gauss消去法程序： clear format I ongeng A1=[0.3e-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; for k=1: n-1 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); end L1=tril(A1,0); for i=1: n L1(i,i)=1； end L=L1 U=triu(A1,0) for j=1: n-1 b1(j)=b1(j)/L(j,j); b1(j+1: n)=b1(j+1: n)-b1(j)*L(j+1: n,j); end b1( n)=b1( n)/L(n,n); for j=n:-1:2 b1(j)=b1(j)/U(j,j); b1(1:j-1)=b1(1:j-1)-b1(j)*U(1:j-1,j); end b1(1)=b1(1)/U(1,1); x1=b1 程序运行结果如下： ■ 1 0 0 0] 15 TOC \o 1-5 \h \z ?17.637X0 1 0 0 L = 15 37.33F0 2.1168 1 0 15 -3.333勺0 0.189 0 1_ 0.^10_15059.1418-1.043 0.^10_15 0 59.14 18 -1.043 10 0 0 3 15 -52.91 10 -16 0 1