afd_butt.m

来自「Matlab环境下带控制面板的低通、高通滤波器设计」· M 代码 · 共 9 行

function [b,a]=afd_butt(Omegap,Omegar,Ap,Ar,N0);
N=ceil((log10((10^(0.1*Ar)-1)/(10^(0.1*Ap)-1)))/(2*log10(Omegar/Omegap)));
if N0>=N
    N=N0;
end
fprintf('\n***巴特沃斯模拟滤波器阶次=%2.0f\n',N);
OmegaC=Omegap/((10^(Ap/10)-1)^(1/(2*N)));
[b,a]=u_buttap(N,OmegaC);