freqz_m.m

来自「是一些关于matlab的实验fir滤波器代码」· M 代码 · 共 14 行

function [db, mag, pha, grd, w]= freqz_m (b, a)
% db = [0 到 pi 弧度]区间内的相对振幅（db）
% mag = [0 到 pi 弧度]区间内的绝对振幅
% pha = [0 到 pi 弧度]区间内的相位响应
% grd = [0 到 pi 弧度]区间内的群延迟
% w = [0 到 pi 弧度]区间内的501个频率样本向量
% b = Ha(z)的分子多项式系数 （FIR: b=h）
% a = Ha(z)的分母多项式系数 （FIR: a=[1]）
[ H, w] =freqz (b, a, 1000, 'whole');  %1000 为所画点数，whole表示整个频率范围
H= (H(1: 1: 501))';  w = (w(1:1:501))';
mag =abs(H);
db = 20* log10((mag +eps ) /max(mag));% eps为MATLAB默认的一很小的数
pha =angle(H);
grd =grpdelay (b, a, w); %群延迟函数