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

1、-PAGE . z.通信系统综合设计与实践题目基于MATLAB的低通滤波器设计院系名称信院通信系专业名称通信工程学生学生*指导教师2013年5月25日目 录摘要21巴特沃斯低通数字滤波器简介3 1.1选择巴特沃斯低通滤波器及双线性变换法的原因4 1.2巴特沃斯低通滤波器的根本原理4 1.2.1巴特沃斯低通滤波器的根本原理4 1.2.2双线性变换法的原理 51.3数字滤波器设计流程图 7 1.4数字滤波器的设计步骤 72.巴特沃斯低通数字滤波器技术指标的设置83.用matlab实现巴特沃斯低通数字滤波器的仿真并分析9 3.1用matlab实现巴特沃斯低通数字滤波器的仿真9 3.2波形图分析 10

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

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

4、硬件电路仿真，观察仿真结果，并通过一个输入波形来验证仿真效果。6通过对用Matlab软件仿真和Simulink硬件仿真的比拟，说明两者的优缺点。7对实验结果进展最后的总结，写出自己的感想。关键字： matlab 低通滤波器 巴特沃斯 巴特沃斯低通数字滤波器简介1.1选择巴特沃斯低通滤波器及双线性变换法的原因1由于低通滤波器是组成其它滤波器的根底，应选用低通滤波器；2在当今社会，数字信号的应用越来越广泛，应选用数字信号；3巴特沃斯滤波器的特点是通频带的频率响应曲线最平滑并且应用围最广，应选巴特沃斯型滤波器；4为了不使数字滤波器在=附近产生频谱混叠，应选用双线性变换法。1.2巴特沃思低通滤波器的根

5、本原理：1.2.1巴特沃思低通滤波器的根本原理：巴特沃斯低通数字滤波器的幅度平方函数用下式表示式中，N称为滤波器的阶数。当=0时，；时，是3dB截止频率。在附近，随加大，幅度迅速下降。幅度特性与与N的关系如图1.1所示。幅度下降的速度与阶数N有关，N愈大，通带愈平坦，过渡带愈窄，过渡带与阻带幅度下降的速度愈快，总的频响特性与理想低通滤波器的误差愈小。图1.1 巴特沃斯低通数字滤波器 图1.2 三阶巴特沃斯滤波器极点幅度特性与与N的关系 分布图以替换，将幅度平方函数写成的函数复变量，此式表示幅度平方函数有2N个极点，极点用下式表示： =0,1,2,3. 2N个极点等间隔分布在半径为的圆上该圆称为

6、巴特沃斯圆，间隔为rad。例如N=3,极点间隔为/3rad，如图1.2所示。 为形成因果稳定的滤波器，2N个极点中只取平面左半平面的的N个极点构成Ha(s), 而右半平面的的N个极点构成Ha(-s),Has的表达式为为使设计公式和图表统一，将频率归一化。巴特沃斯低通数字滤波器采用对3dB截止频率归一化，归一化后的系统函数为令,称为归一化频率， 称为归一化复变量，这样，巴特沃斯低通原型系统函数为1.2.2双线性变换法原理双线性变换法是使数字滤波器的频率响应与模拟滤波器的频率响应相似的一种变换方法。为了克制多值映射的缺点，采用把整个s平面频率压缩方法，将整个频率轴上的频率围压缩到-/T/T之间，再

7、用转换到Z平面上。也就是说，第一步先将整个S平面压缩映射到S1平面的-/T/T一条横带里；第二步再通过标准变换关系将此横带变换到整个Z平面上去。这样就使S平面与Z平面建立了一一对应的单值关系，消除了多值变换性，也就消除了频谱混叠现象。映射关系如图1.3所示。设Has,经过非线性频率压缩后用，表示，这里用正切变换实现频率压缩：图1.3 双线性变换的映射关系式中，T为采样间隔，当从-/T经过0变化到/T时，由-经过0变化到+，实现了s平面上整个虚轴完全压缩到平面上虚轴的+/T之间的转换。即代入,得到再通过从平面转换到z平面，得到 上式是S平面与Z平面之间的单值映射关系，这种变换都是两个线性函数之比

8、，因此称为双线性变换。双线性变换法与冲激响应不变法相比，其主要的优点是防止了频率响应的混叠现象，虽然在线性方面有些欠缺，但是可以通过频率的预畸来加以校正且计算比冲激响应不变法方便，实现起来比拟容易，所以，本设计选择用双线性变换法设计巴特沃斯低通滤波器。1.3数字滤波器设计流程图数字滤波器技术指标指标参数变 换相应的模拟滤波器设计模拟滤波器离散化数字滤波器1.4数字滤波器的设计步骤数字滤波器的设计步骤：根据数字滤波器的技术指标先设计过渡模拟滤波器得到系统函数Ha(s),然后将Ha(s)按*种方法本实验采用双线性变换法转换成数字滤波器的系统函数Hz。具体为：1确定巴特沃斯数字低通滤波器的技术指标：

9、通带边界频率p,阻带截止频率s,通带最大衰减p,阻带最小衰减s。2将数字滤波器的技术指标转换为模拟滤波器的技术指标。这里指p和s的变换而p和s保持不变。此题采用双线性变换法，其转换公式为：3根据技术指标p、s、p和s用下面公式求出滤波器的阶数。根据N由表1.4求出归一化极点和归一化低通原型系统函数Ga(p)。表1.4 巴特沃斯归一化低通滤波器参数 分母因式阶数B(p)=B1(p)B2(p)BN/2(p) N/2表示取大于等于N/2的最小整数 1p2+1 2 (p2+1.4142p+1) 3 (p2+p+1)(p+1)4 (p2+0.7654p+1)(p2+1.8478p+1)5 (p2+0.6

10、180p+1)(p2+1.6180p+1)(p+1)6 (p2+0.5176p+1)(p2+1.4142p+1)(p2+1.9319p+1)7 (p2+0.4450p+1)(p2+1.2470p+1)(p2+1.8019p+1)(p+1)8 (p2+0.3902p+1)(p2+1.1111p+1)(p2+1.6629p+1)(p2+1.9619p+1)9 (p2+0.3473p+1)(p2+p+1)(p2+1.5321p+1)(p2+1.8974p+1)(p+1)5将Gap去归一化，将代入Gap，得到实际的滤波器系统函数：这里c为3dB截止频率。6用双线性变换法将模拟滤波器Ha(s)转换成数字

11、低通滤波器系统函数H(z)。转换公式为：s=巴特沃斯低通数字滤波器技术指标的设置数字低通技术指标为p=30Hz, p=1dBs=60Hz s=30dB采样频率为fs=500设计要求：输入*=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)将对仿真结果进展比拟，从而检验滤波

12、器滤波性能的准确性。3.用matlab实现巴特沃斯低通数字滤波器的仿真并分析3.1用matlab实现巴特沃斯低通数字滤波器的仿真Matlab程序如下：fs=500;t=0:1/fs:1;*=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,*);title(输入信号);N,wc=buttord(wp,ws,Rp,As);B,A=butter(N,wc);H,W=freqz(B,A);y=filter(B,A,*);subplot(312)

13、;plot(W,abs(H);title(低通滤波器);subplot(313);plot(t,y)title(30Hz);用matlab滤波前后的信号波形变化如图3.1所示：图3.1 用matlab滤波前后的信号波形变化3.2波形图分析：由技术指标得：设计的巴特沃斯低通数字滤波器为30Hz以的信号能通过，而高于30Hz的信号将通不过滤波器。因此，我们设计一个输入信号，其为：*=sin(2*20*t)+2*sin(2*100*t)+5*sin(2*200*t)。使其通过设计好的巴特沃斯低通数字滤波器，来验证设计好的巴特沃斯低通数字滤波器是否准确。如图3.1所示，输入信号为*=sin(2*20*

14、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，和输入信号中的sin(2*20*t)分量相差不大，但其频率稍微有些失真，说明了设计的低通滤波器到达了要求。4.用Simulink实现巴特沃斯低通数字滤波器的仿真并分析4.1Simuli

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

16、仿真用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为巴特沃斯低通数字滤波器，其设置如图4.2所示；而Scope为三个信号相加后的波形如图4.3、4.4所示，用来观察输入信号；Scope1为经过巴特沃斯低通数字滤波器后的波形如图4.5所示，用来观察滤波效果。Digital Filter Design的设置如图

17、三所示，其设置的为巴特沃斯低通数字滤波器，Scope为输入信号，由图4.3的，其幅度略大于五，其周期为0.01s，频率为100Hz，与输入信号*=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，说明有噪声干扰，但也根本到达了设计要求。图4.2 Digital Filter Design图4.3 输入信号

18、的波形图4.4 输入信号的波形图4.5 滤波后的信号4.3.2用Matlab与Simulink仿真结果的比拟： 输入信号的比拟用Simulink出来的信号干扰较大，失真较严重，说明存在外界干扰，导致信号衰减，而用matlab出来的信号较好，干扰可以忽略。巴特沃斯低通数字滤波器的比拟由图3.1中的低通滤波器和图4.2比拟可得可知，用matlab设计的滤波器过渡带较窄，而用Simulink设计的滤波器过渡带较大，仿真结果和matlab相比效果较差，用matlab设计的滤波器与理论值差距较小，滤波效果更好。输出信号的比拟用matlab滤波后的信号比理论信号在1秒大了1/4个周期，效果不错，而用Sim

19、ulink滤波后的信号周期没变，幅度大了0.1，但滤波器控制不了幅度，幅度的变化说明存在外界干扰，说明噪声较大，但从滤波效果来看也到达了设计要求。 综述由图4.4和图4.5可知，用Simulink设计的滤波器的输入信号与滤波后信号均存在5s的延迟，说明了加法器的处理速度较慢，而用matlab设计的滤波器不存在此现象。从滤波效果来看，用matlab设计的滤波器效果较好，其虽存在频率误差，但影响不大，而用Simulink设计的滤波器频率虽不失真，但存在较大的外界干扰。个人认为，用matlab设计的巴特沃斯低通数字滤波器滤波效果较好。对设计结果与理论值的比拟，其均存在一定的误差，从侧面说明对于任何实

20、验结果，在一定围的误差是允许的，其与理论值均有一定的误差，这是不可防止的，其中的原因多种多样，如实验设计不合理，外界因素的干扰，信号衰减对于这些，我们只能通过不断改良设计原理、优化设计步骤等措施来减小误差。同时，也告诉我们，任何理论结果要转换成实际应用，均要进展一定的优化，使其到达实际的应用要求。5.总结与体会 数字滤波器是一种对数字信号进展滤波处理以得到期望的响应特性的离散时间系统。它能够有效的对数字信号进展过滤，过滤多数字信号中的干扰和多余成分，能够保证通信过程的正常实现。本次实验设计选择了在通信系统中最常用，也是最根底的巴特沃斯低通数字滤波器。根据给定的数字数字滤波器的技术指标，先设计过渡模拟滤波器得到系统函数Ha(s),然后将Ha(s)采用双线性变换法转换成数字滤波器的系统函数Hz。在这个变换的过程中防止了频率的混叠现象。然后根据数字滤波器的系统函数，编写程序，分别运用matlab和Simulink仿真工具对设计的巴特沃斯低通数字滤波器进展了仿真。经过几次调试和对产生的仿真结果的分析，得到了在误差允许围之的低通数字滤波器。并且设计了验证实验。对一个由多个信号合成的数字