afd_chb2.m

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

function [b,a]=afd_chb2(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
ep=sqrt(10^(Rp/10)-1);
A=10^(Ar/20);
OmegaC=Wp;
OmegaR=Ws/Wp;
g=sqrt(A*A-1)/ep;
N=ceil(log10(g+sqrt(g*g-1))/log10(OmegaR+sqrt(OmegaR*OmegaR-1)));
fprintf('\n 切比雪夫Ⅱ型滤波器阶数 = %2.0f \n',N)
[b,a]=u_chb2ap(N,Ar,Ws);