《课程设计（论文）-基于窗函数的FIR高通数字滤波器设计》由会员分享，可在线阅读，更多相关《课程设计（论文）-基于窗函数的FIR高通数字滤波器设计（13页珍藏版）》请在金锄头文库上搜索。
1、1基于窗函数的 FIR 高通数字滤波器设计摘要无限长脉冲数字滤波器的设计方法只考虑了幅度特性，没有考虑相位特性，所设的滤波器一般是某种确定的非线性相位特性。有限脉冲响应（FIR）滤波器在保证了幅度特性满足技术要求的同时，很容易做到有严格的线性相位特性。本课题利用 MATLAB 软件实现。MATLAB 是“矩阵实验室 ”（MATrix LABoratoy）的缩写，是一种科学计算软件，它使用方便，输入简捷，运算高效，内容丰富，因此利用 MATLAB 软件，通过一系列较为系统的函数法，根据已知的技术指标，就可以设计出满足要求的滤波器。关键字:MATLAB；窗函数；FIR 带阻数字滤波器；线性相位2目
2、录1. FIR 滤波器简介 .31.1 FIR 的特点 .32.2 线性相位 .42主要设计内容 .43.窗函数 .53.1 常用窗函数 .53.2 窗函数的指标 .84 应用窗函数法设计 FIR 数字滤波器的步骤 .94.1 数字高通滤波器的设计： .94.2 Matlab 仿真结果 .10总结 .12致谢 .13参考文献 .1431 FIR 滤波器简介数字滤波器是一种用来过滤时间离散信号的数字系统，通过对抽样数据进行数学处理来达到频域滤波的目的。根据其单位冲激响应函数的时域特性可分为两类：无限冲激响应（IIR）滤波器和有限冲激响应（FIR）滤波器 1。1.1 FIR 的特点FIR 滤波器的
3、主要优点为：系统总是稳定的，FIR 滤波器的系统函数可以表示为(2-1)易知，H(z) 在 Z 平面上有 N-1 个零点，z =0 是 N-1 阶极点，因此 FIR 系统总是稳定的（极点都在单位圆内） 。FIR 滤波器的优点之二：容易实现线性相位。当 FIR 系统的单位冲激响应满足 时，该系统具有线性相位。 （N 为奇数） （2-2）（N 为偶数） （2-3）FIR 滤波器的优点之三：允许设置多通带（或多阻带）滤波器。FIR 滤波器的优点之四：FIR 滤波器可以采用 FFT 方法实现其功能，从而大大提高效率。FIR 滤波器的缺点：由于 FIR 系统只有零点，因此这类系统不像 FIR 滤波器不像
4、 IIR 滤波器那样容易取得比较好的通带与阻带衰减特性。要取得较好的衰减特性，一般要求 H(z) 的阶次较高。综合起来看， FIR 滤波器具有 IIR 滤波器没有的许多特点，得到了越来越广泛的应用。FIR 滤波器的设计方法主要有三种：a.窗函数设计法；b.频率抽样发；c.最小平法抽样法；这里我主要讨论在 MATLAB 环境下通过调用信号分析与处理工具箱的几类窗函数来设计滤波器并分析与比较其性能2.2 线性相位1101)(10 )()()( NnnNnn zfhzhz )1()nhn2/)1()(/4一个单一频率的正弦信号通过一个系统，假设它通过这个系统的时间需要t，则这个信号的输出相位落后原来
5、信号 wt 的相位。从这边可以看出，一个正弦信号通过一个系统落后的相位等于它的 w*t；反过来说，如果一个频率为 w的正弦信号通过系统后，它的相位落后 delta，则该信号被延迟了 delta/w 的时间。在实际系统中，一个输入信号可以分解为多个正弦信号的叠加，为了使得输出信号不会产生相位失真，必须要求它所包含的这些正弦信号通过系统的时间是一样的。因此每一个正弦信号的相位分别落后，w1*t，w2*t，w3*t。因此，落后的相位正比于频率 w，如果超前，超前相位的大小也是正比于频率 w。从系统的频率响应来看，就是要求它的相频特性是一条直线。在 FIR 滤波器的设计中，为了得到线性相位的性质，通常
6、利用实偶对称序列的相频特性为常数 0和实奇对称序列为相频特性为常数 90 度的特点。因此得到的是对称序列，不是因果序列，是不可实现系统，为了称为物理可实现系统，需要将它向右移动半个周期，这就造成了相移特性随时间的变化，同时也是线性变化。单位脉冲响应 h(n)（为实数）具有偶对称或奇对称性，则 FIR 数字滤波器具有严格的线性相位特性。2 主要设计内容利用窗函数法、频率取样法及优化设计方法设计 FIR 滤波器，绘制出滤波器的特性图。利用所设计的滤波器对多个频带叠加的正弦信号进行处理，对比滤波前后的信号时域和频域图，验证滤波器的效果。基本思路：从时域出发设计 h(n)逼近理想 hd(n) 。设理想
7、滤波器的单位响应在时域表达为 hd(n),则 Hd(n) 一般是无限长的，且是非因果的，不能直接作为FIR 滤波器的单位脉冲响应。要想得到一个因果的有限长的滤波器单位抽样响应 h(n)，最直接的方法是先将 hd(n)往右平移，再进行截断，即截取为有限长因果序列：h(n)=hd(n)w(n)，并用合适的窗函数进行加权作为 FIR 滤波器的单位脉冲响应。按照线性相位滤波器的要求，线性相位 FIR 数字低通滤波器的单位抽样响应 h(n)必须是偶对称的。对称中心必须等于滤波器的延时常数，即用矩形窗设计的 FIR 低通滤波器，所设计滤波器的幅度函数在通带和阻带都呈现出振荡现象，且最大波纹大约为幅度的 9
8、%，返个现象称为吉布斯（Gibbs ）效应。5为了消除吉布斯效应，一般采用其他类型的窗函数。MATLAB 设计 FIR 滤波器有多种方法和对应的函数。窗函数设计法不仅在数字滤波器的设计中占有重要的地位，同时可以用于功率谱的估计，从根本上讲，使用窗函数的目的就是消除由无限序列的截短而引起的 Gibbs 现象所带来的影响 2。3. 窗函数在设计 FIR 数字滤波器时，窗函数的频谱应该满足：1 主瓣宽度尽可能的窄，以使过渡带尽量陡峭；2 最大旁瓣相对于主瓣尽可能的小，使能量尽可能集中于主瓣内，这样能够使得波动减小。3.1 常用窗函数1 汉宁（Hanning ）窗 3汉宁窗函数的时域形式表示为：(3-
9、1)频域形式为（3-2）汉宁窗函数的最大旁瓣值比主瓣值低 31dB，但是主瓣宽度比矩形窗函数的主瓣宽度增加了 1 倍，为 8/N。hanning 函数：生成汉宁窗调用方式：(1) w = hanning(n)：输入参数 n 是窗函数的长度；输出参数 w 是由窗函数的值组成的 n 阶向量。注意：此函数不返回是零点的窗函数的首尾两个元素。 1,20),12cos(5.0)( NnNnwL 21)()(25.0)(.)( NjRRR eWWj6(2) w = hanning(n,symmetric)：与上面相类似。(3) w = hanning(n,periodic)：此函数返回包括为零点的窗函数的
10、首尾两个元素。图 1 汉宁窗及其频谱特性注释：汉宁窗又被称为升余弦窗，汉宁窗可以看做三个矩形时间窗的频谱之和，汉宁窗主瓣加宽并降低，旁瓣显著减小，分辨率下降2 汉明（Hamming）窗:函数的时域形式可以表示为 （3-3）12cos46.05)(NkkwN,2L频域形式为： （3-4）123.)(. WWRRR其中， 为矩形窗函数的幅度频率特性函数。)(R汉明窗函数的最大旁瓣值比主瓣值低 41dB，但它和汉宁窗函数的主瓣宽度是一样大的。Hamming 函数：生成海明窗调用方式(1) w = hamming(n)：输入参数 n 是窗函数的长度；输出参数 w 是由窗函数的7值组成的 n 阶向量。(
11、2) w = hamming(n,sflag)：参数 sflag 用来控制窗函数首尾的两个元素值；其取值为 symmetric 或 periodic；默认值为 symmetric。图 2 汉明窗及其频谱特性注释：汉明窗和汉宁窗都是余弦窗，只是加权系数不同，汉明窗旁瓣更小3 布莱克曼窗函数的时域形式可以表示为(3-5)14cos08.12cos5.04)( NkNkkw Nk,21L它的频域特性为（3-6）)(W42.0R5. 1212WRR1404.WRR其中， 为矩形窗函数的幅度频率特性函数。)(布莱克曼窗函数的最大旁瓣值比主瓣值低 57dB，但是主瓣宽度是矩形窗函数的主瓣宽度的 3 倍，为 12/N。Blackman 函数：生成海明窗调用方式(1) w = blackman (n)：输入参数 n 是窗函数的长度；输出参数 w 是由窗函数的值组成的 n 阶向量。8(2) w = blackman (n,sflag)：参数 sflag 用来控制窗函数首尾的两个元素值；其取值为 symmetric 或 pe
