一种机械振动信号的倒频谱优化计算方法和系统技术方案

本发明专利技术实施例公开一种机械振动信号的倒频谱优化计算方法和系统。在一具体实施方式中，该方法包括：对获取的机械振动信号进行补零延展；对补零延展后的机械振动信号进行快速傅里叶变换，获得机械振动信号的单边幅频函数；对单边幅频函数进行去噪声处理，并依次进行取对数操作和补零延展；对补零延展后的去噪声处理的单边频谱函数的对数函数进行快速傅里叶变换，获得机械振动信号的倒频谱。本发明专利技术对倒频谱的优化符合倒频谱的数学定义，结果准确有效，同时本发明专利技术步骤清晰简洁，适合计算机编程实现，可在工程实施中发挥有效作用。可在工程实施中发挥有效作用。可在工程实施中发挥有效作用。

【技术实现步骤摘要】
一种机械振动信号的倒频谱优化计算方法和系统


[0001]本专利技术涉及倒频谱优化计算方法。更具体地，涉及一种机械振动信号的倒频谱优化计算方法和系统。

技术介绍

[0002]目前，倒频谱是一种监测复杂谱图中周期分量的有效工具，在数学上，一个信号函数的倒频谱函数被定义为该信号函数的傅里叶变换的对数的傅里叶变换。倒频谱可以将信号的频谱上的众多边带谱线凝聚为单根谱线，从而使信号分析变得简单。
[0003]倒频谱分析方法在机械故障诊断、地震预报、语音分析等多个领域获得了广泛的应用。在机械故障诊断中，当齿轮或轴承等关键零部件发生故障时，通过振动传感器所敏感的振动信号会出现调制现象，信号调制的结果，在其频谱上的主频周围形成众多的规律难以识别的边频，此时利用倒频谱分析，则可方便提取边频的周期成分从而快速定位故障；另外对于布置在不同位置的振动传感器，由于传递路径不同，其振动信号的频谱也有差异，但倒频谱可以将信号的路径效应和传递效应分离，即针对同一检测对象的倒频谱分析结果对传感器的位置不敏感。以上两个优点使得倒频谱对机械故障识别极为有用。
[0004]机械故障诊断领域的工程应用上，大多机械故障诊断系统以及通用工程软件均集成了倒频谱分析功能。但由于倒频谱分析的具体实现是基于其数学定义进行的，其实现方式涉及信号处理及软件设计等知识，不同的机械故障诊断系统以及通用工程软件所实现的倒频谱分析的效果不尽相同，目前未见到有关倒频谱分析具体实现的通用方法或优化方法。

技术实现思路

[0005]本专利技术的目的在于提供一种机械振动信号的倒频谱优化计算方法和系统，以解决现有技术存在的问题中的至少一个。
[0006]为达到上述目的，本专利技术采用下述技术方案：
[0007]本专利技术第一方面提供一种机械振动信号的倒频谱优化计算方法，包括：
[0008]对获取的机械振动信号进行补零延展；
[0009]对补零延展后的机械振动信号进行快速傅里叶变换，获得机械振动信号的单边幅频函数；
[0010]对单边幅频函数进行去噪声处理，并依次进行取对数操作和补零延展；
[0011]对补零延展后的去噪声处理的单边频谱函数的对数函数进行快速傅里叶变换，获得振动信号的倒频谱。
[0012]可选地，所述获取机械振动信号并对其进行补零延展包括
[0013][0014]其中，x0(t)为机械振动信号，N
x
为机械振动信号数据量，
[0015]K
x
N
x
为补零延展后的机械振动信号数据量，K
x
为整数，且K
x
>1。
[0016]可选地，所述机械振动信号的单边幅频函数为
[0017]y0(f)＝FFT(x(t))
[0018]其中FFT(
·
)为快速傅里叶变换。
[0019]可选地，进行去噪声处理的单边幅频函数为
[0020][0021]其中，单边幅频函数y0(f)的频率自变量f∈[f
a
,f
b
]且幅值不小于y
m
，f
a
>0，f
b
<f
s
/2，y
m
为常数，f
s
为机械振动信号的采样频率，机械振动信号的单边幅频函数的频率自变量f的定义域为f∈[0,f
s
/2]。
[0022]可选地，去噪声处理后的单边频谱函数进行取对数操作得到的对数函数为
[0023]z0(f)＝log(y(f))。
[0024]可选地，对去噪声处理的单边频谱函数的对数函数进行补零延展后为
[0025][0026]其中K
z
N
z
补零延展后的机械振动信号数据量，K
z
为整数，且K
z
>1。
[0027]可选地，所述机械振动信号的倒频谱函数为
[0028]C(q)＝FFT(z(f))
[0029]其中倒频谱函数C(q)的倒频率自变量q的定义域为q∈[0,K
x
T]，T为机械振动信号x0(t)经历时间，T＝N
x
f
s
。
[0030]本专利技术第二方面提供一种机械振动信号的倒频谱优化计算系统，包括
[0031]第一单元，用于对获取的机械振动信号进行补零延展；
[0032]第二单元，用于对补零延展后的机械振动信号进行快速傅里叶变换，获得机械振动信号的单边幅频函数；
[0033]第三单元，用于对单边幅频函数进行去噪声处理；
[0034]第四单元，用于对噪声处理后的单边幅频函数进行取对数操作，得到对数函数；
[0035]第五单元，用于对去噪声处理的单边频谱函数的对数函数进行补零延展；
[0036]第六单元，用于对补零延展后的去噪声处理的单边频谱函数的对数函数进行快速傅里叶变换，得到机械振动信号的倒频谱。
[0037]本专利技术第三方面提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现本专利技术第一方面提供的方法。
[0038]本专利技术第四方面提供一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现本专利技术第一方面提供的方法。
[0039]本专利技术的有益效果如下：
[0040]本专利技术对倒频谱的优化符合倒频谱的数学定义，结果准确有效，同时本专利技术步骤清晰简洁，适合计算机编程实现，可在工程实施中发挥有效作用。
附图说明
[0041]下面结合附图对本专利技术的具体实施方式作进一步详细的说明。
[0042]图1示出本专利技术的一个实施例可以应用于其中的示例性方法流程图。
具体实施方式
[0043]为了更清楚地说明本专利技术，下面结合实施例和附图对本专利技术做进一步的说明。附图中相似的部件以相同的附图标记进行表示。本领域技术人员应当理解，下面所具体描述的内容是说明性的而非限制性的，不应以此限制本专利技术的保护范围。
[0044]具体的一个实施例，
[0045]一种机械振动信号的倒频谱优化计算方法，包括：
[0046]第一步：对获取的机械振动信号进行补零延展；
[0047]一般的，机械设备的振动信号是由振动传感器敏感生成的，振动传感器的后端信号采集单元将振动信号具体化为数字采样的时间序列信号。假设经采样后的机械振动信号为x0(t),设机械振动信号x0(t)数据量为N
x
，记
[0048][0049]对机械振动信号进行补零延展，使得补零扩展后的信号数据量为K
x
N
x
，规定K
x
为整数，且K
x
>1。记补零扩展后的机械振动信号为x(t)，则本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种机械振动信号的倒频谱优化计算方法，其特征在于，包括：对获取的机械振动信号进行补零延展；对补零延展后的机械振动信号进行快速傅里叶变换，获得机械振动信号的单边幅频函数；对单边幅频函数进行去噪声处理，并依次进行取对数操作和补零延展；对补零延展后的去噪声处理的单边频谱函数的对数函数进行快速傅里叶变换，获得机械振动信号的倒频谱。2.根据权利要求1所述的方法，其特征在于，所述补零扩展后的机械振动信号为其中，x0(t)为机械振动信号，N
x
为机械振动信号数据量，K
x
N
x
为补零延展后的机械振动信号数据量，K
x
为整数，且K
x
>1。3.根据权利要求2所述的方法，其特征在于，所述机械振动信号的单边幅频函数为y0(f)＝FFT(x(t))其中FFT(
·
)为快速傅里叶变换。4.根据权利要求3所述的方法，其特征在于，进行去噪声处理的单边幅频函数为其中，单边幅频函数y0(f)的频率自变量f∈[f
a
,f
b
]且幅值不小于y
m
，f
a
>0，f
b
<f
s
/2，y
m
为常数，f
s
为机械振动信号的采样频率，机械振动信号的单边幅频函数的频率自变量f的定义域为f∈[0,f
s
/2]。5.根据权利要求4所述的方法，其特征在于，去噪声处理后的单边频谱函数进行取对数操...

【专利技术属性】
技术研发人员：刘敏，张慧，张龙龙，
申请(专利权)人：北京无线电测量研究所，
类型：发明
国别省市：