freqz_m.m

来自「mb-ofdm-uwb经本人修改」· M 代码 · 共 11 行

function[db,mag,pha,grd,w]=freqz_m(b,a)
%db为相对振幅(dB);mag为绝对振幅;pha为相位响应;grd为群延时;w为频率样本点向量。
[H,w]=freqz(b,a,1024,'whole');
%b和a分别为H(z)分子和分母多项式系数(对FIR而言b=a,a=1).
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);