freqz_m.m

来自「按照给定指标采用模拟低通滤波器（巴特沃思低通）去逼近理想模拟低通滤波器」· M 代码 · 共 8 行

%数字滤波器响应子程序
function [db,mag,pha,grd,w]=freqz_m(b,a);
[H,w]=freqz(b,a,1000,'whole');
H=(H(1:501))';w=(w(1:501))';
mag=abs(H);
db=20*log10((mag+eps)/max(mag));
pha=angle(H);
grd=grpdelay(b,a,w);