afd_chb2.m

数字信号分析处理方面的好教材,有DSP的MATLAB源程序.

function [b,a] = afd_chb2(Wp,Ws,Rp,As);
% 切比雪夫-2型模拟低通滤波器设计
% -----------------------------------------
% [b,a] = afd_chb2(Wp,Ws,Rp,As);
%  b = Ha(s) 分子的系数
%  a = Ha(s) 分母的系数
% Wp = 以弧度/秒为单位的通带边缘频率; Wp > 0
% Ws = 以弧度/秒为单位的阻带边缘频率; Ws > Wp > 0
% Rp = 通带中的振幅波动的+dB数; (Rp > 0)
% As = 阻带衰减的+dB数; (As > 0)
%
if Wp <= 0
        error('通带边缘必须大于 0')
end
if Ws <= Wp
        error('阻带边缘必须大于通带边缘')
end
if (Rp <= 0) | (As < 0)
        error('通带波动或阻带衰减必须大于0')
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+sqrt(OmegaR*OmegaR-1)));
fprintf('\n*** 切比雪夫-2 滤波器阶次 = %2.0f \n',N)
[b,a]=u_chb2ap(N,As,Ws);