[发明专利]一种自适应谱峭度信号处理方法

说明书

本发明公开了一种自适应谱峭度的信号处理方法，属于数字信号处理领域。本发明通过将信号转换到频域，在频域的幅值谱上找出极大值点，根据极大值点设计滤波器，再从时域对信号进行滤波，得到一系列子函数，最后计算子函数的峭度值。根据谱峭度易受噪声和信号幅值影响的特点，本发明将较大幅值的信号成分优先滤除，计算得到的峭度值能够更好反映信号中的瞬态成分。

技术领域

本发明涉及一种非平稳信号处理的方法，具体地涉及一种自适应谱峭度的信号处理方法，属于数字信号处理领域。

背景技术

工程实践中存在着大量的非平稳信号，谱峭度就是一种处理非平稳信号的方法。谱峭度的概念最早是由Dwyer等人在1983年提出，Antoni等人于2006年详细介绍了谱峭度的定义、算法及其应用条件，并先后提出了基于短时傅里叶变换(STFT)的谱峭度法以及谱峭度的快速算法。

谱峭度是在峭度基础上发展而来的。峭度是归一化的四阶中心矩，它是一个反映信号分布特性的数值统计量。峭度指标对冲击信号十分敏感，特别适用于设备早期的故障诊断。但是，峭度是对信号的一个全局性描述，不能准确定位冲击信号。谱峭度是能量归一化的四阶谱累积量，可以度量一个过程在某一频率上的概率密度函数的峰值大小，对非平稳信号非常敏感，能够对时频分解后的各个频带逐个进行分析，找到谱峭度值最大的频带。谱峭度可以实现对冲击信号等瞬态成分的准确定位。

目前，STFT法、滤波器组法和小波变换法及其改进方法是实现谱峭度的主要方法，比如：基于Wigner-Ville分布的谱峭度法、小波变换与EMD相结合的谱峭度法、EEMD谱峭度与EMD相结合的谱峭度法等。但是，这些方法缺乏自适应性。2011年，Wang等人提出了一种以最大谱峭度为准则，不断外扩给定窗函数的自适应谱峭度方法(Wang Y X，LiangM.Identification of multiple transient faults based on the adaptive spectralkurtosis method，Journal of Sound and Vibration，2011，331(2)：470-486)。该方法的不足之处是在瞬态脉冲衰减较慢的情况下，计算得到的峭度值不能很好的与信号瞬态特性相匹配。2013年，Luo等人提出了基于小波变换可调Q的峭度引导自适应解调技术轴承故障检测方法，该方法有时会产生误判断(Luo J S，Yu D J，Liang M.A kurtosis-guidedadaptive demodulation technique for bearing fault detection based on tunable-Q wavelet transform，Measurement Science and Technology，2013，24(2)：1-11)。如何自适应的得到信号的谱峭度，以及如何自适应的用峭度值来反映信号内部隐藏的特定成分，也就是说计算得到的谱峭度值能够很好的反映信号的真实情况，仍然是这一领域面临的一大难题。

发明内容

本发明的目的是为了解决现有谱峭度方法不能自适应反映信号特性的缺点，提供一种自适应的谱峭度信号处理方法。

本发明的自适应谱峭度信号处理方法，具体包括如下步骤：

步骤1，对给定的信号进行快速傅里叶变换(FFT)，将信号从时域变换到频域，得到信号的幅值谱；

步骤2，利用步骤1的幅值谱，在幅值谱上找出极大值点；

步骤3，在步骤2的基础上，对每一个极大值点添加一个窗函数；

步骤4，对步骤3的窗函数进行合并；

步骤5，对步骤4得到的合并后的窗函数，设计带通滤波器；

步骤6，利用步骤5设计的带通滤波器对信号进行滤波处理，得到带通子信号；

步骤7，从滤波前的信号中减去带通滤波得到的带通子信号，得到剩余信号；