数值分析大作业

1、数 值 分 析 大 作 业姓名：黄 晨 晨学号：S 学院：储运与建筑工程学院学院班级：储建研17-2实验3.1 Gauss消去法的数值稳定性实验实验目的：理解高斯消元过程中出现小主元即很小时引起方程组解数值不定性实验内容：求解方程组Ax=b，其中（1）A1=0.310-1559.14315.291-6.130-1211.，b1=59.1746.7812 ;（2）A2=10-701-32.9625-15-10102，b2=85.151 ;实验要求： (1)计算矩阵的条件数，判断系数矩阵是良态的还是病态的 (2)用Gauss列主元消去法求得L和U及解向量x1,x2R4 (3)用不选主元的高斯消去法

2、求得L和U及解向量x1,x2R4 (4)观察小主元并分析对计算结果的影响 （1）计算矩阵的条件数，判断系数矩阵是良态的还是病态的代码：format longeformat compactA1=0.3*10-15,59.14,3,1;5.291,-6.130,-1,2;11.2,9,5,2;1,2,1,1b1=59.17;46.78;1;2n=4C1=cond(A1,1) %C1为A1矩阵1范数下的条件数C2=cond(A1,2) %C2为A1矩阵2范数下的条件数C3=cond(A1,inf) %C3为1矩阵谱范数下的条件数结果：C1 = 1.0448e+02C2 = 6.3409e+01C3

3、= 8.8001e+01显然A1矩阵为病态矩阵将矩阵A2，b2输入上述代码中求得A2矩阵的条件数为：C1 = 1.2894e+01C2 = 8.0119e+00C3 = 1.6072e+01显然A2矩阵也为病态矩阵(2)用Gauss列主元消去法求得L和U及解向量x1,x2R4代码：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

4、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 endendL=tril(A1,0);for i=1:n L(i,i)=1; endLU=triu(A1,0) y1=Lb1x1=Uy1得到如下结果：x1 = 3.1634e+00 1.8522e+00 -1.6655e+01 1.9787e+01将A2=10,-7,0,1;-3,2.9,6,2;5,-1,5,-1;0,1,0,2b2=8;5.1;5;1代入上述代码

5、求得结果如下：X2 = 4.0626e-16 -9.9993e-01 9.9997e-01 1.0000e+00(3)用不选主元的高斯消去法求得L和U及解向量x1,x2R4代码：format longeformat compactA1=0.3*10-15,59.14,3,1;5.291,-6.130,-1,2;11.2,9,5,2;1,2,1,1b1=59.17;46.78;1;2 L,U=lu(A1)y1=Lb1x1=Uy1求得如下结果：x1= 3.1634e+00 1.8522e+00 -1.6655e+01 1.9787e+01将A2=10,-7,0,1;-3,2.9,6,2;5,-1,

6、5,-1;0,1,0,2b2=8;5.1;5;1代入上述代码，求得结果如下：x2= 4.0626e-16 -9.9993e-01 9.9997e-01 9.9999e-01（2）（3）求得结果相同，可验证结果正确。(4)观察小主元并分析对计算结果的影响实验结果表明，小主元的存在对计算结果影响较大。小主元的存在会使得计算结果有较大的误差。实验3.2 方程组的性态和条件数实验实验目的：理解条件数的意义和方程组的性态对解向量的影响。实验内容：已知两个方程组A1x=b和A2x=b，其中：A1=1x0x02x0n1x1x12x1n1x2x22x2n1xnxn2xnn , A2=11/21/31/n1/2

7、1/31/41/(n+1)1/31/41/51/(n+2)1/n1/(n+1)1/(n+2)1/(2n-1)b=j=1na1j,j=1na2j,j=1na1nj,T实验要求：对A1，取xk=1+0.1k,k=0,1,，n，下面均用Matlab函数“x=A/b”计算方程组的解。（1）取n=4,6,8,分别求出A1,A2的条件数，判别它们是否是病态阵？随n的增大，矩阵性态的变化如何？代码：format shorteformat compactA1=1,1,1,1,1;1,1.1,1.12,1.13,1.14;1,1.2,1.22,1.23,1.24;1,1.3,1.32,1.33,1.34;1,1

8、.4,1.42,1.43,1.44C1=cond(A1,1) %C1为A1矩阵1范数下的条件数C2=cond(A1,2) %C2为A1矩阵2范数下的条件数C3=cond(A1,inf) %C3为1矩阵谱范数下的条件数1.对A1矩阵运算n=4时：C1 = 6.5120e+05C2 = 3.5740e+05C3 = 6.2755e+05n=4时，A1为病态矩阵n=6时：A1=1,1,1,1,1,1,1;1,1.1,1.12,1.13,1.14,1.15,1.16;1,1.2,1.22,1.23,1.24,1.25,1.26;1,1.3,1.32,1.33,1.34,1.35,1.36;1,1.4,

9、1.42,1.43,1.44,1.45,1.46;1,1.5,1.52,1.53,1.54,1.55,1.56;1,1.6,1.62,1.63,1.64,1.65,1.66C1 = 1.8531e+08C2 = 8.7385e+07C3 = 1.6574e+08n=6时，A1为病态矩阵n=8时：A1=1,1,1,1,1,1,1,1,1;1,1.1,1.12,1.13,1.14,1.15,1.16,1.17,1.18;1,1.2,1.22,1.23,1.24,1.25,1.26,1.27,1.28;1,1.3,1.32,1.33,1.34,1.35,1.36,1.37,1.38;1,1.4,1.

10、42,1.43,1.44,1.45,1.46,1.47,1.48;1,1.5,1.52,1.53,1.54,1.55,1.56,1.57,1.58;1,1.6,1.62,1.63,1.64,1.65,1.66,1.67,1.68;1,1.7,1.72,1.73,1.74,1.75,1.76,1.77,1.78;1,1.8,1.82,1.83,1.84,1.85,1.86,1.87,1.88C1 = 5.0565e+10C2 = 2.2739e+10C3 = 4.4788e+10n=8时，A1为病态矩阵2.对A2矩阵运算n=4时：A2=1,1/2,1/3,1/4;1/2,1/3,1/4,1/5;

11、1/3,1/4,1/5,1/6;1/4,1/5,1/6,1/7C1 = 2.8375e+04C2 = 1.5514e+04C3 = 2.8375e+04n=4时，A2为病态矩阵n=6时：A2=1,1/2,1/3,1/4,1/5,1/6;1/2,1/3,1/4,1/5,1/6,1/7;1/3,1/4,1/5,1/6,1/7,1/8;1/4,1/5,1/6,1/7,1/8,1/9;1/5,1/6,1/7,1/8,1/9,1/10;1/6,1/7,1/8,1/9,1/10,1/11C1 = 2.9070e+07C2 = 1.4951e+07C3 = 2.9070e+07n=6时，A2为病态矩阵n=8

12、时：A2=1,1/2,1/3,1/4,1/5,1/6,1/7,1/8;1/2,1/3,1/4,1/5,1/6,1/7,1/8,1/9;1/3,1/4,1/5,1/6,1/7,1/8,1/9,1/10;1/4,1/5,1/6,1/7,1/8,1/9,1/10,1/11;1/5,1/6,1/7,1/8,1/9,1/10,1/11,1/12;1/6,1/7,1/8,1/9,1/10,1/11,1/12,1/13;1/7,1/8,1/9,1/10,1/11,1/12,1/13,1/14;1/8,1/9,1/10,1/11,1/12,1/13,1/14,1/15C1 = 3.3873e+10C2 = 1

13、.5258e+10C3 = 3.3873e+10n=8时，A2为病态矩阵观察上述计算可得：随着n的增大，A1、A2矩阵的条件数迅速增大，矩阵病态越来越严重。（2）取n=5，分别 两个方程组的解向量x1,x2.代码：format shorteformat compactA1=1,1,1,1,1,1; 1,1.1,1.12,1.13,1.14,1.15; 1,1.2,1.22,1.23,1.24,1.25; 1,1.3,1.32,1.33,1.34,1.35; 1,1.4,1.42,1.43,1.44,1.45; 1,1.5,1.52,1.53,1.54,1.55;b=1+1+1+1+1+1; 1

14、+1.1+1.12+1.13+1.14+1.15;1+1.2+1.22+1.23+1.24+1.25; 1+1.3+1.32+1.33+1.34+1.35; 1+1.4+1.42+1.43+1.44+1.45; 1+1.5+1.52+1.53+1.54+1.55;X=A1b结果如下：X = 1.0000e+00 1.0000e+00 1.0000e+00 1.0000e+00 1.0000e+00 1.0000e+00A2=1,1/2,1/3,1/4,1/5;1/2,1/3,1/4,1/5,1/6;1/3,1/4,1/5,1/6,1/7;1/4,1/5,1/6,1/7,1/8;1/5,1/6,

15、1/7,1/8,1/9;b=1+1/2+1/3+1/4+1/5;1/2+1/3+1/4+1/5+1/6;1/3+1/4+1/5+1/6+1/7;1/4+1/5+1/6+1/7+1/8;1/5+1/6+1/7+1/8+1/9;X=A2b求得结果如下：X = 1.0000e+00 1.0000e+00 1.0000e+00 1.0000e+00 1.0000e+00（3）取n=5,b不变，对A1的元素a22和a66分别加一个扰动10-14，分别求出第一个方程组的解向量x;若A1不变，对b的元素b6加一个扰动10-4，求出x.1.对a22施加扰动后，矩阵A1，b为：A1=1,1,1,1,1,1; 1

16、,1.1+10-14,1.12,1.13,1.14,1.15; 1,1.2,1.22,1.23,1.24,1.25; 1,1.3,1.32,1.33,1.34,1.35; 1,1.4,1.42,1.43,1.44,1.45; 1,1.5,1.52,1.53,1.54,1.55;b=1+1+1+1+1+1; 1+1.1+1.12+1.13+1.14+1.15;1+1.2+1.22+1.23+1.24+1.25; 1+1.3+1.32+1.33+1.34+1.35; 1+1.4+1.42+1.43+1.44+1.45; 1+1.5+1.52+1.53+1.54+1.55; x=A1b将数值显示改为

17、longe，结果如下：x = 1.5445e+00 9.8064e-01 1.7554e+00 9.7031e-01 1.3978e+00 9.5134e-012.对a66施加扰动后，矩阵A1，b为：A1=1,1,1,1,1,1; 1,1.1,1.12,1.13,1.14,1.15; 1,1.2,1.22,1.23,1.24,1.25; 1,1.3,1.32,1.33,1.34,1.35; 1,1.4,1.42,1.43,1.44,1.45; 1,1.5,1.52,1.53,1.54,1.55+10-14;b=1+1+1+1+1+1; 1+1.1+1.12+1.13+1.14+1.15;1+1

18、.2+1.22+1.23+1.24+1.25; 1+1.3+1.32+1.33+1.34+1.35; 1+1.4+1.42+1.43+1.44+1.45; 1+1.5+1.52+1.53+1.54+1.55; x=A1b求得结果如下：x = 1.6265e+00 9.1665e-01 1.6582e+00 9.6432e-01 1.3006e+00 9.3373e-013. 对矩阵b6施加扰动10-4后，矩阵A1，b为：A1=1,1,1,1,1,1; 1,1.1,1.12,1.13,1.14,1.15; 1,1.2,1.22,1.23,1.24,1.25; 1,1.3,1.32,1.33,1.

19、34,1.35; 1,1.4,1.42,1.43,1.44,1.45; 1,1.5,1.52,1.53,1.54,1.55;b=1+1+1+1+1+1; 1+1.1+1.12+1.13+1.14+1.15;1+1.2+1.22+1.23+1.24+1.25; 1+1.3+1.32+1.33+1.34+1.35; 1+1.4+1.42+1.43+1.44+1.45; 1+1.5+1.52+1.53+1.54+1.55+10-4; x=A1b求得结果为：x = 7.9980e-01 1.8460e+00 -4.2500e-01 2.1958e+00 5.0000e-01 1.0833e+00（4）

20、取n=6，b不变，对A2的元素a22和a66分别加一个扰动10-7，别分求出第二个方程组的解向量x.1.对a22施加扰动10-4后，矩阵A2，b为:A2=1,1/2,1/3,1/4,1/5,1/6;1/2,1/3+10-7,1/4,1/5,1/6,1/7;1/3,1/4,1/5,1/6,1/7,1/8;1/4,1/5,1/6,1/7,1/8,1/9;1/5,1/6,1/7,1/8,1/9,1/10;1/6,1/7,1/8,1/9,1/10,1/11b=1+1/2+1/3+1/4+1/5+1/6;1/2+1/3+1/4+1/5+1/6+1/7;1/3+1/4+1/5+1/6+1/7+1/8;1/

21、4+1/5+1/6+1/7+1/8+1/9;1/5+1/6+1/7+1/8+1/9+1/10;1/6+1/7+1/8+1/9+1/10+1/11x=A2b求得结果如下：x = 1.0001e+00 9.9853e-01 1.0088e+00 9.7886e-01 1.0220e+00 9.9170e-012.对a66施加扰动10-4后，矩阵A2，b为:A2=1,1/2,1/3,1/4,1/5,1/6;1/2,1/3,1/4,1/5,1/6,1/7;1/3,1/4,1/5,1/6,1/7,1/8;1/4,1/5,1/6,1/7,1/8,1/9;1/5,1/6,1/7,1/8,1/9,1/10;1/6,1/7,1/8,1/9,1/10,1/11+10-7b=1+1/2+1/3+1/4+1/5+1/6;1/2+1/3+1/4+1/5+1/6+1/7;1/3+1/4+1/5+1/6+1/7+1/8;1/4+1/5+1/6+1/7+1/8+1/9;1/5+1/6+1/7+1/8+1/9+1/10;1/6+1/7+1/8+1/9+1/10+1/11x=A2b求得结果如下：x = 1.0003e+00 9.9223e-01 1.0544e+00 8.5490e-01 1.1632e+00 9.3471e-01（5）观察和分析A1,A2和b的微小扰动对解向