一种局部超声阵列信号的高精度去噪方法技术

本发明专利技术公开了一种局部超声阵列信号的高精度去噪方法，包含以下步骤：运用快速独立分量FastICA将多通道的阵列信号分离为统计的独立分量，对这些独立分量进行自相关分析，初步消除噪声的影响；运用集合经验模态分解算法EEMD对每个通道的信号进行进一步去噪；对精去噪后的信号进行重组，得到高精度的去噪后的信号。本发明专利技术解决了传统的EMD分解只能对单个传感器、单通道信号去噪、FastICA只能对阵列信号进行粗略去噪的弊端，避免了因现场噪声信号的干扰而造成的测向定位误差太大，从而提高了超声阵列信号去燥效果和测向定位精度。

A High Precision Denoising Method for Local Ultrasound Array Signals

The invention discloses a high-precision denoising method for local ultrasonic array signals, which includes the following steps: separating multi-channel array signals into statistical independent components by fast independent component FastICA, autocorrelation analysis of these independent components, and preliminary elimination of the influence of noise; further denoising of each channel signal by using ensemble empirical mode decomposition algorithm EEMD The denoised signal is reconstructed to obtain the denoised signal with high precision. The invention solves the disadvantage that traditional EMD decomposition can only denoise a single sensor, a single channel signal and FastICA can only denoise the array signal roughly, avoids the large error of direction finding and positioning caused by the interference of field noise signal, thereby improving the effect of ultrasonic array signal dedrying and the precision of direction finding and positioning.

【技术实现步骤摘要】
一种局部超声阵列信号的高精度去噪方法
本专利技术涉及一种去噪方法，特别是一种局部超声阵列信号的高精度去噪方法。
技术介绍
电气设备局放超声阵列检测时往往存在噪声的影响，因此有必要对信号的去噪算法进行研究。现有的去噪算法主要分为两类，第一类为单个信号的去噪，如小波去噪、EMD去噪等等，该类算法仅对单通道信号的去噪效果较好，当运用该方法对阵列信号进行去噪时会影响各个阵列信号之间的相位差，对后续阵列信号的处理上造成了重大影响；第二类为多通道阵列信号的去噪，该类方法虽然不会影响阵列信号的相位差，但是该类算法主要是基于盲源分离原理，限制条件比较多，且去噪效果不明显。
技术实现思路
本专利技术所要解决的技术问题是提供一种局部超声阵列信号的高精度去噪方法，其对阵列信号去燥效果好。为解决上述技术问题，本专利技术所采用的技术方案是：一种局部超声阵列信号的高精度去噪方法，其特征在于包含以下步骤：步骤一：运用快速独立分量FastICA将多通道的阵列信号分离为统计的独立分量，对这些独立分量进行自相关分析，初步消除噪声的影响；步骤二：运用集合经验模态分解算法EEMD对每个通道的信号进行进一步去噪；步骤三：对精去噪后的信号进行重组，得到高精度的去噪后的信号。进一步地，所述快速独立分量FastICA为设局放源信号为S(t)＝[S1(t),S2(t),…SN(t)]T，其由N个相互独立的局放信号组成，阵列传感器采集到的阵列信号表示如下式所示：X(t)＝AS(t)+n(t)(1)其中，t为时间变量，X(t)为M维传感器接收的N个阵列信号，即X(t)＝[x1(t),x2(t),…xN(t)]T,A为M×N阶未知混合矩阵，n(t)为M维观测噪声向量；运用观测到的阵列信号X(t)，寻求分离矩阵W，从而得到源信号的估计信号Y(t)。进一步地，所述步骤一具体为1.1首先对阵列信号X(t)进行预处理，预处理主要包含去均值和白化两个步骤，去均值使观测信号满足零均值，白化消除各个数据之间的相关性，使得各个分量尽可能独立；采用下式去除样本的均值：去均值后的阵列信号为X0(t)，对去均值后的数据进行白化处理，即：X′(t)＝TX0(t)(3)白化处理后的阵列信号为X′(t)，T为线性变化矩阵；1.2设置分离矩阵W的初始值；1.3构建二次函数G1，G2，使得式中a1＝1.5,a2＝2，对G1，G2进行求导，得到导函数为g1，g2；1.4运用牛顿迭代法进行循环迭代，迭代式子为W←E(xg1(WTx′))-E(g2(WTx′))W，直到收敛为止；1.5求分离后的阵列信号，即：Y(t)＝WTX(t)。进一步地，所述集合经验模态分解算法EEMD为运用FastICA去噪后得到的阵列信号，各个通道不具有相关性，对其运用单通道去噪算法去噪后不会影响阵列信号之间的相位差；在Yi(t)中加入不同的白噪声序列fi(t)后得到Y’i(t)，对Y’i(t)进行经验模态分解，得到各阶固有模态函数分量(IMF)，此时其中Ck(t)为各阶固有模态函数分量，rn(t)为余量。进一步地，所述步骤二具体为2.1计算出Y’i(t)的所有局部极值点；2.2将所有的极大值点和极小值点构成的上下包络线分别记为u0i(t)和v0i(t)；2.3上下包络线的均值为单通道局放信号的上下包络线的均值的差为H0i(t)＝Yi′(t)-m0i(t)；2.4判断H0i(t)是否满足极值点数目或过零点数目相等或最多相差一个，且由局部极大值构成的上包络线和由局部极小值构成的下包络线的平均值为零，若满足以上条件，则H0i(t)为固有模态分量，反之重复2.1～2.3步骤，直至找到第一个固有模态分量，记为C1(t)；2.5记r1(t)＝Yi′(t)-C1(t)为新的待分析的量，重复以上步骤得到第二个IMF量，……以此类推得到n个IMF量，最终剩下一个单调信号rn(t)，因此单通道原始信号可以表示为：本专利技术与现有技术相比，具有以下优点和效果：本专利技术将EEMD分解与FastICA相结合应用于局放超声阵列信号的去噪，解决了传统的EMD分解只能对单个传感器、单通道信号去噪、FastICA只能对阵列信号进行粗略去噪的弊端，避免了因现场噪声信号的干扰而造成的测向定位误差太大，从而提高了超声阵列信号去燥效果和测向定位精度。附图说明图1是本专利技术的一种局部超声阵列信号的高精度去燥方法的流程图。图2是本专利技术的FastICA分解机重构原理图。图3是本专利技术的EEMD分解流程图。具体实施方式下面通过实施例对本专利技术作进一步的详细说明，以下实施例是对本专利技术的解释而本专利技术并不局限于以下实施例。如图1所示，本专利技术的一种局部超声阵列信号的高精度去噪方法，包含以下步骤：步骤一：运用快速独立分量FastICA将多通道的阵列信号分离为统计的独立分量，对这些独立分量进行自相关分析，初步消除噪声的影响；快速独立分量FastICA为设局放源信号为S(t)＝[S1(t),S2(t),…SN(t)]T，其由N个相互独立的局放信号组成，阵列传感器采集到的阵列信号表示如下式所示：X(t)＝AS(t)+n(t)(1)其中，t为时间变量，X(t)为M维传感器接收的N个阵列信号，即X(t)＝[x1(t),x2(t),…xN(t)]T,A为M×N阶未知混合矩阵，n(t)为M维观测噪声向量；运用观测到的阵列信号X(t)，寻求分离矩阵W，从而得到源信号的估计信号Y(t)。主要分为以下步骤：1.1首先对阵列信号X(t)进行预处理，预处理主要包含去均值和白化两个步骤，去均值使观测信号满足零均值，白化消除各个数据之间的相关性，使得各个分量尽可能独立；采用下式去除样本的均值：去均值后的阵列信号为X0(t)，对去均值后的数据进行白化处理，即：X′(t)＝TX0(t)(3)白化处理后的阵列信号为X′(t)，T为线性变化矩阵；1.2设置分离矩阵W的初始值；1.3构建二次函数G1，G2，使得式中a1＝1.5,a2＝2，对G1，G2进行求导，得到导函数为g1，g2；1.4运用牛顿迭代法进行循环迭代，迭代式子为W←E(xg1(WTx′))-E(g2(WTx′))W，直到收敛为止；1.5求分离后的阵列信号，即：Y(t)＝WTX(t)。步骤二：运用集合经验模态分解算法EEMD对每个通道的信号进行进一步去噪；集合经验模态分解算法EEMD为运用FastICA去噪后得到的阵列信号，各个通道不具有相关性，对其运用单通道去噪算法去噪后不会影响阵列信号之间的相位差；在Yi(t)中加入不同的白噪声序列fi(t)后得到Y’i(t)，对Y’i(t)进行经验模态分解，得到各阶固有模态函数分量(IMF)，此时其中Ck(t)为各阶固有模态函数分量，rn(t)为余量。具体步骤如下：2.1计算出Y’i(t)的所有局部极值点；2.2将所有的极大值点和极小值点构成的上下包络线分别记为u0i(t)和v0i(t)；2.3上下包络线的均值为单通道局放信号的上下包络线的均值的差为H0i(t)＝Yi′(t)-m0i(t)；2.4判断H0i(t)是否满足极值点数目或过零点数目相等或最多相差一个，且由局部极大值构成的上包络线和由局部极小值构成的下包络线的平均值为零，若满足以上条件，则H0i(t)为固有模态分量，反之重复2.1～2.3步骤，直至找到第一个固有模态本文档来自技高网...

【技术保护点】
1.一种局部超声阵列信号的高精度去噪方法，其特征在于包含以下步骤：步骤一：运用快速独立分量FastICA将多通道的阵列信号分离为统计的独立分量，对这些独立分量进行自相关分析，初步消除噪声的影响；步骤二：运用集合经验模态分解算法EEMD对每个通道的信号进行进一步去噪；步骤三：对精去噪后的信号进行重组，得到高精度的去噪后的信号。

【技术特征摘要】
1.一种局部超声阵列信号的高精度去噪方法，其特征在于包含以下步骤：步骤一：运用快速独立分量FastICA将多通道的阵列信号分离为统计的独立分量，对这些独立分量进行自相关分析，初步消除噪声的影响；步骤二：运用集合经验模态分解算法EEMD对每个通道的信号进行进一步去噪；步骤三：对精去噪后的信号进行重组，得到高精度的去噪后的信号。2.按照权利要求1所述的一种局部超声阵列信号的高精度去噪方法，其特征在于：所述快速独立分量FastICA为设局放源信号为S(t)＝[S1(t),S2(t),…SN(t)]T，其由N个相互独立的局放信号组成，阵列传感器采集到的阵列信号表示如下式所示：X(t)＝AS(t)+n(t)(1)其中，t为时间变量，X(t)为M维传感器接收的N个阵列信号，即X(t)＝[x1(t),x2(t),…xN(t)]T,A为M×N阶未知混合矩阵，n(t)为M维观测噪声向量；运用观测到的阵列信号X(t)，寻求分离矩阵W，从而得到源信号的估计信号Y(t)。3.按照权利要求1或2所述的一种局部超声阵列信号的高精度去噪方法，其特征在于：所述步骤一具体为1.1首先对阵列信号X(t)进行预处理，预处理主要包含去均值和白化两个步骤，去均值使观测信号满足零均值，白化消除各个数据之间的相关性，使得各个分量尽可能独立；采用下式去除样本的均值：去均值后的阵列信号为X0(t)，对去均值后的数据进行白化处理，即：X′(t)＝TX0(t)(3)白化处理后的阵列信号为X′(t)，T为线性变化矩阵；1.2设置分离矩阵W的初始值；1.3构建二次函数G1，G2，使得式中a1＝1.5,a2＝2，对G1，G2进...

【专利技术属性】
技术研发人员：马君鹏，王成亮，亓彦珣，吴晗，刘叙笔，杨贤彪，谢庆，岳贤强，
申请(专利权)人：江苏方天电力技术有限公司，华北电力大学保定，国网江苏省电力有限公司，国家电网公司，
类型：发明
国别省市：江苏,32