freqs_m.m

来自「MATLAB信号处理扩展函数,用于信号处理的扩展函数」· M 代码 · 共 6 行

function [db,mag,pha,w]=freqs_m(b,a,wmax);%freqs修正
         w=[0:1:500]*wmax/500;
         H=freqs(b,a,w);
         mag=abs(H);
         db=20*log10((mag+eps)/max(mag));
         pha=angle(H);