基数字信号处理课程设计-于MATLAB的巴特沃斯滤波器参考.doc

数字信号处理课程设计 2015年 6 月 25 日 目 录 一．设计目的： 2 二．设计要求： 3 三．设计内容： 3 3.1选择巴特涡斯低通数据滤波器及双线性变换法的原因 3 3.2巴特沃思低通滤波器的基本原理 4 3.3双线性变换法原理 5 3.4数字滤波器设计流程图 7 3.5数字滤波器的设计步骤 7 四．用matlab实现巴特沃斯低通数字滤波器的仿真并分析 8 4.1巴特沃斯低通数字滤波器技术指标的设置 8 4.2用matlab实现巴特沃斯低通数字滤波器的仿真 9 4.3波形图分析： 12 五.总结与体会 12 六．附录 参考文献 14 一．设计目的： 该课程设计是测控技术与仪器专业的必修课，开设课程设计的目的使学生掌握数字信号处理的基本概念和基本理论，能够利用辅助工具进行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)的分量信号，来验证设计的滤波器是否达到了设计要求。 三．设计内容： 3.1选择巴特涡斯低通数据滤波器及双线性变换法的原因 （1）由于低通滤波器是组成其它滤波器的基础，故选用低通滤波器； （2）在当今社会，数字信号的应用越来越广泛，故选用数字信号； （3）巴特沃斯滤波器的特点是通频带的频率响应曲线最平滑并且应用范围最广，故选巴特沃斯型滤波器； （4）为了不使数字滤波器在ω=附近产生频谱混叠，故选用双线性变换法。 3.2巴特沃思低通滤波器的基本原理 巴特沃斯低通数字滤波器的幅度平方函数用下式表示 式中，N称为滤波器的阶数。当=0时，；时，，是3dB截止频率。在附近，随加大，幅度迅速下降。幅度特性与与N的关系如图3.1所示。幅度下降的速度与阶数N有关，N愈大，通带愈平坦，过渡带愈窄，过渡带与阻带幅度下降的速度愈快，总的频响特性与理想低通滤波器的误差愈小。 图3.1 巴特沃斯低通数字滤波器 图3.2 三阶巴特沃斯滤波器极点 幅度特性与与N的关系 分布图 以ｓ替换，将幅度平方函数写成ｓ的函数 复变量，此式表示幅度平方函数有2N个极点，极点用下式表示： （=0,1,2,3….） 2N个极点等间隔分布在半径为的圆上（该圆称为巴特沃斯圆），间隔为rad。例如N=3,极点间隔为/3rad，如图3.2所示。 为形成因果稳定的滤波器，2N个极点中只取ｓ平面左半平面的的N个极点构成Ha(s), 而右半平面的的N个极点构成Ha(-s),Ha（s）的表达式为 为使设计公式和图表统一，将频率归一化。巴特沃斯低通数字滤波器采用对3dB截止频率归一化，归一化后的系统函数为 令,,称为归一化频率， 称为归一化复变量，这样，巴特沃斯低通原型系统函数为 3.3双线性变换法原理 双线性变换法是使数字滤波器的频率响应与模拟滤波器的频率响应相似的一种变换方法。为了克服多值映射的缺点，采用把整个s平面频率压缩方法， 将整个频率轴上的频率范围压缩到-π/T～π/T之间，再用转换到Z平面上。也就是说，第一步先将整个S平面压缩映射到S1平面的-π/T～π/T一条横带里；第二步再通过标准变换关系将此横带变换到整个Z平面上去。这样就使S平面与Z平面建立了一一对应的单值关系，消除了多值变换性，也就消除了频谱混叠现象。映射关系如图3.3所示。 设Ha（s）,,经过非