ex4p13a.m

关于数字信号处理的一些matlab例程

w=[0:1:500]*pi/500;
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('Magnitude');
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');