工程振动与测试技术课件-第十四章数字信号分析.pptVIP

每个蝶型有四个数据点，上面两个是参加计算的数据，下面两个是计算的结果，箭头表示参加计算的数与结果之间的联系，蝶型的一边写上“旋转因子”数W。不管N有多长，其蝶型计算流程图是一样的。 工 程 振 动 与 测 试 Testing Techniques 在蝶形计算中，数据是按它的奇偶位置来排列的，每进行一次计算都要排列一次。 以N=16为例来说明排列情况，因为N=24，故共有4排，数据前后共排列了4次，排列情况如表所示。 以上是针对时域离散化后所进行的快速傅立叶变换，对于傅立叶逆变换，同理可进行在频域离散化的快速傅立叶逆变换。 工 程 振 动 与 测 试 Testing Techniques 14.5 泄漏与窗函数 数字信号分析对有限时间长度T的离散时间序列进行离散傅里叶变换（DFT）运算，这意味着首先要对时域信号进行截断。这种截断将导致频谱分析出现误差，其效果是使得本来集中于某一频率的功率（或能量），部分被分散到该频率邻近的频域，这种现象称为“泄漏”效应。 工 程 振 动 与 测 试 Testing Techniques 余弦信号被矩形窗截断形成的泄漏 工 程 振 动 与 测 试 Testing Techniques 以余弦信号x(t)=Acos(2?f0t) 为例说明截断前后的频谱变化的泄漏效应。 a. 有限时间长度T的离散时间序列信号被截断，相当于原来的余弦信号乘以一个矩形窗函数。 b. 无限长度的余弦信号具有一个单一的频率成分，其单边谱是在f0处的单根分布的离散的谱线。 余弦信号被矩形窗截断形成的泄漏 工 程 振 动 与 测 试 Testing Techniques 余弦信号被矩形窗截断形成的泄漏 c. 矩形窗函数的频谱是包含一个主瓣和许多旁瓣的连续谱。 d. 时域中余弦信号乘以矩形窗函数。频域中的频谱等于原信号的频谱与窗函数频谱的卷积。 工 程 振 动 与 测 试 Testing Techniques 余弦信号被矩形窗截断形成的泄漏 e. 卷积的结果将导致离散谱变为在f0处有一主瓣，两旁各有许多旁瓣的连续谱。 f. 这也就是说，原来集中在频率f0处的功率，泄漏到了f0邻近的很宽的频带上。 工 程 振 动 与 测 试 Testing Techniques 泄漏示意图 工 程 振 动 与 测 试 Testing Techniques 余弦信号被矩形窗截断形成的泄漏 为了抑制“泄漏”，需采用特种窗函数来替代矩形窗函数。这一过程，称为窗处理，或者叫加窗。 加窗的目的，是使在时域上截断信号两端的波形突变变为平滑，在频域上尽量压低旁瓣的高度。 工 程 振 动 与 测 试 Testing Techniques 在一般情况下，压低旁瓣通常伴随着主瓣的变宽，但是旁瓣的泄漏是主要考虑因素，然后才考虑主瓣变宽的泄漏问题。 余弦信号被矩形窗截断形成的泄漏 工 程 振 动 与 测 试 Testing Techniques 常用窗函数的时域图象 工 程 振 动 与 测 试 Testing Techniques 常用窗函数的时域图象 在数字信号处理中常用的窗函数有： ①矩形（Rectangular）窗: w(t)=1 0≤t≤T ②汉宁（Hanning）窗: w(t)= 0≤t≤T 工 程 振 动 与 测 试 Testing Techniques ③凯塞—贝塞尔（kaiser-bessel）窗 0≤t≤T ④平顶（Rectangular）窗 0≤t≤T 常用窗函数的时域图象 工 程 振 动 与 测 试 Testing Techniques 常用窗函数的频谱图 工 程 振 动 与 测 试 Testing Techniques 为了保持加窗后的信号能量不变，要求窗函数曲线与时间坐标轴所包围的面积相等。对于矩型窗，该面积为T×1，因此，对于任意窗函数W(t)，必需满足积分关系式 工 程 振 动 与 测 试 Testin