afd_butt.m

巴特沃斯滤波器巴特沃斯滤波器巴特沃斯滤波器巴特沃斯滤波器

function [b,a]=afd_butt(wp,ws,rp,ar)
if wp<=0
    error('通带边缘必须大于0')
end
if ws<=wp
   error('阻带边缘必须大于通带边缘')
end
if (rp<=0)|(ar<0)
    error('通带波动或阻带衰减必须大于0')
end
N=ceil(log10((10^(rp/10)-1)/(10^(ar/10)-1))/(2*log10(wp/ws)));
OmegaC=wp/((10^(rp/10)-1)^(1/(2*N)))
[b,a]=u_buttap(N,OmegaC);