afd_butt.m

原名：数字信号处理基础及MATLAB实现,解压缩密码：[e-mail]

function [b,a]=afd_butt(Wp,Ws,Rp,Ar)
%按给定的技术指标设计模拟巴特沃思低通滤波器
%b直接型滤波器分子系数；a直接型滤波器分母系数
%Wp通带截止频率，单位弧度/秒；Ws阻带截止频率，单位弧度/秒；
%Rp通带中的振幅波动系数（dB）；Ar阻带衰减，（dB）；
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)));
fprintf('\n Butterworth 滤波器阶数 = %2.0f \n',N)
OmegaC=Wp/((10^(Rp/10)-1)^(1/(2*N)));
[b,a]=u_buttap(N,OmegaC);