ex4p11b.m

来自「关于数字信号处理的一些matlab例程」· M 代码 · 共 9 行

w=[0:1:100]*pi/100;
H=freqz(b,a,w);
magH=abs(H);phaH=angle(H);
subplot(2,1,1);plot(w/pi,magH);grid;
xlabel('frequency in units of pi');ylabel('|H|');
title('Magnitude Response');
subplot(2,1,2);plot(w/pi,phaH/pi);grid;
xlabel('frequency in units of pi');ylabel('phase in pi Radians');
title('Phase Response');