FIR带通滤泼器的设计.docVIP

《数字信号处理课程设计报告》 题 目： FIR带通滤泼器的设计 学 院： 沈阳工业大学工程学院 专 业： 测控技术与仪器 班 级： 0901 姓 名： 指导教师： 2012年 6月24日 目 录 1 课程设计要求……………………………………………3 2 基本原理…………………………………………………3 2.1 FIR滤泼器简介………………………………………3 2.2 窗函数法原理…………………………………………4 3 建立模型概述………………………………………………5 3.1 MATLB 常用函数……………………………………5 3.1.1 窗函数………………………………………………5 3.1.2 fir1函数……………………………………………6 3.1.3 freqz函数……………………………………………7 3.1.4 其他函数与命令……………………………………7 3.2 程序流程图……………………………………………8 4 模块功能分析………………………………………………9 5调试过程及结论……………………………………………13 6 心得体会……………………………………………………17 7 参考文献……………………………………………………17 1 课程设计要求 用窗函数法设计FIR带通滤波器。要求低端阻带截止频率ωls=0.2π，低端通带截止频率ωlp=0.35π，高端通带截止频率ωhp=0.65π，高端阻带截止频率ωhs=0.8π。绘出h(n)及其幅频响应特性曲线。 2 基本原理 2.1 FIR滤泼器简介 数字滤波器包括FIR（有限单位脉冲响应）滤波器与IIR（无限单位脉冲响应）滤波器两种。在现代信号处理技术中例如数据传输、雷达接收以及一些要求较高的电子系统都越来越多地要求信道具有线性的相位特性在这方面FIR滤波器具有独到的优点它可以在幅度特性随意设计的同时保证精确、严格的线性相位特性FIR滤波器的单位脉冲响应h（n）是有限长的（0≤n≤N-1），其z变换为的（N-1）阶多项式： 可得FIR滤波器的系统差分方程为： 因此，FIR滤波器又称为卷积滤波器。FIR滤波器的频率响应表达式为： 信号通过FIR滤波器不失真条件是在通带内具有恒定的幅频特性和线性相位特性。理论上可以证明：当FIR滤波器的系数满足下列中心对称条件： 或者 时，滤波器设计在逼近平直幅频特性的同时，还能获得严格的线性相位特性。线性相位FIR滤波器的相位滞后和群延迟在整个频带上是相等且不变的。对于一个 N 阶的线性相位FIR滤波器，群延迟为常数，即滤波后的信号简单地延迟常数个时间步长。这一特性使通带频率内信号通过滤波器后仍保持原有波形形状而无相位失真。 FIR滤波器的设计任务是选择有限长度的h(n)使传输函数满足技术要求。FIR滤波器的设计方法有多种，如窗函数法、频率采样法及其它各种优化设计方法，本窗函数法滤波器。 ，则其对应的单位脉冲响应为： 窗函数设计法的基本原理是用有限长单位脉冲响应序列h(n)逼近hd(n)。由于hd(n)往往是无限长序列，且是非因果的，变成物理上可实现呢？一个自然的想法是只取其中的某些项，即只截取hd(n)中的一部分，比如n=0，…，N-1，N为正整数。这种处理相当于将hd(n)与函数w(n)相乘，w(n)具有下列形式： w(n)相当于一个矩形，我们称之为矩形窗。即我们可采用矩形窗函数w(n)将无限脉冲响应hd(n)截取一段h(n)来近似为hd(n)，这种截取在数学上表示为： 这里应该强调的是，加窗函数不是可有可无的，而是将设计变为物理可实现所必须的。h(n)就作为实际设计的FIR数字滤波器的单位脉冲响应序列，其频率响应函数Hd(ejω)为： 式中，N为所选窗函数w(n)的长度。 经过加矩形窗后所得的滤波器实际频率响应能否很好地逼近理想频率响应呢？图1给出了理想滤波器加矩形窗后的情况。理想低通滤波器的频率响应如图中左上角图，矩形窗的频率响应为左下角图。时间域内的乘积要求实际频率响应为这两个频率响应函数在频域内的卷积（卷积定理），即得到图形为图中右边的结果。 图1 FIR滤波器理想与实际频率响应 用窗函数法设计的滤波器性能取决于窗函数w(n)的类型及窗口长度N的取值。因此，这种方法的重点在于选择某种恰当的窗函数和一种合适的理想滤波器。 3 建立模型概述 3.1 MATLB常用函数 3.1.1 窗函数 矩形窗：w=boxcar(N) 海明（