e313.m

matlab例子

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;
ylabel('|H|');title('振幅特性');
subplot(2,1,2);plot(w/pi,angH/pi);grid;
xlabel('频率(以pi为单位)');ylabel('相位(单位：pi)');
title('相位特性');