% 数字滤波器指标
wp = 0.2*pi;
% digital Passband freq in Hz
ws = 0.3*pi;
% digital Stopband freq in Hz
Rp = 1;
% Passband ripple in dB
As = 15;
% Stopband attenuation in dB
例4:设计一个低通椭圆滤波器,以满足: 通带截止频率:Ωp=0.2π,通带波动:Rp=1dB 阻带截止频率:Ωs=0.3π,阻带波动:As=16dB
wp = 0.2*pi; Ws = 0.3*pi; Rp = 1; As = 16; % 模拟滤波器设计: [b,a] = afd_elip(Wp,Ws,Rp,As); % 计算频率响应: w = [0:1:500]*pi/500; h = freqs(b,a,w); % 画图 subplot(2,2,1); plot(w/pi,abs(h));title('幅度响应');grid; subplot(2,2,3); plot(w/pi,20*log10(abs(h))); title('幅度响应(dB)');grid; subplot(2,2,2); plot(w/pi,angle(h)/pi); title('相位响应');grid; subplot(2,2,4); impulse(b,a);title('脉冲响应');grid;
2、切比雪夫Ⅰ型低通滤波器的设计
MATLAB提供了一个函数[z,p,k]=cheblap(N,Rp), 来设计一个阶数为N,通带波动为Rp的归一化切比雪夫 Ⅰ型原型滤波器。
可以自编写函数u_chblap用于设计一个未归一化的 切比雪夫Ⅰ型原型滤波器,它返回直接形式的Ha(s)。