freqz_m.m

来自「原名：数字信号处理基础及MATLAB实现,解压缩密码：lfjk2002@163.」· M 代码 · 共 11 行

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);