《窗函数设计低通滤波器》由会员分享，可在线阅读，更多相关《窗函数设计低通滤波器（9页珍藏版）》请在金锄头文库上搜索。
1、 目录目录1 1 课题描述课题描述.1.1 2 2 设计原理设计原理.1.12.12.1 滤波器的分类滤波器的分类.1.12.22.2 利用窗函数设计利用窗函数设计 FIRFIR 滤波器滤波器.2.22.2.12.2.1 窗函数的设计原理窗函数的设计原理.3.32.2.22.2.2 典型窗函数的介绍典型窗函数的介绍.4.42.2.32.2.3 用窗函数设计用窗函数设计 FIRFIR 滤波器的步骤滤波器的步骤.5.52.2.42.2.4 窗函数法的窗函数法的 MATLABMATLAB 设计函数简介设计函数简介.5.5 3.3. 设计内容设计内容.6.63.13.1 设计题目设计题目. 6 63.
2、23.2 用用 MATLABMATLAB 编程编程. 7 73.33.3 设计结果分析设计结果分析. 8 84.4. 总结总结.9.91 1 课题描述课题描述数字滤波器是数字信号处理的重要工具之一，它通过数值运算处理改变输入信号所含 频率成分的相对比例或者滤出某些频率成分的数字器件或程序，而数字滤波器处理精度高， 体积小，稳定，重量轻，灵活，不存在阻抗匹配问题，可以实现模拟滤波器无法实现的特 殊功能。故本课题使用 MATLAB 信号处理箱和运用窗函数设计数字滤波器。 2 2 设计原理设计原理2.12.1 滤波器的分类滤波器的分类（1）从功能上分：低通滤波器，高通滤波器，带通滤波器，带阻滤波器。
3、（2）从实现方法上分：FIR 滤波器,IIR 滤波器。（3）从设计方法上分：Chebyshev(切比雪夫)滤波器，Butterworth（巴特沃斯)滤波 器，Ellipse(椭圆)滤波器，Bessel（贝塞尔）滤波器。（4）从处理信号上分：经典滤波器，现代滤波器。2.22.2 利用窗函数设计利用窗函数设计 FIRFIR 滤波器滤波器2.2.12.2.1 窗函数的设计原理窗函数的设计原理设希望逼近的滤波器频率响应函数为 H(w)，其单位脉冲响应是. nhdnjwn djw denheH)()(由已知的求出，经过 Z 变换可以得到滤波器的系统函数。通常以理想)(jw deH)(nhd滤波器作为，其
4、幅频特性逐段恒定，在边界频率处有不连续的点，因而是)(jw deH)(nhd无限时宽的，且是非因果序列的某为了构造一个长度为 N 的线性相位滤波器，只有将截取一段，设截取的那段用 h(n)表示即：)(nhdh(n)= )()(nRnhnd式中是一个矩形序列，长度为 N 。我们实际设计的滤波器的单位脉冲响应为 h(n)。)(nRn长度为 N，其系统函数为 H(Z)= 。这样用一个有限长的序列 h(n)去代替，10)(Nnnznh)(nhd肯定会引起误差，表现在频域就是通常所说的吉布斯效应。该效应引起过度加宽以及通带dweeHnhjwnjwwcwcdd)(21)( 和阻带内的波动，尤其使阻带的衰减
5、小，从而满足不了技术上的要求。吉布斯效应是由于将直接截断引起的，称为截断效应，窗函数法设计 FIE 滤波器就是构造一个窗函数)(nhdw(n)来减少截断效应。 2.2.22.2.2 典型窗函数的介绍典型窗函数的介绍（1）矩形窗： 幅度函数 )()(nRnwNR)2/sin()2/sin()(wwNwwRg（2）三角窗： 12)(NnnwB) 1(210Nn122)(NnnwB1) 1(21NnN幅度函数 22_2/sin()4/sin(2ewwN NwBg频谱函数 wNjjw BewwN NeW21 2)2/sin()4/sin(2)(（3） 汉宁窗 ： )()12cos(1 5 . 0)(n
6、RNnnWNHnwNiHngHniw HnewWnwFTeW21 )()()()2()2(25. 0)(5 . 0)( N1-N1NNwWNwWwWnWRgRgRgHng 时，当（4）哈明窗： )()12cos(46. 054. 0)(nRNnnwNHm频谱函数 )(23. 023. 0-)(54. 0)()12()12(NwjRNwjRiw Wjw HmeWeWeWeW ）（幅度函数 )2(23. 0)2(23. 0)(54. 0)(NwWNwWwWwWRgRgRgHmg（5）贝塞尔窗： 10)()()(00NnIInWk21121）（Nn210)2(!11)(kkkI幅度函数 )1(1)c
7、os()(2)0()(NnkkkgwnnwwwW六种窗函数的基本参数六种窗函数的基本参数过渡带宽度tB窗函数类型旁瓣峰值an/dB 近似值精确值阻带最小衰 减 as/dB矩形窗-13N/4N/8 . 1-21三角窗-25N/8N/1 . 6-25汉宁窗-31N/8N/2 . 6-44哈明窗-41N/8N/6 . 6-53布莱克曼窗-57N/10N/11-74凯赛窗-57N/10-802.2.32.2.3 用窗函数设计用窗函数设计 FIRFIR 滤波器的步骤滤波器的步骤 用用窗函数设计 FIR 滤波器的步骤如下：（1）根据对阻带衰减以及过渡带的指标要求，选择窗函数的类型，并估计窗口长度。N 先按
8、照阻带衰减函数类型。原则是在保证阻带衰减满足要求的情况下，尽量的选择主瓣窄的窗函数。再根据过渡带的宽度估计窗口长度。带球滤波器的过渡带宽近似等于窗函数tB主瓣宽度，且近似与窗口长度成反比,A 取决于窗口类型，例如，矩形窗的NBtAN/ ，哈明窗的等，参数 A 的近似和精确取值参考表所示。4A8A(2)构造希望逼近的频率响应函数即: jwedH 2/1 -j dgj de e）（NHH所谓的“标准窗函数法”，就是选择为线性相位理想滤波器（理想低通，理想高 jwedH通，理想带通，理想带阻）。以低通滤波器为例，应满足： dgH 1wdgHc | 0dgH|c理想滤波器的截止频率近似为最终设计的的过
9、载的中心频率点，幅度函数衰减cFIRDF一半。所以如果设计指标给定带通边界频率和阻带边界频率阻带边界频率和，一般pn取：2sp c(3) 计算，如果给出多的待求滤波器的频响函数，那么单位脉冲响应应用下nhd jwedH式求出: dee21nhnj-j ddH如果较复杂，或者不能用封闭公式表示，这不能够用上式求出。我们可以对 jwedHnhd从到采样 M 点，采样值为进行点,02 k2jddekM MHH MIFFTIDFT得到： MMMHIDFTknhdd根据频域采样理论，与应满足如下关系： nhdM ndh -rMddnrnhnhRMM因此，如果值选的较大，可以保证窗口内有效逼近。M nhd
10、M ndh(4)加窗得到设计结果：。)()()(nwnhnhd2.2.42.2.4 窗函数法的窗函数法的 MATLABMATLAB 设计函数简介设计函数简介 可以调用 MATLAB 工具箱函数 fir1 实现本实验所要求的线性相位 FIR-DF 的设计，调用一维快速傅立叶变换函数 fft 来计算滤波器的频率响应函数。fir1 是用窗函数法设计线性相位 FIRDF 的工具箱函数，调用格式如下：hn=fir1(N, wc, ftype, window)fir1 实现线性相位 FIR 滤波器的标准窗函数法设计。hn=fir1(N,wc)可得到 6 dB 截止频率为 wc 的 N 阶(单位脉冲响应 h
11、(n)长度为 N+1)FIR 低通滤波器，默认(缺省参数 windows)选用 hammiing 窗。其单位脉冲响应 h(n)满足线性相位条件：h(n)=h(N-1-n)其中 wc 为对 归一化的数字频率，0wc1。当 wc=wc1, wc2时，得到的是带通滤波器。hn=fir1(N,wc,ftype)可设计高通和带阻滤波器。当 ftype=high 时，设计高通 FIR 滤波器；当 ftype=stop 时，设计带阻 FIR 滤波器。应当注意，在设计高通和带阻滤波器时，阶数 N 只能取偶数(h(n)长度 N+1 为奇数)。不过，当用户将 N 设置为奇数时，fir1 会自动对 N 加 1。hn=fir1(N,wc,window)可以指定窗函数向量 window。如果缺省 window 参数，则 fir1 默认为 hamming 窗。可用的其他窗函数有 Boxcar, Hanning, Bartlett, Blackman, Kaiser 和Chebwin 窗。例如：hn=fir1(N,wc,bartlett(N+1)使用 Bartlett 窗设计；hn=fir1(N,wc,chebwin(N+1,R)使用 C
