ex0427.asv

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

% Chapter 4: Example 4.27:
%
b = [0 1 1]; a = [1 -0.9 0.81];
[R,p,C] = residuez(b,a)
pause
Mp=abs(p')
Ap=angle(p')/pi
pause
%Frequency Response is below
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')