MATLAB信号处理学习-典型模拟低通滤波器

巴特沃斯滤波器

巴特沃斯滤波器的幅频特性

巴特沃斯滤波器的幅度平方特性定义为

                                ​​​​​​​        ​​​​​​​        

主要参数有2个：N是滤波器的阶次，wc是通带截至频率。

幅频特性有如下特点：

1、直流增益为1。在w=0处，的前2N阶导数均为0，在低频处有最平坦的特性。

2、当w=wc时，=0.707，=0.5，纹波=3dB，且与阶次无关。

3、通带和阻带频率响应具有单调下降的特性。

切比雪夫Ⅰ型/Ⅱ型

巴特沃斯滤波器在通带和阻带上的响应都是单调的，切比雪夫具有波动性，Ⅰ型在通带中具有等纹波响应，Ⅱ型在阻带中具有等纹波响应。等纹波特性比单调特性的滤波器具有较低的阶次。

因此对应相同的指标，切比雪夫比巴特沃斯的阶数低。

椭圆型滤波器

椭圆型低通滤波器的幅频特性在通带和阻带内均具有等纹波特性。

设计-代码

设计一个带通滤波器，带通值为wp1=0.2pai,wp2=0.3pai,带阻值为ws1=0.1pai,ws2=0.4pai,Rp=1,Rs=20。对这些指标分别以巴特沃斯、切比雪夫Ⅰ型Ⅱ型和椭圆型设计4类模拟滤波器。

%设计一个模拟带通滤波器，带通值为wp1=0.2pi,wp2=0.3pi,带阻值为ws1=0.1pi,ws2=0.4pi,Rp=1,Rs=20。

clear all; close all; clc;

wp=[0.2*pi 0.3*pi];              % 设置通带频率
ws=[0.1*pi 0.4*pi];              % 设置阻带频率
Rp=1; Rs=20;                     % 设置波纹系数
% 巴特沃斯滤波器设计
[N,Wn]=buttord(wp,ws,Rp,Rs,'s'); % 求巴特沃斯滤波器阶数
fprintf('巴特沃斯滤波器 N=%4d\n',N) % 显示滤波器阶数
[bb,ab]=butter(N,Wn,'s');        % 求巴特沃斯滤波器系数
W=0:0.01:2;                      % 设置模拟频率
[Hb,wb]=freqs(bb,ab,W);          % 求巴特沃斯滤波器频率响应
plot(wb/pi,20*log10(abs(Hb)),'b')% 作图
hold on

% 切比雪夫I型滤波器设计
[N,Wn]=cheb1ord(wp,ws,Rp,Rs,'s');  % 求切比雪夫I型滤波器阶数
fprintf('切比雪夫I型滤波器 N=%4d\n',N) % 显示滤波器阶数
[bc1,ac1]=cheby1(N,Rp,Wn,'s');     % 求切比雪夫I型滤波器系数
[Hc1,wc1]=freqs(bc1,ac1,W);        % 求切比雪夫I型滤波器频率响应
plot(wc1/pi,20*log10(abs(Hc1)),'k')% 作图

% 切比雪夫II型滤波器设计 
[N,Wn]=cheb2ord(wp,ws,Rp,Rs,'s');  % 求切比雪夫II型滤波器阶数
fprintf('切比雪夫II型滤波器 N=%4d\n',N) % 显示滤波器阶数
[bc2,ac2]=cheby2(N,Rs,Wn,'s');    % 求切比雪夫II型滤波器系数
[Hc2,wc2]=freqs(bc2,ac2,W);       % 求切比雪夫II型滤波器频率响应
plot(wc2/pi,20*log10(abs(Hc2)),'r')% 作图

% 椭圆型滤波器设计
[N,Wn]=ellipord(wp,ws,Rp,Rs,'s');  % 求椭圆型滤波器阶数
fprintf('椭圆型滤波器 N=%4d\n',N) % 显示滤波器阶数
[be,ae]=ellip(N,Rp,Rs,Wn,'s');     % 求椭圆型滤波器系数
[He,we]=freqs(be,ae,W);            % 求椭圆型滤波器频率响应
% 作图
plot(we/pi,20*log10(abs(He)),'g')
axis([0 max(we/pi) -30 2]); %grid;
legend('巴特沃斯滤波器','切比雪夫I型滤波器','切比雪夫II型滤波器','椭圆型滤波器')
xlabel('角频率{\omega}/{\pi}'); ylabel('幅值/dB')
set(gcf,'color','w'); 

line([0 max(we/pi)],[-20 -20],'color','k','linestyle','--');
line([0 max(we/pi)],[-1 -1],'color','k','linestyle','--');
line([0.2 0.2],[-30 2],'color','k','linestyle','--');
line([0.3 0.3],[-30 2],'color','k','linestyle','--');

结果

最后得到4类滤波器阶数分别为4；3；3；2。

巴特沃斯和切比雪夫Ⅱ型在通通带是平坦的，切比雪夫Ⅰ型和椭圆型滤波器在通带有波纹，切比雪夫Ⅱ型和椭圆型滤波器在阻带有波纹。

在满足相同的设计要求的条件下，切比雪夫Ⅱ型滤波器和椭圆型滤波器在过渡带内衰减更陡。