freqz_m2.m

matlab exprience,M文件,实验代码

function[db,mag,pha,w]=freqz_m2(b,a)
%滤波器幅值响应（绝对、相对）、相位响应
%db:相对幅值响应；
%mag:绝对幅值响应；
%pha:相位响应；
%w:采样频率；
%b:系统函数分子式系数向量
%a:系统函数分母式系数向量

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