afd_chb1.m

来自「数字信号处理的MATLAB仿真加实验指导」· M 代码 · 共 23 行

function[b,a]=afd_chb1(Omegap, Omegar,Ap,Ar);


if Omegap<=0
   error('通带边缘必须大于0')
end
if Omegar<=Omegap
   error('阻带边缘必须大于通带边缘')
end
if(Ap<=0)|(Ar<0)
   error('通带波动或阻带衰减必须大于0')
end
   
ep=sqrt(10^(Ap/10)-1);
A=10^(Ar/20);
OmegaC=Omegap;
OmegaR=Omegar/Omegap;
g=sqrt(A*A-1)/ep;
N=ceil(log10(g+sqrt(g*g-1))/log10(OmegaR+sqrt(OmegaR*OmegaR-1)));
fprintf('\n***CB Filter Order=%2.0f\n',N);
[b,a]=u_chb1ap(N,Ap,OmegaC);