关于滤波器的窗函数的应用.docVIP

第5章 窗 函 数 本章对数字信号处理过程中所用到的有关窗函数进行了说明，这些窗函数主要有：矩形窗函数、三角窗函数、巴特利特窗函数、广义余弦窗、汉宁窗函数、海明窗函数、布莱克曼窗函数、凯塞窗函数和切比雪夫窗函数等。 本章的学习目标： 理解窗函数的基本概念以及加窗原理 掌握基本窗函数产生的方法 掌握广义余弦窗函数产生的各种方法 了解凯塞窗函数和切比雪夫窗函数产生的方法 5.1 窗函数基本概念 在实际进行数字信号处理时，往往需要把信号的观察时间限制在一定的时间间隔内，只需要选择一段时间信号对其进行分析。这样，取用有限个数据，即将信号数据截断的过程，就等于将信号进行加窗函数操作。而这样操作以后，常常会发生频谱分量从其正常频谱扩展开来的现象，即所谓的“频谱泄漏”。当进行离散傅立叶变换时，时域中的截断是必需的，因此泄漏效应也是离散傅立叶变换所固有的，必须进行抑制。而要对频谱泄漏进行抑制，可以通过窗函数加权抑制DFT的等效滤波器的振幅特性的副瓣，或用窗函数加权使有限长度的输入信号周期延拓后在边界上尽量减少不连续程度的方法实现。而在后面的FIR滤波器的设计中，为获得有限长单位取样响应，需要用窗函数截断无限长单位取样响应序列。另外，在功率谱估计中也要遇到窗函数加权问题。由此可见，窗函数加权技术在数字信号处理中的重要地位。 下面介绍窗函数的基本概念。设x(n)是一个长序列，w(n)是长度为N的窗函数，用w(n)截断x(n)，得到N点序列xn(n)，即 xn(n) = x(n) w(n) 在频域上则有 由此可见，窗函数w(n)不仅仅会影响原信号x(n)在时域上的波形，而且也会影响到频域内的形状。 MATLAB信号工具箱主要提供了以下几种窗函数，如表5-1所示。 表5-1 MATLAB窗函数 窗 窗 函 数 矩形窗 Boxcar 巴特利特窗 Barlett 三角窗 Triang 布莱克曼窗 Blackman 海明窗 Hamming 汉宁窗 Hanning 凯塞窗 Kaiser 切比雪夫窗 Chebwin 5.2 基本窗函数 数字信号处理领域中所用到的基本窗函数主要有：矩形窗函数、三角窗函数和巴特利特窗函数。下面就对这些窗函数展开介绍。 5.2.1 矩形窗函数 矩形窗(Rectangular Window)函数的时域形式可以表示为： 它的频域特性为 Boxcar函数：生成矩形窗 调用方式 w = boxcar (n)：输入参数n是窗函数的长度；输出参数w是由窗函数的值组成的n阶向量。 从功能上讲，该函数又等价于w = ones(n,1)。 应用说明 例：生成一个长度为50的矩形窗，并观察其频率特性(使用归一化的幅值和频率)，结果如图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 三角窗函数 三角窗是最简单的频谱函数为非负的一种窗函数。三角窗函数的时域形式可以表示为： 当n为奇数时 当n为偶数时 它的频域特性为： 三角窗函数的主瓣宽度为，比矩形窗函数的主瓣宽度增加了一倍，但是它的旁瓣宽度却小得多。 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,2,2); plot(w/pi,20*log(abs(h)/abs(h(1)))); 5.2.3 巴特利特窗函数 巴特利特窗函数的时域形式可以表示为 当n为奇数时 当n为偶数时 Bartlett函数：生成巴特利特窗 调用方式 w = bartlett(n)： (1) 输入参数n是窗函数的长度； (2) 输出参数w是由窗函数的值组成的n阶向量。 (3) 巴特利特窗是两个矩形窗的卷积。 (4) 此函数生成的巴特利特窗的两端总是0。 (5) 该函数与三角窗函数非常类似。当n是奇数时，则有bartlett(n+2)的中间部分等于triang(n)；当n是偶数时，巴特利特窗的傅立叶变换是负数。 应用说明 例：生成一个