freqs_m.m

FIR滤波器设计

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