循环平稳信号EMD去噪分解局限性研究.docVIP

循环平稳信号EMD去噪分解局限性研究 　　摘 要：利用matlab软件进行仿真分析，探讨各类循环平稳信号EMD去噪分解算法的局限性。总结了单独调幅信号、单独调频信号、单独调频调幅信号EMD去噪分解的局限性。 　　关键词：EMD分解;调幅信号;调频信号;调频调幅信号 　　DOI：10.11907/rjdk.143846 　　中图分类号：TP302 　　文献标识码：A 文章编号：1672-7800（2014）012-0033-03 　　作者简介：陆人定（1981-），男，江苏常州人，硕士，常州信息职业技术学院电子与电气工程学院工程师，研究方向为汽车电子技术、嵌入式开发、数字信号处理。 　　0 引言 　　 EMDlt;supgt;[1]lt;/supgt;去噪分解方法应用广泛，在地质勘探lt;supgt;[2]lt;/supgt;、电力检测lt;supgt;[3]lt;/supgt;、音视频处理技术lt;supgt;[4]lt;/supgt;等领域都取得了较好效果。如一些大型发电厂、化工企业的大型齿轮箱有上百个齿轮，发生故障时，技术人员很难直接判断出故障部位，如果逐个部位拆卸检查，不仅影响工作进度，而且会对其它好的齿轮产生影响。在EMD分解与滤波去噪的基础上，通过分析大型齿轮箱振动信号并解调处理lt;supgt;[5]lt;/supgt;，可以准确判断出故障部位。大型故障齿轮箱振动信号表现为故障调幅信号、调频信号、调频调幅信号等lt;supgt;[5]lt;/supgt;。EMD分解算法虽然在一定程度上去除了噪声并分离有用信号，但过度分解会对有用信号产生破坏。 　　1 单个调频、调幅、调频调幅信号EMD分解局限性 　　1.1 单个调幅信号EMD分解局限性 　　 理论上，EMD去噪分解过程中，单个调幅信号可以作为一个IMF分量lt;supgt;[1]lt;/supgt;被分解出，但SD值lt;supgt;[1]lt;/supgt;越小，信号被迭代次数越多，信号的包络越趋于对称化，因此在某一SD下能作为IMF分出的调幅信号在另一更小的SD下将被分解。从式（1）看出，调幅信号可看作3个单频成分的叠加。因此，SD值很小，调幅信号迭代次数过多，原调幅信 　　号被EMD分解完全破坏掉，变成3个单频分量析出。如图1所示，当SD=0.2时，信号被完整析出，当SD=0.002时，信号被破坏分解成两部分（7Hz、8Hz作为一部分，6Hz作为另一部分）。 　　[1+Acos（2πflt;subgt;1lt;/subgt;t）]cos（2πflt;subgt;2lt;/subgt;t）=cos（2πflt;subgt;2lt;/subgt;t）+ 　　A2cos[2π（flt;subgt;1lt;/subgt;+flt;subgt;2lt;/subgt;）t]+A2cos[2π（flt;subgt;1lt;/subgt;-flt;subgt;2lt;/subgt;）t]（1） 　　图1 调幅信号在不同SD值下EMD分解情况 　　1.2 单个调频信号EMD分解局限性 　　单个调频信号与单频信号的包络特性与信号局部均值特性完全一致，因此SD值对信号分解效果影响不大，即使SD值设置很小，迭代次数增加，但其信号仍然能被完整解析出。 　　对比图2与图3，仿真信号x（t）=cos[2π10t+5sin（2π1t）]采样频率60Hz，采样点数1024点。可以看出，取SD=0.2与SD=0.002对信号的分解无太大影响。 　　图2 调频信号在SD=0.2下EMD分解 　　图3 调频信号在SD=0.002下EMD分解 　　1.3 单个调频调幅信号EMD分解局限性 　　调频调幅信号的包络特性类似于调幅信号，又类似于调频信号，如式（2）所示。一个同时调频调幅的信号可看作3个频率成分交叉的纯调频信号叠加。 　　x（t）=[Alt;subgt;1lt;/subgt;+Blt;subgt;1lt;/subgt;cos（2πflt;subgt;1lt;/subgt;t）]cos（2πf11t+Clt;subgt;1lt;/subgt;sin（2πflt;subgt;1lt;/subgt;t）） 　　=Blt;subgt;1lt;/subgt;2cos[2π（f11+flt;subgt;1lt;/subgt;）t+Clt;subgt;1lt;/subgt;sin（2πflt;subgt;1lt;/subgt;t）]+Blt;subgt;1lt;/subgt;2cos[2π（f11 　　-flt;subgt;1lt;/subgt;）t+Clt;subgt;1lt;/subgt;sin（2πflt;subgt;1lt;/subgt;