《[2017年整理]关于滤波器的窗函数的应用》由会员分享，可在线阅读，更多相关《[2017年整理]关于滤波器的窗函数的应用（11页珍藏版）》请在金锄头文库上搜索。
1、第5章 窗 函 数本章对数字信号处理过程中所用到的有关窗函数进行了说明，这些窗函数主要有：矩形窗函数、三角窗函数、巴特利特窗函数、广义余弦窗、汉宁窗函数、海明窗函数、布莱克曼窗函数、凯塞窗函数和切比雪夫窗函数等。本章的学习目标： 理解窗函数的基本概念以及加窗原理 掌握基本窗函数产生的方法 掌握广义余弦窗函数产生的各种方法 了解凯塞窗函数和切比雪夫窗函数产生的方法5.1 窗函数基本概念在实际进行数字信号处理时，往往需要把信号的观察时间限制在一定的时间间隔内，只需要选择一段时间信号对其进行分析。这样，取用有限个数据，即将信号数据截断的过程，就等于将信号进行加窗函数操作。而这样操作以后，常常会发生频
2、谱分量从其正常频谱扩展开来的现象，即所谓的“频谱泄漏”。当进行离散傅立叶变换时，时域中的截断是必需的，因此泄漏效应也是离散傅立叶变换所固有的，必须进行抑制。而要对频谱泄漏进行抑制，可以通过窗函数加权抑制 DFT 的等效滤波器的振幅特性的副瓣，或用窗函数加权使有限长度的输入信号周期延拓后在边界上尽量减少不连续程度的方法实现。而在后面的FIR 滤波器的设计中，为获得有限长单位取样响应，需要用窗函数截断无限长单位取样响应序列。另外，在功率谱估计中也要遇到窗函数加权问题。由此可见，窗函数加权技术在数字信号处理中的重要地位。下面介绍窗函数的基本概念。设 x(n)是一个长序列，w(n )是长度为 N 的窗
3、函数，用w(n)截断 x(n)，得到 N 点序列 xn(n)，即xn(n) = x(n) w(n)在频域上则有 jjj de21eWXNMATLAB 在数字信号处理中的应用 154 由此可见，窗函数 w(n)不仅仅会影响原信号 x(n)在时域上的波形，而且也会影响到频域内的形状。MATLAB 信号工具箱主要提供了以下几种窗函数，如表 5-1 所示。表 5-1 MATLAB 窗函数窗 窗 函 数矩形窗 Boxcar巴特利特窗 Barlett三角窗 Triang布莱克曼窗 Blackman海明窗 Hamming汉宁窗 Hanning凯塞窗 Kaiser切比雪夫窗 Chebwin5.2 基本窗函数数
4、字信号处理领域中所用到的基本窗函数主要有：矩形窗函数、三角窗函数和巴特利特窗函数。下面就对这些窗函数展开介绍。5.2.1 矩形窗函数矩形窗(Rectangular Window)函数的时域形式可以表示为： 其 他,01)(NnRnwN它的频域特性为 2sine21jj WNR第 5 章 窗 函 数 155 Boxcar 函数：生成矩形窗调用方式w = boxcar (n)：输入参数 n 是窗函数的长度；输出参数 w 是由窗函数的值组成的 n 阶向量。从功能上讲，该函数又等价于 w = ones(n,1)。应用说明例：生成一个长度为 50 的矩形窗，并观察其频率特性(使用归一化的幅值和频率)，结
5、果如图 5-1 所示。程序n=50;window=boxcar(n);h,w=freqz(window,1);subplot(2,1,1)stem(window);subplot(2,1,2);plot(w/pi,20*log(abs(h)/abs(h(1);5.2.2 三角窗函数三角窗是最简单的频谱函数 为非负的一种窗函数。三角窗函数的时域形式可以)W(ej表示为：当 n 为奇数时 nknkw21,)()当 n 为偶数时 nknkw2,)1()图 5-1 矩形窗及其频谱特性MATLAB 在数字信号处理中的应用 156 它的频域特性为： 221jj sin41eNWNR三角窗函数的主瓣宽度为
6、，比矩形窗函数的主瓣宽度增加了一倍，但是它的旁瓣8宽 度却小得多。Triang 函数：生成三角窗调用方式w = triang(n)：输入参数 n 是窗函数的长度；输出参数 w 是由窗函数的值组成的 n 阶向量。三角窗也是两个矩形窗的卷积。三角窗函数的首尾两个数值通常是不为零的。当 n 是偶数时，三角窗的傅立叶变换总是非负数。应用说明例：生成一个长度为 50 的三角窗，并观察其频率特性(使用归一化的幅值和频率) ，结果如图 5-2 所示。程序n=31;window=triang(n);h,w=freqz(window,1);subplot(1,2,1)stem(window);subplot(1
7、,2,2);plot(w/pi,20*log(abs(h)/abs(h(1);5.2.3 巴特利特窗函数巴特利特窗函数的时域形式可以表示为当 n 为奇数时图 5-2 三角窗及其频谱特性第 5 章 窗 函 数 157 nknkw21,)()(当 n 为偶数时 nknkw2,1)(Bartlett 函数：生成巴特利特窗调用方式w = bartlett(n)：(1) 输入参数 n 是窗函数的长度；(2) 输出参数 w 是由窗函数的值组成的 n 阶向量。(3) 巴特利特窗是两个矩形窗的卷积。(4) 此函数生成的巴特利特窗的两端总是 0。(5) 该函数与三角窗函数非常类似。当 n 是奇数时，则有 bar
8、tlett(n+2)的中间部分等于triang(n)；当 n 是偶数时，巴特利特窗的傅立叶变换是负数。应用说明例：生成一个长度为 50 的巴特利特窗，并观察其频率特性(使用归一化的幅值和频率) ，结果如图 5-3 所示。程序n=51;window=bartlett(n);h,w=freqz(window,1);subplot(1,2,1)stem(window);subplot(1,2,2);plot(w/pi,20*log(abs(h)/abs(h(1);5.3 广义余弦窗汉宁窗、海明窗和布莱克曼窗，都可以用一种通用的形式表示，这就是广义余弦窗。图 5-3 巴特利特窗及其频谱特性MATLAB
9、 在数字信号处理中的应用 158 这些窗都是广义余弦窗的特例，汉宁窗又被称为余弦平方窗或升余弦窗，海明窗又被称为改进的升余弦窗，而布莱克曼窗又被称为二阶升余弦窗。采用这些窗可以有效地降低旁瓣的高度，但是同时会增加主瓣的宽度。这些窗都是频率为 0、2/(N 1)和 4/(N1)的余弦曲线的合成，其中 为窗的长度。N通常采用下面的命令来生成这些窗： 12:0ind ind)*cos(Ccs(ind)BAwo其中，A、B、C 适用于自己定义的常数。根据它们取值的不同，可以形成不同的窗函数，分别是： 汉宁窗 A=0.5，B=0.5，C=0； 海明窗 A=0.54，B=0.54，C=0； 布莱克曼窗 A
10、=0.5，B=0.5，C=0.08；5.3.1 汉宁窗函数汉宁窗函数的时域形式可以表示为：12cos5.0)(nkkwNk,21它的频域特性为： 21je125.0. NRRRWNW其中， 为矩形窗函数的幅度频率特性函数。)(R汉宁窗函数的最大旁瓣值比主瓣值低 31dB，但是主瓣宽度比矩形窗函数的主瓣宽度增加了 1 倍，为 8/N。hanning 函数：生成汉宁窗调用方式(1) w = hanning(n)：输入参数 n 是窗函数的长度；输出参数 w 是由窗函数的值组成的n 阶向量。注意：此函数不返回是零点的窗函数的首尾两个元素。 (2) w = hanning(n,symmetric)：与上
11、面相类似。图 5-4 汉宁窗及其频谱特性第 5 章 窗 函 数 159 (3) w = hanning(n,periodic)：此函数返回包括为零点的窗函数的首尾两个元素。应用说明例：生成一个长度为 50 的汉宁窗，并观察其频率特性(使用归一化的幅值和频率)，结果如图 5-4 所示。程序n=51;window=hanning(n);h,w=freqz(window,1);subplot(1,2,1)stem(window);subplot(1,2,2);plot(w/pi,20*log(abs(h)/abs(h(1);5.3.2 海明窗函数海明窗函数的时域形式可以表示为12cos46.05)(
12、NkkwNk,21它的频域特性为 1123.0)(54.)( WWRRR其中， 为矩形窗函数的幅度频率特性函数。R海明窗函数的最大旁瓣值比主瓣值低 41dB，但它和汉宁窗函数的主瓣宽度是一样大的。Hamming 函数：生成海明窗调用方式(1) w = hamming(n)：输入参数 n 是窗函数的长度；输出参数 w 是由窗函数的值组成的n 阶向量。(2) w = hamming(n,sflag)：参数 sflag 用来控制窗函数首尾的两个元素值；其取值为symmetric 或 periodic；默认值为 symmetric。应用说明例：生成一个长度为 50 的海明窗，并观察图 5-5 海明窗及
13、其频谱特性MATLAB 在数字信号处理中的应用 160 其频率特性(使用归一化的幅值和频率 )，结果如图 5-5 所示。程序n=51;window=hamming(n);h,w=freqz(window,1);subplot(1,2,1)stem(window);subplot(1,2,2);plot(w/pi,20*log(abs(h)/abs(h(1);5.3.3 布莱克曼窗函数布莱克曼窗函数的时域形式可以表示为14cos08.12cos5.04)( NkNkkw Nk,21它的频域特性为 )(W42.0R5. 112WRR4. N其中， 为矩形窗函数的幅度频率特性函数。)(R布莱克曼窗函
14、数的最大旁瓣值比主瓣值低 57dB，但是主瓣宽度是矩形窗函数的主瓣宽度的 3 倍，为 12/N。Blackman 函数：生成海明窗调用方式(1) w = blackman (n)：输入参数 n 是窗函数的长度；输出参数 w 是由窗函数的值组成的 n 阶向量。(2) w = blackman (n,sflag)：参数 sflag 用来控制窗函数首尾的两个元素值；其取值为symmetric 或 periodic；默认值为 symmetric。应用说明例：生成一个长度为 50 的布莱克曼窗，并观察其频率特性(使用归一化的幅值和频率 )，结果如图 5-6 所示。图 5-6 布莱克曼窗及其频谱特性第 5
15、 章 窗 函 数 161 程序n=51;window=blackman(n);h,w=freqz(window,1);subplot(1,2,1)stem(window);subplot(1,2,2);plot(w/pi,20*log(abs(h)/abs(h(1);5.4 凯 塞 窗上面所讨论的几种窗函数，在获得旁瓣抑制的同时却增加了主瓣的宽度。而凯塞窗定义了一组可调的窗函数，它是由零阶贝塞尔函数构成的，其主瓣能量和旁瓣能量的比例是近乎最大的。而且，这种窗函数可以在主瓣宽度和旁瓣高度之间自由选择它们的比重，使用户的设计变得非常灵活。凯塞窗函数的时域形式可表示为)(12)(0INkkw10Nk其中， 是第 1 类变形零阶贝塞尔函数， 是窗函数的形状参数，由下式确定：)(0I21, 50),(786.)21(548.,.40其中， 为凯塞窗函数的主瓣值和旁瓣
