数值分析作业课题.doc

数 值 分 析 作 业 姓名：王森 学号： 学院：化学工程学院 班级：化工8班 实验3.1 Gauss消去法的数值稳定性实验 实验目的： 理解高斯消元过程中出现小主元即很小时引起方程组解数值不稳定性 实验题： 求解线性方程组 实验要求 1 计算矩阵的条件数判断系数矩阵是良态的还是病态的 2 用高斯列主元消去法求得L和U及解向量 3 用不选主元的高斯消去法求得L和U及解向量 4 观察小主元并分析对计算结果的影响 解1 1.1判断矩阵A1 是否病态 程序如下: A1=[0.3*10^-15,59.14,3,1;5.291,-6.130,-1,2;11.2,9,5,2;1,2,1,1]; cond(A1 =,1) 结果如下: ans = 136.2945 A1=[0.3*10^-15,59.14,3,1;5.291,-6.130,-1,2;11.2,9,5,2;1,2,1,1]; cond(A1 ,inf) 结果如下 ans = 84.3115 A1=[0.3*10^-15,59.14,3,1;5.291,-6.130,-1,2;11.2,9,5,2;1,2,1,1]; cond(A 1,2) 结果如下 ans = 68.4296 因为cond(A1 ,1) =136.2945》1；cond(A1 ,2) = 68.4296》1；cond(A1 ,inf) =84.3115》1所以该矩阵A1 是病态矩阵。 1.2判断矩阵A2是否病态 程序如下: A2=[10,-7,0,1;-3,2.099999999999,6,2;5,-1,5,-1;0,1,0,2]; cond(A2,1) ans = 19.2832 A2 =[10,-7,0,1;-3,2.099999999999,6,2;5,-1,5,-1;0,1,0,2]; cond(A2,2) ans = 8.9939 A2=[10,-7,0,1;-3,2.099999999999,6,2;5,-1,5,-1;0,1,0,2]; cond(A2,inf) ans = 18.3564 因为cond(A2,1) =19.2832 》1；cond(A2,2) = 8.9939》1；cond(A2,inf) = 18.3564 》1所以该矩阵A2是病态矩阵。 2高斯列主元消去法程序如下： function x=gauss(A,b) %高斯求解方程组 %x=gauss(A,b) n=length(A); a=[A,b]; for k=1:n-1 maxa=max(abs(a(k:n,k))); if maxa==0 return; end for i=k:n if abs(a(i,k))==maxa y=a(i,k:n+1);a(i,k:n+1)=a(k,k:n+1);a(k,k:n+1)=y; break; end end for i=k+1:n l(i,k)=a(i,k)/a(k,k); a(i,k+1:n+1)=a(i,k+1:n+1)-l(i,k).*a(k,k+1:n+1); end end %回代 if a(n,n)==0 return end x(n)=a(n,n+1)/a(n,n); for i=n-1:-1:1 x(i)=(a(i,n+1)-sum(a(i,i+1:n).*x(i+1:n)))/a(i,i); end 执行下列程序： A=[0.3*10^(-15),59.14,3,1;5.291,-6.130,-1,2;11.2,9,5,2;1,2,1,1];b=[59.17;46.78;1;2];x=gauss(A,b) 得到结果如下： x = 3.8457 1.6095 -15.4761 10.4113 A=[10,-7,0,1;-3,2.099999999999,6,2;5,-1,5,-1;0,1,0,2];b=[8;5.900000000001;5;1];x=gauss(A,b) 得到结果如下： x = 0.0000 -1.0000 1.0000 1.0000 2.1 选列主元的分解如下： 2.1.1对A1进行列主元分解程序如下： A1=[0.3*10^-15,59.14,3,1;5.291,-6.130,-1,2;11.2,9,5,2;1,2,1,1]; [L,U,P]=lu(A1) 结果如下： L = 1.0000 0 0 0 0.0000 1.0000 0