matlab实习报告二.docxVIP

MATLAB实习报告(2) 实验二 MATLAB矩阵分析与处理 王夏 实验目的 掌握生成特殊矩阵的方法。 掌握矩阵分析的方法。 用矩阵求逆发解线性方程组。 实验内容 设有分块矩阵A=[E3×3 R3×2 ;O2×3 S2×2],其中E、R、O、S分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证A2=[E R+RS;O S2]。 程序清单：E=eye(3); R=rand(3,2); O=zeros(2,3); S=diag([4,5]); A=[E R ;O S] ; A2=A^2; C=[E R+R*S;O S^2]; length(find(A2==C))==25 运行结果：ans =1 产生5阶希尔伯特矩阵H和5阶帕斯卡矩阵P，且求其行列式的值Hh和Hp以及他们的条件数Th和Tp，判断哪个矩阵的性能更好，为什么？ 程序清单：format rat H=hilb(5) ;format short P=pascal(5); Hh=det(H); Hp=det(P); Th=cond(A) ;Tp=cond(P); 运行结果：Hh =3.7493e-012 Hp =1 Th =5.5228 Tp =8.5175e+003 实验收获：会建立希尔伯特矩阵和帕斯卡矩阵，知道怎么求矩阵行列式的值以及条件数。希尔伯特矩阵的性能更好，条件数越接近1的矩阵性能越好。 建立一个5×5矩阵，求它的行列式值、迹、秩和范数。 程序清单：A=[1:5;1:0.1:1.5;2 5 7 3 9;2:6;3:0.4:4.6] Ha=det(A);Ra=rank(A) ;Ta=trace(A); Na=norm(A); 运行结果：Ha =1.4421e-031 Ra = 3 Ta =18.7000 Na =19.4966 已知向量A，求A的特征值及特征向量，并分析其数学意义。 程序清单：A=[-29 6 18 ;20 5 12;-8 8 5] [V,D]=eig(A) 运行结果：V = 0.7130 0.2803 0.2733 -0.6084 -0.7867 0.8725 0.3487 0.5501 0.4050 D = -25.3169 0 0 0 -10.5182 0 0 0 16.8351 求解下列的线性方程组： [1/2 1/3 1/4;1/3 1/4 1/5;1/4 1/5 1/6]*[x1;x2;x3]= [0.95;0.67;0.52] （1）求方程的解 程序清单：format rat A=[1/2 1/3 1/4; 1/3 1/4 1/5;1/4 1/5 1/6] format short B=[0.95;0.67;0.52]; x=inv(A)*B 运行结果：x =1.2000 0.6000 0.6000 （2）将方程右边向量元素b3改为0.53，在求解，并比较b3的变化和解的相对变化。 程序清单：B=[0.95;0.67;0.53]; x=inv(A)*B 运行结果：x = 3.0000 -6.6000 6.6000 （3）计算系数矩阵的条件数并分析结论。 程序清单：cond(A) 运行结果：ans = 1.3533e+003 建立A矩阵，试比较sqrtm(A)和sqrt(A)，并分析他们的区别。 程序清单：a=rand(5) ；a1=sqrtm(a)； a2=sqrt(a)； 运行结果：a = 0.4565 0.7919 0.9355 0.3529 0.1987 0.0185 0.9218 0.9169 0.8132 0.6038 0.8214 0.7382 0.4103 0.0099 0.2722 0.4447 0.1763 0.8936 0.1389 0.1988 0.6154 0.4057 0.0579 0.2028 0.0153 a1 = 0.5983 0.4165 0.5202 0.1609 -0.0723 -0.3534 0.9522 0.4754 0.6539 0.3489 0.6282 0.3061 0.4941 -0.2370 0.2398 -0.0864 -0.0324 0.8597 0.614