海南大学MATLAB第1章习题解答 下载本文

内容发布更新时间 : 2024/5/20 13:00:39星期一 下面是文章的全部内容请认真阅读。

9 7 6 12 4 14 15 1 b =

19 5 6 16 8 14 13 11 12 10 9 15 7 17 18 4 c =

35 7 9 29 13 25 23 19 21 17 15 27 11 31 33 5

12.矩阵a为2×3的均匀分布随机矩阵,b为3×2的均匀分布随机矩阵,求a与b矩阵相乘,生成矩阵c。

解:>> a=rand(2,3),b=rand(3,2),c=a*b a =

0.8147 0.1270 0.6324 0.9058 0.9134 0.0975 b =

0.2785 0.9649 0.5469 0.1576 0.9575 0.9706 c =

0.9018 1.4199 0.8452 1.1126

13.已知A为3×3均匀分布随机矩阵,B为3×2均匀分布随机矩阵,C为2×3均匀分布随机矩阵,求Q= c*a^2*b。

解:>> a=rand(3),b=rand(3,2),c=rand(2,3),Q= c*a^2*b a =

0.0462 0.6948 0.0344 0.0971 0.3171 0.4387 0.8235 0.9502 0.3816 b =

0.7655 0.4898 0.7952 0.4456 0.1869 0.6463 c =

0.7094 0.2760 0.6551 0.7547 0.6797 0.1626 Q =

1.4237 1.2051 1.0678 0.9204

14.已知A为4×4均匀分布随机矩阵,B为全是“1”的4元列向量,分别利用求逆矩阵和矩阵除法解方程组,比较计算后的结果。

解:>> a=rand(4) %产生(4×4)均匀分布随机矩阵

a =

0.1190 0.5853 0.5060 0.5472 0.4984 0.2238 0.6991 0.1386 0.9597 0.7513 0.8909 0.1493 0.3404 0.2551 0.9593 0.2575

>> b=ones(4,1) % 产生全为1的4元列向量 b = 1 1 1 1

>> x1=inv(a)*b % 矩阵求逆解方程组 x1 =

2.0257 -1.7991 -0.1154 3.4178

>> x2=a\\b % 矩阵除法解方程组

x2 =

2.0257 -1.7991 -0.1154 3.4178

15.指出下列矩阵函数所实现的具体运算。

A=rand(5), B=rank(A),C= eig (A),D= sqrtm(A),E= det(A) 解:A=rand(5) %生成矩阵A,A为5×5均匀分布随机矩阵 B=rank(A) %求矩阵A的秩

C= eig (A) %求方阵A的特征值和特征向量 D= sqrtm(A) %求矩阵A的平方根 E= det(A) %求方阵A的行列式

16.利用MATLAB逻辑运算算符实现a?b&b运算,假设a=1,b=0。 解:>> a=1,b=0,c=a|b&b a = 1 b = 0 c = 1

5432f(x)?x?4x?10x?16x?17x?12的根。 17.利用MATLAB的roots函数求

解:>> p=[1,4,10,16,17,12]

p =

1 4 10 16 17 12 >> x=roots(p) x =

-1.6506

-1.0000 + 1.4142i -1.0000 - 1.4142i -0.1747 + 1.5469i -0.1747 - 1.5469i

cos(3t)在t区间为[0,4π]内的曲线。 18.画出函数y?e解:解题步骤为:

1)在MATLAB主界面下点击file?new?blank M-file, 2)在新建的M文件里输入如下程序段:

t=0:pi/50:4*pi; y0=exp(-t/6);

y=exp(-t/6).*cos(3*t);

plot(t,y,'-r',t,y0,':b',t,-y0,':b')

3)程序编写完成以后,在主菜单栏中选择Debug?Save File and Run,如果编译过程中出现错误,则在MATLAB的command window中显示错误行;若没有,即可显示运行结果,如图1.16所示。

?t6

图1.16 第1-18题运行结果示意图

?19.画出一个幅度为2,频率为4Hz,相位为3的正弦信号。

解:解题步骤为:

1)创建新的M文件(操作与题18同), 2)在M文件中输入如下程序段: A=2; f=4; phi=pi/3; w0=2*pi*f; t=0:0.01:1;

x=A*sin(w0*t+phi); plot(t,x);

3)编译、显示结果,如图1.17所示。

图1.17第1-19题运行结果示意图