afd_butt.m

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

function [b,a]=afd_butt(Omegap,Omegar,Ap,Ar);
if Omegap<=0
    error('通带边缘必须大于0')
end
if (Omegar<=Omegap)
    error('阻带边缘必须大于通带边缘必须')
end
N=ceil((log10((10^(0.1*Ar)-1)/(10^(0.1*Ap)-1)))/(2*log10(Omegar/Omegap)));
fprintf('\n***巴特沃斯模拟低通滤波器阶数=%2.0f\n',N);
OmegaC=Omegap/((10^(0.1*Ap)-1)^(1/(2*N)));
[b,a]=u_buttap(N,OmegaC);