freqz_m.asv

来自「FIR滤波器设计」· ASV 代码 · 共 15 行

function [db,mag,pha,w]=freqz_m(b,a);
%db=[0-pi弧度]区间内的相对振幅
%mag=[0-pi弧度]区间内的绝对振幅
%pha=[0-pi弧度]区间内的相位响应
%grd=[0-pi弧度]区间内的群延迟
%w=[0-pi弧度]区间内的501个频率样本向量
%b=Ha(z)的分子多项式系数（对FIR:b=h）
%a=Ha(z)的分子多项式系数（对FIR:b=[1]）
[H,w]=freqz(b,a,1000,'whole');
H=(H(1:1:501))';
w=(w(1:1:501))';
mag=abs(H)
db=20*log10((mag+eps)/max(mag))
pha=angle(H)