• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 10
    • 下载费用:30 金币  

    重庆时时彩漏洞破解: 自适应总体平均经验模式分解EEMD协助噪声大小确定方法.pdf

    摘要
    申请专利号:

    重庆时时彩单双窍门 www.4mum.com.cn CN201410178920.2

    申请日:

    2014.04.29

    公开号:

    CN103954443A

    公开日:

    2014.07.30

    当前法律状态:

    授权

    有效性:

    有权

    法律详情: 专利权人的姓名或者名称、地址的变更IPC(主分类):G01M 13/00变更事项:专利权人变更前:华电电力科学研究院变更后:华电电力科学研究院有限公司变更事项:地址变更前:310030 浙江省杭州市西湖区西湖科技经济园西园一路10号变更后:310030 浙江省杭州市西湖区西湖科技经济园西园一路10号|||授权|||实质审查的生效IPC(主分类):G01M 13/00申请日:20140429|||公开
    IPC分类号: G01M13/00; G06F19/00(2011.01)I 主分类号: G01M13/00
    申请人: 华电电力科学研究院
    发明人: 孔德同; 吴峥峰; 郝建刚; 曾庆猛; 柯愈龙
    地址: 310030 浙江省杭州市西湖区西湖科技经济园西园一路10号
    优先权:
    专利代理机构: 杭州天欣专利事务所 33209 代理人: 张狄峰
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201410178920.2

    授权公告号:

    |||||||||

    法律状态公告日:

    2018.04.24|||2017.01.04|||2014.08.27|||2014.07.30

    法律状态类型:

    专利权人的姓名或者名称、地址的变更|||授权|||实质审查的生效|||公开

    摘要

    本发明涉及一种自适应总体平均经验模式分解EEMD协助噪声大小确定方法。本发明包括如下步骤:(1)初始化噪声加入最大值emax及网格数N,每个网格i(i=1,2…N)对应噪声大小为ei=i/N*emax;(2)对信号x(t)依次加入幅值大小为ei的高斯白噪声ni(t);(3)求取信号xi(t)的极大值点所在位置Max(j)、极小值点所在位置Min(k),得到新序列ai(j)和bi(j);计算新序列ai、bi的标准差Sai、Sbi,得到噪声分布特性Qi与ei的对应曲线关系;(4)Qi随ei迅速下降的点对应的er,即为自适应总体平均经验模式分解EEMD协助噪声大小的最优值。本发明减少了人为选择参数的盲目性。

    权利要求书

    权利要求书
    1.  一种自适应总体平均经验模式分解EEMD协助噪声大小确定方法,其特征在于:包括如下步骤:
    (1)初始化噪声加入最大值emax及网格数N,每个网格i(i=1,2…N)对应噪声大小为ei=i/N*emax;
    (2)对信号x(t)依次加入幅值大小为ei的高斯白噪声ni(t),得到xi(t)=x(t)+ni(t);
    其中ni(t)=ei*mi(t),mi(t)是均值为0,且标准差为1的高斯白噪声序列;
    (3)求取信号xi(t)的极大值点所在位置Max(j)、极小值点所在位置Min(k),得到新序列ai(j)=Max(j+1)-Max(j),bi(j)=Min(k+1)-Min(k);计算新序列ai、bi的标准差Sai、Sbi,得到构造信号xi(t)的噪声分布特性得到Qi与ei的对应曲线关系;
    (4)根据Qi与ei的对应曲线关系,Qi随ei迅速下降的点对应的er,即为自适应总体平均经验模式分解EEMD协助噪声大小的最优值。

    2.  根据权利要求1所述的自适应总体平均经验模式分解EEMD协助噪声大小确定方法,其特征在于:所述步骤(4)中:
    ①对Qi与ei用最小二乘法进行拟合得到Pi,Pi的最小值为Pmin,提取所有Py≥1.2*Pmin的y值构成集合Y;
    ②er根据所有满足Pi≥Pi+L(i∈Y,i+L∈Y)所确定的最小i值imin对应的噪声大小er=imin/N*emax。

    说明书

    说明书自适应总体平均经验模式分解EEMD协助噪声大小确定方法
    技术领域
    本发明涉及一种自适应总体平均经验模式分解EEMD协助噪声大小确定方法,属于机械设备故障诊断领域,可以自适应地选择EEMD方法中协助噪声幅值大小最优值,实现机械设备早期故障的有效诊断。
    背景技术
    随着现代工业和科学技术的飞速发展,工业已成为衡量一个国家科学技术发展的标准?;瞪璞甘枪ひ捣⒄沟脑靥?,为工业发展提供关键技术,在国民经济中发挥着越来越重要的作用。同时,机电设备也越来越朝着大型化、复杂化、精密化发展,设备的功能越来越多,性能指标越来越高,其组成与结构越来越复杂,这样势必会使得故障出现的概率大大增加。
    由于机电设备工况的复杂多样化,机械故障也越来越复杂,设备的故障特征往往是非平稳、非线性的,传统的机械设备故障诊断方法往往只能针对平稳信号才能得到较好的分解结果,对非平稳、非线性信号束手无策。经验模式分解EMD是针对非线性、非平稳信号而提出的一种基于信号局部极值点的处理方法。但是当信号极值点分布不均时,EMD分解结果往往出现模式混淆现象,在其基础上提出了总体平均经验模式分解EEMD,如公开日为2011年11月23日,公开号为CN102254103A的中国专利中,公开了一种自适应总体平均经验模式分解EEMD筛选次数确定方法,以及公开日为2011年11月23日,公开号为CN102254103A的中国专利中,公开了一种自适应总体平均经验模式分解EEMD筛选次数确定方法。通过对信号加入高斯白噪声,改善信号极值点分布,减少模式混淆不足,已被成功应用于转子、齿轮、轴承等旋转机械故障诊断中。
    然而,EEMD分解结果较大程度的依赖于高斯白噪声幅值大小,目前,噪声幅值大小较大程度的依赖于人为选择。为了减少参数人为选择导致的分解结果盲目性,从极值点角度出发,考虑加入噪声后信号的极值点分布情况,在不改变原始信号极值点分布的情况下,最大程度的改善信号极值点,使其分布最为均匀??悸堑降奔尤氲母咚拱自肷到洗笫?,噪声振动剧烈,破坏了信号高频成分的极值点分布,EEMD分解结果会产生多余成分,出现模式混淆;当加入的高斯白噪声幅值较小时,噪声振动微弱,无法改善信号低频成分的极值点分布,EEMD分解结果与EMD类似,无法克服模式混淆不足。因此,加入的噪声大小应在未破坏原始信号高频成分极值点分布特性基础上尽可能改善低频成分极值点分布。
    发明内容
    本发明的目的在于克服现有技术中存在的上述不足,而提供一种减少了参数人为选择带来盲目性的自适应总体平均经验模式分解EEMD协助噪声大小确定方法。
    本发明解决上述问题所采用的技术方案是:该自适应总体平均经验模式分解EEMD协助噪声大小确定方法,其特征在于:包括如下步骤:
    (1)初始化噪声加入最大值emax及网格数N,每个网格i(i=1,2…N)对应噪声大小为ei=i/N*emax;
    (2)对信号x(t)依次加入幅值大小为ei的高斯白噪声ni(t),得到xi(t)=x(t)+ni(t);
    其中ni(t)=ei*mi(t),mi(t)是均值为0,且标准差为1的高斯白噪声序列;
    (3)求取信号xi(t)的极大值点所在位置Max(j)、极小值点所在位置Min(k),得到新序列ai(j)=Max(j+1)-Max(j),bi(j)=Min(k+1)-Min(k);计算新序列ai、bi的标准差Sai、Sbi,得到构造信号xi(t)的噪声分布特性得到Qi与ei的对应曲线关系;
    (4)根据Qi与ei的对应曲线关系,Qi随ei迅速下降的点对应的er,即为自适应总体平均经验模式分解EEMD协助噪声大小的最优值。
    作为优选,本发明所述步骤(4)中:
    ①对Qi与ei用最小二乘法进行拟合得到Pi,Pi的最小值为Pmin,提取所有Py≥1.2*Pmin的y值构成集合Y;
    ②er根据所有满足Pi≥Pi+L(i∈Y,i+L∈Y)所确定的最小i值imin对应的噪声大小er=imin/N*emax。
    本发明与现有技术相比,具有以下优点和效果:核心是实现了EEMD方法协助噪声大小的自适应选取,减少了参数人为选择带来的盲目性。以极值点分布特性为评价参数,既保证了噪声幅值不会过小导致的无法改善信号极值点分布,又避免了噪声幅值不会过大导致的破坏信号高频成分的极值点分布。
    附图说明
    图1为自适应总体平均经验模式分解EEMD协助噪声大小确定方法的流程图。
    图2(a)为仿真信号及其各个组成部分,(b)为仿真信号极值点与噪声大小的对应关系,(c)为仿真信号极值点与噪声大小的对应关系的最小二乘拟合曲线,(d)为自适应EEMD分解结果,(e)为EMD分解结果。
    图3(a)某转子实际振动时域信号,(b)是该转子实际振动信号的频谱,(c)是该实际信号极值点与噪声大小的对应关系,(d)是该实际信号极值点与噪声大小的对应关系的 最小二乘拟合曲线,(e)为实际信号的自适应EEMD分解结果,(f)为实际信号的EMD分解结果。
    具体实施方式
    下面结合附图并通过实施例对本发明作进一步的详细说明,以下实施例是对本发明的解释而本发明并不局限于以下实施例。
    实施例。
    参见图1至图3,本实施例中的自适应总体平均经验模式分解EEMD协助噪声大小确定方法包括如下步骤:
    (1)初始化噪声加入最大值emax及网格数N,每个网格i(i=1,2…N)对应噪声大小为ei=i/N*emax。
    (2)对信号x(t)依次加入幅值大小为ei的高斯白噪声ni(t),得到xi(t)=x(t)+ni(t);
    其中ni(t)=ei*mi(t),mi(t)为均值为0,标准差为1的高斯白噪声序列。
    (3)求取信号xi(t)的极大值点所在位置Max(j)、极小值点所在位置Min(k),得到新序列ai(j)=Max(j+1)-Max(j),bi(j)=Min(k+1)-Min(k);计算新序列ai、bi的标准差Sai、Sbi,得到构造信号xi(t)的噪声分布特性得到Qi与ei的对应曲线关系。
    (4)根据Qi与ei的对应曲线关系,对Qi与ei用最小二乘法进行拟合得到Pi,Pi的最小值为Pmin,提取所有Py≥1.2*Pmin的y值构成集合Y。
    (5)er根据所有满足Pi≥Pi+L(i∈Y,i+L∈Y)所确定的最小i值imin对应的噪声大小er=imin/N*emax。
    根据上述发明内容和图1的自适应选择总体平均经验模式分解(EEMD)噪声大小技术的流程图,首先初始化噪声最大值及网格数,得到各网格对应的噪声大小,然后依次将各网格对应幅值大小的高斯白噪声加入信号中,评估噪声协助信号的极值点分布特性,得到极值点分布与噪声大小的对应关系,最后,对Qi与ei用最小二乘法进行拟合得到Pi,Pi的最小值为Pmin,提取所有Py≥1.2*Pmin的y值构成集合Y;er根据所有满足Pi≥Pi+L(i∈Y,i+L∈Y)所确定的最小i值imin对应的噪声大小er=imin/N*emax。
    按照上面选择的最优噪声大小对信号进行EEMD分解,得到的分解结果就是自适应求得噪声大小的总体平均经验模式分解(EEMD)结果。
    为了证明上述发明内容的有效性,先仿真一组信号,仿真信号s如图2中(a)所示,由冲击信号c1、调制信号c2、谐波c3组成。对合成信号s利用本发明提出的自适应EEMD方法分 解。初始化噪声最大值e=0.5,网格数1000,得到信号的极值点分布与噪声大小的曲线如图2(b)所示,最小二乘拟合得到光滑曲线如图2(c)。从图2(c)中可以看出,当加入的噪声较小时(ei<0.008),噪声振动微弱,极值点分布的离散性随噪声的增加而增大,这是由于加入的噪声振动微弱,仅仅改善了信号中某些低频成分的极值点分布,导致极值点分布呈现较大的离散型;当噪声大小进一步增大,达到ei=0.008时,极值点分布特性Fi达到最大值,信号的极值点离散型最大。随后,当加入的噪声较大时(ei>0.008),极值点分布的不均匀性随噪声的增加而减小,最后趋于稳定,主要是由于随着噪声较大,振动剧烈,剧烈振动的噪声破坏了高频成分的振动特性,噪声协助信号的极值点全部是由加入噪声引起。因此,根据图2(c)可以确定,最优噪声大小为0.02,将确定大小的高斯白噪声加入仿真信号s进行EEMD分解,分解得到的IMF如图2(d)所示??梢钥闯?,分解得到的前三个IMF分别是冲击、调制、谐波分量,依次对应原始信号中的各组成成分,每个IMF的幅值与构成原始信号的各成分几乎相等,各IMF之间不存在模式混淆现象。对仿真信号进行EMD分解,分解结果如图2(e)所示。分解得到的第一个IMF中,既包括了冲击成分,也含有部分调制成分,很明显分解结果发生了模式混淆现象。此外,第三个IMF相比原始信号的谐波成分,存在严重失真。
    同时将上述方法应用在实际数据的分析中。图3(a)为用位移传感器测得的某炼油厂重催机组振动信号,图3(b)为其对应的频谱,频谱中的三个主要频率为f1=25.35Hz,f2=97.65Hz与f3=193.36Hz,分别对应齿轮箱低速轴的转频、机组转频和二次谐波。因此从信号的时域波形及频谱中均无法判别振动过大的原因。利用自适应EEMD对信号进行分析,初始化噪声最大值e=0.5,网格数1000,得到极值点分布特性与噪声大小的对应关系如图3(c),最小二乘拟合得到光滑曲线如图3(d),从图3(d)中可以看出,随着加入噪声的不断增加,极值点的分布不均匀性逐渐减少。在噪声大小为0.02时,信号极值点分布不均匀特性趋于平稳,之后,随着噪声的不断增大,信号极值点分布不均匀性急剧下降,最后趋于平稳,这是由于信号的极值点特性已被加入的噪声破坏。因此,最佳噪声大小为0.02。将确定大小的高斯白噪声加入实际信号进行EEMD分解,分解得到的前三个IMF及趋势项如图3(e)所示。从图中可以看出,IMF1、IMF2、IMF3分别为分解得到的前三个分量,R为趋势项。IMF2对应机组的工频振动,频率为f2=97.65Hz,IMF3代表齿轮箱低速轴的转频,频率为f1=25.35Hz。仔细观察IMF1,其中含有周期为0.031s的微弱冲击成分,接近机组旋转周期的3倍(1/97.65Hz=0.01024s),表明冲击信号的频率是1/3倍的机组工频,表明轴与轴瓦之间产生了早期碰摩故障。为了进行对比,利用EMD方法对信号进行分解,分解得到的前三个IMF及趋势项如图3(f)所示,可以看出,IMF1中周期性冲击不明显,且IMF2与IMF3已出现了模式混淆,无法判别故障原因。
    通过以上仿真信号和实际信号的分解结果对比,可以得到本发明的自适应选取协助噪声大小的EEMD方法能够在一定程度上减少模式混叠,能够分解出物理意义明确的本征模式函数,分解结果比较好、操作方便简单,说明本发明的自适应总体平均经验模式分解EEMD协助噪声大小确定方法可以更好的实现机械设备的故障诊断。
    虽然本发明已以实施例公开如上,但其并非用以限定本发明的?;し段?,任何熟悉该项技术的技术人员,在不脱离本发明的构思和范围内所作的更动与润饰,均应属于本发明的?;し段?。

    关 键 词:
    自适应 总体 平均 经验 模式 分解 EEMD 协助 噪声 大小 确定 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    关于本文
    本文标题:自适应总体平均经验模式分解EEMD协助噪声大小确定方法.pdf
    链接地址://www.4mum.com.cn/p-6142890.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03