freqs_m.m

巴特沃斯滤波器设计

%计算系统函数的幅度响应和相位响应
function [db,mag,pha,w]=freqs_m(b,a)
w=[0:1:500]*pi/500;
H=freqs(b,a,w);
mag=abs(H);
db=20*log10((mag+eps)/max(mag));
pha=angle(H);