本发明涉及电缆局部放电在线监测技术领域，更具体地，涉及一种基于emd分解的局部放电信号去噪方法。
背景技术：
电缆局部放电在线监测中，在电缆局部放电在线检测过程中，由于设备处于运行状态，现场存在严重的电磁干扰，这些干扰会和原始局部放电信号进行叠加，造成原始局部放电信号的严重畸变，给依靠局部放电信号波形进行缺陷类型判断和缺陷位置定位的运维人员带来极大困难。因此，对检测到的局部放电信号进行去噪预处理具有重要意义。
在局部放电信号去噪方面，小波分析法由于具有良好的局部时频特性，已广泛应用于局部放电信号的去噪处理，但其容易收到噪声的影响，且基小波难于选择。而经验模式分解(empiricalmodedecomposition，emd)方法基于信号的局部特征时间尺度，能把非平稳信号分解为有限的基本模态分量(intrinsicmodefunction,imf)之和，是一种自适应的信号处理方法，适合非线性和非平稳过程。
技术实现要素：
本发明提供一种基于emd分解的局部放电信号去噪方法，该方法将小波分析方法和emd方法进行有效结合，综合了两者的优势，能够取得更好的去噪效果。
为解决上述技术问题，本发明采用的技术方案是：一种基于emd分解的局部放电信号去噪方法，其中，包括以下步骤：
s1.输入待处理的局部放电信号；
s2.对局部放电信号进行小波分解，得到最高分解尺度的尺度系数和各分解尺度的小波系数；
s3.保留最高分解尺度的尺度系数，对各分解尺度的小波系数进行emd分解，并进行相应的去噪处理，得到处理后的小波系数；
s4.对最高分解尺度的尺度系数和处理后的各分解尺度的小波系数进行小波重构，得到去噪后的局部放电信号。
进一步的，所述步骤s2中，对局部放电信号进行小波分解，得到最高分解尺度的尺度系数和各分解尺度的小波系数的步骤为：
采用db8小波对局部放电信号进行尺度为4的小波分解，得到最高分解尺度的尺度系数和各分解尺度的小波系数。其中a表示尺度系数，dj(j＝1,2,3,4)表示小波系数。
进一步的，所述步骤s3中，保留最高分解尺度的尺度系数，对各分解尺度的小波系数进行emd分解，并进行相应的去噪处理，得到处理后的小波系数的步骤为：
s31.对小波系数dj(t)，经过经验模态分解后，其可表示为：
式中，imfi(t)表示固有模态函数分量，rn(t)表示余项。
s32.对于dj(t)，计算其固有模态函数分量的能量：
式中，ei为固有模态函数分量imfi(t)对应的能量，t为imfi(t)的长度。
s33.对能量ei进行归一化处理：
ei＝ei/emax
式中，emax＝max(ei)，ei为能量值ei归一化后的值。
s34.根据ei的数量级分布，设定阈值eλ，将ei划分为两类，即将大于或等于eλ的ei对应的固有模态函数分量定义为有效分量，将小于eλ的ei对应的固有模态函数分量定义为无效分量。
s35.则进行去噪后的dj(t)可以表示为：
式中，
进一步的，所述步骤s4中，对最高分解尺度的尺度系数和处理后的各分解尺度的小波系数进行小波重构，得到去噪后的局部放电信号的步骤为：
对最高分解尺度的尺度系数a和去噪处理后的各分解尺度的小波系数
本发明与现有技术相比，其有益效果是：
本发明提供的方法将小波分析方法和emd方法进行有效结合，综合了两者的优势，能够取得更好的去噪效果。
附图说明
图1是本发明的流程示意图。
具体实施方式
附图仅用于示例性说明，不能理解为对本专利的限制；为了更好说明本实施例，附图某些部件会有省略、放大或缩小，并不代表实际产品的尺寸；对于本领域技术人员来说，附图中某些公知结构及其说明可能省略是可以理解的。附图中描述位置关系仅用于示例性说明，不能理解为对本专利的限制。
如图1所示，一种基于emd分解的局部放电信号去噪方法，其中，包括以下步骤：
s1.输入待处理的局部放电信号；
s2.对局部放电信号进行小波分解，得到最高分解尺度的尺度系数和各分解尺度的小波系数；具体的，包括如下步骤：
采用db8小波对局部放电信号进行尺度为4的小波分解，得到最高分解尺度的尺度系数和各分解尺度的小波系数。其中a表示尺度系数，dj(j＝1,2,3,4)表示小波系数。
s3.保留最高分解尺度的尺度系数，对各分解尺度的小波系数进行emd分解，并进行相应的去噪处理，得到处理后的小波系数；具体的，包括如下步骤：
s31.对小波系数dj(t)，经过经验模态分解后，其可表示为：
式中，imfi(t)表示固有模态函数分量，rn(t)表示余项。
s32.对于dj(t)，计算其固有模态函数分量的能量：
式中，ei为固有模态函数分量imfi(t)对应的能量，t为imfi(t)的长度。
s33.对能量ei进行归一化处理：
ei＝ei/emax
式中，emax＝max(ei)，ei为能量值ei归一化后的值。
s34.根据ei的数量级分布，设定阈值eλ，将ei划分为两类，即将大于或等于eλ的ei对应的固有模态函数分量定义为有效分量，将小于eλ的ei对应的固有模态函数分量定义为无效分量。
s35.则进行去噪后的dj(t)可以表示为：
式中，
s4.对最高分解尺度的尺度系数和处理后的各分解尺度的小波系数进行小波重构，得到去噪后的局部放电信号；具体的，包括如下步骤：
对最高分解尺度的尺度系数a和去噪处理后的各分解尺度的小波系数
显然，本发明的上述实施例仅仅是为了清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明权利要求的保护范围之内。