function.m

本代码是我自己写的

function [b,a]=afd_chb1(wp,ws,rp,as);
if wp<=0
    error
end
if ws<=wp
    error
end
if(rp<=0)|(as<0)
    error
end
ep=sqrt(10^(rp/10)-1);
a=10^(as/20);
OmegaC=wp;
OmegaR=ws/wp;
g=sqrt(a*a-1)/ep;
N=ceil(log10(g+sqrt(g*g-1))/log10(OmegaR+OmegaR*OmegaR-1)));
fprintf('\n*** 切比雪夫-1滤波器阶次=%2.0f\n',N)
[b,a]=u_chb1ap(N,rp,OmegaC);