测试技术信号分析与处理

1、测试技术信号分析与处理 第第2 2章章 信号的分析与处理信号的分析与处理 Signal Analysis and Processing 2.0 序序(Introduction) 2.1 信号的时域分析信号的时域分析(Signal Analysis in Time Domain) 2.2 信号的相关分析信号的相关分析(Signal Correlation) 2.3 信号的频域分析信号的频域分析(Signal Analysis in Frequency Domain) 2.4 数字信号处理基础数字信号处理基础(Basic of Digital Signal Processing) 测试技术信号分析

2、与处理 第第2 2章章 信号的分析与处理信号的分析与处理 信号分析与处理的目的：信号分析与处理的目的： 1）剔除信号中的噪声和干扰，即提高信噪比；）剔除信号中的噪声和干扰，即提高信噪比； 2）消除测量系统的误差，修正畸变的波形；）消除测量系统的误差，修正畸变的波形； 3）强化、突出有用信息，削弱无用部分；）强化、突出有用信息，削弱无用部分； 4）将信号加工、处理、变换，以便更容易识别和分析）将信号加工、处理、变换，以便更容易识别和分析 信号的特征，解释被测对象所表现的各种物理现象。信号的特征，解释被测对象所表现的各种物理现象。 序序 (Introduction) 测试技术信号分析与处理 信号分

3、析和信号处理是密切相关的，二者并没有明确信号分析和信号处理是密切相关的，二者并没有明确 的界限。的界限。 本章重点讨论频域分析。本章重点讨论频域分析。 信号分析和处理的方法主要有模拟分析方法和数字处信号分析和处理的方法主要有模拟分析方法和数字处 理分析方法。理分析方法。 数字信号处理可以在专用计算机上进行，也可以在通数字信号处理可以在专用计算机上进行，也可以在通 用计算机上实现。用计算机上实现。 序序 测试技术信号分析与处理 2.1 2.1 信号的时域分析信号的时域分析 ( (Signal Analysis in Time Domain) 离散时间序列统计参数离散时间序列统计参数 x N n

4、n N x x N 1 1 lim 2.1.1 特征值分析特征值分析 离散信号的离散信号的绝对平均值绝对平均值(absolute mean) 离散信号的离散信号的均值均值(mean) N 为离散点数为离散点数 x N n n N x x N 1 1 lim 测试技术信号分析与处理 离散信号的均方值离散信号的均方值(mean square) 2 x N n N x n x N 1 22 1 lim 信号的均方根值信号的均方根值(root of mean square) ，即为有效值，即为有效值 2 smrx x 离散信号的方差离散信号的方差(variance) 2 x n n xn N x x

5、N 1 22 1 lim 信号的时域分析信号的时域分析 测试技术信号分析与处理 特征值分析的应用特征值分析的应用 信号的时域分析信号的时域分析 旋转机械振动标准旋转机械振动标准 测试技术信号分析与处理 2.1.2 概率密度概率密度(probability density)函数分析函数分析 2 1 d)()( x x xxpxP 正弦信号正弦信号 正弦加随机噪声正弦加随机噪声 窄带随机信号窄带随机信号 宽带随机信号宽带随机信号 概率密度函数概率密度函数 常见信号的概率密度函数：常见信号的概率密度函数： 信号的时域分析信号的时域分析 测试技术信号分析与处理 正态分布随机信号的概率密度函数正态分布随

6、机信号的概率密度函数 正态分布又叫高斯分布，正态分布又叫高斯分布， 是概率密度函数中最重要是概率密度函数中最重要 的一种分布。的一种分布。 2 2 2 )( exp 2 1 )( x x x xp 因此，因此， 信号的时域分析信号的时域分析 997. 033 95. 022 68. 0 xxxx xxxx xxxx xP xP xP 测试技术信号分析与处理 2.2 2.2 信号的相关分析信号的相关分析 (Signal Correlation Analysis) 2.2.1 相关系数相关系数 x与与y变量的相关性变量的相关性 x x y y x x y y x x y y 不相关不相关相关相关

7、000 测试技术信号分析与处理 22 )()( )( yx yx yx yx yx yExE yxE 变量变量x和和y之间的相关程度常用相关系数表示：之间的相关程度常用相关系数表示： 由柯西由柯西-许瓦兹不等式许瓦兹不等式 所以，所以， 信号的相关分析信号的相关分析 )()()( 222 yxyx yExEyxE 1 xy 测试技术信号分析与处理 2.2.2 自相关自相关(self-correlation)分析分析 2 2 0 2 0 d)()( 1 lim d)()( 1 lim )( x x T T x x T x T x ttxtx T ttxtx T 相关系数相关系数 信号的相关分析信

8、号的相关分析 测试技术信号分析与处理 自相关函数定义自相关函数定义 T T x ttxtx T R 0 d)()( 1 lim)( 2 2 )( )( x xx x R T x ttxtx T R 0 d)()( 1 )(周期信号：周期信号： 非周期信号：非周期信号： ttxtxRxd)()()( 进一步，对于周期信号和非周期信号有：进一步，对于周期信号和非周期信号有： 信号的相关分析信号的相关分析 测试技术信号分析与处理 自相关函数的性质自相关函数的性质 自相关函数为实偶函数自相关函数为实偶函数 )()( xx RR )( )(d)()( 1 lim d)()( 1 lim)( 0 0 x

9、T T T T x R ttxtx T ttxtx T R 证明：证明： T xxx T x ttx T R 0 2222 d)( 1 lim)0( 信号的相关分析信号的相关分析 2222 )( xxxxx R 测试技术信号分析与处理 自相关函数的性质自相关函数的性质 0)( x 2 )( xx R 信号的相关分析信号的相关分析 0 )( x R 22 xx 22 xx 2 x 测试技术信号分析与处理 周期函数的自相关函数仍为同频率的周期函数周期函数的自相关函数仍为同频率的周期函数 )()(d)()( 1 lim )(d)()( 1 lim)( 0 0 x T T T T x Rttxtx T

10、 nTtnTtxnTtx T nTR 例例2.1 求正弦函数求正弦函数 的自相关函数。的自相关函数。 )sin()( 0 txtx T x ttxtx T R 0 d)()( 1 )( T tttx T 0 0 d )(sin)sin( 1 2Tt 把把 解：解： 代入代入 2 0 2 0 2 0 cos 2 d)sin(sin 2 )( xx Rx 信号的相关分析信号的相关分析 ， 测试技术信号分析与处理 2.2.3 互相关互相关(cross-correlation)分析分析 互相关函数的概念互相关函数的概念 T T xy ttytx T R 0 d)()( 1 lim)( yx yxxy

11、yx T yx T yx T xx T xy R ttytx T ttytx T )( d)()( 1 lim d )()( 1 lim )( 0 0 互相关系数互相关系数 互相关函数互相关函数 信号的相关分析信号的相关分析 测试技术信号分析与处理 互相关函数的性质互相关函数的性质 1）互相关函数是可正）互相关函数是可正 、可负的实函数。、可负的实函数。 2）互相关函数非偶函数、亦非奇函数，具有关系）互相关函数非偶函数、亦非奇函数，具有关系 )()( yxxy RR )(d)()( 1 lim d)()( 1 limd)()( 1 lim)( 0 00 yx T T T T T T xy Rt

12、txty T ttytx T ttytx T R 因为：因为： 信号的相关分析信号的相关分析 测试技术信号分析与处理 3） 的峰值不在的峰值不在 处，其峰值偏离原点的位置处，其峰值偏离原点的位置 反映了两信号时移的大小，相关程度最高。反映了两信号时移的大小，相关程度最高。 )( xy R 0 互相关函数的性质互相关函数的性质 信号的相关分析信号的相关分析 0 )( xy R yxyx yx yxyx 0 测试技术信号分析与处理 5）两个统计独立的随机信号，当均值为零时，则）两个统计独立的随机信号，当均值为零时，则 0)( xy R 信号的相关分析信号的相关分析 4）互相关函数的限制范围为）互相

13、关函数的限制范围为 0 )( d)()( 1 lim d)()( 1 lim)( 0 0 yx yx yx T yx T T T xy R ttytx T ttytx T R 证明证明 有上述结论。有上述结论。 yxyxxyyxyx R)( 测试技术信号分析与处理 6) 两个不同频率的周期信号，其互相关为零。两个不同频率的周期信号，其互相关为零。 T T T T xy tttyx T ttytx T R 0 221100 0 d )(sin()sin( 1 lim d)()( 1 lim)( 0 7）周期信号与随机信号的互相关函数为零。）周期信号与随机信号的互相关函数为零。 信号的相关分析信号

14、的相关分析 测试技术信号分析与处理 例例2.2 求两个同频率的正弦函数求两个同频率的正弦函数 和和 的互相关函数。的互相关函数。 )sin()( 0 txtx )sin()( 0 tyty 解：因为信号是周期函数，可以用一个共同周期内解：因为信号是周期函数，可以用一个共同周期内 的平均值代替其整个历程的平均值，故的平均值代替其整个历程的平均值，故 )cos( 2 1 d)(sin)sin( 1 d)()( 1 lim)( 00 0 00 0 yx ttytx T ttytx T R T T T xy 信号的相关分析信号的相关分析 测试技术信号分析与处理 d 速度v 透镜 光电池 可调延迟 相关

15、器 钢带 )( xy R d dv )( xy R )(tx )(ty 0 d 信号的相关分析信号的相关分析 钢带运动速度的非接触测量钢带运动速度的非接触测量 测试技术信号分析与处理 相关分析在故障诊断中的应用相关分析在故障诊断中的应用 m 2 1 vS 信号的相关分析信号的相关分析 x1(t)x2(t) t s 测试技术信号分析与处理 2.3 2.3 信号的频域分析信号的频域分析 ( (Signal Analysis in Frequency Domain) 信号的时域描述反映了信号幅值随时间变化的特征；信号的时域描述反映了信号幅值随时间变化的特征； 相关分析从时域为在噪声背景下提取有用信息

16、提供了相关分析从时域为在噪声背景下提取有用信息提供了 手段；手段； 信号的频域描述反映了信号的频率结构和各频率成分信号的频域描述反映了信号的频率结构和各频率成分 的幅值大小；的幅值大小； 功率谱密度函数、相干函数、倒谱分析则从频域为研功率谱密度函数、相干函数、倒谱分析则从频域为研 究平稳随机过程提供了重要方法。究平稳随机过程提供了重要方法。 测试技术信号分析与处理 2.3.1 巴塞伐尔（巴塞伐尔（Paseval）定理）定理 信号在时域中的总能量与信号在频域中的总能量相等信号在时域中的总能量与信号在频域中的总能量相等 ffXttxd| )(|d)( 22 )(*)()()( 2121 fXfXt

17、xtx fffXfXttxtx tf d)()(de )()( 021 j2 21 0 0 0 f ffXfXttxtxd)()(d)()( 2121 )()()( 21 txtxtx ffXffXfXttxd| )(|d)()(d)( 2*2 由卷积定理由卷积定理 即即 令令 令令 信号的频域分析信号的频域分析 ，则则 测试技术信号分析与处理 2.3.2功率谱功率谱(power spectrum)分析及其应用分析及其应用 de )()( j2f xx RfS ffSR f xx de )()( 2j 定义随机信号的自功率谱密度函数（自谱）为定义随机信号的自功率谱密度函数（自谱）为 其逆变换为

18、其逆变换为 定义两随机信号的互功率谱密度函数（互谱）为定义两随机信号的互功率谱密度函数（互谱）为 de )()( 2jf xyxy RfS ffSR f yxyx de )()( 2j 其逆变换为其逆变换为 信号的频域分析信号的频域分析 测试技术信号分析与处理 功率谱密度函数的物理意义功率谱密度函数的物理意义 表示信号的功率密度沿频率轴的分布，故又称表示信号的功率密度沿频率轴的分布，故又称 为功率谱密度函数。为功率谱密度函数。 )( fSx 信号的频域分析信号的频域分析 测试技术信号分析与处理 自功率谱密度函数自功率谱密度函数 和幅值谱和幅值谱 及能谱及能谱 之间的关系之间的关系 )( fSx

19、)( fX 2 |)(|fX ffX T ttx T P T T T av d 1 limd 1 lim 2 0 2 由巴塞伐尔定理：由巴塞伐尔定理： T x T x ffSttx T R 0 2 d)(d)( 1 lim)0( 由功率谱定义：由功率谱定义： 21 limfX T S T x 因此，有因此，有 信号的频域分析信号的频域分析 测试技术信号分析与处理 自功率谱密度函数是偶函数自功率谱密度函数是偶函数,它的频率范围它的频率范围 ， 又称双边自功率谱密度函数。它在频率范围又称双边自功率谱密度函数。它在频率范围 的函数值是其在的函数值是其在 频率范围函数值的对称映射，频率范围函数值的对称

20、映射， 因此因此 。 ),( ),0( )(2)(fSfG xx 单边谱和双边谱单边谱和双边谱 信号的频域分析信号的频域分析 0 f )(),(fGfS xx )( fGx )( fSx )0 ,( 测试技术信号分析与处理 功率谱的应用功率谱的应用 1 1）自）自功率谱密度功率谱密度 与幅值谱与幅值谱 及系统频率响应函数及系统频率响应函数 的关系的关系 Sf x ( )| )(|fX )( fH )( )( )( fX fY fH 信号的频域分析信号的频域分析 )( )( fH th )( )( fY ty )( )( fX tx 测试技术信号分析与处理 )( )( )( )( )( )( )

21、( )( )( fG fG fS fS fX fX fX fY fH xx xy xx xy )(j)()(fXfXfX IR )(j)()(fXfXfX IR 若 2 22 )()()()()(fXfXfXfXfX IR 2 * )( )( )( )( )( )( )( )()(fH fS fS fX fY fX fY fHfH x y )( )( )( fS fS fH x y 信号的频域分析信号的频域分析 测试技术信号分析与处理 )(| )(|)( 2 fSfHfS xy )(| )(|)( 2 fGfHfG xy 输入输入/输出自功率谱密度函数与系统频率响应函数关系输出自功率谱密度函数

22、与系统频率响应函数关系 通过输入、输出自谱的分析，就能得出系统的幅频特性。通过输入、输出自谱的分析，就能得出系统的幅频特性。 但这样的谱分析丢失了相位信息，不能得出系统的相频但这样的谱分析丢失了相位信息，不能得出系统的相频 特性。特性。 单输入、单输出的理想线性系统单输入、单输出的理想线性系统 )()()(fSfHfS xxy 信号的频域分析信号的频域分析 测试技术信号分析与处理 2）互谱排除噪声影响）互谱排除噪声影响 )()()()()( 321 tntntntxty )()()()()( 321 xnnxnxxxxy RRRRR 由于输入和噪声是独立无关的，由于输入和噪声是独立无关的， )

23、()( xxxy RR )()()()( fSfHfSfS xxxxy 信号的频域分析信号的频域分析 )( )( 1 1 fH th )( )( fY ty )( )( fX tx + )( )( 2 2 fH th + + + + )( )( 1 1 fN tn )( )( 2 2 fN tn )( )( 3 3 fN tn 测试技术信号分析与处理 3）功率谱在设备诊）功率谱在设备诊 断中的应用断中的应用 汽车变速箱上加速汽车变速箱上加速 度信号的功率谱图度信号的功率谱图 正常正常 异常异常 故障频率故障频率 信号的频域分析信号的频域分析 （a） （b） 测试技术信号分析与处理 压缩机振动瀑

24、布图压缩机振动瀑布图 频率 f /Hz 转速 r /min 振幅 A /um 0 100200300 4）瀑布）瀑布( (water fall) )图图 信号的频域分析信号的频域分析 测试技术信号分析与处理 5）坎贝尔（）坎贝尔（Canbel）图）图 信号的频域分析信号的频域分析 测试技术信号分析与处理 2.3.3 相干函数相干函数(coherence function) ) 1)(0( )()( | )(| )( 2 2 2 f fSfS fS f xy yx xy xy 相干函数为零相干函数为零 - 输出信号与输入信号不相干。输出信号与输入信号不相干。 相干函数为相干函数为1 - 输出与输

25、入信号完全相干。输出与输入信号完全相干。 相干函数在相干函数在01之间之间 - 有如下三种可能：有如下三种可能：测试中有外测试中有外 界噪声干扰；界噪声干扰；输出是输入和其他输入的综合输出；输出是输入和其他输入的综合输出； 系统是非线性的。系统是非线性的。 1 )()( )()( )()( )()( )()( | )(| )( 2 2 2 fSfS fSfS fSfS fSfH fSfS fS f yx xy yx x yx xy xy 对于线性系统对于线性系统 信号的频域分析信号的频域分析 测试技术信号分析与处理 油压脉动与油管振动的相干分析油压脉动与油管振动的相干分析 压油管压力脉动的基频

26、压油管压力脉动的基频 润滑油泵转速为润滑油泵转速为n=781 r/min, 油泵齿轮的齿数为油泵齿轮的齿数为z=14 信号的频域分析信号的频域分析 （a）信号）信号x(t)的自谱的自谱 （b）信号）信号y(t)的自谱的自谱 （c）相干函数）相干函数 测试技术信号分析与处理 2.3.4 倒谱倒谱(cepstrum)分析分析 倒频谱分析亦称为二次频谱分析倒频谱分析亦称为二次频谱分析 检测复杂信号频谱上的周期结构，分离和提取在密集泛检测复杂信号频谱上的周期结构，分离和提取在密集泛 频谱信号中周期成分频谱信号中周期成分 功率倒频谱函数功率倒频谱函数 )cepstrumpower)(qC p |)(lo

27、g|)()( 0 fSqCqC xp F fSqC yy log 1 F fSR y 1 F 与自相关函数比较与自相关函数比较 信号的频域分析信号的频域分析 幅值倒频谱函数幅值倒频谱函数 也可以定义为也可以定义为 2 |)(log|)(fSqC xp F 测试技术信号分析与处理 倒频谱的应用倒频谱的应用 （1）分离信息通道对信号的影响）分离信息通道对信号的影响 在机械状态监测和故障诊断中所测得的信号往往是由故在机械状态监测和故障诊断中所测得的信号往往是由故 障源经系统路径的传输而得到的响应，如欲得到该源信障源经系统路径的传输而得到的响应，如欲得到该源信 号，必须消除传递通道的影响。号，必须消除

28、传递通道的影响。 2 | )(| )()(fHfSfS xy 信号的频域分析信号的频域分析 )( )( fH th )( )( fY ty )( )( fX tx 测试技术信号分析与处理 | )(|log)(log)(log 2 fHfSfS xy FFF qCqCqC xyh 信号的频域分析信号的频域分析 测试技术信号分析与处理 （2）用倒频谱诊断齿轮故障）用倒频谱诊断齿轮故障 )sin()cos1 ()( 0m ttmAty t mA t mA tAty mm000 sin 2 sin 2 sin 齿轮的振动齿轮的振动 转轴频率转轴频率 啮合频率啮合频率 功率谱功率谱倒谱倒谱 信号的频域分

29、析信号的频域分析 Sy Cy f /Hz q /ms 测试技术信号分析与处理 信号的频域分析信号的频域分析 频谱分析仪 振动 传感器 不对中 轴承异常 不平衡 轴弯曲 合成 分解 频谱分析仪测得振动波形 频谱分析 轴弯曲 不平衡 轴承问题 不对中 频率 振 幅 利用利用FFT频谱频谱 分析，将复杂分析，将复杂 的波形转换成的波形转换成 频谱，以便进频谱，以便进 一步了解振动一步了解振动 的构成原因。的构成原因。 0 测试技术信号分析与处理 2.4 2.4 数字信号处理基础数字信号处理基础 (Basic of Digital Signail Processing) 2.4.1 数字信号处理的基本

30、步骤数字信号处理的基本步骤 信号预处理信号预处理：幅值调理、滤波、隔离直流分量、解调等。：幅值调理、滤波、隔离直流分量、解调等。 A/D转换转换：采样、量化为数字量。：采样、量化为数字量。 数字信号处理器或计算机数字信号处理器或计算机：信号分析与处理（数据截断、加窗、奇：信号分析与处理（数据截断、加窗、奇 异点剔除、趋势分离、数字滤波、时域分析、频域分析等）。异点剔除、趋势分离、数字滤波、时域分析、频域分析等）。 结果显示结果显示：数据或图形显示、：数据或图形显示、D/A、记录、打印等。、记录、打印等。 预处理A/D转换 x(t) 数字信号 处理器 或 计算机预处理A/D转换 x(t) 结果

31、显示 数字信号处理系统简图数字信号处理系统简图 测试技术信号分析与处理 2.4.2 时域采样和采样定理时域采样和采样定理 采样采样(sampling)：连续时间信号离散化的过程。：连续时间信号离散化的过程。 采样时间间隔为采样时间间隔为Ts，则，则x(t)经采样后的离散序列经采样后的离散序列x(n)为为 nn nTttxnTttxtstxnTxnx)( )()( )()()()()( sss x(n)与与x(t)是局部与整体的关系。是局部与整体的关系。 能否由能否由x(n)唯一确定或恢复出唯一确定或恢复出x(t)，或能否通过对，或能否通过对x(n)的分的分 析获得析获得x(t)的全部信息是采样

32、最关心的问题。的全部信息是采样最关心的问题。 数字信号处理基础数字信号处理基础 测试技术信号分析与处理 混叠混叠(aliasing)和采样定理和采样定理 时域采样间隔过长，造成频域周期化间隔不够大时，在重时域采样间隔过长，造成频域周期化间隔不够大时，在重 复频谱交界处出现的局部互相重叠现象，称为频率混叠复频谱交界处出现的局部互相重叠现象，称为频率混叠 （如主教材图（如主教材图5 - 2c）。）。 rn T r f T fSnTtts ss s 1 )()()( )()()()(fSfXtstx r T r fX T fSfXnx ss 1 )()()(F所以所以 数字信号处理基础数字信号处理基

33、础 测试技术信号分析与处理 A、B、C 被误认为是一条曲线被误认为是一条曲线 高频正弦信号高频正弦信号 被误认为是低频正弦信号被误认为是低频正弦信号 数字信号处理基础数字信号处理基础 x(t) 0 1 2 3 t A B C Ts 混叠现象混叠现象 4 测试技术信号分析与处理 混叠的后果是原来的高频信号将被误认为是某种相应的低频混叠的后果是原来的高频信号将被误认为是某种相应的低频 信号。信号。 发生混叠的高频成分（大于频率）发生混叠的高频成分（大于频率）f1和低频成分和低频成分f2（低于频（低于频 率）之间满足：率）之间满足： 22)( s21 fff 即即f1和和f2以以 为轴对称，可以将混

34、叠视为以为轴对称，可以将混叠视为以 为轴将为轴将 高频分量高频分量f1折叠至低频分量折叠至低频分量f2处。因此，处。因此， 称为折叠频率。称为折叠频率。 也称为奈奎斯特频率也称为奈奎斯特频率(Nyquist frequency) 。 数字信号处理基础数字信号处理基础 2 s f2 s f 2 s f 测试技术信号分析与处理 若原始信号是带限信号，则采样后信号频谱不发生重叠的条若原始信号是带限信号，则采样后信号频谱不发生重叠的条 件为件为 fs2fh 。其中。其中fh为信号中的最高频率。此即为采样定理。为信号中的最高频率。此即为采样定理。 实际工作中，实际工作中，fs常取为信号最高频率的常取为信

35、号最高频率的2.56倍以上。倍以上。 数字信号处理基础数字信号处理基础 |X(f)S (f )| f 不产生混叠的条件不产生混叠的条件 0-fhfh fs fh 测试技术信号分析与处理 消除混叠的措施消除混叠的措施 提高采样频率。但提高采样频率将导致在同样信号提高采样频率。但提高采样频率将导致在同样信号 长度下采样点数随之提高，增加计算负担。长度下采样点数随之提高，增加计算负担。 应用抗混滤波器降低信号中的最高频率。从理论上应用抗混滤波器降低信号中的最高频率。从理论上 讲，由于抗混滤波器的非理想特性，信号中高频分讲，由于抗混滤波器的非理想特性，信号中高频分 量不可能完全衰减，因此不可能彻底消除

36、混叠。量不可能完全衰减，因此不可能彻底消除混叠。 数字信号处理基础数字信号处理基础 测试技术信号分析与处理 2.4.3 截断截断(truncation)、泄漏、泄漏(leakage)和窗函数和窗函数(window) 计算机处理的数据长度是有限的，进行数字信号处理必须计算机处理的数据长度是有限的，进行数字信号处理必须 对过长时间历程的信号进行截断处理。截断相当于对信号对过长时间历程的信号进行截断处理。截断相当于对信号 进行加窗处理，截断即是将信号乘以时域的有限宽矩形窗进行加窗处理，截断即是将信号乘以时域的有限宽矩形窗 函数：函数： )2(0 )2(1 )( Tt Tt tw 即：采样后信号即：采

37、样后信号x(t)s(t)经截断成为经截断成为x(t)s(t)w(t)。 数字信号处理基础数字信号处理基础 测试技术信号分析与处理 矩形窗函数的频谱为无限带宽的矩形窗函数的频谱为无限带宽的sinc函数，即使函数，即使x(t)为带限为带限 信号，经截断后必然成为无限带宽信号，这种信号的能量信号，经截断后必然成为无限带宽信号，这种信号的能量 在频率轴分布扩展的现象称为泄漏。在频率轴分布扩展的现象称为泄漏。 无论采样频率多高，信号不可避免地出现混叠。无论采样频率多高，信号不可避免地出现混叠。 减小泄漏的措施：减小泄漏的措施： 提高截断信号长度，即提高矩形窗宽度，此时提高截断信号长度，即提高矩形窗宽度，

38、此时sinc函数主函数主 瓣变窄，旁瓣向主瓣密集，由于旁瓣衰减较快，故可减小瓣变窄，旁瓣向主瓣密集，由于旁瓣衰减较快，故可减小 泄漏，但显然采样点数随之提高，增加计算负担。泄漏，但显然采样点数随之提高，增加计算负担。 数字信号处理基础数字信号处理基础 测试技术信号分析与处理 采用其他窗函数。一个好的窗函数应当：主瓣尽可能窄采用其他窗函数。一个好的窗函数应当：主瓣尽可能窄 （提高频率分辨力）、旁瓣相对于主瓣尽可能小，且衰减（提高频率分辨力）、旁瓣相对于主瓣尽可能小，且衰减 快（减小泄漏）。快（减小泄漏）。 常用窗函数常用窗函数 矩形窗矩形窗(rectangle window) )2(0 )2(1

39、 )( Tt Tt tw fTTfWsinc)( 数字信号处理基础数字信号处理基础 测试技术信号分析与处理 三角窗三角窗( (triangle windsow) ) 20 2 2 1 )( T Tt Ttt Ttw 2 sinc 2 )( 2 T fTT fW 数字信号处理基础数字信号处理基础 W(f) -2/T w(t) 1 0 T/2T/2t T/2 0 2/T f 测试技术信号分析与处理 汉宁窗汉宁窗( (Hanning window) )（余弦窗）（余弦窗） )2(0 )2( 2 cos 2 1 2 1 )( Tt Tt T t tw )1()1( 4 1 )( 2 1 )( RRR

40、TfWTfWfWfW ftTfWRsinc)(其中其中 数字信号处理基础数字信号处理基础 2/T w(t) 1 0 T/2T/2t W(f) T/2 -2/T 0 f 测试技术信号分析与处理 指数窗指数窗( (exponent window) ) )0(0 )0, 0(e )( t at tw at 22 )2( 1 )( fa fW 数字信号处理基础数字信号处理基础 w(t) 1 0 t W(f) 1/a 0 f 测试技术信号分析与处理 几种典型窗函数的技术指标几种典型窗函数的技术指标 数字信号处理基础数字信号处理基础 窗函数类型 主瓣宽度 最大旁瓣幅度 旁瓣衰减速度 矩形窗 2/T 13

41、dB 6dB/oct 三角形窗 4/T 26 dB 12dB/oct 汉宁窗 4/T 32 dB 18dB/oct 测试技术信号分析与处理 2.4.4 频域采样与栅栏效应频域采样与栅栏效应 频域采样与时域采样类似，频域采样导致对时域截断信号频域采样与时域采样类似，频域采样导致对时域截断信号 进行周期延拓，将时域截断信号进行周期延拓，将时域截断信号“改造改造”为周期信号。为周期信号。 数字信号处理基础数字信号处理基础 x(t) w(t) 0 t -f0f00 T T-T s2(t) 0 S2(f) 0 f f X(f)*W(f)S2(f)x(t) w(t)*s2(t) T 频域采样频域采样 -f

42、0f0 0 f 0 测试技术信号分析与处理 经频域采样后的频谱仅在各采样点上存在，而非采样点的经频域采样后的频谱仅在各采样点上存在，而非采样点的 频谱则被频谱则被“挡住挡住”无法显示（视为无法显示（视为0），这种现象称为栅栏），这种现象称为栅栏 效应。显然，采样必然带来栅栏效应。效应。显然，采样必然带来栅栏效应。 在时域，只要满足采样定理，栅栏效应不会丢失信号信息在时域，只要满足采样定理，栅栏效应不会丢失信号信息 在频域，则有可能丢失重要的或具有特征的频率成分（由在频域，则有可能丢失重要的或具有特征的频率成分（由 于泄漏，丢失频率成分附近的频率有可能存在），导致谱于泄漏，丢失频率成分附近的频率

43、有可能存在），导致谱 分析结果失去意义。分析结果失去意义。 数字信号处理基础数字信号处理基础 测试技术信号分析与处理 X(f)S0(f) x(t)s0(t) s0(t) x(t) S0(f) X(f) w(t) W(f) x(t)s0(t)w(t) X(f)S0(f)W(f) s1(t) S1 (f) X(f)S0(f)W(f) S1 (f ) x(t)s0(t)w(t) * s1(t) 离散傅里叶变换图解说明离散傅里叶变换图解说明 (a) (b) (c) (d) (e) (f) (g) 测试技术信号分析与处理 频率分辨力、整周期截断频率分辨力、整周期截断 频率采样间隔频率采样间隔 f 决定了频率分辨力。决定了频率分辨力。 f 越小，分辨力越越小，分辨力越 高，被挡住的频率成分越少。高，被挡住的频率成分越少。 由于由于DFT在频域的一个周期内（周期为：在频域的一个周期内（周期为：1/Ts）输出）输出N个有个有 效谱值，故频率间隔为：效谱值，故频率间隔为： TN f N T f 11 ss 显然，可以通过降低显然，可以通过降低 fs 或提高或提高N 以提高以提高 f。但前者受采样。但前者受采样 定理的限制，不可能随意降低，后者必然增加计算量。定理的限制，不可能随意降低，后者必然增加计