function [db,mag,pha,grd,w]=freqz_m(b,a)
%z域數字頻率響應計算freqz的改進版本
%b直接型濾波器分子系數；a直接型濾波器分母系數
%db為[0:pi]之間的相對振幅；mag為[0 pi]區間的絕對值；pha為[0 pi]區間的相位響應；
%grd為[0 pi]區間的群延遲；w為[0 pi]區間內的501個頻率樣本數組
[H,w]=freqz(b,a,1000,'whole');
   H=(H(1:1:501))';w=(w(1:501))';
mag=abs(H);
db=20*log10((mag+eps)/max(mag));
pha=angle(H);
grd=grpdelay(b,a,w);