MATLAB程序设计及应用(第二版)课后实验答案.docVIP

MATLAB程序设计及应用(第二版)课后实验答案 MATLAB程序设计及应用(第二版)课后实验答案 PAGE / NUMPAGES MATLAB程序设计及应用(第二版)课后实验答案 Matlab 课后实验题答案 实验一 MATLAB 运算基础 先求下列表达式的值，然后显示 MATLAB工作空间的使用情况并保存全部变量。 2sin85 0 (1) z1 1 e2 (2) z2 1 ln( x 2 1 2i 1 x2 ) ，其中 x 5 2 0.45 (3) z3 e0.3a e 0.3a sin(a 0.3) ln 0.3 a , a 3.0, 2.9, , 2.9, 3.0 2 2 t2 0 t 1 (4) z4 t2 1 1 t 2 ，其中 t =0:: t2 2t 1 2 t 3 解： M文件 : z1=2*sin(85*pi/180)/(1+exp(2)) x=[2 1+2*i; 5]; z2=1/2*log(x+sqrt(1+x^2)) a=::; z3=(exp.*a)-exp.*a))./2.*sin(a++log(+a)./2) t=0::; z4=(t=0t1).*(t.^2)+(t=1t2).*(t.^2-1)+(t=2t3) .*(t. ^2-2*t+1) 已知： 12 34 4 1 3 1 A 34 7 87,B 2 0 3 3 65 7 3 2 7 求下列表达式的值： A+6*B 和 A-B+I （其中 I 为单位矩阵） A*B 和 A.*B A^3 和 A.^3 A/B 及 B\A [A,B] 和[A([1,3],:);B^2] 解： 文件： A=[12 34 -4;34 7 87;3 65 7];B=[1 3 -1;2 0 3;3 -2 7]; A+6.*B A-B+eye(3) A*B A.*B A^3 A.^3 A/B B\A [A,B] [A([1,3],:);B^2] 设有矩阵 A和 B 1 2 3 4 5 3 0 16 6 7 8 9 10 17 6 9 A 11 12 13 14 15,B 0 23 4 16 17 18 19 20 9 7 0 21 22 23 24 25 4 13 11 求它们的乘积 C。 将矩阵 C的右下角 3×2 子矩阵赋给 D。 查看 MATLAB工作空间的使用情况。 解： . 运算结果： E=(reshape(1:1:25,5,5));F=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11]; C= E*F H=C(3:5,2:3) C = 93 150 77 258 335 237 423 520 397 588 705 557 753 890 717 H = 397 557 717 完成下列操作： 求 [100,999] 之间能被 21整除的数的个数。 建立一个字符串向量，删除其中的大写字母。 解： (1) 结果： m=100:999; n=find(mod(m,21)==0); length(n) ans = 43 (2). 建立一个字符串向量 例如： ch=ABC123d4e56Fg9; 则要求结果是： ch=ABC123d4e56Fg9; k=find(ch=Ach=Z); ch(k)=[] ch = 123d4e56g9 实验二 MATLAB矩阵分析与处理 1. 设有分块矩阵 A E3 3 R3 2 ，其中 E、R、 O、S 分别为单位矩阵、 O2 3 S2 2 随机矩阵、零矩阵和对角阵，试通过数值计算验证A2 E R RS 。 O S2 解: M 文件如下； 由ans, 所以 A2 E R RS O S2 2. 产生 5 阶希尔伯特矩阵 H 和 5 阶帕斯卡矩阵 P，且求其行列式的 值 Hh和 Hp以及它们的条件数 Th 和 Tp，判断哪个矩阵性能更好。为什么 解： M文件如下： 因为它们的条件数 ThTp,所以 pascal 矩阵性能更好。 建立一个 5×5 矩阵，求它的行列式值、迹、秩和范数。解： M 文件如下： 已知 29 6 18 20512 8 8 5 求 A 的特征值及特征向量，并分析其数学意义。解： M文件如图： 数学意义： V 的 3 个列向量是 A 的特征向量， D的主对角线上 3 个是 A 的特征值，特别的， V 的 3 个列向量分别是 D 的 3 个特征值的特征 向量。 下面是一个线性方程组： 1 1 1 2 3 4 x1 0.95 1 1 1 x2 0.67 3 4 5 x3 0.52 1 1 1 4 5 6 求方程的解。 将方程右边向量元素 b3 改为再求解，并比较 b3 的变化和解的相对变化。 计算系数矩阵