freqs_m.m

来自「数字信号处理的比较全面的程序代码」· M 代码 · 共 10 行

function [db,mag,pha,w]=freqs_m(b,a,wmax);
w=[0:1:500]*wmax/500;
H=freqs(b,a,w);
mag=abs(H);
db=20*log10((mag+eps)/max(mag));
pha=angle(H);
subplot(2,2,1);plot(w,mag);
subplot(2,2,2);plot(w,db);
subplot(2,2,3);plot(w,pha);