[发明专利]针刺机传动机构振动信号特征提取方法

摘要

本发明提出了一种基于MEEMD与指数小波阈值结合的针刺机传动机构振动信号处理方法。将采集的针刺机传动机构振动信号通过波形匹配的自适应端点延拓法进行端点延拓处理，然后进行小波分解，利用指数小波阈值函数对采集的信号进行降噪并重构，得到降噪后的信号，用三次B样条法对降噪后的振动信号进行包络线拟合再对信号进行MEEMD分解，得到若干IMF分量与降噪后信号的互相关系数和峭度值，对IMF和原信号进行相关分析，选择有用的IMF，最后对IMF分量进行希尔伯特变换，得到信号的希尔伯特谱，为针刺机传动机构振动信号的信号特征。本方法适用于工业领域针刺机传动机构的故障诊断。

主权项

1.一种针刺机传动机构振动信号特征提取方法，其特征在于将采集的针刺机传动机构振动信号通过波形匹配的自适应端点延拓法进行端点延拓处理，然后进行小波分解，利用指数小波阈值函数对采集的信号进行降噪并重构，得到降噪后的信号，用三次B样条法对降噪后的振动信号进行包络线拟合再对信号进行MEEMD分解，得到若干IMF分量与降噪后信号的互相关系数和峭度值，对IMF和原信号进行相关分析，选择有用的IMF，最后，对IMF分量进行希尔伯特变换，得到信号的希尔伯特谱，即为针刺机传动机构振动信号的信号特征；/n针刺机传动机构振动信号特征提取方法的具体步骤如下：/n(1) 合理布置传感器，采集针刺机传动机构运行时的振动信号；/n(2) 振动信号两端采用自相关镜像法对其进行延拓；/n1)设信号长度为，对应的时间依次为；/n其中的左端点是（不是极值点），极值点为（从左到右），对应的时间为，过这四个点的波形组成一个特征波形，记为A，长度为L；/n其中包含至少一个极大值、极小值和过零点；/n2) 将A进行自相关计算，公式如下：/n /n上式中的N为信号长度，n为时延数；/n3) 取的镜像，延拓至信号的左端点，延拓后信号的左端点为一个极值点；/n4) 右端点同上进行延拓；/n(3) 指数小波阈值降噪/n /n其中且为正数/n1)选择一个恰当的小波并确定小波分解的层数N，然后对含噪信号进行N层小波分解，从而得到小波的分解系数；/n2）对分解后的小波系数进行阈值量化，选择指数阈值函数对含噪信号进行降噪处理，得到估计小波系数；/n3）对降噪后的信号进行重构即对处理过的小波系数进行小波逆变换，得到估计信号即降噪后的信号；/n(4) 对降噪的针刺机传动机构振动信号进行MEEMD分解之前，采用三次B样条插值法拟合包络线，具体步骤如下：/n由各控制点连接构造成一个控制多边形，设控制多边形的边长为，则多边形周长为，那么次B样条曲线的节点矢量可写为：/n /n为了满足B样条曲线的局部性质要求，计算时控制多边形的边由其相应条边的和来代替，然后对其规范化；/n计算出定义域内节点区间长度：/n /n以给定信号极值点为控制顶点，次B样条插值曲线控制多边形相应条边的和为：/n /n 次B样条插值曲线的节点矢量规范化分母为：/n /n由以上分析可得次B样条拟合极值曲线所有节点值为：/n /n则参数化后的节点矢量可以表示为：/n /n由三次B样条曲线定义并将节点值代入方程可知拟合包络方程为：/n /n这样便求得包络线方程；/n(5) 进行MEEMD分解，步骤如下：/nMEEMD分解是在原始信号中加入一组幅值和标准差一致、方向相反的白噪声，再重复进行EEMD和EMD分解，以减少EMD分解的模态混淆和EEMD分解的计算量与残余分量；具体步骤如下：/n1）：向原始信号中所加入的白噪声均方根值应接近待分解的内部噪声，或不超过待分解信号均方根值的 0.3 倍；/n2）：设向待分解信号中所加入的绝对值相等、方向相反的一组白噪声为和，分别对其进行EEMD分解：/n /n上式中的和为EEMD分解过程中得到的 2 组 IMFs分量；/n3）：取和的平均值：/n /n4）：将以下列的EMD分解方式进行分解：/n /n上式中表示经过EMD分解后得到的IMF分量；表示剩余残余分量的叠加；/n5）：经过MEEMD的分解后，原始信号可表示为：/n /n上式中的表示最终得到的IMF分量；表示最终得到的残余分量；/n (6) 将MEEMD分解得到的IMF分量进行Hilbert变换，得到信号的Hilbert谱，即为针刺机传动机构振动信号的信号特征。/n