Chapter02故障诊断的信号处理方法

《Chapter02故障诊断的信号处理方法》由会员分享，可在线阅读，更多相关《Chapter02故障诊断的信号处理方法（138页珍藏版）》请在装配图网上搜索。

1、机械电子工程学院机械电子工程学院第二章第二章 故障诊断的信号处理方法故障诊断的信号处理方法本章内容本章内容1、信号的定义与分类、信号的定义与分类2、信号特征的时域提取方法、信号特征的时域提取方法3、信号特征的频域提取方法、信号特征的频域提取方法4、信号特征的图像表示、信号特征的图像表示5、希尔伯特变换与解调分析、希尔伯特变换与解调分析6、全息谱理论和方法、全息谱理论和方法机械电子工程学院机械电子工程学院本章学习要求本章学习要求1、了解转子振型、轴颈涡动中心位置、波特图、奈奎斯、了解转子振型、轴颈涡动中心位置、波特图、奈奎斯 特图、三维坐标图、阶比谱分析。特图、三维坐标图、阶比谱分析。2、理解信

2、号的功率谱、细化谱、倒频谱、希尔伯特变换原、理解信号的功率谱、细化谱、倒频谱、希尔伯特变换原理及结果的物理意义。理及结果的物理意义。3、理解机械信号处理技术的物理意义、轴心轨迹图技术、理解机械信号处理技术的物理意义、轴心轨迹图技术、全息谱技术。全息谱技术。4、掌握振动监测的基本参数、时域指标、频域分析结果的、掌握振动监测的基本参数、时域指标、频域分析结果的 物理意义等。物理意义等。第二章第二章 故障诊断的信号处理方法故障诊断的信号处理方法机械电子工程学院机械电子工程学院2.1 信号的定义和分类信号的定义和分类 信号信号是表征客观事物状态或行为信息的载体。是表征客观事物状态或行为信息的载体。信号

3、具有信号具有能量能量，它描述了物理量的变化过程，在数学上可表示为一个，它描述了物理量的变化过程，在数学上可表示为一个或几个独立变量的函数，也可以取为随时间或空间变化的图形。或几个独立变量的函数，也可以取为随时间或空间变化的图形。1 1、按信号随时间的变化规律分、按信号随时间的变化规律分机械电子工程学院机械电子工程学院1、按信号随时间的变化规律分、按信号随时间的变化规律分2.1 信号的定义和分类信号的定义和分类质量质量-弹簧振动系统（弹簧振动系统（无阻尼无阻尼）余弦信号的波形图余弦信号的波形图机械电子工程学院机械电子工程学院2.1 信号的定义和分类信号的定义和分类2、按信号幅值随时间变化的连续性

4、分、按信号幅值随时间变化的连续性分汽车速度汽车速度(连续信号连续信号)(a)含第一类间断点的信号含第一类间断点的信号(b)锯齿波锯齿波机械电子工程学院机械电子工程学院2.1 信号的定义和分类信号的定义和分类2、按信号幅值随时间变化的连续性分、按信号幅值随时间变化的连续性分(c)矩形脉冲矩形脉冲(d)截断信号截断信号连续信号连续信号每日股市的指数变化每日股市的指数变化(离散信号离散信号)每隔每隔2us对正弦信号采样获得的离散信号对正弦信号采样获得的离散信号机械电子工程学院机械电子工程学院 当信号当信号 在在 内满足下式（即平方可积）时：内满足下式（即平方可积）时：则该信号的能量是有限的，称为则该

5、信号的能量是有限的，称为能量（有限）信号能量（有限）信号。2.1 信号的定义和分类信号的定义和分类3、按信号的能量特征分、按信号的能量特征分 若信号若信号 在在 内内 ，而在有限区间，而在有限区间 内的平均功率是有限的，即：内的平均功率是有限的，即：则信号称为则信号称为功率信号功率信号。机械电子工程学院机械电子工程学院 频域有限信号频域有限信号是指信号经过傅立叶变换，在频域内占据一是指信号经过傅立叶变换，在频域内占据一定带宽，在带宽外恒等于定带宽，在带宽外恒等于0。例如，正弦信号、。例如，正弦信号、sinc(t)函数、带函数、带限白噪声等为限白噪声等为时域无限、频域有限时域无限、频域有限信号。

6、信号。函数、白噪声、理函数、白噪声、理想采样信号等，则为频域无限信号。想采样信号等，则为频域无限信号。2.1 信号的定义和分类信号的定义和分类4、按信号的持续范围分、按信号的持续范围分 时域有限信号时域有限信号是在有限时间区间内有定义，而在区间外恒是在有限时间区间内有定义，而在区间外恒等于等于0。例如，矩形脉冲、三角脉冲、余弦脉冲等。而周期信。例如，矩形脉冲、三角脉冲、余弦脉冲等。而周期信号、指数衰减信号、随机过程等，则称为时域无限信号。号、指数衰减信号、随机过程等，则称为时域无限信号。时域有限信号的频谱，在频率轴上可以延伸至无限远。而时域有限信号的频谱，在频率轴上可以延伸至无限远。而一个在频

7、域上具有有限带宽的信号，必然在时间轴上延伸至无一个在频域上具有有限带宽的信号，必然在时间轴上延伸至无限远处。限远处。一个信号不能够在时域和频域上都是有限的一个信号不能够在时域和频域上都是有限的。机械电子工程学院机械电子工程学院2.2.1 时域分解时域分解1、直流分量和交流分量、直流分量和交流分量 2.2 信号特征的时域提取方法信号特征的时域提取方法信号信号 可以分解为直流分量可以分解为直流分量 与交流分量与交流分量 ，即：，即：信号分解为直流信号分解为直流分量和交流分量分量和交流分量信号分解为趋势信号分解为趋势项和交流分量项和交流分量机械电子工程学院机械电子工程学院2、脉冲分量、脉冲分量 2.

8、2.1 时域分解时域分解信号分解为矩形窄脉冲之和信号分解为矩形窄脉冲之和机械电子工程学院机械电子工程学院3、实部分量和虚部分量、实部分量和虚部分量 2.2.1 时域分解时域分解信号的实数表示法信号的实数表示法 信号的复数表示法信号的复数表示法信号的实数和复数表示法及其对应关系信号的实数和复数表示法及其对应关系机械电子工程学院机械电子工程学院4、正交函数分量、正交函数分量 2.2.1 时域分解时域分解 信号信号 可以用正交函数集可以用正交函数集 来表示，来表示，即：即：正交条件为：正交条件为：即在区间即在区间 内分量乘积的积分为零，任一分量在此区内分量乘积的积分为零，任一分量在此区间内能量为有限

9、值。间内能量为有限值。分量系数分量系数 代表了该正交函数分量的大小，可在满足最代表了该正交函数分量的大小，可在满足最小均方差条件下求得：小均方差条件下求得：机械电子工程学院机械电子工程学院2.2.2 时域相关分析时域相关分析1、相关的概念、相关的概念 2.2 信号特征的时域提取方法信号特征的时域提取方法相关是指客观事物变化量之间的相互依赖关系。相关是指客观事物变化量之间的相互依赖关系。变量变量x x和和y y之间的不同相关情况之间的不同相关情况 两个随机变量两个随机变量x x和和y y之之间的线性相关程度可用相间的线性相关程度可用相关系数来描述，即：关系数来描述，即：机械电子工程学院机械电子工

10、程学院1、相关的概念、相关的概念 2.2.2 时域相关分析时域相关分析 相关系数可以定量地描述两个变量相关系数可以定量地描述两个变量x x和和y y之间的相似或相之间的相似或相依关系，但它也有局限性。依关系，但它也有局限性。信号信号 和它的时延信号和它的时延信号机械电子工程学院机械电子工程学院2、相关函数、相关函数 2.2.2 时域相关分析时域相关分析 设随机变量设随机变量x、y是一个与时间有关的函数，令两个信号之是一个与时间有关的函数，令两个信号之间产生时差间产生时差（即令某个信号在时间轴上平移，平移量为（即令某个信号在时间轴上平移，平移量为），），互相关函数的定义为：互相关函数的定义为：如

11、果如果x和和y为同一函数，则成为自相关函数：为同一函数，则成为自相关函数：机械电子工程学院机械电子工程学院2、相关函数、相关函数 2.2.2 时域相关分析时域相关分析对于功率信号，相关函数的定义为：对于功率信号，相关函数的定义为：机械电子工程学院机械电子工程学院3、相关分析的工程应用（测距）、相关分析的工程应用（测距）2.2.2 时域相关分析时域相关分析两传感器中点至泄漏点的距离为：两传感器中点至泄漏点的距离为：机械电子工程学院机械电子工程学院3、相关分析的工程应用（消除噪声求相位）、相关分析的工程应用（消除噪声求相位）2.2.2 时域相关分析时域相关分析基准正弦信号：基准正弦信号：基准余弦信

12、号：基准余弦信号：转轴振动信号：转轴振动信号：由此可直接获得同频振由此可直接获得同频振动信号的幅值及其相对动信号的幅值及其相对于基准信号的相位：于基准信号的相位：机械电子工程学院机械电子工程学院2.2.3.1 有量纲指标有量纲指标1、平均值、平均值 平均值描述信号的稳定分量，又称平均值描述信号的稳定分量，又称直流分量直流分量。指信号在观测。指信号在观测时间时间T内取值的时间平均，即：内取值的时间平均，即：式中式中T为信号的观测区间。为信号的观测区间。均值的离散形式为：均值的离散形式为：2.2 信号特征的时域提取方法信号特征的时域提取方法2.2.3 时域统计指标时域统计指标机械电子工程学院机械电

13、子工程学院位移传感器测得的振动信号位移传感器测得的振动信号2.2.3.1 有量纲指标有量纲指标 在不存在摩擦碰撞的情况下，测量在不存在摩擦碰撞的情况下，测量加速度、速度时，平均值反映了测量系加速度、速度时，平均值反映了测量系统的温漂、时漂等参数变化；测量位移统的温漂、时漂等参数变化；测量位移时，平均值反映磨损量的变化。时，平均值反映磨损量的变化。机械电子工程学院机械电子工程学院2、均方值、均方值用于描述振动信号的用于描述振动信号的能量能量(功率功率)。3、均方根值（有效值）、均方根值（有效值）有效值是机械故障诊断系统中用于判别运转状态是否正常的有效值是机械故障诊断系统中用于判别运转状态是否正常

14、的重要指标。有效值也描述振动信号的重要指标。有效值也描述振动信号的能量（功率）能量（功率），稳定性、，稳定性、重复性好，当这项指标超出正常值（故障判定限）较多时，通重复性好，当这项指标超出正常值（故障判定限）较多时，通常表示机械设备存在常表示机械设备存在故障隐患或故障故障隐患或故障。若有效值的物理参数是速度（单位：若有效值的物理参数是速度（单位：mm/s），则有效值就），则有效值就成为用于判定机械状态等级的成为用于判定机械状态等级的振动烈度振动烈度指标。指标。2.2.3.1 有量纲指标有量纲指标机械电子工程学院机械电子工程学院4、方均根值、方均根值5、方差、方差 方差反映信号中的动态部分（方差

15、反映信号中的动态部分（波动程度波动程度）。）。方差的平方根称为标准差。方差的平方根称为标准差。若信号的均值为零，则均方值等于方差。若信号的均值为零，则均方值等于方差。6、峰值、峰值 2.2.3.1 有量纲指标有量纲指标机械电子工程学院机械电子工程学院2.2.3.2 无量纲指标无量纲指标1、波形指标、波形指标2、峰值指标、峰值指标峰值指标是用来检测信号中是否存在峰值指标是用来检测信号中是否存在冲击冲击的一个统计指标。的一个统计指标。3、脉冲指标、脉冲指标也是用来检测信号中是否存在也是用来检测信号中是否存在冲击冲击的一个统计指标。的一个统计指标。机械电子工程学院机械电子工程学院4、裕度指标、裕度指

16、标 裕度指标用于检测机械设备的裕度指标用于检测机械设备的磨损磨损情况。情况。若偏度指标变化不大，峰值与若偏度指标变化不大，峰值与方均根值方均根值的比值增大，说明的比值增大，说明由于磨损导致间隙增大，因而振动的峰值比方均根值增加快，由于磨损导致间隙增大，因而振动的峰值比方均根值增加快，其裕度指标也增大了。其裕度指标也增大了。2.2.3.2 无量纲指标无量纲指标机械电子工程学院机械电子工程学院 偏度指标反映振动信号的偏度指标反映振动信号的不对称性不对称性。表示信号概率密度。表示信号概率密度函数的函数的中心偏离标准正态分布中心偏离标准正态分布的程度的程度，反映信号幅值分布相，反映信号幅值分布相对其对

17、其理想均值理想均值的不对称性。的不对称性。除有急回特性的机械设备除有急回特性的机械设备外，如果存在着某一方向的外，如果存在着某一方向的摩摩擦或碰撞擦或碰撞，就会造成振动波性，就会造成振动波性的不对称，使偏度指标增大。的不对称，使偏度指标增大。5、偏度指标、偏度指标2.2.3.2 无量纲指标无量纲指标机械电子工程学院机械电子工程学院 峭度指标表示信号概率密度峭度指标表示信号概率密度函数峰顶的函数峰顶的陡峭程度，陡峭程度，反映振反映振动信号中的动信号中的冲击特征（冲击特征（波形中波形中的冲击分量的大小）。的冲击分量的大小）。峭度指标对信号中的冲击特峭度指标对信号中的冲击特征很敏感，征很敏感，正常情

18、况下其值应正常情况下其值应该在该在3左右左右，如果这个值接近，如果这个值接近4或超过或超过4，则说明机械的运动状，则说明机械的运动状况中存在况中存在冲击性振动冲击性振动。一般情。一般情况下是间隙过大、滑动副表面况下是间隙过大、滑动副表面存在破碎等原因。存在破碎等原因。6、峭度指标、峭度指标2.2.3.2 无量纲指标无量纲指标机械电子工程学院机械电子工程学院 参数指标诊断是参数指标诊断是使用较早且比较有效使用较早且比较有效的诊断方法。诊断参数的诊断方法。诊断参数指标一般应满足如下要求：指标一般应满足如下要求：1、易于测量和计算，所需计算机存储量小。、易于测量和计算，所需计算机存储量小。2、能够敏

19、锐地反映和预报机器的早期故障。、能够敏锐地反映和预报机器的早期故障。3、不受机器运行状态不受机器运行状态，如负载、转速等变化的影响。，如负载、转速等变化的影响。4、能够指示故障的存在，以便及时排查故障。、能够指示故障的存在，以便及时排查故障。在流程生产工业中，往往有这样的情况，当发现设备的情况在流程生产工业中，往往有这样的情况，当发现设备的情况不好，某项或多项特征指标上升，但设备不能停产检修，只能不好，某项或多项特征指标上升，但设备不能停产检修，只能让设备带故障运行。当这些让设备带故障运行。当这些指标从峰值跌落指标从峰值跌落时，往往预示某个时，往往预示某个零件已经损坏零件已经损坏，若这些指标（

20、含其它指标）再次上升，则预示，若这些指标（含其它指标）再次上升，则预示大的设备故障将要发生，此时大的设备故障将要发生，此时需要格外注意需要格外注意。2.2.2.3 运用统计指标的注意之处运用统计指标的注意之处 机械电子工程学院机械电子工程学院 时域统计特征指标只能反映机械设备的时域统计特征指标只能反映机械设备的总体运转状态总体运转状态是否是否正常，因而在设备故障诊断系统中用于故障正常，因而在设备故障诊断系统中用于故障监测监测，趋势，趋势预报预报。要识别机械运动状态，要识别机械运动状态，知道故障的知道故障的部位部位、故障的、故障的类型，类型，就就需要进一步做精密分析，把反映故障部位和类型的需要进

21、一步做精密分析，把反映故障部位和类型的相关信号从相关信号从传感器测得的合成信号中传感器测得的合成信号中分离分离出来。出来。频谱是信号在频域上的重要特征，它反映了信号的频谱是信号在频域上的重要特征，它反映了信号的频率成频率成分以及分布情况分以及分布情况。2.3 信号特征的频域提取方法信号特征的频域提取方法 2.3.1 频域信号与时域信号的关系频域信号与时域信号的关系 机械电子工程学院机械电子工程学院 法国数学家，物理学家。傅立叶出身平民，是一个裁缝的法国数学家，物理学家。傅立叶出身平民，是一个裁缝的儿子，早在小学时就对数学产生浓厚的兴趣。后来他也曾在儿子，早在小学时就对数学产生浓厚的兴趣。后来他

22、也曾在他的母校担任数学教师。法国革命的浪潮中，他投身于政治，他的母校担任数学教师。法国革命的浪潮中，他投身于政治，从此以后，它的生活一直充满了冒险。从此以后，它的生活一直充满了冒险。1798年，傅立叶和其他队员一起，陪同拿破仑远征埃及。年，傅立叶和其他队员一起，陪同拿破仑远征埃及。并在拿破仑建立的并在拿破仑建立的Cairo研究所担任三年秘书，在工程技术以研究所担任三年秘书，在工程技术以及外交任务方面都提出了许多意见。及外交任务方面都提出了许多意见。1801年，他开始着手大年，他开始着手大范围研究埃及古迹，回国后，他被任命出版了大量的有关埃范围研究埃及古迹，回国后，他被任命出版了大量的有关埃及的

23、刊物。及的刊物。1809年拿破仑封他为男爵。年拿破仑封他为男爵。1815年，拿破仑垮台，年，拿破仑垮台，此后傅立叶在巴黎过了一段平静的学术研究生活。此后傅立叶在巴黎过了一段平静的学术研究生活。1817年，年，他被选为科学院院士，他被选为科学院院士，1822年，担任科学院常任秘书。年，担任科学院常任秘书。傅立叶于傅立叶于1807年开始他的学术论文写作，并提出求解偏微分方程的分离变量法和可年开始他的学术论文写作，并提出求解偏微分方程的分离变量法和可以将解表示成一系列任意函数的概念。于以将解表示成一系列任意函数的概念。于1822年完成论文，发表了著名论著年完成论文，发表了著名论著热的解热的解析理论析

24、理论，解决了热在非均匀加热的固体中分布传播问题，成为分析学在物理中应用，解决了热在非均匀加热的固体中分布传播问题，成为分析学在物理中应用的最早例证之一，对的最早例证之一，对19世纪数学和理论物理学的发展产生深远影响。这一著作奠定了世纪数学和理论物理学的发展产生深远影响。这一著作奠定了导热的理论基础，描述导热的定律就是以他的名字命名的。并在求解该方程时发现导热的理论基础，描述导热的定律就是以他的名字命名的。并在求解该方程时发现解解函数可以由三角函数构成的级数形式表示，从而提出任一函数都可以展成三角函数的函数可以由三角函数构成的级数形式表示，从而提出任一函数都可以展成三角函数的无穷级数无穷级数。傅

25、立叶被公认为导热理论的奠基人。傅立叶被公认为导热理论的奠基人。其他贡献有：其他贡献有：最早使用定积分符号最早使用定积分符号，改进了代数方程符号法则的证法和实根个数的，改进了代数方程符号法则的证法和实根个数的判别法等。判别法等。傅立叶傅立叶（Fourier 1768-1830）2.3.1 频域信号与时域信号的关系频域信号与时域信号的关系机械电子工程学院机械电子工程学院2.3.1 频域信号与时域信号的关系频域信号与时域信号的关系 图中左侧所示为一图中左侧所示为一组在时间坐标轴上表组在时间坐标轴上表示简谐运动信号的时示简谐运动信号的时域波形曲线。域波形曲线。图中右侧表示简谐图中右侧表示简谐运动的频率

26、和幅值。运动的频率和幅值。由于它们的幅值相同，由于它们的幅值相同，所以各个所以各个幅值谱都相幅值谱都相同同，只有相位谱中的，只有相位谱中的初相位各不相同初相位各不相同一组时域波形曲线（余弦）及其幅值谱和相位谱一组时域波形曲线（余弦）及其幅值谱和相位谱机械电子工程学院机械电子工程学院2.3.1 频域信号与时域信号的关系频域信号与时域信号的关系时域波形的分解及其频域表示时域波形的分解及其频域表示机械电子工程学院机械电子工程学院 信号是由多个正弦波组成，信号是由多个正弦波组成，频率比为：频率比为：1:3:5:7:，幅，幅值比为：值比为：1:1/3:1/5:1/7:，信号之间无相位差。信号之间无相位差

27、。需要注意的是，如果在需要注意的是，如果在频频率比率比、幅值比幅值比、相位差相位差这三这三个方面有任一个不满足以上个方面有任一个不满足以上条件，其叠加的波形便不是条件，其叠加的波形便不是方波。即使所有信号都是周方波。即使所有信号都是周期信号，只有当各信号的频期信号，只有当各信号的频率比是整数，其叠加合成信率比是整数，其叠加合成信号才表现出周期性特征。号才表现出周期性特征。2.3.1 频域信号与时域信号的关系频域信号与时域信号的关系信号的时域和频域关系信号的时域和频域关系机械电子工程学院机械电子工程学院2.3.2 周期信号的频谱周期信号的频谱 如果正弦信号的周期为如果正弦信号的周期为T，则周期，

28、则周期T与频率与频率f和角频率和角频率之间之间的关系为：的关系为：根据傅里叶级数理论，满足根据傅里叶级数理论，满足狄里赫利条件狄里赫利条件的周期信号，可以的周期信号，可以表示为若干表示为若干正弦函数的叠加正弦函数的叠加（三角函数展开式）。（三角函数展开式）。狄里赫利（狄里赫利（Dirichlet）条件：）条件：1、函数在任意有限区间内连续，或只有有限个第一类间断、函数在任意有限区间内连续，或只有有限个第一类间断点（当点（当t从左或右走向于该间断点时，函数存在有限的左极限从左或右走向于该间断点时，函数存在有限的左极限或右极限）；或右极限）；2、在一个周期内，函数存在有限个极大值或极小值。、在一个

29、周期内，函数存在有限个极大值或极小值。机械电子工程学院机械电子工程学院 在机械故障诊断的信号中，在机械故障诊断的信号中，常数分量常数分量 是直流分量，代表是直流分量，代表某个变动缓慢的物理因素，如某个间隙。基频和它的某个变动缓慢的物理因素，如某个间隙。基频和它的n n次谐波在次谐波在机械故障诊断领域都有明确的物理意义。机械故障诊断领域都有明确的物理意义。2.3.2 周期信号的频谱周期信号的频谱 机械电子工程学院机械电子工程学院傅里叶级数也可以写成复指数函数的形式。根据欧拉公式：傅里叶级数也可以写成复指数函数的形式。根据欧拉公式：2.3.2 周期信号的频谱周期信号的频谱 机械电子工程学院机械电子

30、工程学院周期性方波信号周期性方波信号2.3.2 周期信号的频谱周期信号的频谱 机械电子工程学院机械电子工程学院周期信号的频谱具有下列三个特征：周期信号的频谱具有下列三个特征：1、离散性离散性 周期信号的频谱是离散谱。周期信号的频谱是离散谱。2、谐波性谐波性 周期信号的谱线仅出现在基频及各次谐波频率处。周期信号的谱线仅出现在基频及各次谐波频率处。3、收敛性收敛性 周期信号的幅值谱中各频率分量的幅值随着频率的升高而周期信号的幅值谱中各频率分量的幅值随着频率的升高而减小，频率越高，幅值越小。减小，频率越高，幅值越小。2.3.2 周期信号的频谱周期信号的频谱 机械电子工程学院机械电子工程学院2.3.3

31、 非周期信号的频谱非周期信号的频谱 非周期信号分为准周期信号和瞬变信号。非周期信号分为准周期信号和瞬变信号。当周期信号的当周期信号的周期趋向于无穷大周期趋向于无穷大时，原来的周期信号便可时，原来的周期信号便可当作非周期信号来处理。此时，信号的相邻谱线间隔趋向于当作非周期信号来处理。此时，信号的相邻谱线间隔趋向于无穷小，谱线变得越来越密集，最终成为一条连续的频谱。无穷小，谱线变得越来越密集，最终成为一条连续的频谱。各频率分量的幅值尽管也相应地趋向于无穷小，但这些分量各频率分量的幅值尽管也相应地趋向于无穷小，但这些分量间仍保持着一定的间仍保持着一定的比例关系比例关系。对于非周期信号，需要用对于非周

32、期信号，需要用傅里叶变换傅里叶变换来求其频谱。来求其频谱。机械电子工程学院机械电子工程学院 非周期函数非周期函数 存在傅里叶变换的充分条件是存在傅里叶变换的充分条件是 在区间在区间 上上绝对可积绝对可积，即：，即：信号信号 的傅里叶变换的傅里叶变换 定义为：定义为：对应的傅里叶逆变换为：对应的傅里叶逆变换为：2.3.3 非周期信号的频谱非周期信号的频谱 机械电子工程学院机械电子工程学院 一个非周期函数可分解成频率一个非周期函数可分解成频率f f连续变化的谐波的叠加，连续变化的谐波的叠加，式中式中 是谐波是谐波 的系数，决定着信号的振幅和相位。的系数，决定着信号的振幅和相位。由于不同的频率由于不

33、同的频率f，项中的项中的 是相同的，而只有是相同的，而只有 才反映不同谐波分量的振幅与相位的变化情况，因此，称才反映不同谐波分量的振幅与相位的变化情况，因此，称 为为 的的连续频谱连续频谱。由于。由于 一般为实变量一般为实变量f的的复函数复函数，故，故可写为：可写为：式中的式中的 称为非周期信号的幅值谱，称为非周期信号的幅值谱，称为相位谱。称为相位谱。2.3.3 非周期信号的频谱非周期信号的频谱 机械电子工程学院机械电子工程学院 需要注意的是，尽管非周期信号的幅值谱需要注意的是，尽管非周期信号的幅值谱 与周期与周期信号的幅值谱信号的幅值谱 在名称上相同，但在名称上相同，但 是是连续连续的，而的

34、，而 是是离散离散的。的。此外，两者在此外，两者在量纲量纲上也不一样。上也不一样。与信号幅值量纲一致，与信号幅值量纲一致，而而 的量纲与信号量纲不一致。的量纲与信号量纲不一致。与与 的量纲的量纲一致，一致，是单位频宽上的幅值。因此，严格地说，是单位频宽上的幅值。因此，严格地说，是是频谱密度函数频谱密度函数。2.3.3 非周期信号的频谱非周期信号的频谱 机械电子工程学院机械电子工程学院矩形窗函数的时域表达式为：矩形窗函数的时域表达式为：幅幅-频谱频谱相相-频谱频谱 1、谱线是连续的谱线是连续的，这是，这是瞬变信号与周期信号在谱瞬变信号与周期信号在谱图上的显著区别。图上的显著区别。2、矩形窗的时间

35、长度、矩形窗的时间长度T越长，幅频图中主瓣越高越长，幅频图中主瓣越高而窄，意味着而窄，意味着能量越集中能量越集中于主瓣于主瓣。2.3.3 非周期信号的频谱非周期信号的频谱 机械电子工程学院机械电子工程学院2.3.4 离散傅立叶变换离散傅立叶变换 在计算机上实现在计算机上实现Fourier变换，必须做到：变换，必须做到：1）把连续信号（包括时域和频域）改造为离散数据；）把连续信号（包括时域和频域）改造为离散数据；2）把计算范围收缩到一个有限区间；）把计算范围收缩到一个有限区间；3）实现正、逆傅立叶变换。）实现正、逆傅立叶变换。在这种条件下所构成的变换对，在在这种条件下所构成的变换对，在时域时域和

36、和频域频域都只取有限都只取有限个离散数据，这些数据分别构成个离散数据，这些数据分别构成周期周期性的性的离散时间离散时间函数和函数和离离散频率散频率函数。函数。离散信号离散信号 的傅立叶变换表达式为：的傅立叶变换表达式为：机械电子工程学院机械电子工程学院2.3.4.1 截断、泄露与窗函数截断、泄露与窗函数 截断就是对无限长的信号进行截取，也就是对截断就是对无限长的信号进行截取，也就是对x(t)信号乘以信号乘以矩形窗函数矩形窗函数w(t)。当。当w(t)=0时，乘积的结果时，乘积的结果y(t)=0；当；当w(t)=1时，乘积的结果时，乘积的结果y(t)=x(t)。两个信号在两个信号在时域时域内的内

37、的乘积乘积，对应于这两个信号在，对应于这两个信号在频域频域内的内的卷积卷积（）。）。由于由于w(t)对应的频域函数对应的频域函数W(f)是一个是一个无限带宽无限带宽的的sinc函数，它函数，它与信号与信号x(t)对应的频域函数对应的频域函数X(f)在频域的卷积，必然造成在频域的卷积，必然造成x(t)信号信号的能量沿频率轴扩展开来，这就是所谓的的能量沿频率轴扩展开来，这就是所谓的谱泄漏谱泄漏。频域卷积的结果，将使得在频谱图中出现不属于频域卷积的结果，将使得在频谱图中出现不属于x(t)信号的谱信号的谱线（它们是线（它们是w(t)的谱线）。的谱线）。机械电子工程学院机械电子工程学院2.3.4.1 截

38、断、泄露与窗函数截断、泄露与窗函数 信号时域加矩形窗及其频域变化信号时域加矩形窗及其频域变化机械电子工程学院机械电子工程学院 为了抑制或减小泄漏效应，需要选择性能更好的特殊窗来为了抑制或减小泄漏效应，需要选择性能更好的特殊窗来代替矩形窗，这种处理称为代替矩形窗，这种处理称为加窗加窗处理。处理。加窗的目的，在时域是加窗的目的，在时域是平滑平滑截断信号两端的波形突变，而截断信号两端的波形突变，而在频域则是尽可能地在频域则是尽可能地压低旁瓣的高度压低旁瓣的高度。一般来说，一个好的窗函数其频谱的一般来说，一个好的窗函数其频谱的主瓣应窄，旁瓣应小主瓣应窄，旁瓣应小。主瓣窄意味着能量集中，分辨率高；旁瓣小

39、意味着能量泄漏主瓣窄意味着能量集中，分辨率高；旁瓣小意味着能量泄漏少。少。第一种措施，加大矩形窗的第一种措施，加大矩形窗的时间长度时间长度，即增大采样的样本点，即增大采样的样本点数。也就是使数。也就是使W(f)的主瓣尽量地高而窄，能量最大限度地集中的主瓣尽量地高而窄，能量最大限度地集中于主瓣，将旁瓣尽量压缩。于主瓣，将旁瓣尽量压缩。第二种措施，采用第二种措施，采用旁瓣较低的函数作为采样窗函数旁瓣较低的函数作为采样窗函数，如汉宁，如汉宁窗、海明窗等等。窗、海明窗等等。2.3.4.1 截断、泄露与窗函数截断、泄露与窗函数 机械电子工程学院机械电子工程学院汉宁窗及其幅频特性汉宁窗及其幅频特性2.3.

40、4.1 截断、泄露与窗函数截断、泄露与窗函数 三角窗及其幅频特性三角窗及其幅频特性矩形窗及其幅频特性矩形窗及其幅频特性机械电子工程学院机械电子工程学院2.3.4.1 截断、泄露与窗函数截断、泄露与窗函数 汉宁窗的时域、汉宁窗的时域、频域波形图频域波形图海明窗的时域、频海明窗的时域、频域波形图域波形图矩形窗函数的时矩形窗函数的时域、频域波形图域、频域波形图机械电子工程学院机械电子工程学院 矩形窗矩形窗主瓣最窄，旁瓣则较高，泄漏较大，适合于要获得主瓣最窄，旁瓣则较高，泄漏较大，适合于要获得精确精确主峰的频率主峰的频率、而幅值精度要求不高的场合。、而幅值精度要求不高的场合。汉宁窗汉宁窗旁瓣明显降低，

41、具有抑制泄漏的作用，但主瓣较宽，致旁瓣明显降低，具有抑制泄漏的作用，但主瓣较宽，致使频率分辨能力较差，在截断随机信号，或对周期信号进行非整使频率分辨能力较差，在截断随机信号，或对周期信号进行非整周期截断时，为了平滑或削弱截取信号的两端，减小泄漏，宜加周期截断时，为了平滑或削弱截取信号的两端，减小泄漏，宜加汉宁窗。汉宁窗。指数窗指数窗无旁瓣，主瓣很宽，频率分辨力低，对脉冲响应类信号无旁瓣，主瓣很宽，频率分辨力低，对脉冲响应类信号宜加指数窗，若适当选择衰减系数，可起到抑制噪声的作用。宜加指数窗，若适当选择衰减系数，可起到抑制噪声的作用。2.3.4.1 截断、泄露与窗函数截断、泄露与窗函数 除矩形窗

42、之外，其它的窗函数存在如下的不足：除矩形窗之外，其它的窗函数存在如下的不足：1、初相位信息消失初相位信息消失。采用它们的频谱结果没有相频谱图。采用它们的频谱结果没有相频谱图。2、谱图中的、谱图中的幅值幅值相对实际信号该频率成份的幅值存在着相对实际信号该频率成份的幅值存在着失真失真。失真度的大小与所取的修正值相关。失真度的大小与所取的修正值相关。机械电子工程学院机械电子工程学院2.3.4.2 采样、频混和采样定理采样、频混和采样定理 数字信号处理时，首先要将一个模拟信号转换为一个数数字信号处理时，首先要将一个模拟信号转换为一个数字信号。信号的字信号。信号的采样由模采样由模/数转换电路来实施数转换

43、电路来实施的。的。如果以如果以 代表原始的连续时间信号，代表原始的连续时间信号，代表采样后代表采样后获得的离散信号，则采样信号获得的离散信号，则采样信号 可以看成是原始信号可以看成是原始信号 与与周期脉冲序列周期脉冲序列 的乘积。的乘积。脉冲序列脉冲序列 是一系列周期为是一系列周期为T T的脉冲函数：的脉冲函数：时域采样的数学表达式为：时域采样的数学表达式为：机械电子工程学院机械电子工程学院采样过程在时域和频域的表示采样过程在时域和频域的表示 (a)原函数原函数(b)原函数频谱原函数频谱(c)采样冲击函数采样冲击函数(d)采样冲击函数的频谱采样冲击函数的频谱(e)离散时间信号离散时间信号(d)

44、采样序列的频谱采样序列的频谱2.3.4.2 采样、频混和采样定理采样、频混和采样定理 机械电子工程学院机械电子工程学院 对一个一定长度的模拟信号，若对它的采样间隔小，亦即对一个一定长度的模拟信号，若对它的采样间隔小，亦即采样率高采样率高，则采样的数据量大，要求计算机具有较大内存及，则采样的数据量大，要求计算机具有较大内存及较长的处理时间。较长的处理时间。若采样率过低，即采样间隔大，则系列的离散时间序列不若采样率过低，即采样间隔大，则系列的离散时间序列不能真正反映原始信号的波形特征，在频域处理时会出现能真正反映原始信号的波形特征，在频域处理时会出现频率频率混淆混淆的现象，又称混叠（的现象，又称混

45、叠（aliasing）。）。采样后得到间隔为采样后得到间隔为T的等距脉冲序列，这个序列的包络线应的等距脉冲序列，这个序列的包络线应与原始信号一致。即与原始信号一致。即采样后的信号应能恢复原信号采样后的信号应能恢复原信号，不发生，不发生失真。这主要取决采样间隔失真。这主要取决采样间隔T。2.3.4.2 采样、频混和采样定理采样、频混和采样定理 机械电子工程学院机械电子工程学院采样序列及还原曲线采样序列及还原曲线2.3.4.2 采样、频混和采样定理采样、频混和采样定理 上面两个图的原信上面两个图的原信号的频率较高，采样间号的频率较高，采样间隔隔T T过小，因此采样序列过小，因此采样序列不能复原原信

46、号（实线不能复原原信号（实线表示原始信号曲线，虚表示原始信号曲线，虚线表示采样点描述的信线表示采样点描述的信号曲线）。号曲线）。最下面的信号因为最下面的信号因为频率低，采样信号就能频率低，采样信号就能复原原信号曲线。复原原信号曲线。机械电子工程学院机械电子工程学院 设设 所包含的各信号成分中最高频率为所包含的各信号成分中最高频率为 ，当采，当采样频率样频率低于低于 时，采样得到的离散信号时，采样得到的离散信号频率频率不等于原信号不等于原信号频率。因此，对数字信号处理来说，当一个信号包含多个频率频率。因此，对数字信号处理来说，当一个信号包含多个频率成分时，为避免混叠产生，要求采样频率成分时，为避

47、免混叠产生，要求采样频率 必须高于信号必须高于信号频频率成分中最高频率率成分中最高频率 的的2倍倍，即：，即：这就是这就是采样定理采样定理，也称香农（，也称香农（Shannon）采样定理。）采样定理。在给定的采样频率在给定的采样频率 条件下，信号中条件下，信号中能被分辨的最高频能被分辨的最高频率率称为奈奎斯特（称为奈奎斯特（Nyquist）频率：）频率：2.3.4.2 采样、频混和采样定理采样、频混和采样定理 机械电子工程学院机械电子工程学院 实际进行信号处理时，不可能无限制地提高采样频率，因实际进行信号处理时，不可能无限制地提高采样频率，因此往往在信号进入此往往在信号进入A/D转换器之前先通

48、过一个转换器之前先通过一个模拟低通滤波器模拟低通滤波器，滤除信号中不加以考虑的高频成分，降低信号中的最高频率，滤除信号中不加以考虑的高频成分，降低信号中的最高频率，从而可以降低采样频率。这种用途的滤波器称为从而可以降低采样频率。这种用途的滤波器称为抗混低通滤抗混低通滤波器波器。抗混滤波器的截止频率通常选择为抗混滤波器的截止频率通常选择为等于等于信号中分析的最高信号中分析的最高频率。不管何种形式的滤波器，均不可能具有理想的滤波特频率。不管何种形式的滤波器，均不可能具有理想的滤波特性，在其截止频率之外总还有一段过渡带。因此，在实际中性，在其截止频率之外总还有一段过渡带。因此，在实际中常将采样频率选

49、择为抗混滤波器截止频率的常将采样频率选择为抗混滤波器截止频率的34倍倍。2.3.4.2 采样、频混和采样定理采样、频混和采样定理 机械电子工程学院机械电子工程学院2.3.4.3 量化误差量化误差 采样所得的离散序列值，需用采样所得的离散序列值，需用有限字长的二进制码来表示，这有限字长的二进制码来表示，这一过程称为一过程称为量化量化。模拟信号的幅值是连续的，而模拟信号的幅值是连续的，而数字信号受到位数的限制，其值数字信号受到位数的限制，其值是跳跃的。模拟信号在量化过程是跳跃的。模拟信号在量化过程中，若采样点的幅值落在两相邻中，若采样点的幅值落在两相邻的量化值之间，就要舍入到邻近的量化值之间，就要

50、舍入到邻近的一个量化值上，这造成了的一个量化值上，这造成了量化量化误差误差。减小量化误差只能选用位数更减小量化误差只能选用位数更多的多的A/D转换装置。转换装置。量化误差示意图量化误差示意图机械电子工程学院机械电子工程学院 一个离散傅里叶变换的过程可分为一个离散傅里叶变换的过程可分为时域采样时域采样、时域截断时域截断和和频频域采样域采样3个步骤。个步骤。信号信号 的频谱的频谱 经离散傅里叶变换计算之后，得到经离散傅里叶变换计算之后，得到的的N N根谱线的位置是在根谱线的位置是在 （k=0,1,2,）的地方，的地方，即仅在基频的整数倍的频率点上才有其各个频率成分，所有那些即仅在基频的整数倍的频率

51、点上才有其各个频率成分，所有那些位于离散谱线之间的频谱图形都得不到显示，不能知道其精确的位于离散谱线之间的频谱图形都得不到显示，不能知道其精确的值。值。若信号中某频率成分的频率若信号中某频率成分的频率 等于等于K/T，即它与输出的频，即它与输出的频率采样点相重合，那么该谱线便可被精确地显示出来；若率采样点相重合，那么该谱线便可被精确地显示出来；若 与与频率采样点不重合，便得不到显示，所得的频谱便会产生误差。频率采样点不重合，便得不到显示，所得的频谱便会产生误差。这种现象称为这种现象称为栅栏效应栅栏效应。2.3.4.4 栅栏效应栅栏效应 机械电子工程学院机械电子工程学院信号时域采样及其频域变化信

52、号时域采样及其频域变化2.3.4.4 栅栏效应栅栏效应 机械电子工程学院机械电子工程学院采样后信号时域加矩形窗及其频域变化采样后信号时域加矩形窗及其频域变化2.3.4.4 栅栏效应栅栏效应 机械电子工程学院机械电子工程学院信号频域采样及其时域变化信号频域采样及其时域变化2.3.4.4 栅栏效应栅栏效应 机械电子工程学院机械电子工程学院 在离散傅里叶变换中，将两条谱线间的距离称为在离散傅里叶变换中，将两条谱线间的距离称为频率分辨频率分辨率率 ，谱线间隔越小，频率分辨率便越高，被栅栏效应所漏，谱线间隔越小，频率分辨率便越高，被栅栏效应所漏掉的频率成分便越少。掉的频率成分便越少。当被分析的时域信号长

53、度当被分析的时域信号长度T（即窗宽（即窗宽 ）和采样频）和采样频率率 被确定之后，频率分辨率被确定之后，频率分辨率 也被确定：也被确定：2.3.4.4 栅栏效应栅栏效应 机械电子工程学院机械电子工程学院 对于工程信号来说，一旦根据其分析的频带确定对它的最低对于工程信号来说，一旦根据其分析的频带确定对它的最低采样频率采样频率 之后，为获得足够的频率分辨率，便必须要之后，为获得足够的频率分辨率，便必须要增加数增加数据点数据点数N，由此使计算机的计算机急剧增加。为解决这一问题，由此使计算机的计算机急剧增加。为解决这一问题，通常有不同的途径加以选择，如频率细化（通常有不同的途径加以选择，如频率细化（Z

54、oom）技术、）技术、Z变变换换及现代谱分析等方法。及现代谱分析等方法。对于对于周期信号，周期信号，作作整周期截取整周期截取是获得正确频谱的先决条件。是获得正确频谱的先决条件。对信号做离散傅里叶变换的结果是将用窗函数截取的时域信对信号做离散傅里叶变换的结果是将用窗函数截取的时域信号作号作周期性延拓周期性延拓。如果实施整周期截取，则截取的整周期信号经。如果实施整周期截取，则截取的整周期信号经延拓之后仍为周期信号，没有产生任何畸变。但若不是整周期截延拓之后仍为周期信号，没有产生任何畸变。但若不是整周期截取，被截取的信号经延拓之后将在原先连续的波形上产生间断点，取，被截取的信号经延拓之后将在原先连续

55、的波形上产生间断点，从而造成波形畸变，不能再复现原来的信号，而对应的频谱亦将从而造成波形畸变，不能再复现原来的信号，而对应的频谱亦将发生畸变。发生畸变。2.3.4.4 栅栏效应栅栏效应 机械电子工程学院机械电子工程学院2.3.5 随机信号的频谱随机信号的频谱 2.3.5.1 自功率谱密度自功率谱密度功率信号功率信号 的平均功率可用均方值来表示，即：的平均功率可用均方值来表示，即：如果令：如果令：具有单位频率的平均功率量纲，故称为具有单位频率的平均功率量纲，故称为功率谱密度函数功率谱密度函数。描述信号的平均功率相对于频率的分布情况。描述信号的平均功率相对于频率的分布情况。则平均功率：则平均功率：

56、机械电子工程学院机械电子工程学院2.3.5.1 自功率谱密度自功率谱密度 根据维纳根据维纳-辛钦辛钦（Wiener-Khintchine）公式）公式，平稳随机过，平稳随机过程的程的功率谱密度功率谱密度 与与自相关函数自相关函数 是一对傅立叶变换，是一对傅立叶变换，即：即：通过自相关函数的傅立叶变换就可以得到功率谱密度函数。通过自相关函数的傅立叶变换就可以得到功率谱密度函数。机械电子工程学院机械电子工程学院2.3.5.2 互功率谱密度互功率谱密度 两个随机信号两个随机信号 和和 之间的之间的互谱密度函数互谱密度函数 与与互互相关函数相关函数 构成一对傅立叶变换，即：构成一对傅立叶变换，即：单边互

57、谱密度函数为：单边互谱密度函数为：因为互相关函数为非因为互相关函数为非偶函数，所以互谱函数偶函数，所以互谱函数是一个是一个复数复数。在实际应用中，常用在实际应用中，常用谱密度的幅值和相位来谱密度的幅值和相位来表示，即：表示，即：机械电子工程学院机械电子工程学院2.3.5.3 相干函数和频率响应函数相干函数和频率响应函数 利用互谱密度函数可以定义利用互谱密度函数可以定义相干函数相干函数 及系统的及系统的频率频率响应函数响应函数 ，即：，即：相干函数相干函数（Coherence Function）又称又称凝聚函数凝聚函数，它是，它是在频域内描述两个信号因果关系的一种无因次比例系数，是用在频域内描述

58、两个信号因果关系的一种无因次比例系数，是用来说明两个信号在频域内是否相关的一个判别指标。来说明两个信号在频域内是否相关的一个判别指标。它把两个测点信号之间的它把两个测点信号之间的互谱互谱与各自的与各自的自谱自谱联系起来，用联系起来，用来确定输出信号来确定输出信号 中有哪些频率成分、多大程度上来自输入中有哪些频率成分、多大程度上来自输入信号信号 ，可以了解到输入与输出信号之间的影响程度，这在，可以了解到输入与输出信号之间的影响程度，这在故障原因的识别故障原因的识别方面是很有用的。方面是很有用的。机械电子工程学院机械电子工程学院频谱细化示意图频谱细化示意图2.3.6 频率细化分析频率细化分析机械电

59、子工程学院机械电子工程学院2.3.6 频率细化分析频率细化分析 若若 ，在时域中信号沿时间轴平移一常值，在时域中信号沿时间轴平移一常值 ，则：，则：在频域中信号沿频率轴平移一常值在频域中信号沿频率轴平移一常值 ，则：，则：时移特性时移特性表明：如果信号在时域中表明：如果信号在时域中延迟延迟了时间了时间 ，则其，则其频谱幅值不会改变，而相频谱中各次谐波的频谱幅值不会改变，而相频谱中各次谐波的相移相移为为 ，与，与频率成正比。频率成正比。频移特性频移特性表明：如果频谱函数在频率坐标上表明：如果频谱函数在频率坐标上平移平移了了 ，则其代表的信号波形将与频率为则其代表的信号波形将与频率为 的正、余弦信

60、号相乘，即的正、余弦信号相乘，即进行了进行了调制调制。机械电子工程学院机械电子工程学院 频率细化分析的基本思想是利用频率细化分析的基本思想是利用频移定理频移定理，对被分析信号在，对被分析信号在时域进行时域进行复调制复调制，再，再重新采样重新采样作傅里叶变换，从而在频域得到更作傅里叶变换，从而在频域得到更高的频率分辨率。主要计算步骤如下。高的频率分辨率。主要计算步骤如下。2.3.6 频率细化分析频率细化分析机械电子工程学院机械电子工程学院2.3.6 频率细化分析频率细化分析 1、选用采样频率、选用采样频率 进行采样，得到进行采样，得到N N点离散序列点离散序列 。假设需要细化的频带是中心频率为。

61、假设需要细化的频带是中心频率为 的一个窄带的一个窄带 ，这里的，这里的 和和 分别是以分别是以 为中心频率的窄带的左、右端点为中心频率的窄带的左、右端点频率。频率。机械电子工程学院机械电子工程学院 2、用一个复序列、用一个复序列 乘以乘以 进行复调制，得到进行复调制，得到N N点点新离散复序列新离散复序列 。根据傅里叶变换的频移定理，复调制将。根据傅里叶变换的频移定理，复调制将频率频率原点移到了频率原点移到了频率 处处，即，即 成为新的频率坐标原点。相应的成为新的频率坐标原点。相应的正、负采样频率正、负采样频率 也同样移动了一个量也同样移动了一个量 。2.3.6 频率细化分析频率细化分析机械电

62、子工程学院机械电子工程学院 3、对、对 进行低通滤波得到离散复序列进行低通滤波得到离散复序列 。设。设 为原来为原来信号抗混滤波的截止频率，由于新的序列信号抗混滤波的截止频率，由于新的序列 的频率上限的频率上限 可能高于原序列可能高于原序列 的奈奎斯特频率的奈奎斯特频率 ，由此产生频率混淆。，由此产生频率混淆。因此，需进一步进行因此，需进一步进行低通滤波低通滤波，把围绕，把围绕 的一个窄带的一个窄带 以外以外的所有频率分量都的所有频率分量都滤掉滤掉，消除可能出现的混叠频率成分。，消除可能出现的混叠频率成分。2.3.6 频率细化分析频率细化分析机械电子工程学院机械电子工程学院 4、对、对 进行进

63、行重新采样重新采样，得到离散复序列，得到离散复序列 。若滤波后。若滤波后的总带宽的总带宽 是原采样频率是原采样频率 的的1/D倍，则就有可能把新序倍，则就有可能把新序列列 的采样频率降低到的采样频率降低到1/D，而不会在新的奈奎斯特频率附近，而不会在新的奈奎斯特频率附近产生混叠。产生混叠。实际进行细化分析时，首先必须保证原始信号采样时有足够实际进行细化分析时，首先必须保证原始信号采样时有足够的长度。如果要对原始信号进行的长度。如果要对原始信号进行D倍的细化分析，就得保证原始倍的细化分析，就得保证原始信号的采样长度为信号的采样长度为DN。这样对滤波后的复序列以采样频率。这样对滤波后的复序列以采样

64、频率 进行重抽样，即每隔进行重抽样，即每隔D个点抽取一个数据，得到个点抽取一个数据，得到新的长度为新的长度为N的的复序列。这时新采样序列的时间跨度增长复序列。这时新采样序列的时间跨度增长D倍，频率分辨率也将倍，频率分辨率也将提高提高D倍倍。2.3.6 频率细化分析频率细化分析机械电子工程学院机械电子工程学院 5、对重抽样后的、对重抽样后的复序列复序列 进行复数进行复数FFT变换，即可得到变换，即可得到细化后中心频率为细化后中心频率为 带宽为带宽为 的细化谱。的细化谱。由于由于 是是复序列复序列，变换后的全部数据都是有用的信息，变换后的全部数据都是有用的信息，且以新频率零点（即调制频率且以新频率

65、零点（即调制频率 ）为基准。）为基准。频谱不存在对称性，频谱不存在对称性，频谱上的负频率和正频率成分实质上分别是原始频率低于和高于频谱上的负频率和正频率成分实质上分别是原始频率低于和高于 的分量，的分量，应将它移到原来的正确位置。应将它移到原来的正确位置。复调制细化包括复调制细化包括幅值细化幅值细化与与相位细化相位细化。由于复调制过程中需。由于复调制过程中需通过数字滤波器，会产生附加相移，所以一般要按滤波器的通过数字滤波器，会产生附加相移，所以一般要按滤波器的相位相位特性予以修正特性予以修正，才能得到真实的细化相位谱。，才能得到真实的细化相位谱。2.3.6 频率细化分析频率细化分析机械电子工程

66、学院机械电子工程学院2.3.6 频率细化分析频率细化分析仿真信号的细化分析结果（采样频率仿真信号的细化分析结果（采样频率200Hz，采样点数采样点数8192）机械电子工程学院机械电子工程学院卫星天线传动机构的频谱细化图（中心频率卫星天线传动机构的频谱细化图（中心频率750Hz）（a）原始信号的频谱，（）原始信号的频谱，（b）细化）细化2倍后的频谱图，倍后的频谱图，（c）细化）细化4倍后的频谱图，（倍后的频谱图，（d）细化）细化8倍后的频谱图倍后的频谱图2.3.6 频率细化分析频率细化分析机械电子工程学院机械电子工程学院2.3.7 倒频谱分析倒频谱分析 倒频谱分析也称为二次频谱分析，是检测复杂倒频谱分析也称为二次频谱分析，是检测复杂谱图中周期分谱图中周期分量的有用工具量的有用工具。在语音中语音音调的测定、机械振动中故障监测。在语音中语音音调的测定、机械振动中故障监测和诊断以及排除回波（反射波）等方面均得到广泛的应用。和诊断以及排除回波（反射波）等方面均得到广泛的应用。设信号设信号 的傅里叶变换为的傅里叶变换为 ，功率谱密度函数为，功率谱密度函数为 。所谓倒频谱，就是对功率谱。所谓倒频谱，