基于频域相似度的局部均值分解端点效应改进方法
【技术领域】
[0001] 本发明属于自适应信号处理与分析技术领域，具体涉及一种基于频域相似度的局 部均值分解端点效应改进方法的设计。
【背景技术】
[0002] 旋转机械在发生故障时，其振动观测信号包含多种频率成分，且其幅值和频率受 到周期性冲击力的调制，整体呈现非线性和非平稳性。为了实现准确的故障诊断，如何从振 动观测信号中提取出故障诊断所需的特征信号是一个重点，也是一个难点。局部均值分解 (Local Mean Decomposition，LMD)是依据信号本身的信息进行自适应分解，产生的乘积函 数(Product Function，PF)具有真实的物理意义，由此得到的时频分布能够准确地反映出 信号能量在空间各尺度上的分布规律，因此LMD是一种具有自适应性的先进信号处理方法， 适用于处理非线性、非平稳信号，该方法已经在很多领域取得了不错的结果。但是，该方法 的一个重要局限性在于分解过程中存在的端点效应。
[0003] 局部均值分解方法与经验模态分解方法都存在端点效应问题。当采用局部均值分 解方法从待分析信号中提取乘积函数分量时，需要不断地平滑由信号局部极值点构成的局 部均值函数线与包络估计函数线。若信号的端点是极值点，则对信号分解没有不良影响。若 信号的端点既不是局部极大值也不是局部极小值，这样便会造成平滑过程中局部均值函数 线与包络估计函数线在数据两端出现发散现象。这种发散现象会随着信号筛分和迭代过程 逐渐向内"污染"整个分析信号，使得分解结果中的误差不断累积，最终导致分解结果出现 失真。这一现象就是端点效应。待分析信号频率越低，端点效应的误差影响范围越大。随着 LMD分解过程的进行，高频成分不断被筛选出来，余量函数的极值点数目不断减少，端点与 相邻极值点之间的距离进一步增大，误差由信号两端不断向信号中部扩散，导致信号分解 结果出现偏差或失真，这是导致信号分解性能下降的重要原因。
[0004] 为了解决局部均值分解方法中的端点效应问题，由端点延拓出的信号应尽可能符 合信号本身的时频特征，才能避免由端点处误差引发的分解结果失真，提高信号分解方法 的性能。
【发明内容】
[0005] 本发明的目的是为了解决现有技术中局部均值分解方法的端点效应问题，减少由 于信号处理方法引入信号中的误差，提高信号分解方法的性能，提出了一种基于频域相似 度的局部均值分解端点效应改进方法。
[0006] 本发明的技术方案为:一种基于频域相似度的局部均值分解端点效应改进方法， 包括以下步骤：
[0007] S1、信号分段:提取信号的零点信息，按照零点数相等原则将待分析信号分成若干 段，每段信号近似为分析信号的一个周期或部分周期；
[0008] S2、波形重复:对每段信号进行时域波形重复；
[0009] S3、频谱比较:对每段信号进行傅里叶变换，然后计算各端点处信号段与其他信号 段的频域相似度；
[0010] S4、波形左延拓:对比起始端信号段与其他信号段的频域相似度，标记具有最大频 域相似度对应的信号段，利用其相邻左侧信号段对信号进行左延拓，即对信号起始端向前 延拓；
[0011] S5、波形右延拓:对比终止端信号段与其他信号段的频域相似度，标记具有最大频 域相似度对应的信号段，利用其相邻右侧信号段对信号进行右延拓，即对信号终止端向后 延拓；
[0012] S6、信号分解:采用局部均值分解方法处理S4和S5中延拓后的信号，得到初步分解 结果；
[0013] S7、信号截取:截去初步分解结果延拓的部分，得到最终信号分解结果。
[0014] 进一步地，S2中时域波形重复次数为分析信号段数目减1。
[0015] 进一步地，S3中在进行傅里叶变换时分析点数选择按照数据点数较多的信号段进 行设置。
[0016] 进一步地，S3中频域相似度的计算公式为：
[0017]
[0018] 其中Cj(F)代表信号段Cj(t)的傅里叶变换;频域相似度的取值范围为[0，1];频域 相似度数值为零说明两个信号段的频谱完全不同;频域相似度数值越接近1，表明两个信号 段Cj(t)和Cj+l(t)的频谱相似度越高。
[0019] 进一步地，S3中计算频域相似度时，对于起始端，应计算待分析信号的起始端信号 段与其余信号段的频域相似度，包括待分析信号的终止端信号段;对于终止端，应计算待分 析信号的终止端信号段与其余信号段的频域相似度，包括待分析信号的起始端信号段。
[0020] 进一步地，S4和S5中自动搜索与起始端或终止端具有最大频域相似度的信号，若 有多于一个信号段与比较信号段具有相同的频域相似度数值，则任意选择其中一个信号段 进行后续步骤。
[0021] 进一步地，S6具体为：
[0022]确定信号所有的局部极值点，计算相邻两个极值点的平均值和局部幅值，采用滑 动平均法进行平滑处理可得局部均值函数和包络估计函数;从待分析信号中提取乘积函 数，并重复本步骤，直至余量信号为一个单调函数停止;所得每一个乘积函数为一个单分量 的调幅调频信号。
[0023] 进一步地，S7中截去的数据点数与S4和S5中延拓的数据点数相同。
[0024] 本发明的有益效果是:本发明提出的基于频域相似度的局部均值分解端点效应改 进方法，利用信号内部的波形进行延拓，使得延拓后的信号保持原信号的波动规律和时频 特征。通过信号延拓移除局部均值分解方法中存在的端点效应问题。对于具有周期性和准 周期性的信号，采用该方法可实现快速自适应波形延拓。该方法可用于旋转机械装备振动 信号处理与故障诊断领域，亦可用于其他具有类似问题的信号处理中。
【附图说明】
[0025] 图1为本发明提供的基于频域相似度的局部均值分解端点效应改进方法流程图。
[0026] 图2为本发明实施例仿真信号波形与延拓信号波形对比图。
[0027] 图3为本发明实施例含有轴承外圈故障的牵引电机示意图。
[0028] 图4为本发明实施例由牵引电机采集到的轴承振动观测信号时域波形图。
[0029]图5为本发明实施例轴承振动信号的起始端信号段与其他信号段的频谱及频域相 似度数值示意图。
[0030] 图6为本发明实施例轴承振动信号的延拓结果及分解出的信号分量示意图。
[0031] 图7为本发明实施例原始振动观测信号和延拓信号的频谱图。
[0032] 图8为本发明实施例轴承振动信号分解出的乘积函数分量的包络谱低频细节图。
【具体实施方式】
[0033]下面结合附图对本发明的实施例作进一步的说明。
[0034]本发明提供了一种基于频域相似度的局部均值分解端点效应改进方法，包括以下 步骤：
[0035] S1、信号分段:提取信号的零点信息，按照零点数相等原则将待分析信号分成若干 段，每段信号近似为分析信号的一个周期或部分周期。
[0036] 提取信号x(t)的零点信息，按照零点数相等原则将待分析信号x(t)分成若干段， 每段信号(^(〇(1 = 1，…，N，N为段数)近似为x(t)的一个周期或部分。以图2(a)中的仿真信 号为例，起始端为ts = 0s，终止端为te = 0. ls，该信号分为10段。
[0037] S2、波形重复:对每段信号进行时域波形重复。
[0038] 由于S1中分出的每段信号仅为待分析信号的一小部分，数据点数过少，在该步骤 中需将每段信号进行时域波形重复，重复的次数为N-1。
[0039] S3、频谱比较:对每段信号进行傅里叶变换，然后计算各端点处信号段与其他信号 段的频域相似度。
[0040] 由于每个信号段的数据点数不同，在进行傅里叶变换时分析点数选择按照数据点 数较多的信号段进行设置。
[0041] 频域相似度的计算公式为：
[0042]
[0043] 其中Cj(F)代表信号段Cj(t)的傅里叶变换;频域相似度的取值范围为[0，1];频域 相似度数值为零说明两个信号段的频谱完全不同;频域相似度数值越接近1，表明两个信号 段Cj(t)和Cj+l(t)的频谱相似度越高。
[0044] 计算频域相似度时，对于起始端，应计算待分析信号的起始端信号段与其余信号 段的频域相似度，包括待分析信号的终止端信号段;对于终止端，应计算待分析信号的终止 端信号段与其余信号段的频域相似度，包括待分析信号的起始端信号段