基于MATLAB的低通滤波器的设计要点

通信系统综合设计与实践

题目基于MATLAB的低通滤波器设计院（系）名称信院通信系

专业名称通信工程

学生姓名

学生学号

指导教师

2013年 5 月25 日

目录

通信系统综合设计与实践第 1 页摘要 (2)

1．巴特沃斯低通数字滤波器简介 (3)

1.1选择巴特沃斯低通滤波器及双线性变换法的原因 (4)

1.2巴特沃斯低通滤波器的基本原理 (4)

1.2.1巴特沃斯低通滤波器的基本原理 (4)

1.2.2双线性变换法的原理............... . (5)

1.3数字滤波器设计流程图............... (7)

1.4数字滤波器的设计步骤............... (7)

2.巴特沃斯低通数字滤波器技术指标的设置 (8)

3.用matlab实现巴特沃斯低通数字滤波器的仿真并分析 (9)

3.1用matlab实现巴特沃斯低通数字滤波器的仿真 (9)

3.2波形图分析............... ............... (10)

4.用Simulink实现巴特沃斯低通数字滤波器的仿真并分析 (11)

4.1 Simulink简介............... ............... . (11)

4.2用Simulink实现巴特沃斯低通数字滤波器的仿真 (11)

4.3波形图分析............... ............... . (11)

4.3.1 Simulink波形图分析............... ....................... .. (11)

4.3.2与matlab波形的比较............... .. (14)

5.总结与体会.............................. (15)

6.附录 (16)

摘要

通信系统综合设计与实践第 2 页低通滤波器是让规定频率以下的信号分量通过，而对该频率以上的信号分量抑制的电容、电感与电阻等器件的组合装置。巴特沃斯滤波器是电子滤波器的一种，特点是通频带内的频率响应曲线最大限度平坦，没有起伏，而在阻频带则逐渐下降为零。可以构成低通、高通、带通和带阻四种组态，是目前最为流行的一类数字滤波器 ,经过离散化可以作为数字巴特沃思滤波器较模拟滤波器具有精度高、稳定、灵活、不要求阻抗匹配等众多优点因而在自动控制、语音、图像、通信、雷达等众多领域得到了广泛的应用，是一种具有最大平坦幅度响应的低通滤波器。本文将介绍其中最常用的一种——巴特沃斯低通数字滤波器。本文侧重于理论分析、matlab编程和结果分析。

（1）本文将先概述巴特沃斯低通数字滤波器的工作原理和特点。

（2）设置技术指标;

（3）用matlab进行软件编程，将仿真波形与理论值进行比较，分析其中的异同，并通过一个输入波形来验证设计的巴特沃斯低通数字滤波器的准确性。

（4）对实验结果和理论结果进行比较，分析它们的异同点并进行总体分析。

（5）用Simulink进行硬件电路仿真，观察仿真结果，并通过一个输入波形来验证仿真效果。

（6）通过对用Matlab软件仿真和Simulink硬件仿真的比较，说明两者的优缺点。

（7）对实验结果进行最后的总结，写出自己的感想。

关键字： matlab低通滤波器巴特沃斯

1.巴特沃斯低通数字滤波器简介

通信系统综合设计与实践 第 3 页

1.1选择巴特沃斯低通滤波器及双线性变换法的原因

（1）由于低通滤波器是组成其它滤波器的基础，故选用低通滤波器；

（2）在当今社会，数字信号的应用越来越广泛，故选用数字信号；

（3）巴特沃斯滤波器的特点是通频带的频率响应曲线最平滑并且应用范围最广，故选巴特沃斯型滤波器；

（4）为了不使数字滤波器在ω=π附近产生频谱混叠，故选用双线性变换法。

1.2巴特沃思低通滤波器的基本原理：

1.2.1巴特沃思低通滤波器的基本原理：

N a c

22)(11)j (H ΩΩ+=Ω 式中，N 称为滤波器的阶数。当Ω=0

c Ω=

c Ω是3dB 截止频率。在c Ω=Ω附近，随Ω加大，幅度迅速下降。幅度特性与Ω与N 的关系如图1.1所示。幅度下降的速度与阶数N 有关，N 愈大，通带愈平坦，过渡带愈窄，过渡带与阻带幅度下降的速度愈快，总的频响特性与理想低通滤波器的误差愈小。

图1.1 巴特沃斯低通数字滤波器 图1.2 三阶巴特沃斯滤波器极点

幅度特性与Ω与N 的关系 分布图

通信系统综合设计与实践 第 4 页

N a a c

j s s s 2)(11)(H )(H Ω+=- 复变量Ω+=j s σ，此式表示幅度平方函数有2N 个极点，极点k s 用下式表示：)21221(21

()1(N k j c N k e c j s ++Ω=Ω-=π） （k =0,1,2,3….）

2N 个极点等间隔分布在半径为c Ω的圆上（该圆称为巴特沃斯圆），间隔为N /πrad 。例如N=3,极点间隔为π/3rad ，如图1.2所示。

为形成因果稳定的滤波器，2N 个极点中只取ｓ平面左半平面的的N 个极点构成Ha(s), 而右半平面的的N 个极点构成Ha(-s),Ha （s ）的表达式为)

()(Ha 10k N k N c

s s s -Ω=∏-= 为使设计公式和图表统一，将频率归一化。巴特沃斯低通数字滤波器采用对3dB 截止频率c Ω归一化，归一化后的系统函数为

∏=ΩΩ=Ω1-0k )c

s -(1)(H N k c c s s a 令ληj p +=,c ΩΩ=/λ,λ称为归一化频率，p 称为归一化复变量，这样，巴特沃斯低通原型系统函数为

∏==1-0)

-(1

Ga(p)N k k p

p

1.2.2双线性变换法原理

双线性变换法是使数字滤波器的频率响应与模拟滤波器的频率响应相似的一种变换方法。为了克服多值映射的缺点，采用把整个s 平面频率压缩方法，

将整个频率轴上的频率范围压缩到-π/T ～π/T 之间，再用sT e Z =转换到Z 平面上。也就是说，第一步先将整个S 平面压缩映射到S1平面的-π/T ～π/T 一条横带里；第二步再通过标准变换关系sT e Z =将此横带变换到整个Z 平面上去。这样就使S 平面与Z 平面

通信系统综合设计与实践 第 5 页

建立了一一对应的单值关系，消除了多值变换性，也就消除了频谱混叠现象。映射关系如图1.3所示。

设Ha （s ）,Ω=j s ,经过非线性频率压缩后用^

1)(s H a ，11Ω=j s 表示，这里用正切变换实现频率压缩： )2

1tan(21T T Ω=Ω

图1.3 双线性变换的映射关系

式中，T 为采样间隔，当1Ω从-π/T 经过0变化到π/T 时，Ω由-∞经过0变化到+∞，实现了s 平面上整个虚轴完全压缩到1s 平面上虚轴的+π/T 之间的转换。即 T j T j T j T j tT j T j e e T e e e e T j 11111111222/2/2/2/Ω-Ω-ΩΩΩ-Ω+-=+-=Ω

代入Ω=j s ,11Ω=j s ,得到 T

s T

s e e T s 11112--+-= 再通过T s e z 1=从1s 平面转换到z 平面，得到 1

1

112--+-=z z T s s T

s T z -+=22 上式是S 平面与Z 平面之间的单值映射关系，这种变换都是两个线性函数之比，因此称为双线性变换。

通信系统综合设计与实践 第 6 页

双线性变换法与冲激响应不变法相比，其主要的优点是避免了频率响应的混叠现象，虽然在线性方面有些欠缺，但是可以通过频率的预畸来加以校正且计算比冲激响应不变法方便，实现起来比较容易，所以，本设计选择用双线性变换法设计巴特沃斯低通滤波器。

1.3数字滤波器设计流程图

1.4数字滤波器的设计步骤

数字滤波器的设计步骤：根据数字滤波器的技术指标先设计过渡模拟滤波器得到系统函数Ha(s),然后将Ha(s)按某种方法（本实验采用双线性变换法）转换成数字滤波器的系统函数H （z ）。具体为：

（1）确定巴特沃斯数字低通滤波器的技术指标：通带边界频率ωp,阻带截止频率ωs,通带最大衰减аp,阻带最小衰减аs 。

（2）将数字滤波器的技术指标转换为模拟滤波器的技术指标。这里指ωp 和ωs 的变换而аp 和аs 保持不变。本题采用双线性变换法，其转换公式为：

2tan 2

p T p ω=

Ω 2

tan 2

s T s ω=

Ω （3）根据技术指标Ωp 、Ωs 、ωp 和ωs 用下面公式求出滤波器的阶数。

p s

ΩΩ=

sp λ 1

1011010

/10

/--=p s sp k αα sp g sp g l k l λ=N

（4）根据N 由表1.4求出归一化极点k p 和归一化低通原型系统函数Ga(p)。

通信系统综合设计与实践 第 7 页

表1.4 巴特沃斯归一化低通滤波器参数

（5）将Ga （p ）去归一化，将c s p Ω=代入Ga （p ），得到实际的滤波器系统函数： c s

p p G s Ha Ω==)()(

这里Ωc 为3dB 截止频率。

（6）用双线性变换法将模拟滤波器Ha(s)转换成数字低通滤波器系统函数H(z)。转换公式为：

)()(s Ha z H =s=11

112--+-z

z T 2. 巴特沃斯低通数字滤波器技术指标的设置

数字低通技术指标为

ωp=30Hz, αp =1dB

ωs=60Hz αs =30dB

通信系统综合设计与实践第8 页采样频率为fs=500

设计要求：

输入x=sin(2*π*20*t)+2*sin(2*π*100*t)+5*sin(2*π*200*t)合成信号，经过滤波器后滤除30Hz以上的分量，即只保留sin(2*π*20*t)分量信号，来验证设计的滤波器是否达到了设计要求。

设计目的：

(1)理解低通滤波器的过滤方法。

(2)进一步熟悉低通滤波器的基本应用。

(3)用仿真工具matlab和Simulink分别对设计的滤波器进行软件和硬件仿真。

(6)将对仿真结果进行比较，从而检验滤波器滤波性能的准确性。

3.用matlab实现巴特沃斯低通数字滤波器的仿真并分析

3.1用matlab实现巴特沃斯低通数字滤波器的仿真

Matlab程序如下：

fs=500;

t=0:1/fs:1;

x=sin(2*pi*20*t)+2*sin(2*pi*100*t)+5*sin(2*pi*200*t);

wp=2*30/fs;

ws=2*60/fs;

Rp=1;

As=30;

subplot(311);

plot(t,x);

title('输入信号');

[N,wc]=buttord(wp,ws,Rp,As);

[B,A]=butter(N,wc);

[H,W]=freqz(B,A);

y=filter(B,A,x);

subplot(312);

plot(W,abs(H));

通信系统综合设计与实践 第 9 页

title('低通滤波器');

subplot(313);

plot(t,y)

title('30Hz');

用matlab 滤波前后的信号波形变化如图3.1所示：

10

输入信号

00.5

1 1.5

2 2.5

3 3.50

1

2

低通滤波器

00.20.40.60.81

2

30Hz

图3.1 用matlab 滤波前后的信号波形变化

3.2波形图分析：

由技术指标得：设计的巴特沃斯低通数字滤波器为30Hz 以内的信号能通过，而高于30Hz 的信号将通不过滤波器。因此，我们设计一个输入信号，其为：x=sin(2*π*20*t)+2*sin(2*π*100*t)+5*sin(2*π*200*t)。使其通过设计好的巴特沃斯低通数字滤波器，来验证设计好的巴特沃斯低通数字滤波器是否准确。

如图3.1所示，输入信号为x=sin(2*π*20*t)+2*sin(2*π*100*t)+5*sin(2*π*200*t)，其图形在1秒内重复了20次，故其频率为20Hz ，幅度为8，由sin(2*π*20*t)、2*sin(2*π*100*t)和5*sin(2*π*200*t)三个正弦波合成，从图中可知，输入信号的幅度与频率与输入信号的理论波形几乎完全相同。而通过低通滤波器后信号只剩下一个，

通信系统综合设计与实践第10 页且在1秒内波形重复了20.25次，即周期为1/20.25s，其频率为20.25Hz, 幅度为1，和输入信号中的sin(2*π*20*t)分量相差不大，但其频率稍微有些失真，说明了设计的低通滤波器达到了要求。

4.用Simulink实现巴特沃斯低通数字滤波器的仿真并分析

4.1Simulink简介

Simulink是Matlab最重要的组件之一，是MATLAB中的一种可视化仿真工具，是一种基于MATLAB的框图设计环境，是实现动态系统建模、仿真和分析的一个软件包它提供一个动态系统建模、仿真和综合分析的集成环境。在该环境中，无需大量书写程序，而只需要通过简单直观的鼠标操作，就可构造出复杂的系统。Simulink具有适应面广、结构和流程清晰及仿真精细、贴近实际、效率高、灵活等优点，并基于以上优点Simulink已被广泛应用于控制理论和数字信号处理的复杂仿真和设计。同时有大量的第三方软件和硬件可应用于或被要求应用于Simulink。

4.2用Simulink实现巴特沃斯低通数字滤波器的仿真

用Simulink硬件仿真的电路图4.1如下所示

图4.1 用Simulink硬件仿真的原理图

4.3波形图分析

4.3.1 Simulink波形图分析

通信系统综合设计与实践第11 页图中Sine Wave为正弦波，函数分别为sin(2*π*20*t)、2*sin(2*π*100*t)和5*sin(2*π*200*t)；三个加号的为求和函数；Digital Filter Design为巴特沃斯低通数字滤波器，其设置如图4.2所示；而Scope为三个信号相加后的波形（如图4.3、4.4所示），用来观察输入信号；Scope1为经过巴特沃斯低通数字滤波器后的波形（如图4.5所示），用来观察滤波效果。

Digital Filter Design的设置如图三所示，其设置的为巴特沃斯低通数字滤波器，Scope为输入信号，由图4.3的，其幅度略大于五，其周期为0.01s，频率为100Hz，与输入信号x=sin(2*π*20*t)+2*sin(2*π*100*t)+5*sin(2*π*200*t)的幅度8，频率200Hz有较大差距，说明失真较严重存在不小的信号衰减，而Scope1（如图4.5所示）为滤波后的信号，其幅度为1.1，在5s到5.5s内波形重复了10次，故其周期为0.05s，频率为20Hz，与输入信号中的sin(2*π*20*t)相当，幅度变大了0.1，说明有噪声干扰，但也基本达到了设计要求。

通信系统综合设计与实践第12 页

图4.2 Digital Filter Design

通信系统综合设计与实践第13 页

图4.3 输入信号的波形

图4.4 输入信号的波形

通信系统综合设计与实践第14 页

图4.5 滤波后的信号

4.3.2用Matlab与Simulink仿真结果的比较：

输入信号的比较

用Simulink出来的信号干扰较大，失真较严重，说明存在外界干扰，导致信号衰减，而用matlab出来的信号较好，干扰可以忽略。

巴特沃斯低通数字滤波器的比较

由图3.1中的低通滤波器和图4.2比较可得可知，用matlab设计的滤波器过渡带较窄，而用Simulink设计的滤波器过渡带较大，仿真结果和matlab相比效果较差，用matlab 设计的滤波器与理论值差距较小，滤波效果更好。

输出信号的比较

用matlab滤波后的信号比理论信号在1秒内大了1/4个周期，效果不错，而用Simulink滤波后的信号周期没变，幅度大了0.1，但滤波器控制不了幅度，幅度的变化说明存在外界干扰，说明噪声较大，但从滤波效果来看也达到了设计要求。

综述

由图4.4和图4.5可知，用Simulink设计的滤波器的输入信号与滤波后信号均存在5s的延迟，说明了加法器的处理速度较慢，而用matlab设计的滤波器不存在此现象。

本文来源：https://www.bwwdw.com/article/490q.html