[数字信号处理实验三.docVIP

实验三：离散LSI系统的频域分析 一、实验内容 2、求以下各序列的z变换： 程序清单如下： syms w0 n z a; x1=n*a^n;X1=ztrans(x1) x2=sin(w0*n);X2=ztrans(x2) x3= exp(-a*n)*sin(w0*n);X3=ztrans(x3) 程序运行结果如下： X1 =z/(a*(z/a - 1)^2) X2 =(z*sin(w0))/(z^2 - 2*cos(w0)*z + 1) X3 =(z*exp(a)*sin(w0))/(exp(2*a)*z^2 - 2*exp(a)*cos(w0)*z + 1) 3、求下列函数的逆z变换 程序清单如下： syms w0 n z a; X1=z/(z-a);x1=iztrans(X1) X2= z/(a-z)^2;x2=iztrans(X2) X3=z/ z-exp(j*w0);x3=iztrans(X3) X4=(1-z^-3)/(1-z^-1);x4=iztrans(X4) 程序运行结果如下： x1 =a^n x2 =n*a^n/a x3 =charfcn[0](n)-iztrans(exp(i*w0),w0,n) x4 =charfcn[2](n)+charfcn[1](n)+charfcn[0](n) 4、求一下系统函数所描述的离散系统的零极点分布图，并判断系统的稳定性 （1） z1=[0,0.3];p1=[-1+j,-1-j];k=1; [b1,a1]=zp2tf(z1,p1,k); subplot(1,2,1);zplane(z1,p1); title(极点在单位圆外); subplot(1,2,2);impz(b1,a1,20); 由图可见:当极点位于单位圆内，系统的单位序列响应随着频率的增大而收敛；当极点位于单位圆上，系统的单位序列响应为等幅振荡；当极点位于单位圆外，系统的单位序列响应随着频率的增大而发散。由此可知系统为不稳定系统。 （2） b=[4,-1.6,-1.6,4]; a=[1,0.4,0.35,-0.4]; rz=roots(b) rp=roots(a) subplot(2,1,1);zplane(b,a); title(系统的零极点分布图); subplot(2,1,2);impz(b,a,20); title(系统的单位序列响应); xlabel(n);ylabel(h(n)); rz = -1.0000 0.7000 + 0.7141i 0.7000 - 0.7141i rp = -0.4500 + 0.7730i -0.4500 - 0.7730i 0.5000 由零极点分布图可见，该系统的所有极点均在单位圆内，因此该系统是一个因果稳定系统。 5、已知某离散时间系统的系统函数为 求该系统在0～П频率范围内的绝对幅频响应与相频响应、相对幅频响应与相频响应及群时延。 b=[0.187632,0,-0.241242,0,0.241242,0,-0.187632]; a=[1,0,0.602012,0,0.495684,0,0.035924]; n=(0:500)*pi/500; [h,w]=freqz(b,a,n); db=20*log10(abs(h)); subplot(2,2,1);plot(w/pi,abs(h));grid; axis([0,1,1.1*min(abs(h)),1.1*max(abs(h))]); title(幅频特性（V）); xlabel(\omega/\pi);ylabel(幅度(V)); subplot(2,2,2);plot(w/pi,angle(h));grid; axis([0,1,1.1*min(angle(h)),1.1*max(angle(h))]); xlabel(\omega/\pi);ylabel(相位); title(相频特性); subplot(2,2,3);plot(w/pi,db);grid axis([0,1,-100,5]); title(幅频特性（dB）); subplot(2,2,4);zplane(b,a); title(零极点分布); 二、思考题答案 ① 系统函数零极点的位置与系统单位序列响应有何关系？ 当极点位于单位圆内，系统的单位序列响应随着频率的增大而收敛；当极点位于单位圆上，系统的单位序列响应为等幅振荡；当极点位于单位圆外，系统的单位序列响应随着频率的增大而发散。 ② 离散系统的零极点对系统幅频响应有何影响？ ① 在原点（z=0）处的零点或极点至单位圆的距离始终保持不变，其值|ejω|=1，所以，对幅度响应不起