数字信号处理上机题32.docxVIP

32*.下面四个二阶网络的系统函数具有一样的极点分布: Hi⑵二 1 T.6z‘ 0.9425Z2 H2(z) H2(z) 1 - 1.6Z1 0.9425Z2 H3(z) H3(z) 1 - 1.6z_v 0.9425z「2 %程序 ex232.m A=[1,-1.6,0.9425]; %系统函数的分母多项式系数 B1=1;B2=[1,-0.3];B3=[1,-0.8];B4=[1,-1.6,0.8]; % 系统函数的分子多项式系数 b仁[1,0,0];b2=[1,-0.3,0];b3=[1,-0.8,0];b4=[1,-1.6,0.8]; % 系统函数的正次幕分子多项式系数 p=roots(A) z1=roots(b1) z2=roots(b2) z3=roots(b3) z4=roots(b4) %求系统函数的极点 %求H1(z)的零点 %求H2(z)的零点 %求H3(z)的零点 %求H4(z)的零点 [h1 n,n]=impz(B1,A,100); %计算单位脉冲响应 h1(n)的100个样值 [h2n,n]=impz(B2,A,100); %计算单位脉冲响应 h2(n)的100个样值 [h3n,n]=impz(B3,A,100); %计算单位脉冲响应 h3(n)的100个样值 [h4n,n]=impz(B4,A,100); %计算单位脉冲响应 h4(n)的100个样值 %================================================================== %以下是绘图部分 figure(1); subplot(221); zplane(B1,A); %绘制H1(z)的零极点图 subplot(222); stem(n,h1 n,.); %绘制 h1(n)的波形图 lin e([0,100],[0,0]) xlabel( n);ylabel(h1( n)) subplot(223); zplane(B2,A); %绘制H2(z)的零极点图 subplot(224); stem(n,h2n,.); %绘制 h2(n)的波形图 lin e([0,100],[0,0]) xlabel( n);ylabel(h2( n)) figure(2); subplot(221); zplane(B3,A); %绘制H3(z)的零极点图 subplot(222); stem(n,h3n,.); %绘制 h3(n)的波形图 lin e([0,100],[0,0]) xlabel(n );ylabel(h3( n)) subplot(223); zplane(B4,A); %绘制H4(z)的零极点图 subplot(224); stem(n,h4n,.); %绘制 h4(n)的波形图 lin e([0,100],[0,0]) xlabel(n );ylabel(h4 (n)) 程序运行结果如题 32*解图所示。 thap w^anroa 卩 1 0.5 0 -0.5 -1 - * -i 2 - 2 _ 1 r P 0.5 y a 0 g m -0.5 -1 -1 -0.5 0 0.5 1 Real Part -1 -0.5 0 0.5 1 2 -1 -2 L 0 50 100 Real Part 结果分析：四种系统函数的极点分布一样，只是零点不同， 第一种零点在原点，不影 响系统的频率特性，也不影响单位脉冲响应。 第二种的零点在实轴上，但离极点较远。 第三种的零点靠近极点。第四种的零点非常靠近极点，比较它们的单位脉冲响应， 会发现 零点愈靠近极点， 单位脉冲响应的变化愈缓慢， 因此零点对极点的作用起抵消作用； 同时, 第四种有两个零点，抵消作用更明显。