2010级电子科学与技术专业学生MATLAB控制与系统仿真期末实验考试试题

黄淮学院电子科学与工程系2012-2013学年第二学期期末考试MATLAB控制与系统仿真实验试题适用对象：2010级电技本科1、先建立自己的工作目录，再将自己的工作目录设置到MATLAB的搜索路径下，再试验用help命令能否查询到自己的工作目录。 helpHELP topics桌面陈强 - (No table of contents file)2、已知， 求下列表达式的值：1）和（其中为单位矩阵）。2）、和。3）及4)A,B和A（1，3,:);B2。A=-1,5,-4;0,7,8;3,61,7 A = -1 5 -4 0 7 8 3 61 7 B=8,3,-1;2,5,3;-3,2,0 B = 8 3 -1 2 5 3 -3 2 0（1） A+6*B ans = 47 23 -10 12 37 26 -15 73 7 I=1,0,0;0,1,0;0,0,1 I = 1 0 0 0 1 0 0 0 1 A2-B+Ians = -18 -217 17 22 533 109 21 867 526（2） A*Bans = 14 14 16 -10 51 21 125 328 180 A.*Bans = -8 15 4 0 35 24 -9 122 0 B*Aans = -11 0 -15 7 228 53 3 -1 28（3） A/Bans = 1.2234 -0.9255 2.9787 -0.9468 2.3511 -0.9574 4.6170 3.8723 13.8936 B/Aans = -5.3173 -3.5661 0.8942 -1.2768 -0.4743 0.2411 2.1229 1.3173 -0.2924（4） A,Bans = -1 5 -4 8 3 -1 0 7 8 2 5 3 3 61 7 -3 2 0 A(1,3,:);B2ans = -1 5 -4 3 61 7 73 37 1 17 37 13 -20 1 93、已知 完成下列操作：1）输出A在10.25范围内的全部元素。2）取出A前3行构成矩阵B，前两列构成矩阵C，右下角32子矩阵构成矩阵D，B与C的乘积构成矩阵E。3）分别求表达式E=10)&(A A(F)ans = 23 10B=A(1:3,:)B = 23.0000 10.0000 -0.7780 0 41.0000 -45.0000 65.0000 5.0000 32.0000 5.0000 0 32.0000 C=A(:,1:2)C = 23.0000 10.0000 41.0000 -45.0000 32.0000 5.0000 6.0000 -9.5400 D=A(2:4,3:4)D = 65.0000 5.0000 0 32.0000 54.0000 3.1400 E=B*CE = 1.0e+003 * 0.9141 -0.2239 1.2080 2.7123 1.1330 -0.2103 E E&Dans = 1 1 0 1 1 1 E|Dans = 1 1 1 1 1 1 E|Dans = 0 0 1 0 0 04、产生5阶希尔伯特矩阵H和5阶帕斯卡矩阵P，且求其行列式的值Hh和Hp以及他们的条件数Th和Tp，判断哪个矩阵性能更好，为什么？format rat h=hilb(5)h = 1 1/2 1/3 1/4 1/5 1/2 1/3 1/4 1/5 1/6 1/3 1/4 1/5 1/6 1/7 1/4 1/5 1/6 1/7 1/8 1/5 1/6 1/7 1/8 1/9 p=pascal(5)p = 1 1 1 1 1 1 2 3 4 5 1 3 6 10 15 1 4 10 20 35 1 5 15 35 70 Hh=det(h)Hh = 1/0 Hp=det(p)Hp = 1 th=cond(h)th = tp=cond(p)tp = /21 5阶帕斯卡矩阵p比5阶希尔伯特矩阵H矩阵性能更好。因为thtp,所以值更小的更稳定。5、求二重定积分 syms x ya=int(int(exp(-(x2+y2),x,0,1),y,0,1);b=simple(a);c=vpa(b,4) c =0.55786、求函数的符号导数1），求。（1）syms x;y=(sin(1/x);dx1=diff(y),dx2=diff(y,2) dx1 = -cos(1/x)/x2 dx2 = (2*cos(1/x)/x3 - sin(1/x)/x4（2）syms x;y=(1-cos(2*x)/x;dx1=diff(y),dx2=diff(y,2) dx1 = (2*sin(2*x)/x + (cos(2*x) - 1)/x2dx2 = (4*cos(2*x)/x - (4*sin(2*x)/x2 - (2*(cos(2*x) - 1)/x37、设 在同一图形窗口采用子图的形式绘制条形图、阶梯图、杆图和对数坐标图等不同图形。t=-pi:pi/10:pi;y=1./(1+exp(-t);subplot(2,2,1);bar(t,y,g);title(bar(t,y,g);axis(-3.5,3.5,-1,1；subplot(2,2,2);stairs(t,y,b);title(stairs(t,y,b);axis(-3.5,3.5,-1,1);subplot(2,2,3);stem(t,y,k);title(stem(t,y,k);axis(-3.5,3.5,-1,1);subplot(2,2,4);semilogx(t,y);title(loglog(t,y);grid on;8、计算分段函数 。x=input(请输入一个数);if x=0 f(x)=(x+sqrt(pi)/exp(2);else f(x)=(1/2)*log(x+sqrt(1+x2); endf(x) 请输入一个数8ans = 1.38829、系统闭环传递函数为，试绘制系统在单位脉冲、单位阶跃、单位斜坡及正弦函数作用下的响应曲线。10、