FIR带通滤波器设计成熟版.docVIP

目 录 1 技术要求 1 2 基本原理 1 2.1 FIR滤波器简介 1 2.2 窗函数法原理 2 3 建立模型描述 4 3.1 MATLAB常用函数 4 3.1.1 矩形窗函数 4 3.1.2 三角窗函数 5 3.1.3 广义余弦窗 5 3.1.4 汉宁窗（升余弦窗） 6 3.1.5 fir1函数 7 3.1.6 freqz函数 7 3.1.7 其他函数与命令 8 3.2 方案设计与论证 8 3.3 程序流程图 9 4 模块功能分析或源程序代码 10 5 调试过程及结论 13 5.1 实验结果 13 5.2 结果分析 15 6 思考题 15 7 心得体会 16 8 参考文献 17 FIR带通滤波器的设计 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滤波器的设计方法有多种，如窗函数法、频率采样法及其它各种优化设计方法，本窗函数法滤波器。数字滤波器可以理解为是一个计算程序或算法，将代表输入信号的数字时间序列转化为代表输出信号的数字时间序列，并在转化过程中，使信号按预定的形式变化。数字滤波器有多种分类，根据数字滤波器冲激响应的时域特征，可将数字滤波器分为两种，即无限长冲激响应（）滤波器和有限长冲激响应（）滤波器。数字滤波器具有无限宽的冲激响应，与模拟滤波器相匹配。所以滤波器的设计可以采取在模拟滤波器设计的基础上进一步变换的方法。数字滤波器的单位脉冲响应是有限长序列。它的设计问题实质上是确定能满足所要求的转移序列或脉冲响应的常数问题，设计方法主要有窗函数法、频率采样法和等波纹最佳逼近法等 因此设计FIR滤波器的方法之一可以从时域出发，截取有限长的一段冲击响应作为H(z)的系数，冲击响应长度N就是系统函数H(z)的阶数。只要N足够长，截取的方法合理，总能满足频域的要求。一般这种时域设计、频域检验的方法要反复几个回合才能成功。 要设计一个线性相位的FIR数字滤波器，首先要求理想频率响应。是w的周期函数，周期为，可以展开成傅氏级数： （公式1-1），一个理想的频率响应的傅立叶反变换： （公式1-2）是与理想频响对应的理想单位抽样响应序列。但不能用来作为设计FIR DF用的h(n)，因为一般都是无限长、非因果的，物理上无法实现。为了设计出频响类似于理想频响的滤波器，可以考虑用来近似。 窗函数的基本思想：先选取一个理想滤波器（它的单位抽样响应是非因果、无限长的），再截取（或加窗）它的单位抽样响应得到线性相位因果FIR滤波器。这种方法的重点是选择一个合适的窗函数和理想滤波器。 设是一个长序列，是长度为N的窗函数，用截断，得到N点序列，即 （公式1-3） 在频域上则有 （公式1-4） 由此可见，窗函数不仅仅会影响原信号在时域上的波形，而且也会影响到频域内的形状。 MATLAB信号工具箱主要提供了以下几种窗函数，如表1 表1 MATLAB窗函数 窗 窗 函 数 矩形窗 Boxcar 三角窗 Triang 海明窗 Hamming 汉宁窗 Hanning 布莱克曼 Blackman 切比雪夫窗 Chebyshev 凯塞窗 Kaiser