2013东南大学信息学院短学期MATLAB试卷答案.docx

《2013东南大学信息学院短学期MATLAB试卷答案.docx》由会员分享，可在线阅读，更多相关《2013东南大学信息学院短学期MATLAB试卷答案.docx（12页珍藏版）》请在第一文库网上搜索。

1、MATLAB上机测验题（考试时间：2： 10-4:10）姓名 学号考试要求：1、要求独立完成不得与他人共享，答卷雷同将做不及格处理。2、答卷用Word文件递交，文件名为学号+姓名.doc,试卷写上姓名及学号。3、答卷内容包括：（1）程序；（2）运行结果及其分析；（3）图也要粘贴在文档中。上机考题：一、已知系统传递函数为：-1130-，画出系统的零极点54 + 9s3 + 45/ + 97s + 50图，判断稳定性。给定频率范围为0,10,步长为0.1,画出其频率响应，同时求系统的单位脉冲响应曲线图。%第一问画出极零点判断稳定性程序 clear; b=l, 11,30, 0; a=l,9, 45

2、,97,50; sys=tf（b, a）; pzmap（sys）;极零图如下，无在Y轴右边极点，所以系统稳定。Pole-Zero Map_ _ _ _ i G _ _ _543210123- - s_xv 一6ee-兴兴Y)-4-3-2Real Axis-1%第二问画出频率响应图 w=0:0. 1:10;freqs (b, a, w);-101101010Frequency (rads)%第三间单位脉冲响应曲线图 impulse(sys);Impulse Response9pr-d 专234Time (sec)二、系统函数H(z) =l + 0.8z, 0.5z2,画出系统的幅频响应和相位响应

3、图。当系统输入为 x(h) = 5 cos(0.2) + 2sin(0.7),0 n title。时域波形(f=0.0625)时域波形（印0625） %DTFT 方法;f=0.0625 w=-pi:0.1:pi; for i=l:length(w);for k=l:129;ff(k)=x(k)*exp( 1 j*(k-1 )*w);endf(i)=sum(ff);end subplot(2,l,l); plot(w,abs(f); title(幅频特性(f=0.0625)X subplot(2,l,2); plot(w,angle(f); titleC相频特性(f=0.0625)R幅频特性户口

4、 0625)相频特性户0.0625)%f=0.4375时域波形 n=0:128; a=0.1; f=0.4375; x=exp(-l*a*n).*sin(2*pi*fX stem(n,x); title(时域波形(f=0.4375)% DTFT 方法 f=0.4375 w=-pi:0.1:pi; for i=l:length(w);for k=l:129;ff(k)=x(k)*exp( 1 j*(k-1 )*w)；endf(i)=sum(ff);end subplot(2,l,l); plot(w,abs(f)i subplot(2,l,l); plot(w,abs(f); titleC幅频特性(f=0.4375)R subplot(2,l,2); plot(w,angle(f); title(湘频特性(f=0.43756)幅频特性户。.4375)相频特性产。.43756)由上图可知f=0.0625时的频谱图比f=0.4375时的频谱图波动更大，f=0.0625频谱图比较陡，f=().4375的频谱图比较光滑。分析原因：f=0.0625发生截断效应