freqs_m.m

原名：数字信号处理基础及MATLAB实现,解压缩密码：[e-mail]

function [db,mag,pha,w]=freqs_m(b,a,wmax)
%s域频率响应计算freqs的改进版本
%wmax希望获得的最大频率响应；b直接型滤波器分子系数；a直接型滤波器分母系数
%mag为[0 wmax]区间的绝对值；pha为[0 wmax]区间的相位响应；w为[0 wmax]区间内的500个频率样本数组
w=[0:1:500]*wmax/500;
H=freqs(b,a,w);
mag=abs(H);
db=20*log10((mag+eps)/max(mag));
pha=angle(H);