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

1、通信系统综合设计与实践题目院（系）名称专业名称学生姓名学生学号指导教师基于MATLAB的低通滤波器设计信院通信系通信工程2013 年 5 月 25 日通信系统综合设计与实践第 11 页摘要1巴特沃斯低通数字滤波器简介31.1 选择巴特沃斯低通滤波器及双线性变换法的原因41.2 巴特沃斯低通滤波器的基本原理41.2.1 巴特沃斯低通滤波器的基本原理41.2.2 双线性变换法的原理51.3 数字滤波器设计流程图71.4 数字滤波器的设计步骤72. 巴特沃斯低通数字滤波器技术指标的设置83. 用 matlab 实现巴特沃斯低通数字滤波器的仿真并分析93.1 用 matlab 实现巴特沃斯低通数字滤波

2、器的仿真93.2 波形图分析104. 用 Simulink 实现巴特沃斯低通数字滤波器的仿真并分析114.1 Simulink 简介 114.2 用 Simulink 实现巴特沃斯低通数字滤波器的仿真114.3 波形图分析114.3.1 Simulink波形图分析114.3.2 与 matlab 波形的比较145. 总结与体会156. 附录 16摘要低通滤波器是让规定频率以下的信号分量通过，而对该频率以上的信号分量抑制的电容、电感与电阻等器件的组合装置。巴特沃斯滤波器是电子滤波器的一种，特点是通频带内的频率响应曲线最大限度平坦，没有起伏，而在阻频带则逐渐下降为零。可以构成低通、高通、带通和带阻

3、四种组态，是目前最为流行的一类数字滤波器, 经过离散化可以作为数字巴特沃思滤波器较模拟滤波器具有精度高、稳定、灵活、不要求阻抗匹配等众多优点因而在自动控制、语音、图像、通信、雷达等众多领域得到了广泛的应用，是一种具有最大平坦幅度响应的低通滤波器。本文将介绍其中最常用的一种巴特沃斯低通数字滤波器。本文侧重于理论分析、matlab 编程和结果分析。（ 1）本文将先概述巴特沃斯低通数字滤波器的工作原理和特点。（ 2）设置技术指标;（ 3）用matlab 进行软件编程，将仿真波形与理论值进行比较，分析其中的异同，并通过一个输入波形来验证设计的巴特沃斯低通数字滤波器的准确性。（ 4）对实验结果和理论结果

4、进行比较，分析它们的异同点并进行总体分析。（ 5）用Simulink 进行硬件电路仿真，观察仿真结果，并通过一个输入波形来验证仿真效果。（ 6）通过对用Matlab 软件仿真和Simulink 硬件仿真的比较，说明两者的优缺点。（ 7）对实验结果进行最后的总结，写出自己的感想。关键字：matlab 低通滤波器巴特沃斯1. 巴 特 沃 斯 低 通 数 字 滤 波 器 简 介1.1 选择巴特沃斯低通滤波器及双线性变换法的原因1）由于低通滤波器是组成其它滤波器的基础，故选用低通滤波器；2）在当今社会，数字信号的应用越来越广泛，故选用数字信号；3）巴特沃斯滤波器的特点是通频带的频率响应曲线最平滑并且应

5、用范围最广，故选巴特沃斯型滤波器；4）为了不使数字滤波器在 = 附近产生频谱混叠，故选用双线性变换法。1.2 巴特沃思低通滤波器的基本原理：1.2.1 巴特沃思低通滤波器的基本原理：巴特沃斯低通数字滤波器的幅度平方函数H（a j ） 2用下式表示Ha(j )21(2Nc式中，N称为滤波器的阶数。当=0时，H（a j ）1 ； c 时，H(a j )1/ 2，c 是 3dB 截止频率。在c 附近，随加大，幅度迅速下降。幅度特性与与 N 的关系如图1.1 所示。幅度下降的速度与阶数N有关， N愈大，通带愈平坦，过渡带愈窄，幅度特性与与 N的关系分布图2j ，将幅度平方函数H （ a j ） 写成的

6、函数1H a (s) H a( s)1 ( s )2Njc复变量 s j ，此式表示幅度平方函数有2N 个极点，极点sk 用下式表示：sk( 1)2N(j c)ce1 2k 1j (2 2N )k =0,1,2,32N 个极点等间隔分布在半径为c 的圆上(该圆称为巴特沃斯圆)，间隔为/ N rad。例如 N=3,极点间隔为/3rad ，如图 1.2 所示。为形成因果稳定的滤波器，2N 个极点中只取平面左半平面的的N 个极点构成NHa(s), 而右半平面的的N个极点构成Ha(-s),Ha ( s)的表达式为Ha(s) N 1 c(s sk)k0为使设计公式和图表统一，将频率归一化。巴特沃斯低通数

7、字滤波器采用对3dB截止频率 c归一化，归一化后的系统函数为Ha( s )N-11sskc (s- k)k0 c c令 p j ,/ c , 称为归一化频率，p 称为归一化复变量，这样，巴特 沃斯低通原型系统函数为1Ga(p) N-1(p-pk)k01.2.2 双线性变换法原理双线性变换法是使数字滤波器的频率响应与模拟滤波器的频率响应相似的一种变换方法。为了克服多值映射的缺点，采用把整个s 平面频率压缩方法，将整个频率轴上的频率范围压缩到- /T /T 之间，再用Z esT转换到Z平面上。也就是说， 第一步先将整个S 平面压缩映射到S1 平面的 - /T /T 一条横带里；第二步再通过标准变换

8、关系Z esT将此横带变换到整个Z平面上去。这样就使S平面与 Z平面建立了一一对应的单值关系，消除了多值变换性，也就消除了频谱混叠现象。映射关系 如图 1.3 所示。设 Ha( s) , s j , 经过非线性频率压缩后用Ha(s1)， s1 j 1表示，这里用正切变换实现频率压缩：21 tan(1T)T 21图 1.3 双线性变换的映射关系式中， T为采样间隔，当 1从 - /T 经过 0变化到 /T 时， 由 -经过 0变化到 +，实现了 s 平面上整个虚轴完全压缩到s1 平面上虚轴的+ /T 之间的转换。即2 ej 1T/2 e j 1tT/22 1 e j 1TjT ej 1T/2ej

9、 1T/2T 1 e j 1Teee代入 s j , s1j 1 , 得到2 1 e s1T sT 1 e s1T2s再通过 zes1T 从s1 平面转换到z 平面，得到21 z1s1T1 z1上式是 S 平面与 Z 平面之间的单值映射关系，这种变换都是两个线性函数之比，因 此称为双线性变换。双线性变换法与冲激响应不变法相比，其主要的优点是避免了频率响应的混叠现象，虽然在线性方面有些欠缺，但是可以通过频率的预畸来加以校正且计算比冲激响应不变法方便，实现起来比较容易，所以，本设计选择用双线性变换法设计巴特沃斯低通滤波器。1.3 数字滤波器设计流程图1.4 数字滤波器的设计步骤数字滤波器的设计步骤

10、：根据数字滤波器的技术指标先设计过渡模拟滤波器得到系统函数Ha（s）, 然后将 Ha（s）按某种方法（本实验采用双线性变换法）转换成数字滤波器的系统函数H（ z） 。具体为：（ 1）确定巴特沃斯数字低通滤波器的技术指标：通带边界频率 p, 阻带截止频率 s,通带最大衰减 p, 阻带最小衰减 s。（ 2）将数字滤波器的技术指标转换为模拟滤波器的技术指标。这里指 p 和 s 的变换而 p和 s 保持不变。本题采用双线性变换法，其转换公式为：2p p tanT2 2ss tanT23）根据技术指标 p、 s、p和 s 用下面公式求出滤波器的阶数。spksp10 s/10110 p/101lkg sp

11、g sp4）根据N由表 1.4 求出归一化极点pk和归一化低通原型系统函数Ga（p）通信系统综合设计与实践第 19 页表 1.4 巴特沃斯归一化低通滤波器参数分母因 式阶数B(p)=B1(p)B2(p) BN/2 (p) N/2 表示取大于等于N/2 的最小整数1( p2+1)2(p 2+1.4142p+1)3(p 2+p+1)(p+1)4(p 2+0.7654p+1)(p 2+1.8478p+1)5(p 2+0.6180p+1)(p 2+1.6180p+1)(p+1)6(p 2+0.5176p+1)(p 2+1.4142p+1)(p 2+1.9319p+1)7(p 2+0.4450p+1)(

12、p 2+1.2470p+1)(p 2+1.8019p+1)(p+1)8(p 2+0.3902p+1)(p 2+1.1111p+1)(p 2+1.6629p+1)(p 2+1.9619p+1)9(p 2+0.3473p+1)(p 2+p+1)(p 2+1.5321p+1)(p 2+1.8974p+1)(p+1)5)将Ga( p)去归一化，将p s 代入Ga( p) ，得到实际的滤波器系统函数：cHa(s) G(p) spc这里 c 为 3dB截止频率。( 6)用双线性变换法将模拟滤波器Ha(s)转换成数字低通滤波器系统函数H(z) 。转换公式为：H (z) Ha(s)21 z1s=T1 z2.

13、巴特沃斯低通数字滤波器技术指标的设置数字低通技术指标为 p=30Hz, p=1dBs=30dBs=60Hz采样频率为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) 将对仿真结果进行比较，从而检验滤波器滤波性

14、能的准确性。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

15、,x);subplot(312);plot(W,abs(H);title( ' 低通滤波器' );subplot(313);plot(t,y)title( '30Hz');matlab 滤波前后的信号波形变化如图3.1 所示：10输入信号0-1000.20.40.6低通滤波器0.8121000.511.522.533.530Hz20-200.20.40.60.81图 3.1 用 matlab滤波前后的信号波形变化3.2 波形图分析：由技术指标得：设计的巴特沃斯低通数字滤波器为30Hz以内的信号能通过，而高30Hz 的信号将通不过滤波器。因此，我们设计一个输入信号

16、，其为：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) 三个正弦波合成，从图中可知，输入信号的幅且在 1 秒内波形重复了20.25 次，即周期为1/20.25s ，其频率为20.25Hz, 幅度为1，和输

17、入信号中的sin(2* *20*t) 分量相差不大，但其频率稍微有些失真，说明了设计的低通滤波器达到了要求。4. 用 Simulink 实现巴特沃斯低通数字滤波器的仿真并分析4.1 mulink 简介Simulink 是 Matlab 最重要的组件之一，是 MATLAB中的一种可视化仿真工具，是一种基于MATLAB的框图设计环境，是实现动态系统建模、仿真和分析的一个软件包它提供一个动态系统建模、仿真和综合分析的集成环境。在该环境中，无需大量书写程序，而只需要通过简单直观的鼠标操作，就可构造出复杂的系统。Simulink 具有适应面广、结构和流程清晰及仿真精细、贴近实际、效率高、灵活等优点，并基

18、于以上优点Simulink 已被广泛应用于控制理论和数字信号处理的复杂仿真和设计。同时有大量的第三方软件和硬件可应用于或被要求应用于 Simulink 。4.2 用 Simulink 实现巴特沃斯低通数字滤波器的仿真用 Simulink 硬件仿真的电路图4.1 如下所示4.1 用 Simulink 硬件仿真的原理图4.3 波形图分析4.3.1 Simulink 波形图分析图中 Sine Wave 为正弦波，函数分别为sin(2* *20*t) 、 2*sin(2* *100*t) 和5*sin(2* *200*t) ；三个加号的为求和函数；Digital Filter Design 为巴特沃斯

19、低通数字滤波器，其设置如图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(如图

20、 4.5 所示)为滤波后的信号，其幅度为1.1 ， 在 5s到 5.5s 内波形重复了10 次， 故其周期为0.05s，频率为20Hz， 与输入信号中的sin(2* *20*t) 相当， 幅度变大了0.1 ， 说明有噪声干扰，但也基本达到了设计要求。4.2 Digital Filter Design图 4.3 输入信号的波形4.4 输入信号的波形通信系统综合设计与实践第 23 页4.5 滤波后的信号4.3.2 用 Matlab 与 Simulink 仿真结果的比较：输入信号的比较用 Simulink 出来的信号干扰较大，失真较严重，说明存在外界干扰，导致信号衰减，而用 matlab 出来的信号

21、较好，干扰可以忽略。巴特沃斯低通数字滤波器的比较由图 3.1 中的低通滤波器和图4.2 比较可得可知，用matlab设计的滤波器过渡带较窄， 而用 Simulink 设计的滤波器过渡带较大，仿真结果和matlab相比效果较差，用 matlab设计的滤波器与理论值差距较小，滤波效果更好。输出信号的比较用 matlab 滤波后的信号比理论信号在1 秒内大了1/4 个周期，效果不错，而用Simulink 滤波后的信号周期没变，幅度大了0.1，但滤波器控制不了幅度，幅度的变化说明存在外界干扰，说明噪声较大，但从滤波效果来看也达到了设计要求。综述由图 4.4和图 4.5可知，用Simulink 设计的滤

22、波器的输入信号与滤波后信号均存在5s 的延迟，说明了加法器的处理速度较慢，而用matlab 设计的滤波器不存在此现象。从滤波效果来看，用matlab 设计的滤波器效果较好，其虽存在频率误差，但影响不大，而用 Simulink 设计的滤波器频率虽不失真，但存在较大的外界干扰。个人认为，用 matlab设计的巴特沃斯低通数字滤波器滤波效果较好。对设计结果与理论值的比较，其均存在一定的误差，从侧面说明对于任何实验结果，在一定范围内的误差是允许的，其与理论值均有一定的误差，这是不可避免的，其中的原因多种多样，如实验设计不合理，外界因素的干扰，信号衰减 对于这些，我们只能通过不断改进设计原理、优化设计步

23、骤等措施来减小误差。同时，也告诉我们，任何理论结果要转换成实际应用，均要进行一定的优化，使其达到实际的应用要求。5. 总结与体会数字滤波器是一种对数字信号进行滤波处理以得到期望的响应特性的离散时间系统。它能够有效的对数字信号进行过滤，过滤多数字信号中的干扰和多余成分，能够保证通信过程的正常实现。本次实验设计选择了在通信系统中最常用，也是最基础的巴特沃斯低通数字滤波器。 根据给定的数字数字滤波器的技术指标，先设计过渡模拟滤波器得到系统函数Ha(s),然后将 Ha(s) 采用双线性变换法转换成数字滤波器的系统函数H( z) 。在这个变换的过程中避免了频率的混叠现象。然后根据数字滤波器的系统函数，编写程序，分别运用matlab 和 Simulink 仿真工具对设计的巴特沃斯低通数字滤波器进行了仿真。经过几次调试