一种鲁棒的振动信号特征值计算方法与流程

技术特征：
1.一种鲁棒的振动信号特征值计算方法，其特征在于：包括以下步骤：s1、短时傅里叶变换：对振动信号x进行短时傅里叶变换得到短时傅里叶变换频谱s2、获得均值频谱：对所述短时傅里叶变换频谱按列计算获得均值频谱y；s3、获得均值频谱的均值和标准差：计算获得所述均值频谱y的均值m
y
和标准差s
y
；s4、获得局部极大值点：根据所述均值m
y
和所述标准差s
y
设定阈值t，利用局部极大值法获得所述均值频谱y中大于所述阈值t的局部极大值点p；s5、获得异常突变信号区：所述振动信号x中以所述局部极大值点p为中心，向左、右分别延拓异常点区间范围长度参数tol的长度的区域为异常突变信号区，所述异常突变信号区中的信号为异常突变信号，将所述异常突变信号中的振动信号设定为所述振动信号x的均值，得到去除异常突变信号后振动信号x
new
；s6、获得去除异常突变信号后特征值：根据有效值计算公式估计所述去除异常突变信号后振动信号x
new
的特征值x
rms
。2.根据权利要求1所述的一种鲁棒的振动信号特征值计算方法，其特征在于：步骤s1中，振动信号x＝[x1,x2,x3...,x
l
]，其中，l为所述振动信号x的长度。3.根据权利要求2所述的一种鲁棒的振动信号特征值计算方法，其特征在于：步骤s1中，进行短时傅里叶变换前设定如下参数：移动窗大小window、相邻两帧信号重叠长度overlap、短时傅里叶变换数据长度nfft和异常点区间范围长度参数tol。4.根据权利要求3所述的一种鲁棒的振动信号特征值计算方法，其特征在于：步骤s1中，m＝fix(nfft/2) 1，其中，fix表示向下取整。5.根据权利要求4所述的一种鲁棒的振动信号特征值计算方法，其特征在于：步骤s1中，n＝fix((l-overlap)/(window-overlap))。6.根据权利要求5所述的一种鲁棒的振动信号特征值计算方法，其特征在于：步骤s1中，短时傅里叶变换的方法为：所述振动信号x某一时刻t的振动信号为x(t)，选择一个时频局部化的窗函数g(t)，假定所述窗函数g(t)在一个短时间间隔内是平稳的，移动所述窗函数g(t)，按照如下公式计算出各个不同时刻的频谱数g(t)，按照如下公式计算出各个不同时刻的频谱7.根据权利要求6所述的一种鲁棒的振动信号特征值计算方法，其特征在于：所述各个不同时刻的频谱组成所述短时傅里叶变换频谱8.根据权利要求1所述的一种鲁棒的振动信号特征值计算方法，其特征在于：步骤s2中，所述均值频谱y为：y＝[y1,y2,y3...,y
n
]。9.根据权利要求1所述的一种鲁棒的振动信号特征值计算方法，其特征在于：步骤s4中，所述阈值t为：t＝m
y
a*s
y
，其中a为异常点监测阈值。10.根据权利要求1所述的一种鲁棒的振动信号特征值计算方法，其特征在于：步骤s6中，所述有效值计算公式为：

技术总结
本发明提供一种鲁棒的振动信号特征值计算方法，对振动信号X进行短时傅里叶变换得到短时傅里叶变换频谱按列计算获得均值频谱Y，根据均值频谱Y均值m


技术研发人员：胡勇 彭六保 曾志生
受保护的技术使用者：航天智控（北京）监测技术有限公司
技术研发日：2021.10.22
技术公布日：2022/2/24