ms10_6.m

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

p=[0]';k=1;
z1=[0.2]';
n=(0:500)*pi/500;
[b1,a1]=zp2tf(z1,p,k);
[h1,w]=freqz(b1,a1,n);
subplot(2,3,1),zplane(b1,a1);
title('零点q1= 0.2');
z2=[0.5]';
[b2,a2]=zp2tf(z2,p,k);
[h2,w]=freqz(b2,a2,n);
subplot(2,3,2),zplane(b2,a2);
title('零点q1= 0.5');
z3=[0.8]';
[b3,a3]=zp2tf(z3,p,k);
[h3,w]=freqz(b3,a3,n);
subplot(2,3,3),zplane(b3,a3);
title('零点q1= 0.8');
subplot(2,1,2),plot(w/pi,abs(h1),w/pi,abs(h2),w/pi,abs(h3));
text(0.1,1,'q1= 0.2');text(0.04,0.6,'q1= 0.5');
text(0.1,0.3,'q1= 0.8');title('幅频特性');