实验4 Z变换和系统频域特性的MATLAB实现

《实验4 Z变换和系统频域特性的MATLAB实现》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、信息工程学院实验报告成绩：指导老师(签名)：课程名称：数字信号处理实验项目名称：实验4Z变换和系统频域特性的MATLAB实现实验时间：2016.11.11班级：姓名：学号：一、实验目的:学习通过Z变换来分析离散系统的频率响应，并用MATLAB实现。加深对系统的零、极点分布概念的理解。二、实验设备与器件MATLAB2008软件三、实验内容与结果分析3.1求系统的零极点和频率响应。指令语句如下：a=[1,-1.8107,2.494,-1.8801,0.9537,-0.2336,0,0,0];b=[0.0528,0.0797,0.1295,

2、0.1295,0.797,0.0528,0,0,0];subplot(2,2,1);zplane(b,a);[H,w]=freqz(b,a,100);magH=abs(H);phaH=angle(H);subplot(2,2,3);plot(w/pi,magH);gridonxlabel('FrequencyinPiuints');ylabel('Magnitude');title('Magnituderesponse');subplot(2,2,4);plot(w/pi,phaH/pi);xlabel('FrequencyinPiu

3、ints');ylabel('Phase');title('Phaseresponse');[h,n]=impz(b,a,40);subplot(2,2,2);stem(n,h,'.');ylabel('h[n]');gridon执行结果如图3-1所示：第11页共11页图3-13.2求出矩形序列的零极点和频率响应。指令语句如下：a=[1,0,0,0,0,0,0,0];b=[1,1,1,1,1,1,1,1];subplot(2,2,1);zplane(b,a);[H,w]=freqz(b,a,100);magH=abs(H);phaH=

4、angle(H);subplot(2,2,3);plot(w/pi,magH);gridonxlabel('FrequencyinPiuints');ylabel('Magnitude');title('Magnituderesponse');subplot(2,2,4);plot(w/pi,phaH/pi);xlabel('FrequencyinPiuints');ylabel('Phase');title('Phaseresponse');[h,n]=impz(b,a,15);subplot(2,2,2);stem(n,h,'.'

5、);ylabel('h[n]');gridon执行结果如图3-2所示：第11页共11页图3-23.3已知离散系统差分方程为：求其Z变换，画出零极点图，并判断系统的稳定性。指令语句如下：a=[1,0.75,0.125];b=[1,-0.5,0];subplot;zplane(b,a);[H,w]=freqz(b,a,100);magH=abs(H);phaH=angle(H);执行结果如图3-3所示：第11页共11页图3-3极点分布在单位圆内，所以系统是稳定的。3.4教材P74-32题指令语句如下：%H1a=[1,-1.6,0.9425

6、];b=[1,0,0];subplot(2,1,1);zplane(b,a);title('H1零点-极点图');[h,n]=impz(b,a,40);subplot(2,1,2);stem(n,h,'.');title('H1单位脉冲响应');ylabel('h[n]');gridon%%%H2a=[1,-1.6,0.9425];b=[1,-0.3,0];subplot(2,1,1);zplane(b,a);title('H2零点-极点图');[h,n]=impz(b,a,40);subplot(2,1,2);stem(n,h,'.

7、');title('H2单位脉冲响应');ylabel('h[n]');gridon%%%H3a=[1,-1.6,0.9425];第11页共11页b=[1,-0.8,0];subplot(2,1,1);zplane(b,a);title('H3零点-极点图');[h,n]=impz(b,a,40);subplot(2,1,2);stem(n,h,'.');title('H3单位脉冲响应');ylabel('h[n]');gridon%%%H4b=[1,-1.6,0.8];a=[1,-1.6,0.9425];subplot(2,1,1)

8、;zplane(b,a);title('H4零点-极点图');[h,n]=impz(b,a,40);subplot(2,1,2);stem(n,h,'.');title('H4单位脉冲响应');ylabel('h[n]'