freqz_m.m

求取系统的绝对幅度响应、相对的db值幅度响应、相位响应和群延时响应的函数 ％ ％db为相对振幅(dB) ％mag为绝对振幅 ％pha为相位响应 ％grd为群延时 ％w为频率样本点向量 ％

function [db,mag,pha,grd,w]=freqz_m(b,a)
%求取系统的绝对幅度响应、相对的db值幅度响应、相位响应和群延时响应的函数
%
%db为相对振幅(dB)
%mag为绝对振幅
%pha为相位响应
%grd为群延时
%w为频率样本点向量
%b为Ha(z)分子多项式系数(对FIR而言，b=h)
%a为Hz(z)分母多项式系数(对FIR而言，a=1)
%


[H,w]=freqz(b,a,1000,'whole');
%freqz显示数字滤波器频域中的图形
%[H,W] = FREQZ(B,A,N,'whole') uses N points around the whole unit circle.
H=(H(1:501))';
w=(w(1:501))';
mag=abs(H);
db=20*log10((mag+eps)/max(mag));
%dB=-20log10|H(ejw)|/|H(ejw)|max>=0
pha=angle(H);
grd=grpdelay(b,a,w);