离散LSI系统的频域分析.docxVIP

院系 信息工程学院 院系 信息工程学院 班级 13普本测控 学号 姓名 日期 2016.4.18 nBKnwMKniffliEi实验报告 nBKnwMKniffliEi 课程名称数字信号处理 实验成绩 指导教师 实验3 离散LSI系统的频域分析 一、实验目的： 1、 加深对离散系统变换域分析一一 z变换的理解，掌握使用 MATLAB进行z变换和逆 z变换的常用函数的用法。 2、 了解离散系统的零极点与系统因果性和稳定性的关系，熟悉使用 MATLAB进行离散 系统的零极点分析的常用函数的用法。 3、 加深对离散系统的频率响应特性基本概念的理解，掌握使用 MATLAB进行离散系统 幅频响应和相频响应特性分析的常用方法。 二、实验原理 1、z变换和逆z变换 用ztrans函数求无限长序列的 z变换。该函数只给出 z变换的表达式，而没有给 出收敛域。另外，由于这一函数还不尽完善，有的序列的 z变换还不能求出，逆 z变换也存 在同样的问题。 例7-1求以下各序列的z变换 n j ?on X1(n )=a X2( n)=n X3(n )=n(n _1)/2 X4(n )=e x5( n)=1/[ n(n _1)] 程序清单如下： syms w0 n z a; x1=0;X 1= ztra ns(x1) x2=si n(w0* n);X2=ztra ns(x2) x3=exp(-a* n)*si n(w0* n);X3=ztra ns(x3) 程序运行结果如下： X1 =z/a/(z/a-1) X2 =z/(z-1)A2 X3 =1/2*z*(z+1)/(z-1)A3-1/2*z/(z-1)A2 X4 =z/exp(i*w0)/(z/exp(i*w0)-1) X5 =z/(z-1)-ztra ns(1/n ,n,z) 用iztrans函数求无限长序列的逆 z变换。 例3-2求下列函数的逆 z变换。 X(z)z-1/ 、 az x‘z)=(a-z)X 3( X(z) z-1 / 、 az x‘z)=(a-z) X 3(z) (z-1) 3 X 4(z) -n 1-z= 1-z -1 程序清单如下: syms n z a; X1= z/(z-a);x1=iztra ns(X1) X2= z/(z-a)A2;x2=iztra ns(X2) X3=z/[z-exp(j*w0)];x3=iztra ns(X3) X4=(1-zA-3)/(1-zA-1)；x4=iztra ns(X4) 程序运行结果如下： x1 =1 x2 =aAn*n x3 =1/2* nA2-1/2*n x4 =iztra ns((1-zA(-n))/(1-1 /z) ,z, n) 2、离散系统的零极点分析(系统极点位置对系统响应的影响) 例3-3研究z右半平面的实数极点对系统的影响。 已知系统的零极点增益模型分别为： H 1(z)= z z-0.85 H Q)=z-1 H 3(z) z-1.5 求这些系统的零极点分布图以及系统的单位序列响应，判断系统的稳定性。 程序清单如下: z1=[0];p 仁[0.85];k=1; [b1,a1]=zp2tf(z1,p1,k); subplot(3,2,1);zpla ne(z1,p1); title(极点在单位圆内’)； subplot(3,2,2);impz(b1,a1,20); z2=[0];p2=[1]; [b2,a2]=zp2tf(z2,p2,k); subplot(3,2,3);zpla ne(z2,p2); title(极点在单位圆上’)； subplot(3,2,4);impz(b2,a2,20); z3=[0];p3=[1.5]; [b3,a3]=zp2tf(z3,p3,k);subplot(3,2,5);zpla ne(z3,p3); title(极点在单位圆外); subplot(3,2,6);impz(b3,a3,20); 程序运行结果如图3-1所示。由图可见，这三个系统的极点均为实数且处于 z平面的右 半平面。由图可知，当极点位于单位圆内，系统的单位序列响应随着频率的增大而收敛；当 极点位于单位圆上，系统的单位序列响应为等幅振荡； 当极点位于单位圆外，系统的单位序 列响应随着频率的增大而发散。由此可知系统 1、2为稳定系统。 trdp w^anLaa 卩■■ar i：Real Part极点在单位圆内0-1-20 20.51J极点在单位圆上2o2- o 1 trdp w^anLaa 卩 ■ ■a r i： Real Part 极点在单位圆内 0 -1 -2 0 2 0.5 1 J 极点在单位圆上 2 o 2 - o 1 - yapcw ml Impulse Resp onse 10 15 n (sa