freqz_m.m

matla窗口制作

%函数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);
%__________________________________________________________________________