matlab解线性方程组 线性方程组及MATLAB应用.docVIP

matlab解线性方程组 线性方程组及MATLAB应用 数值实验 线性方程组与MATLAB应用 王 1．实验目的：理解矩阵的范数与条件数。 实验内容：已知矩阵 11?1? ??11?1 A?? ?1?11? ?1?1?1?1??1? 求A1，A2，A?和cond2(A)。 ?1?1?? 解：编写了一个M文件来求矩阵A的范数与条件数： test3_1.m如下： A=[1 1 1 1;-1 1 -1 1;-1 -1 1 1;1 -1 -1 1]; norm(A,1) norm(A,2) norm(A,inf) cond(A,2) 计算结果依次是： 4 2 4 1.0000 2．实验目的：研究高斯消去法的数值稳定性（出现小主元）。 实验内容：设方程组Ax?b，其中两个矩阵如下，分别对以上两个方程组 ?0.3?10?15? 5.291 （1）A1?? ?11.2? 1??1??59.17? ??46.78? ?6.130?12?? ，b1?? ?1?952? ??? 2211????59.14 3 ?7?10 ??32.09999999999999 （2）A2?? ?5?1? 1?001?8?? ?5.9000000000?000162?? ?，b2????55?1? ??? 102??? （1）计算矩阵的条件数，判断系数矩阵是良态的还是病态的？ 解： 本题编写了一个test3_21的M文件如下： A1=[0.3*1e-15 59.14 3 1;5.291 -6.130 -1 2;11.2 9 5 2;1 2 1 1]; A2=[10 -7 0 1;-3 2.099999999999999 6 2;5 -1 5 -1;0 1 0 2]; cond(A1) cond(A2) 求得两个矩阵的条件数分别为68.4296和8.9939，易知这矩阵A1的条件数远远大于1，而矩阵A2的条件数刚大于1，故这,矩阵A1为病态矩阵，矩阵A2为良态矩阵。 （2）用列主元消去法求得L和U及解向量x1,x2?R； 解：本题利用Matlab的列主元三角分解函数lu()；具体求解如下： A1=[0.3*1e-15 59.14 3 1;5.291 -6.130 -1 2;11.2 9 5 2;1 2 1 1]; A2=[10 -7 0 1;-3 2.099999999999999 6 2;5 -1 5 -1;0 1 0 2]; 4 b1=[59.17;46.78;1;2]; b2=[8;5.0000000000001;5;1]; [L1,U1]=lu(A1) L1 = 0.0000 1.0000 0 0 0.4724 -0.1755 1.0000 0 1.0000 0 0 0 0.0893 0.0202 -0.1738 1.0000 U1 = 11.2000 9.0000 5.0000 2.0000 0 59.1400 3.0000 1.0000 0 0 -2.8354 1.2307 0 0 0 1.0151 [L2,U2]=lu(A2) L2 =1.0000 0 0 0 -0.3000 -0.0000 1.0000 0 0.5000 1.0000 0 0 0 0.4000 -0.3333 1.0000 U2 =10.0000 -7.0000 0 1.0000 0 2.5000 5.0000 -1.5000 0 0 6.0000 2.3000 0 0 0 3.3667 y1=L1\b1; x1=U1\y1 x1 =3.8457 1.6095 -15.4761 10.4113 y2=L2\b2; x2=U2\y2 x2 =0.1337 -0.8218 0.8842 0.9109 ?1,x?2?R4x用不选主元的高斯消去法求得L和U及解向量； 解：编写一个LU_Fact的M文件储存不选主元的LU分解法然后调用求解。具