基于CEEMDAN和HT的谐波检测新方法.pdf
《基于CEEMDAN和HT的谐波检测新方法.pdf》由会员分享,可在线阅读,更多相关《基于CEEMDAN和HT的谐波检测新方法.pdf(6页珍藏版)》请在咨信网上搜索。
1、第6 0 卷第6 期2023年6 月15日电测与仪 表Electrical Measurement&InstrumentationVol.60 No.6Jun.15,2023基于CEEMDAN和HT的谐波检测新方法张乐乐,王海云,王维庆(新疆大学可再生能源发电与并网技术教育部工程研究中心,乌鲁木齐8 3 0 0 47)摘要:经验模态分解(EMD)作为希尔伯特-黄变换(HHT)的重要组成部分,为了克服其在谐波检测中出现的模态混叠、端点效应问题,提出采用自适应噪声完备集合经验模态分解(CEEMDAN)和希尔伯特变换(HT)相结合的谐波检测新方法。文章首先在理论上对比分析了EMD、EEM D 以及C
2、EEMDAN算法,研究CEEMDAN算法的特性。再用CEEMDAN算法对原始信号进行分解,得到固有模态函数(IMF)。最后用HT算法对每阶IMF分量进行分析,检测到谐波中包含的瞬时幅频信息。算例仿真结果表明,相对于HHT算法对信号的处理能力,文中提出的方法在谐波检测中有效地克服了EMD算法的弊端,提高了信号分解精度。关键词:经验模态分解;希尔伯特-黄变换;自适应噪声完备集合经验模态分解;希尔伯特变换;谐波D0I:10.19753/j.issn1001-1390.2023.06.021中图分类号:TM935A novel harmonic detection method based on CE
3、EMDAN and HT(Engineering Research Center of Ministry of Education for Renewable Energy Power Generation and Grid Technology,Abstract:Empirical mode decomposition(EMD)is an important part of Hilbert-Huang transform(HHT),in order to o-vercome its modal aliasing and endpoint effect in harmonic detectio
4、n,a novel harmonic detection method combing adaptivenoise complete integration of empirical mode decomposition(CEEMDAN)and Hilbert transform(HT)is proposed in thispaper.This method firstly compares and analyzes EMD,EEMD and CEEMDAN algorithm in theory,and studies the char-acteristics of CEEMDAN algo
5、rithm,Then,CEEMDAN algorithm is adopted to decompose the original signals to obtain theintrinsic modal function(IMF).Finally,the HT algorithm is used to analyze the IMF component of each order,and theinstantaneous amplitude-frequency information contained in the harmonic is detected.The simulation r
6、esults show thatcompared with the signal processing capability of HHT algorithm,the proposed method overcomes the disadvantages ofEMD algorithm in harmonic detection and improves the accuracy of signal decomposition.Keywords:empirical mode decomposition,Hilbert-Huang transform,complete ensemble empi
7、rical mode decompositionwith adaptive noise,Hilbert transform,harmonic0引言为了更好的利用我国的风、光、水资源,以及各种电子设备的完善和跨距离远路程直流输电技术的发展,致使电网中新能源的占比日益趋增。另一方面,随着用户端精密设备的投人使用,电力用户对电能质量基金项目:国家自然科学基金资助项目(516 6 7 0 2 0);新疆维吾尔自治区教育厅重点项目(XJEDU2019I009);新疆维吾尔自治区实验室开放课题(2 0 18 D04005);教育部创新团队滚动项目(IRT-16R633)文献标识码:AZhang Lele,
8、Wang Haiyun,Wang WeiqingXinjiang University,Urumqi 830047,China)文章编号:10 0 1-13 90(2 0 2 3)0 6-0 147-0 6的要求也随之增加。电能质量 1能否恒定达标成为永恒的话题,而谐波检测是其中重要的组成部分。经验模态分解(EmpiricalModeDecomposition,EMD)是将任意目标信号分解成一系列有限个具有不同特征尺度的固有模态函数(Intrinsic Mode Functions,IMFs)和残余分量,分解后的各IMF分量按照瞬时频率由高频到低频的顺序进行排列,这个强大的自适应算法通常被用于
9、对非线性、非平稳信号进行线性化、平稳化的处理。但EMD算法本身也存在一些问题:模态混叠一147 一第6 0 卷第6 期2023年6 月15日现象。模态混叠是同一IMF中出现不同幅值或频率的信号,或者同一幅值或频率的信号出现在多个不同的IMF中。这种混叠现象会增大检测误差,影响检测结果。EMD相对于小波变换2 4 系列的信号处理方法而言,不需要考虑小波分解时需要合理选择小波基函数和分解尺度的问题。文献5 以异常信号为例,分析了模态混叠现象产生的原因,并通过加人高频谐波的方法解决模态混叠现象,该方法在一定程度上削弱了模态混叠现象,但加入的多余高频谐波在分解过程中无法根除,且EMD对间歇信号的处理并
10、不理想。关于EMD分解过程中出现的模态混叠问题,文献6 采用总体平均经验模态分解(Ensemble Empirical Mode Decomposition,EEMD)。该方法在一定程度上抑制了模态混叠现象,但集成次数的多少严重影响重构误差的大小,导致计算量大,且因无法消除加人的辅助噪声,加大了重建后的信号与原信号的误差。为了克服EEMD算法不能精确重构的问题、减少计算量,文献7 采用了自适应噪声完备集合经验模态分解(CompleteEnsembleEmpiricalMode Decomposition with Adaptive Noise,CEEMDAN),该方法实现了近似于完美的信号重构
11、的同时,改善了模态混叠现象,也是对EEMD算法的重要改进。尽管希尔伯特-黄变换8 能够克服快速傅里叶变换9(FastFourierTransform,FFT)无法将信号频域和时域相结合的缺点,也不需要考虑根据经验选择的基函数会造成误差,更不需要通过大量的数据收集、训练样本得到误差标准,但HHT算法本身的组成部分EMD算法10 1存在一些不足。因此,文章提出采用CEEM-DAN和希尔伯特变换(Hilbert Transform,HT)相结合的谐波检测新方法。该方法通过处理目标信号,得到时频谱和幅频谱,通过分析对比,仿真结果表明CEEM-DAN-HT谐波分析法的优越性。1算法原理步骤与分析1.1C
12、EEMDAN算法在CEEMDAN算法I中,定义s(t)是原始信号,8。为噪声标准差,U;(t)是第i次加到原信号中的零均值单位方差的白噪声,E()是由CEEMDAN产生的第k阶模态分量算子,IMF(t)是信号经过分解后得到的第j阶模态分量,t为时间变量,则CEEMDAN算法的实现过程如下:(1)将噪声8.U;(t)加到原始信号s(t)中,每添加一次噪声都通过CEEMDAN算法对信号进行一次加噪处理后,得到加人I次噪声后分解出的CEEMDAN一一148 一电测与仪表Electrical Measurement&Instrumentation阶模态分量:Mr()=1MP:(C)一(2)计算步骤(1
13、)中的第一个余量信号ri(t):ri(t)=s(t)-IMFl(t)(3)按照步骤(1)向余量信号ri(t)中添加噪声分量8,E,u(t),然后对添加噪声的余量信号进行分解,直到零点个数和极值点个数相差不超过1时,分解结束,得到第一个IMF分量,则得到CEEMDAN二阶模态分量为:(3)(4)重复步骤(2)和步骤(3),计算第j个余量信号和第i+1模态分量:r;(t)=Tj-I(t)-IMF)(t)MF)m(C)=+2 E(c(0)+,(0)(5)执行步骤(4)直至CEEMDAN分解的余量信号无法进行EMD分解,即残余分量满足以下两个条件:(a)包络线上的零点个数和极值点个数相差不超过1;(b
14、)每一点处局部极大值和局部极小值的均值为零。满足以上两个条件后算法终止,得到分解的模态分量个数为J,原始信号可表示为:Js(t)=IMF(t)+R(t)j=1余量信号R(t)可表示为:R(t)=s(t)-IMF)(t)1.2对比分析CEEMDAN设定一组目标信号如式(8)所示,分别用EMD、EEMD以及CEEMDAN算法对这组目标信号进行分解,通过计算对比以上三种算法的完备性、正交性、处理时间以及分量重构个数四个数量指标,以此验证CEEMDAN算法处理谐波的特性。设目标信号为:x(t)=220sin(2 50 t)+100sin(2 150 t)h,(t)+50sin(2 250 t)hz(t
15、)(8)式中hi(t)、h(t)是两个门信号,分别是h(t-0.2)-h(t-0.4)、h(t-0.6)-h(t-0.8),构成了含有零信号的间歇信号,仿真结果如图1所示。图1中的目标信号是由三组信号合成:50 Hz、有Vol.60 No.6Jun.15,2023(1)(2)(4)(5)(6)(7)第6 0 卷第6 期2023年6 月15日效幅值2 2 0 V的信号,0.2 0.4时间段有150 Hz、有效幅值10 0 V的信号,0.6 0.8 时间段含有2 50 Hz、有效幅值50 V的信号。500/0AM-5000.10.20.30.40.50.60.70.80.901150Hz有效幅值1
16、0 01000-10000.10.20.30.40.50.60.70.80.91250Hz有效幅值50500-5000.10.20.30.40.50.60.70.80.91目标信号5000AWW-50000.10.2 0.30.40.50.60.70.80.91时间/s图1目标信号Fig.1Target signal用EMD、EEM D 以及CEEMDAN算法分别对目标信号进行处理,根据三个算法理论,目标信号被分解成一系列有限个不同幅值的IMF分量,各分量按照瞬时频率由高频到低频的顺序进行排列,仿真结果如图2 所示。由图2 仿真结果可以看出,EMD的各IMF分量的模态混叠现象比较严重,无法准确
17、的判别构成原信号的各分量信号信息。因此,处理这种含有间歇信号的复杂信号,EMD算法远不及EEMD和CEEMDAN算法。而EEMD和CEEMDAN算法对目标信号的分解结果比较相近,能够清晰的看出各IMF分量中含有的原信号信息,以及含有门信号的时间段。同时,通过MATLAB分别用EMD、EEM D 以及CEEMDAN算法多次处理图1中的目标信号,计算得到以上三个算法的完备性、正交性、处理时间以及分量重构个数四个数量指标值的平均值,计算结果如表1所示。正交性体现各IMF分量间的相关性,其值越小说明各IMF分量间的模态混叠程度越低,求解公式如下:nC,(t)C;(t)dtj=1K+iW(t)dt式中C
18、,(t)、C(t)是两个不同的IMF分量;s(t)为原始信号。定义原始信号的能量为:E,=s(t)dt电测与仪 表Electrical Measurement&Instrumentation500LNI050000.10.2 0.30.40.50.60.70.80.91500050050Hz有效幅值2 2 000.10.20.30.40.5 0.60.7 0.80.91200020000.10.20.30.40.50.60.70.8 0.91500-5000.10.2 0.30.4 0.5 0.6 0.70.80.91500-5000.10.2 0.30.40.50.60.7 0.80.91时
19、间/s(a)对目标信号进行EMD分解202000.10.2 0.30.40.5 0.60.70.80.91000-100500.10.20.30.40.50.60.70.8 0.91500EHWI0ww50000.10.20.30.40.50.60.70.80.91200020000.10.20.30.40.50.60.70.80.9120-200.10.20.30.40.50.60.70.80.91时间/s(b)对目标信号进行EEMD分解20LHWI0-2000.10.20.30.40.50.60.70.81000-10000.10.2 0.30.40.5 0.60.70.80.91500E
20、HI0-50000.10.20.30.40.50.60.70.80.912000-20000.10.20.30.40.50.60.70.80.9150-500.10.20.30.40.50.60.70.80.91时间/s(c)对目标信号进行CEEMDAN分解图2对目标信号进行EMD、EEM D、CEEM D A N分解Fig.2 EMD,EEMD,CEEMDAN decompositionof the target signal各IMF分量的能量为:(9)E,=Jc(t)dt各IMF分量的能量总和为:E,=E,(10)E,与E,的差体现了完备性,即原始信号的能量与一149 一Vol.60 No
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 CEEMDAN HT 谐波 检测 新方法
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。