内容发布更新时间 : 2024/11/19 23:24:30星期一 下面是文章的全部内容请认真阅读。
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题运行结果示意图