freqz_m.m

本matlab函数

function[db,mag,pha,grd,w]=freqz_m(b,a);
%db是0到π内的相对振幅，mag是0到π内的绝对振幅，pha是0到π内的相位响应
%grd是0到π内的群延迟，w是0到π内的501个频率样本向量
%b是直接型滤波器分子多项式系数，a是直接型滤波器分母多项式系数
[H,w]=freqz(b,a,1000,'whole');
H=(H(1:1:501))';
w=(w(1:1:501))';
mag=abs(H);
db=20*log10((mag+eps)/max(mag));    %num2str(eps) 的数值是2.22e-16
pha=angle(H);
grd=grpdelay(b,a,w);