基于MATLAB的巴特沃斯滤波器

1数字信号处理课程设计2015年6月25日2目录一．设计目的：.........................................................................................................................................................................3二．设计要求：.........................................................................................................................................................................3三．设计内容：.........................................................................................................................................................................43.1选择巴特涡斯低通数据滤波器及双线性变换法的原因............................................................................................43.2巴特沃思低通滤波器的基本原理...............................................................................................................................43.3双线性变换法原理.......................................................................................................................................................53.4数字滤波器设计流程图...............................................................................................................................................73.5数字滤波器的设计步骤...............................................................................................................................................7四．用matlab实现巴特沃斯低通数字滤波器的仿真并分析.................................................................................................94.1巴特沃斯低通数字滤波器技术指标的设置................................................................................................................94.2用matlab实现巴特沃斯低通数字滤波器的仿真.......................................................................................................94.3波形图分析：.............................................................................................................................................................12五.总结与体会..........................................................................................................................................................................13六．附录参考文献.................................................................................................................................................................143一．设计目的：该课程设计是测控技术与仪器专业的必修课，开设课程设计的目的使学生掌握数字信号处理的基本概念和基本理论，能够利用辅助工具进行FIR和IIR数字滤波器的设计，进行一维信号的频谱分析，并进行仿真验证。加强实践教学环节，加强学生独立分析、解决问题的能力，培养学生动手能力和解决实际问题的能力，实现宽口径教育。(1)理解低通滤波器的过滤方法。(2)进一步熟悉低通滤波器的基本应用。(3)用仿真工具matlab软件对设计的滤波器进行软件和硬件仿真。(6)将对仿真结果进行比较，从而检验滤波器滤波性能的准确性。二．设计要求：地震发生时，除了会产生地震波，还会由地层岩石在断裂、碰撞过程中所发生的震动产生次声波。它的频率大约在每秒十赫兹到二十赫兹之间（可以用11Hz和15Hz的两个信号的和进行仿真，幅度可以分别设定为1、2）。大气对次声波的吸收系数很小，因此它可以传播的很远，而且穿透性很强。通过监测次声波信号可以监测地震的发生、强度等信息，因为自然界中广泛存在着各种次声波，这就对地震产生的次声波产生了干扰（可以用白噪声模拟，方差为5），需要采取一定的处理方法，才能检测到该信号，要求设计检测方案；并处理方法给出具体的软件（可以以51系列单片机、STM32F407、TMS320F28335或TMS320F6745为例）。假设地震次声波信号为x,输入x=sin(2*π*11*t)+2*sin(2*π*15*t)和伴有白噪声的合成信号，经过滤波器后滤除15Hz以上的分量，即只保留x=sin(2*π*11*t)+2*sin(2*π*15*t)的分量信号，来验证设计的滤波器是否达到了设计要求。4三．设计内容：3.1选择巴特涡斯低通数据滤波器及双线性变换法的原因（1）由于低通滤波器是组成其它滤波器的基础，故选用低通滤波器；（2）在当今社会，数字信号的应用越来越广泛，故选用数字信号；（3）巴特沃斯滤波器的特点是通频带的频率响应曲线最平滑并且应用范围最广，故选巴特沃斯型滤波器；（4）为了不使数字滤波器在ω=附近产生频谱混叠，故选用双线性变换法。3.2巴特沃思低通滤波器的基本原理巴特沃斯低通数字滤波器的幅度平方函数2ajH）（用下式表示Nac22)(11)j(H式中，N称为滤波器的阶数。当=0时，1jHa）（；c时，2/1jHa）（，c是3dB截止频率。在c附近，随加大，幅度迅速下降。幅度特性与与N的关系如图3.1所示。幅度下降的速度与阶数N有关，N愈大，通带愈平坦，过渡带愈窄，过渡带与阻带幅度下降的速度愈快，总的频响特性与理想低通滤波器的误差愈小。图3.1巴特沃斯低通数字滤波器图3.2三阶巴特沃斯滤波器极点5幅度特性与与N的关系分布图以ｓ替换j，将幅度平方函数2ajH）（写成ｓ的函数Naacjsss2)(11)(H)(H复变量js，此式表示幅度平方函数有2N个极点，极点ks用下式表示：)21221(21()1(NkjcNkecjs）（k=0,1,2,3….）2N个极点等间隔分布在半径为c的圆上（该圆称为巴特沃斯圆），间隔为N/rad。例如N=3,极点间隔为/3rad，如图3.2所示。为形成因果稳定的滤波器，2N个极点中只取ｓ平面左半平面的的N个极点构成Ha(s),而右半平面的的N个极点构成Ha(-s),Ha（s）的表达式为)()(Ha10kNkNcsss为使设计公式和图表统一，将频率归一化。巴特沃斯低通数字滤波器采用对3dB截止频率c归一化，归一化后的系统函数为1-0k)cs-(1)(HNkccssa令jp,c/,称为归一化频率，p称为归一化复变量，这样，巴特沃斯低通原型系统函数为1-0)-(1Ga(p)Nkkpp3.3双线性变换法原理双线性变换法是使数字滤波器的频率响应与模拟滤波器的频率响应相似的一种变换方法。为了克服多值映射的缺点，采用把整个s平面频率压缩方法，6将整个频率轴上的频率范围压缩到-π/T～π/T之间，再用sTeZ转换到Z平面上。也就是说，第一步先将整个S平面压缩映射到S1平面的-π/T～π/T一条横带里；第二步再通过标准变换关系sTeZ将此横带变换到整个Z平面上去。这样就使S平面与Z平面建立了一一对应的单值关系，消除了多值变换性，也就消除了频谱混叠现象。映射关系如图3.3所示。设Ha（s）,js,经过非线性频率压缩后用^1)(sHa，11js表示，这里用正切变换实现频率压缩：)21tan(21TT图3.3双线性变换的映射关系式中，T为采样间隔，当1从-π/T经过0变化到π/T时，由-∞经过0变化到+∞，实现了s平面上整个虚轴完全压缩到1s平面上虚轴的+π/T之间的转换。即TjTjTjTjtTjTjeeTeeeeTj11111111222/2/2/2/代入js,11js,得到TsTseeTs11112再通过Tsez1从1s平面转换到z平面，得到711112zzTssTsTz22上式是S平面与Z平面之间的单值映射关系，这种变换都是两个线性函数之比，因此称为双线性变换。双线性变换法与冲激响应不变法相比，其主要的优点是避免了频率响应的混叠现象，虽然在线性方面有些欠缺，但是可以通过频率的预畸来加以校正且计算比冲激响应不变法方便，实现起来比较容易，所以，本设计选择用双线性变换法设计巴特沃斯低通滤波器。3.4数字滤波器设计流程图3.5数字滤波器的设计步骤数字滤波器的设计步骤：根据数字滤波器的技术指标先设计过渡模拟滤波器得到系统函数Ha(s),然后将Ha(s)按某种方法（本实验采用双线性变换法）转换成数字滤波器的系统函数H（z）。具体为：（1）确定巴特沃斯数字低通滤波器的技术指标：通带边界频率ωp,阻带截止频率ωs,通带最大衰减аp,阻带最小衰减аs。（2）将数字滤波器的技术指标转换为模拟滤波器的技术指标。这里指ωp和ωs的变换而аp和аs保持不变。本题采用双线性变换法，其转换公式为：2tan2pTp数字滤波器技术指标指标参数变换相应的模拟滤波器设计模拟滤波器离散化数字滤波器82tan2sTs（3）根据技术指标Ωp、Ωs、ωp和ωs用下面公式求出滤波器的阶数。pssp11011010/10/psspkspgspglklN（4）根据N由表3.1求出归一化极点kp和归一化低通原型系统函数Ga(p)。表3.1巴特沃斯归一化低通滤波器参数分母因式阶数B(p)=B1(p)B2(p)…B[N/2](p)[N/2表示取大于等于N/2的最小整数]1（p2+1）2(p2+1