本发明专利技术公开了一种自适应总体平均经验模式分解EEMD筛选次数确定方法，该方法对每一个本征模式函数IMF采用不同的筛选停止准则，首先计算振动信号的数据点数，根据数据点数求出本征模式函数IMF的个数，针对每一个本征模式函数IMF，自适应的选择出筛选次数。本发明专利技术的优点在于能自适应的对每一个本征模式函数IMF采用不同的筛选停止准则，克服了传统的对每一个本征模式函数IMF采用相同固定筛选次数而带来的频率混叠现象，从而分解出物理意义明确的本征模式函数IMF，能够对非线性非平稳信号进行分析，实现机械设备故障的有效诊断。
【技术实现步骤摘要】
本专利技术属于机械设备故障诊断领域，涉及一种总体平均经验模式分解筛选次数的选择方法，可以自适应地选择最优的筛选次数，实现机械设备故障的有效诊断。
技术介绍
随着现代工业和科学技术的飞速发展，工业已成为衡量一个国家科学技术发展的标准。机械设备是工业发展的载体，为工业发展提供关键技术，在国民经济中发挥着越来越重要的作用。同时，机电设备也越来越朝着大型化、复杂化、精密化发展，设备的功能越来越多，性能指标越来越高，其组成与结构越来越复杂，这样势必会使得故障出现的概率大大增加。由于机电设备工况的复杂多样化，机械故障也越来越复杂，设备的故障特征往往是非平稳、非线性的，传统的机械设备故障诊断方法往往只能针对平稳信号才能得到较好的分解结果，对非平稳、非线性信号束手无策。经验模式分解EMD是针对非线性、非平稳信号而提出的一种信号处理方法，它是基于信号局部极值点的一种分解方法用三次样条函数根据信号的局部极值点拟合出上下包络，求取上下包络的均值，再将信号减去所求取的均值，重复上述步骤，直到筛选出来的函数是本征模式函数；再从信号中减去筛选出来的本征模式函数继续筛选，如此循环筛选，直到信号的极值点数目少于3。在经验模式分解中， 存在模式混叠的问题，在其基础上提出了总体平均经验模式分解EEMD。总体平均经验模式分解能够减少模式混叠问题。在总体平均经验模式分解方法中，采用的筛选次数iter是固定不变的10次，也就是对高频与低频分量采用相同的筛选次数iter，在总体平均经验模式分解方法中，高频分量出现频率混叠的概率大于低频分量，因此对于高频分量即频率高的本征模式函数IMF需要采用多的筛选次数iter，对高频分量严格筛选，避免不同模式的函数分在一个本征模式函数IMF中；对于低频分量即频率低的本征模式函数IMF需要采用少的筛选次数iter，对低频分量不严格筛选，以免将一个模式函数分解到两个本征模式函数 IMF 中。
技术实现思路
本专利技术的目的在于克服上述现有技术的缺点，提供一种自适应选择经验模式分解 EEMD筛选次数iter的方法。该方法对不同的本征模式函数采取不同的筛选次数iter，其本质是改进了 EMD方法对应的滤波器特性曲线，调整滤波器的中心频率与通带范围，减少各通带的频率折叠，增大了 EMD所对应的相邻高频部分滤波器中心频率的距离，消除频率混叠。得到物理意义明确的本征模式函数IMF，实现机械设备故障的有效诊断。本专利技术的技术方案是按照如下步骤进行的(1)先计算振动信号的数据点数L，根据数据点数计算信号中本征模式函数IMF的个数N，本征模式函数的个数N= Iog2 L-I ；(2)信号中本征模式函数由高频到低频序号是n(n = 1，2. . . N)，对每一个本征模式函数自适应确定筛选次数iter，最后得到总体平均经验模式分解的筛选次数。所述步骤⑵中本征模式函数自适应筛选次数iter是由公式iter(n) =2(wV) +2确定的。对于序号η的本征模式函数，采用的筛选次数自适应确定为/terw=2(wV)+2，对于n= 1也就是第一个本征模式函数而言，对应的它的筛选次数/ter⑴=2(λμ2)+2=2(λμ)+2， 对于η = 2. . . N的本征模式函数而言，得到对应的筛选次数iter(2). . . iterw。本专利技术的核心是实现了本征模式函IMF数筛选次数iter的自适应选取，使得对于高频分量采用多的筛选次数，严格筛选，避免不同模式的函数分在一个IMF中；对于低频分量(例如低频的谐波)采用少的筛选次数，避免将一个模式分量分解到不同的IMF中，克服了由于采用相同的筛选次数iter而带来的模式混淆问题，能够达到自适应选择筛选次数 iter的功能，实现机械设备故障的有效诊断，分解结果相对比较精确。附图说明图1为自适应总体平均经验模式分解EEMD筛选次数确定方法的流程图；图2(a)为改进筛选次数后经验模式分解EMD所对应的滤波器特性曲线，(b)为固定筛选次数iter = 10经验模式分解EMD所对应的滤波器特性曲线；图3(a)是仿真数据及其各个合成部分，(b)是改进前即筛选次数固定为iter = 10的总体平均经验模式分解EEMD分解结果，(c)是改进后即筛选次数/terw=2(wV)+2的总体平均经验模式分解EEMD分解结果；图4(a)是实际的振动信号，(b)是改进前即筛选次数固定为iter = 10的总体平均经验模式分解EEMD对实际的振动信号的分解结果，(c)是改进后即筛选次数 iter(n) =2(wV) +2的总体平均经验模式分解EEMD对实际的振动信号的分解结果。具体实施例方式下面结合附图对本专利技术做进一步详细描述1)先计算信号的数据长度L，得到本征模式函数的个数N = Iog2 L-I ；2)对于不同的IMF，选择不同的筛选次数/terw=2(wV)+2，N为基于数据的一个量，η为IMF的序号，因此最后得到的就是自适应求得的总体平均经验模式分解(EEMD)的筛选次数。根据上述
技术实现思路
和图1的自适应选择总体平均经验模式分解(EEMD)筛选次数技术的流程图，首先计算数据长度L，然后按照公式N= 10 L-I求出数据中本征模式函数的个数，对于不同的本征模式函数按照由高频到低频的顺序为n(n = 1，2...N)，按照 iter(n) =2(n-"2) +2对序号为η的不同的本征模式函数选择不同的筛选次数。按照上面选择的筛选次数对信号进行EEMD分解，得到的分解结果就是自适应求得的筛选次数的总体平均经验模式分解(EEMD)结果。按照上述
技术实现思路
，分别用基于iter = 10和iter⑷:fn、+2的EMD方法对512 个高斯白噪声点进行分解，对分解后的每个本征模式函数进行傅里叶变换，重复上述过程10000次，得到本征模式函数标准化的频率分布曲线如图2所示，由此可以得出EMD本质是类似于小波变换的二进滤波器，图2中(a)改进筛选次数为/terw=2(w_"2)+2后的滤波器特性曲线，(b)为改进前筛选次数iter = 10的滤波器特性曲线，通过对比，就可以看出筛选次数汝rw=2(wV)+2的EMD方法所对应的滤波器特性曲线中心频率，在高频处比iter = 10的EMD方法所对应的滤波器特性曲线中心频率向高频移动，并且相邻中心频率的带宽变大了，这样就能改善本征模式函数的频率混叠。为了证明上述
技术实现思路
的有效性，先仿真一组数据，仿真信号如图3中(a)所示， 数据长度是1024，采样频率为2000Hz。合成信号s由冲击信号C1、调制信号C2、高频谐波c3、 低频谐波C4和趋势项C5组成。对合成信号s分别用基于iter = 10和Zterw =2(wV) +2的 EEMD方法分解，二者的噪声大小同样选取0. 05，平均次数同样选择100次，两种方法的分解结果分别如图3(b)与(c)所示，图3(b)为筛选次数/terw=2(w_"2)+2的EEMD分解结果，图 3(c)为筛选次数iter= 10的EEMD分解结果。从图中可以看出，筛选次数Zterw =2(wV)+2时，能够将信号的各组成部分Cl、c2、c3, C4与C5较好的分解出来，各个模式函数没有发生频率混叠；iter = 10时，分解结果发生了频率混叠，因此专利技术的此方法能够避本文档来自技高网...
【技术保护点】
１．一种自适应总体平均经验模式分解ＥＥＭＤ筛选次数确定方法，包括：（１）先计算振动信号的数据点数Ｌ，根据数据点数计算振动信号中本征模式函数ＩＭＦ的个数Ｎ，本征模式函数的个数Ｎ＝ｌｏｇ２ Ｌ－１；（２）信号中本征模式函数ＩＭＦ由高频到低频序号是ｎ（ｎ＝１，２．．．Ｎ），对每一个本征模式函数ＩＭＦ自适应确定筛选次数ｉｔｅｒ，最后得到总体平均经验模式分解ＥＥＭＤ的筛选次数。
【技术特征摘要】
【专利技术属性】
技术研发人员：雷亚国，林京，孔德同，王琇峰，廖与禾，
申请(专利权)人：西安交通大学，
类型：发明
国别省市：87
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。