freqz_m.m

有关多相滤波的matlab的程序实现

%p185，西电“信号处理”
%freqz_m函数可获得滤波器的幅值响应（绝对和相对），相位响应及群迟延响应
%this function  is the modified version of freqz subroutine
%db is relative magnitude in db computed over 0 to pi radians
%mag is absolute magnitude computed over 0 to pi radians
%pha is phase response in radians over 0 to pi radians
%grd is group delay over 0 to pi radians
%w is 501 frequency samples between 0 and pi radians
%b is numerater polynomial of H(z) (for FIR:b=h)
%a is denominator polynomial of H(z)  (for FIR:a=[1])
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);