一种声音辅助振动微弱信号增强检测的随机共振滤波器的制作方法

本发明涉及振动微弱信号增强检测技术领域，具体涉及一种声音辅助振动微弱信号增强检测的随机共振滤波器。

背景技术：

轴承的状态监测和故障诊断是保障机械关键设备安全运行和减少停机损失的重要途径。基于振动信号分析的轴承故障诊断方法具有准确度高，易于操作，能够连续监测等优点而被广泛应用。轴承的故障诊断通常包含三个连续的步骤：信号采集和预处理，特征提取，故障模式识别。由于传感器采集的振动信号通常包含较强的背景噪声，因此信号的滤波和增强对于保证故障诊断的精度具有重要的意义。

随机共振是一种能够利用噪声增强微弱周期信号的非线性滤波器。通过调整合适的系统参数，噪声和微弱周期信号能够实现最优匹配，从而实现微弱信号的增强检测。由于随机共振滤波器独特的优点，其被应用于轴承振动微弱信号的增强检测和故障诊断中。振动信号通过安装在轴承座上的加速度传感器采集，其噪声强度和分布在信号采集完成之后是确定的，因此只能通过调整合适的随机共振滤波器模型和势阱参数去实现最优的信号增强。但是，如果噪声的强度远高于随机共振系统的触发阈值，则无论如何调整系统参数都难以实现有效的微弱信号增强。

由此可知，如果采用随机共振滤波器对单通道的振动信号进行处理，则振动微弱信号的增强效果和故障诊断的精确度有限，为解决此问题，研究基于多通道信号的随机共振滤波器以提高振动微弱信号的检测效果值得进一步的探讨。

技术实现要素：

本发明的目的在于提供一种声音辅助振动微弱信号增强检测的随机共振滤波器，并将其用于轴承故障诊断。通过采用振动和声音双通道信号信息融合技术，提高振动微弱信号增强检测的效果和轴承故障诊断的精度。

本发明技术解决方案：一种声音辅助振动微弱信号增强检测的随机共振滤波器，该滤波器的实现步骤如下：

步骤1、分别采用加速度计和麦克风传感器采集轴承的振动和声音信号，将振动和声音信号进行归一化预处理，归一化之后的离散振动和声音信号分别表示为yv[i],i＝1,2,...,mv和ys[j],j＝1,2,...,ms,mv和ms分别为振动和声音信号的采样点数，随后通过共振解调技术得到振动包络信号eyv[i],i＝1,2,...,mv,以及声音包络信号eys[j],j＝1,2,...,ms；

步骤2、建立与振动信号长度相等的滑动窗wk[i],i＝1,2,...,mv，wk[i]沿声音包络信号eys[j]的时间轴进行滑动，即wk[i]＝eys[k+i-1],i＝1,2,...,mv,k＝1,2,…,ms-mv+1，随后构造一个振动和声音包络信号的融合信号fk[i]＝wk[i]+eyv[i],i＝1,2,...,mv；

步骤3、为了获得最优的融合信号fk[i]opt，构造一个评价指标rpk，移动滑动窗并不断计算rpk，当rpk的值最大时，则获得最优的融合信号，rpk的计算公式为：

rpk＝rms(fk[i])×psnr(fk[i],eyv[i]),i＝1,2,...,mv,

式中rms为融合信号fk[i]的均方根值，由下式计算：

psnr为融合信号fk[i]与原振动包络信号eyv[i]的峰值信噪比，由下式计算：

其中fk[i]max为融合信号fk[i]的最大幅值；

步骤4、初始化随机共振滤波器参数，将最优的融合信号fk[i]opt输入随机共振滤波器，调整滤波器的参数得到最优输出x[i]opt，随后计算x[i]opt的频谱，从频谱中诊断轴承是否有故障。

本发明的优点和积极效果为：

(1)本发明所提供的声音辅助振动微弱信号增强检测的随机共振滤波器利用声音信号增强振动微弱信号，提供了一种新的多传感器信息融合的方法。

(2)本发明所提供的声音辅助振动微弱信号增强检测的随机共振滤波器能够在振动信号低信噪比的情况下利用声音的能量对振动信号进行有效增强。

(3)本发明所提供的声音辅助振动微弱信号增强检测的随机共振滤波器用于轴承故障诊断，能够实现轴承振动微弱故障特征频率在强背景噪声下的增强，提高轴承故障诊断准确率。

附图说明

图1为本发明滤波器的实现流程图；

图2为轴承外圈故障振动信号及频谱；

图3为轴承外圈故障振动包络信号及频谱；

图4为滑动窗移动时计算得到的融合信号的rpk指标值趋势与信噪比(snr)趋势对比图；

图5为利用本发明随机共振滤波器通过声音增强轴承外圈故障振动微弱信号的最优输出及频谱；

图6为轴承内圈故障振动信号及频谱；

图7为轴承内圈故障振动包络信号及频谱；

图8为利用本发明随机共振滤波器通过声音增强轴承内圈故障振动微弱信号的最优输出及频谱。

具体实施方式

下面结合附图及具体实施方式对本发明进行进一步说明。

实施例一：

根据本发明附图1所示的实现流程图，滤波器的实施步骤如下：

步骤1、分别采用加速度计和麦克风传感器采集轴承的振动和声音信号，将振动和声音信号进行归一化预处理，归一化之后的离散振动和声音信号分别表示为yv[i],i＝1,2,...,mv和ys[j],j＝1,2,...,ms,mv和ms分别为振动和声音信号的采样点数，随后通过共振解调技术得到振动包络信号eyv[i],i＝1,2,...,mv,以及声音包络信号eys[j],j＝1,2,...,ms；

步骤2、建立与振动信号长度相等的滑动窗wk[i],i＝1,2,...,mv，wk[i]沿声音包络信号eys[j]的时间轴进行滑动，即wk[i]＝eys[k+i-1],i＝1,2,...,mv,k＝1,2,…,ms-mv+1，随后构造一个振动和声音包络信号的融合信号fk[i]＝wk[i]+eyv[i],i＝1,2,...,mv；

步骤3、为了获得最优的融合信号fk[i]opt，构造一个评价指标rpk，移动滑动窗并不断计算rpk，当rpk的值最大时，则获得最优的融合信号，rpk的计算公式为：

rpk＝rms(fk[i])×psnr(fk[i],eyv[i]),i＝1,2,...,mv,

式中rms为融合信号fk[i]的均方根值，由下式计算：

psnr为融合信号fk[i]与原振动包络信号eyv[i]的峰值信噪比，由下式计算：

其中fk[i]max为融合信号fk[i]的最大幅值；

步骤4、初始化随机共振滤波器参数，将最优的融合信号fk[i]opt输入随机共振滤波器，调整滤波器的参数得到最优输出x[i]opt，随后计算x[i]opt的频谱，从频谱中诊断轴承是否有故障。

实施例二：

为了验证本发明提出的滤波器的有效性，对轴承故障信号进行分析。实验轴承型号为6205轴承，轴承的振动与声音信号分别由安装在轴承座上的加速度传感器和麦克风传感器进行采集，采样频率设置为20khz，轴承转速为1300rpm。根据轴承的参数和转速可以算出当轴承发生外圈故障时其故障特征频率fbpfo为78.1hz，当轴承发生内圈故障时其故障特征频率fbpfi为117.2hz。首先对轴承外圈故障振动信号进行分析，轴承振动信号及其频谱如附图2所示，从波形图中可见轴承振动信号特征不明显，从频谱图中也难以看出和故障有关的频率成分。对振动信号进行共振解调分析，得到的包络信号和包络谱如图3所示，从包络谱中可以看到78.1hz的故障频率，如图中箭头所指出。但是fbpfo周围的噪声很明显，影响对轴承故障的准确判断。

根据本发明滤波器的实施步骤2，建立一个与振动信号长度相等的滑动窗并沿声音包络信号的时间轴进行滑动，随后构造一个振动和声音包络信号的融合信号。根据本发明滤波器的实施步骤3，构造一个评价指标rpk，移动滑动窗并不断计算rpk，结果如图4上所示。当t＝0.375秒时，rpk达到最大值，则此刻的融合信号为最优融合信号。为了进一步验证rpk指标的有效性，将该指标与信噪比(snr)指标进行对比，snr由下式计算：

式中p[j]为融合后信号的功率谱，fcf为故障特征频率，δf为频率分辨率，nfft为快速傅里叶变换的数据长度。信噪比随着滑动窗移动的趋势如图4下所示。对比图4上下可见，本发明构造的rpk指标与snr指标具有类似的变化趋势，值得注意的是，snr的计算需要知道精确的轴承故障频率(fcf)信息，而本发明构造的rpk指标完全是自适应的。因此通过rpk指标可以在故障频率未知的情况下构造出最优的声音和振动融合信号。

随后，根据本发明滤波器实现步骤4，将最优的融合信号输入随机共振系统，调整系统参数得到最优的滤波信号输出及频谱如图5所示。从频谱中可见，故障特征频率fbpfo＝78.1hz十分突出，其他噪声分量跟fbpfo相比能量很低，因而能够准确地判断出轴承存在外圈故障，以上实施例验证了本发明的随机共振滤波器的效果。

实施例三：

为进一步验证本发明滤波器的有效性，对轴承内圈故障信号进行处理。内圈故障振动信号的波形和频谱如图6所示，从波形和频谱中都可以发现信号的噪声干扰很大，无法从信号中得到有用的信息来诊断轴承故障。对振动信号进行包络解调，其包络信号和包络谱如图7所示，从频谱中可以发现故障特征频率fbpfi＝117.2hz几乎被淹没在了背景噪声中。采用本发明的随机共振滤波器对包络信号进行处理，在轴承的声音信号的辅助下，振动微弱信号得到增强，结果如图8所示。

实施例四：

为了验证本发明的声音辅助振动微弱信号增强检测滤波器的优越性，将本发明的滤波器和传统的只处理振动信号的随机共振滤波器进行对比试验。试验共进行200组，其中外圈和内圈故障信号各测试100组。若经过随机共振滤波器处理后信号的信噪比得到提升，则认为该组试验成功，反之，若滤波后信号信噪比不变或者降低，则认为该组试验失败。试验的结果如表1所示，可见，本发明的滤波器的试验成功次数明显高于传统的随机共振滤波器的成功次数。该结果表明通过合适的声音信号辅助，振动微弱信号在随机共振滤波器框架下能得到更好地增强，从而提高了轴承故障诊断的准确性。

表1定量分析对比试验结果

本发明说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。尽管为说明目的公开了本发明的实施例和附图，但是本领域的技术人员可以理解：在不脱离本发明及所附的权利要求的精神和范围内，各种替换、变化和修改都是可能的。因此，本发明不应局限于本发明实施例和附图所公开的内容。