ex3p13.m

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

w=[0:1:500]*pi/500;
H=exp(j*w)./(exp(j*w)-0.9*ones(1,501));
magH=abs(H);angH=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,angH/pi);grid;
xlabel('frequency in units of pi');ylabel('phase in pi Radians');
title('Phase Response');