ms10r5.m

《数字信号处理实验(MATLAB版)》这

z=[-1]';k=1;
n=(0:500)*pi/500;
p1=[-0.2]';
[b1,a1]=zp2tf(z,p1,k);
[h1,w]=freqz(b1,a1,n);
subplot(2,3,1),zplane(b1,a1);
title('极点p1=0.2');
p2=[-0.5]';
[b2,a2]=zp2tf(z,p2,k);
[h2,w]=freqz(b2,a2,n);
subplot(2,3,2),zplane(b2,a2);
title('极点p1=0.5');
p3=[-0.8]';
[b3,a3]=zp2tf(z,p3,k);
[h3,w]=freqz(b3,a3,n);
subplot(2,3,3),zplane(b3,a3);
title('极点p1=0.8');
subplot(2,1,2),plot(w/pi,abs(h1),'b'),hold on;
subplot(2,1,2),plot(w/pi,abs(h2),'g'),hold on;
subplot(2,1,2),plot(w/pi,abs(h3),'r');
text(0.08,1.7,'p1=0.2');
text(0.05,1.35,'p1=0.5');
text(0.08,1,'p1=0.8');title('幅频特性');