%函數freqz_m
function[dp,mag,pha,grd,w]=freqz_m(b,a);
%[db,mag,pha,grd,w]=freqz_m[b,a]
%求某傳輸函數幅頻特性和相頻特性
%db為0到pi弧度區間的相對振幅
%mag為0到pi弧度區間的絕對振幅
%grd為0到pi弧度區間的群延遲
%H為包含了離散系統頻響在 0——pi范圍內N個頻率等分點的值，
%向量w則包含范圍內N個頻率等分點
%b為Ha(z)的分子多項式的系數（對FIR：b=h)
%a為Ha(z)的分子的多項式系數（對FIR：a=[1])
%__________________________________________________________________________
[H,w]=freqz(b,a,1000,'whole');
H=(H(1:1:501))';
w=(w(1:1:501))';
mag=abs(H);
dp=20*log((mag+eps)/max(mag));
pha=angle(H);
grd=grpdelay(b,a,w);
%__________________________________________________________________________