freqs_m.m

来自「FIR滤波器设计」· M 代码 · 共 7 行

function [dB,mag,pha,Omega]=freqs_m(b,a,pi);
Omega=[0:1:500]*pi/500;
H=freqs(b,a,Omega);
mag=abs(H);
dB=20*log10((mag+eps)/max(mag));
pha=angle(H);