多个经验模态分解对振动信号作用的对比_赵楠.pdf

信 息 与 智 能 2023 NO.8 SCIENCE&TECHNOLOGY INFORMATION科技资讯SCIENCE&TECHNOLOGY INFORMATION科技资讯多个经验模态分解对振动信号作用的对比赵楠 石振刚*(沈阳理工大学信息科学与工程学院 辽宁沈阳 110158)摘要：经振动传感器采集到的信号是非线性、非稳定的，这种信号无论是在时域还是频域上都不易分析。所以通过经验模态分解将原始信号分解成为多个本质模态函数（Intrinsic Mode Function，IMF），之后对其进行特征提取等进一步处理。但是经验模态分解存在模态混叠与端点效应的问题，所以文章采用互补集合经验模态分解（Complementary Ensemble Empirical Mode Decomposision，CEEMD）。CEEMD是在进行经验模态分解之前加入多组符号相反的白噪声，这不仅减少了模态混叠，分解出的IMF分量还更精进。这种互补集合经验模态分解有效地处理了所采集的非线性、非稳定性的振动信号。关键词：振动信号 经验模态分解 本质模态函数 互补集合经验模态分解中图分类号：TN911文献标识码：A 文章编号：1672-3791(2023)08-0017-04Comparison of the Effects of Multiple Empirical Mode Decomposition on Vibration SignalsZHAO Nan SHI Zhengang*(School of Information Science and Engineering,Shenyang Ligong University,Shenyang,Liaoning Province,110158 China)Abstract:The signal collected by vibration sensors is nonlinear and unstable,and this kind of signal is difficult to analyze in either time domain or frequency domain.Therefore,the original signal is decomposed into multiple IMFs through empirical mode decomposition and then is further carried out processing such as feature extraction.However,modal aliasing and endpoint effects exist in empirical mode decomposition,so this paper adopts CEEMD.CEEMD is to add multiple groups of white noises with opposite signs before empirical mode decomposition,which not only reduces mode aliasing,but also makes the decomposed IMF component more refined.The complementary ensemble empirical mode decomposition effectively deals with the collected nonlinear and unstable vibration signals.Key Words:Vibration signal;Empirical mode decomposition;Intrinsic mode function;Complementary ensemble empirical mode decomposition经振动传感器采集得到的振动信号，是一个非线性、非平稳的信号。通常于信号的时域或频域上对信号进行分析与特征提取。而振动信号在时域上无法观测出信号的特征，在频域上，能看出是多个信号进行了叠加，这对于信号的分析造成了干扰。经验模态分解（Empirical Mode Decomposision，EMD）是一种适用于非线性非平稳信号的时频分析方法。王若平等人1采用了一种用EMD改进Mel倒谱系DOI：10.16661/ki.1672-3791.2209-5042-2405作者简介：赵楠（1997），女，硕士在读，研究方向为通信工程。通信作者：石振刚（1971），男，博士，副教授，研究方向为图像处理与模式识别、网络与通信，E-mail：szg_。17 2023 NO.8科技资讯SCIENCE&TECHNOLOGY INFORMATION 信 息 与 智 能SCIENCE&TECHNOLOGY INFORMATION科技资讯数的方法，来对低信噪比的交通环境进行特征提取。该方法在一定程度上提高了特征参数的抗噪性，同时提高了后续的识别率。邹瑛珂等人2在改进EMD基础上，对野外的人和车振动信号进行识别。对改进之后的经验模态分解分解出的IMF分量进行特征提取，提高了后续的分类器分类的准确率。曾成3也在地震信号研究中，运用到了CEEMD方法。张恒4在研究无人值守这一课题中，提到EMD可以自适应对振动信号进行分解，且有较高的准确率。CEEMD也被经常应用于轴承故障的识别，NIU D P、WANG J Q5、汪一飞6用CEEMD来判断滚动轴承的故障，涉及的信号为振动信号，这种方法能准确判断出故障原因及位置。孟繁雯等人7在研究矿井主通风机轴承时，也将 EMD 与CEEMD进行对比，发现CEEMD分解效果优于EMD。经验模态分具有自适应性，无须预先对信号进行分析，就会将信号分解成为多个本质模态函数（IMF）。但是经验模态分解有两个缺点，一个是模态混叠，另一个是端点效应。这两个缺点一个会影响信号的特征提取，另一个会造成信号两端的包络严重扭曲，使信号信息遭受损害，影响判断。为了解决模态混叠，在信号进行分解前加入随机的白噪声，也就是EMD的衍生方法集合经验模态分解（EEMD）。但是添加的白噪声存在个体差异，在下一步进行经验模态分解时，会造成新的混叠，同时可能会有无法消除的信号累积在原始信号中，影响信号的分析及特征提取，在对IMF分量成分进行集合平均时，还会造成运算时间的增加，所以此次实验使用互补集合经验模态分解（CEEMD）。互补集合经验模态与集合经验模态的区别在于互补集合经验模态分解加入的白噪声是一对互为相反数的正负白噪声，可有效地降低运算律与消除信号重构后的残余辅助白噪声。1 经验模态分解经验模态分解的原理是通过寻找信号的极大值点雨帘云栋极小值点，分别形成上、下两个将极值点都包含在内的包络，然后取两个包络间的平均值，用原始信号减去平均值，得到一个新的信号，对此信号重复上诉步骤，直到满足IMF的条件，即可得到一个IMF分量，具体运算过程如下。（1）通过三次样条函数拟合出极大值包络线()maxet，以及极小值包络线()minet。（2）计算上、下包络的平均值。()()()maxmin2etetm t+=（1）（3）构造一个新信号()11p t，该信号是原始信号去掉了一个低频信号。()()()11ptx tm t=-（2）但此时得到的信号并不是稳定的，所以继续重复上述步骤，直到满足IMF的定义为止，才算是找到原信号的一阶IMF分量。()()()111kctIMF tpt=（3）用原始信号减去上述得到的IMF分量，得到了一新的去掉高频成分的信号()1r t。()()()11r tx tct=-（4）重复上述步骤，直到EMD停止分解。满足n阶的IMF分量或残余分量()nr t小于预设值，或者残余分量为单调函数或者常量二者条件之一即可。（4）原信号与IMF分量和残余分量的关系如下。()()()1ninix tctrt=+（5）本质模态函数（IMF）有两个约束条件，一个是在整个数据段内，极值点的个数和过零点的个数必须相等或相差最多不能超过一个。另一个条件是在任意时刻，由局部极大值点形成的上包络线和由局部极小值点形成的下包络线的平均值为零，即上、下包络线相对于时间轴局部对称。2 互补集合经验模态分解为了缓解经验模态分解的混叠问题，法国的EMD算法研究小组与Huang共同研究发现在原始信号中加入白噪声，使原始信号的缺失的时间尺度被填补，也就是EEMD。但是由于造成新的混叠以及运算时间的增加，提出了新的运算方法互补集合经验模态分解（CEEMD）。CEEMD中所加入的白噪声数值一致，但符号相反，这种方法消除了白噪声在原始信号中的累积。CEEMD的算法过程如下。（1）在原始信号x(t)进行经验模态分解之前，先添加多组符号相反的白噪声。s+i(t)=x(t)+a+m,i(i=1,2,3)s-i(t)=x(t)+a-m,i(i=1,2,3)（6）（2）对所添加了白噪声的信号进行经验模态分解处理，将得到的多个IMF分量，分别进行平均集合。18信 息 与 智 能 2023 NO.8 SCIENCE&TECHNOLOGY INFORMATIONSCIENCE&TECHNOLOGY INFORMATION科技资讯科技资讯 11immjjaai+=(i=j=1,2,3)11immjjaai-=(i=j=1,2,3)（7）（3）将所得的正负平均集合再次进行平均集合，即得到总体平均的IMF分量。2mmmaaa+-+=（8）经过互补集合经验模态分解的信号的IMF，要比经验模态分解得更加准确，大大地减少了模态混叠，CEEMD的流程图如图1所示。3 振动信号仿真实验的数据有两种，一种是人体脚步的振动信号，另一种是车辆行驶的振动信号。首先，振动传感器采集数据以后，将数据封装在txt文件中，通过Numpy函数中的loadtxt语句对数据文件进行调用读取。之后，在对数据进行模态分解之前，用pywt小波函数，对原始信号进行去噪，将去噪之后的信号再进行模态分解。最后，引用PyEMD函数，对去噪之后的函数分别进行经验模态分解、集合经验模态分解与互补集合经验模态分解。使用matplotlib函数绘画原始信号图、去噪信号图、各个经验模态分解之后的IMF分量图，直观观测哪个算法更适合处理本实验数据。根据图2所示，图2（a）为经验模态分解图，图2（b）为集合经验模态分解图，图2（c）为互补集合经验模态分解图。对比图中这3种分解的前几个IMF分量，能明显看出CEEMD的分解结果要更精准，EMD与EEMD的模态混叠现象要严重。对比这3种分解的后几个低频IMF分量，无明显区别，但是可以看出EMD和CEEMD的IMF数量要少，从而证明了EMD和CEEMD的运算次数要少，也就是程序的运算速度高一些。（a）EMD分解图（b）EEMD分解图（c）CEEMD分解图图2 人体脚步信号分析 原始信号原始信号X X(t t)添加多组相反的白噪声添加多组相反的白噪声第一组正第一组正白噪声白噪声第一组负第一组负白噪声白噪声第二组正第二组正白噪声白噪声第二组负第二组负白噪声白噪声第 第n组正白组正白噪声噪声第 第n组负白组负白噪声噪声EMDEMD第一组正白噪第一组正白噪声 声IMF分量分量第二组正白噪第二组正白噪声 声IMF分量分量第 第n组正白噪组正白噪声 声IMF分量分量第一组负白噪第一组负白噪声 声IMF分量分量第二组负白噪第二组负白噪声 声IMF分量分量第 第n组负白噪组负白噪声 声IMF分量分量正白噪声集合平均正白噪声集合平均负白噪声集合平均负白噪声集合平均+总体平均的总体平均的IMFIMF分量分量图1 CEEMD处理流程图19 2023 NO.8科技资讯SCIENCE&TECHNOLOGY INFORMATION 信 息 与 智 能SCIENCE&TECHNOLOGY INFORMATION科技资讯综上所述，CEEMD算法在人体脚步信号分析中要好一些。图3是