ex0429.m

数字信号处理（机械工业出版社）的源码

b = [1 0 -1]; a = [1 0 -0.81];
[R,p,C] = residuez(b,a)
w=[0:500]*pi/500;
H=freqz(b,a,w);
magH=abs(H);phaH=angle(H);
subplot(211);plot(w/pi,magH);grid
xlabel('frequency in pi units');ylabel('Magnitude');
title('Magnitude Response')
subplot(212);plot(w/pi,phaH/pi);grid
xlabel('frequency in pi units');ylabel('Phase in pi units');
title('Phase Response')