freqs_m.m

来自「原名：数字信号处理基础及MATLAB实现,解压缩密码：lfjk2002@163.」· M 代码 · 共 9 行

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