freqz_m2.m

来自「<Matlab 7辅助信号处理技术与应用》的源码」· M 代码 · 共 16 行

function[db,mag,pha,w]=freqz_m2(b,a)
% 滤波器幅值响应(绝对、相对)、相位响应
% db: 相对幅值响应；
% mag: 绝对幅值响应；
% pha: 相位响应；
% w: 采样频率；
% b: 系统函数H(z)的分子项(对FIR，b＝h)
% a: 系统函数H(z)的分母项(对FIR，a＝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);