[发明专利]基于多小波阶比双谱分析的旋转机械故障诊断方法

权利要求书

1.基于多小波阶比双谱分析的旋转机械故障诊断方法，其特征在于，包括如下步骤：

步骤一：通过把原始信号输入到塔式算法中进行多小波变换；

步骤二：根据信号分析理论，寻找使分解信号能量最大的基元函数组合；

步骤三：通过连接算法使得所有基元函数信号的总能量最大，进而求得转速信号；

步骤四：对转速信号进行等角度重采样，对重采样信号进行双谱分析得到振动信号的阶比谱；

步骤五：对振动信号的阶比谱进行分析，从振动信号中进行旋转设备的故障诊断。

2.根据权利要求1所述的基于多小波阶比双谱分析的旋转机械故障诊断方法，其特征在于，步骤一的具体方法如下：

步骤101：对采集到的旋转设备转速信号进行过样预处理，将原始信号f(n)变为r个矢量信号，

vi,j=vi,j0vi,jl=f(n)0---(1)]]>

式中：n表示时间节点；i表示子空间序列的行号，i＝1,2,3,4....r；j表示子空间序列的列号，j＝1,2,3,4....r；r代表矢量信号个数，r∈n,且r≥2；v_i,j表示对原始信号进行小波变换后的第i行第j列的子空间序列；表示对原始信号进行小波变换后的第0个输出块；表示对原始信号进行小波变换后的第l个输出块；

步骤102，将r个矢量信号输入到塔式算法中进行多小波变换，多小波变换包括分解和重构，具体公式如下：

分解公式为：H(w)＝H(w)Q(w)，重构公式为：G(w)＝G(w)P(w)； (2)

其中：w代表角频率，H_k，G_k为r×r维系数矩阵；0≤k≤L；L为滤波器的长度；此处i表示虚数单位，i2＝-1；Q(w)为多小波分解算法之前的预滤波器，P(w)为多小波重构算法之后的滤波器。