接地网干扰信号的去噪方法、装置及计算机设备与流程

1.本发明涉及接地网干扰信号降噪技术领域，具体涉及一种接地网干扰信号的去噪方法、一种接地网干扰信号的去噪装置和一种计算机设备。


背景技术：

2.接地网干扰场易受周围环境干扰，产生随机误差，从而影响保护设备稳定性，严重时可能会导致闸刀误动。由于接地网电位差的干扰，智能变电站的数据采集设备可能会得到异常数据，该数据不是一次电气信号的正确反映，但一些情况下其数据帧中的品质因数位是正常的，测量和保护装置视其为正常数据进行处理，会对结果造成较大影响，严重时会造成保护设备的误动。接地网干扰信号主要由周围高压环境引起，由于这类干扰没有规律可循，难以用建模的方式进行消除或弱化。对高频干扰采用滤波、屏蔽等手段加以抑制，但对于频率与接地网干扰场接近的干扰成分，传统的滤波方法难以起到明显的去噪效果，导致随机噪声对接地网干扰信号的影响较大，进一步导致接地网干扰匹配精度较低。


技术实现要素：

3.本发明为解决上述技术问题，提供了一种接地网干扰信号的去噪方法，可以有效降低随机噪声对接地网干扰信号的影响，提高接地网干扰匹配精度。
4.本发明采用的技术方案如下：
5.本发明的第一方面实施例提出了一种接地网干扰信号的去噪方法，包括以下步骤：
6.s1，采集接地网干扰信号b
s
(t)；
7.s2，利用形态学滤波器对所述接地网干扰信号b
s
(t)进行交替混合形态滤波处理，以得到滤波信号b
s
(t)
′
；
8.s3，对所述滤波信号b
s
(t)
′
进行经验模态分解处理，以得到多个本征模态函数imf分量和一个余项；
9.s4，确定所述接地网干扰信号b
s
(t)分别与各个imf分量间的互相关系数，并确定所述接地网干扰信号b
s
(t)分别与各个imf分量间的欧氏距离；
10.s5，选取各个所述互相关系数中的极小值，作为第一分界点，并选取各个所述欧氏距离中的极大值点的下一个点，作为第二分界点，根据所述第一分界点及所述第二分界点，确定所述多个imf分量中的噪声imf分量、混合及信号imf分量，其中，所述第二分界点与所述第一分界点不同，所述噪声imf分量中仅包含噪声，所述混合imf分量中包含噪声和信号，所述信号imf分量中仅包含信号；
11.s6，将每个所述混合imf分量进行归一化最小均方差自适应降噪处理；
12.s7，将处理后的各个混合imf分量和所述信号imf分量进行重构，以得到去除噪声后的接地网干扰信号。
13.另外，根据本发明上述实施例提出的接地网干扰信号的去噪方法，还可以具有如
下附加的技术特征：
14.根据本发明的一个实施例，步骤s2包括以下步骤：
15.s21，利用形态学滤波器对接地网干扰信号进行腐蚀运算及膨胀运算，运算公式分别为；
16.(b
s
θg)(n)＝min{b
s
(n m)
‑
g(m):m∈d
g
}，
[0017][0018]
其中，θ为腐蚀运算，为膨胀运算，b
s
(n)为接地网干扰信号，b
s
(n)的定义域为d
f
＝{0,1,2,3,
……
，n}，g(n)为一维结构元素序列，g(n)的定义域为d
g
＝{0,1,2,
……
，p}，其中p和n均是整数，且n≥p，n为接地网干扰信号的长度，g(m)为形态学滤波器的结构元素，m为结构元素长度；
[0019]
s22，对所述接地网干扰信号进行形态学开运算，并对接地网干扰信号b
s
(n)进行形态学闭运算，运算公式分别为：
[0020][0021][0022]
其中，ο为开运算，
·
为闭运算；
[0023]
s23，利用交替混合滤波器对所述接地网干扰信号进行交替混合滤波运算，运算公式为：
[0024]
[(b
s
)altmix(g)](n)＝[(b
s
)oc(g) (b
s
)co(g)](n)/2，
[0025]
其中，altmix为交替混合滤波运算，oc为先开后闭的交替滤波运算，co为先闭后开的交替滤波运算，
[0026][0027][0028]
根据本发明的一个实施例，步骤s3包括以下步骤：
[0029]
s31，令k＝1，i＝1，选取所述滤波信号b
s
(t)
′
中的所有极大值点和所有极小值点，并将所述所有极大值点连接起来形成上包络线、将所述所有极小值点连接起来形成下包络线，确定所述上包络线与所述下包络线间的均值信号m(t)；
[0030]
s32，将所述滤波信号b
s
(t)
′
减去所述均值信号m(t)，得到中间信号h
i
(t)，判断所述中间信号h
i
(t)是否满足emd分解的条件：
[0031][0032]
且所述中间信号h
i
(t)的极值点个数和过零点个数间的差值是否满足不大于1，若满足，则记第i个本征模态函数imf
i
(t)为：
[0033]
imf
i
(t)＝h
i
(t)，
[0034]
若不满足，则令k＝k 1，i＝1，返回步骤s31中；
[0035]
s33，将所述滤波信号b
s
(t)
′
减去所述第i个本征模态函数imf
i
(t)，得到一个余项信号r
i
(t)，然后判断r
i
(t)是否满足余项条件，若满足，令emd分解的最后一个余项为r
l
(t)；
若不满足，则令i＝i 1，i＝1，b
s
(t)'＝r
i
(t)返回步骤s32继续分解；最终分解后的信号为：
[0036][0037]
其中，h
i
(t)代表第i个imf分量，l表示所有imf分量的个数，r
l
(t)为余项。
[0038]
根据本发明的一个实施例，步骤s4包括以下步骤：
[0039]
s41，计算所述接地网干扰信号b
s
(t)分别与各个imf分量间的互相关系数，分别记为ρ1,ρ2,
……
ρ
l
；
[0040]
s42，计算所述接地网干扰信号b
s
(t)及各个imf分量的概率密度函数，分别记为pdf(b
s
(t))和pdf(h
i
(t))，计算pdf(b
s
(t))分别和各个pdf(h
i
(t))之间的欧氏距离，分别记为s(1),s(2),
……
，s(l)。
[0041]
根据本发明的一个实施例，步骤s5包括以下步骤：
[0042]
s51，确定各个所述互相关系数中的极小值ρ
min
＝[ρ1,ρ2,
……
ρ
l
]，及确定各个所述欧氏距离中的极大值s
max
＝[s(1),s(2),
……
，s(l)]；
[0043]
s52，令第一分界点m1＝ρ
min
、第二分界点m2＝s
max
1，确定互相关系数小于m1时对应的imf分量为所述噪声imf分量，确定互相关系数大于或者等于m1、且欧氏距离大于或者等于m2时对应的imf分量为所述混合imf分量，确定欧氏距离小于m2时对应的imf分量为所述信号imf分量。
[0044]
根据本发明的一个实施例，步骤s6包括以下步骤：
[0045]
通过变步长nlms滤波算法对每个所述混合imf分量进行自适应降噪处理，处理过程为：
[0046]
设任意一个混合imf分量为信号x
imfi
(n)，对应的滤波输出信号为：
[0047]
y
imfi
(n)＝w
t
(n)x
imfi
(n)，
[0048]
误差信号为：
[0049]
e(n)＝d(n)
‑
y
imfi
(n)，
[0050]
权系数更新信号为：
[0051][0052]
其中，α(0＜α＜1)为收敛因子，β(0＜β≤1)是正常数，||x(n)||2为抽头输入向量x
imf
(n)的欧氏范数平方。
[0053]
根据本发明的一个实施例，去除噪声后的接地网干扰信号为：
[0054][0055]
其中，m1为第一分界点，m2为第二分界点，y
imfi
(t)为将混合imf分量进行归一化最小均方差自适应降噪处理后得到的第i个信号，h
i
(t)代表第i个imf分量，为信号imf分量之和，l表示所有imf分量的个数，r
l
(t)为余项。
[0056]
本发明第二方面实施例提出了一种接地网干扰信号的去噪装置，包括：
[0057]
采集模块，用于采集接地网干扰信号b
s
(t)；
[0058]
滤波模块，用于利用形态学滤波器对所述接地网干扰信号b
s
(t)进行交替混合形态滤波处理，以得到滤波信号b
s
(t)
′
；
[0059]
分解模块，用于对所述滤波信号b
s
(t)
′
进行经验模态分解处理，以得到多个本征模态函数imf分量和一个余项；
[0060]
确定模块，用于确定所述接地网干扰信号b
s
(t)分别与各个imf分量间的互相关系数，并确定所述接地网干扰信号b
s
(t)分别与各个imf分量间的欧氏距离；
[0061]
选取模块，用于选取各个所述互相关系数中的极小值，作为第一分界点，并选取各个所述欧氏距离中的极大值点的下一个点，作为第二分界点，根据所述第一分界点及所述第二分界点，确定所述多个imf分量中的噪声imf分量、混合及信号imf分量，其中，所述第二分界点与所述第一分界点不同，所述噪声imf分量中仅包含噪声，所述混合imf分量中包含噪声和信号，所述信号imf分量中仅包含信号；
[0062]
降噪模块，用于将每个所述混合imf分量进行归一化最小均方差自适应降噪处理；
[0063]
重构模块，用于将处理后的各个混合imf分量和所述信号imf分量进行重构，以得到去除噪声后的接地网干扰信号。
[0064]
本发明第三方面实施例提出了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时，实现本发明上述实施例提出的接地网干扰信号的去噪方法。
[0065]
本发明第三方面实施例提出了一种非临时性计算机可读存储介质，其上存储有计算机程序，其该计算机程序被处理器执行时实现本发明上述实施例提出的接地网干扰信号的去噪方法。
[0066]
本发明实施例的技术方案，首先通过形态学滤波器对接地网干扰信号进行形态学滤波处理，以消除其中的脉冲干扰，之后，对滤波处理后的干扰信号进行emd分解得到多个imf分量和一个余项，利用互相关系数和欧氏距离准则将emd分解的imf分量分解为噪声imf分量、混合imf分量和信号imf分量，将混合imf分量进行nlms自适应滤波处理，之后，将处理后的混合imf分量和信号imf分量进行重构，即可得到去噪后的接地网干扰信号。由此，本发明实施例的接地网干扰信号的去噪方法，可以有效地去除瞬时脉冲干扰，降低随机噪声对接地网干扰信号的影响，提高接地网干扰匹配精度；且可以通过互相关系数和欧氏距离准则有效确定emd分解中模态分界点，从而可以提高去噪可靠性。
附图说明
[0067]
图1为本发明实施例的接地网干扰信号的去噪方法的流程图。
[0068]
图2为本发明一个实施例的接地网干扰信号的去噪示意图。
[0069]
图3为本发明实施例的接地网干扰信号的去噪装置的方框示意图。
具体实施方式
[0070]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于
本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0071]
相关技术中，通过小波变换进行干扰信号的变换分析方法，然而，在相同的条件下，不同的小波基和分解层数将会产生不同的降噪效果，不利于小波去噪的应用。而经验模态分解(empirical mode decomposition，简称emd)是目前处理非平稳信号的一个新的重要发展方向。基于此，本发明实施例在emd分解的基础上提出一种接地网干扰信号的去噪方法。
[0072]
图1为本发明实施例的接地网干扰信号的去噪方法的流程图。
[0073]
如图1所示，该接地网干扰信号的去噪方法包括以下步骤：
[0074]
s1，采集接地网干扰信号b
s
(t)。
[0075]
在实际应用中，通过接地网干扰传感器测量接地网干扰信号，进而通过采集接地网干扰传感器输出信号或者模拟输出信号即可采集到的接地网干扰信号b
s
(t)。
[0076]
s2，利用形态学滤波器对接地网干扰信号b
s
(t)进行交替混合形态滤波处理，以得到滤波信号b
s
(t)
′
。
[0077]
其中，形态学滤波器可包括交替滤波器和交替混合滤波器。
[0078]
具体地，在采集到的接地网干扰信号b
s
(t)之后，利用形态学滤波器对b
s
(t)进行交替混合形态滤波处理，以消除脉冲干扰，得到滤波信号b
s
(t)
′
。
[0079]
s3，对滤波信号b
s
(t)
′
进行经验模态分解处理，以得到多个本征模态函数imf分量和一个余项。
[0080]
具体地，在得到滤波信号b
s
(t)
′
之后，可根据emd算法对b
s
(t)
′
进行emd分解，分解后得到多个imf分量和一个余项。
[0081]
s4，确定接地网干扰信号b
s
(t)分别与各个imf分量间的互相关系数，并确定接地网干扰信号b
s
(t)分别与各个imf分量间的欧氏距离。
[0082]
需要说明的是，两个信号之间的互相关系数可以用来表征该两个信号间的关联程度，互相关系数的绝对值越接近于1，表示两个信号的相关性越强；当互相关系数为1，表示两个信号完全相关；互相关系数为0，表示两个信号完全不相关，其中，当互相关系数为正数时，表示两个信号为正相关；当互相关系数为负数时，表示两个信号为负相关。两个信号之间的欧氏距离可以用来表征该两个信号间的相似度，欧氏距离越大，相似度越小；欧氏距离越小，相似度越大。
[0083]
具体地，在得到多个imf分量之后，可根据互相关系数算法确定每个imf分量与接地网干扰信号b
s
(t)间的互相关系数，以得到多个互相关系数，该互相关系数即可表征每个分量与接地网干扰信号间的相关程度。并确定每个imf分量与接地网干扰信号b
s
(t)间的欧氏距离，该欧氏距离即可表征每个分量与接地网干扰信号间的相似度。
[0084]
s5，选取各个互相关系数中的极小值，作为第一分界点，并选取各个欧氏距离中的极大值点的下一个点，作为第二分界点，根据第一分界点及第二分界点，确定多个imf分量中的噪声imf分量、混合及信号imf分量，其中，第二分界点与第一分界点不同，噪声imf分量中仅包含噪声，混合imf分量中包含噪声和信号，信号imf分量中仅包含信号。
[0085]
其中，噪声imf分量中仅包含噪声，混合imf分量中包含噪声和信号，信号imf分量中仅包含信号。即：噪声imf分量是指多个imf分量中仅包含噪声的分量，混合imf分量是指
多个imf分量中包含噪声和信号(即没有噪声的信号)的分量，信号imf分量是指多个imf分量中仅包含信号(即没有噪声的信号)的分量。
[0086]
具体地，在得到多个互相关系数及多个欧氏距离之后，可选取出多个互相关系数中的极小值，并选取出多个欧氏距离中的极大值点的下一个点，将多个互相关系数中的极小值作为第一分界点，将多个欧氏距离中的极大值点的下一个点作为第二分界点，该两个分界点将经emd分解得到的多个imf分量划分成三种模态分量：噪声imf分量(仅噪声)、混合imf分量(噪声 信号)及信号imf分量(仅信号)，其中，混合imf分量中的信号分量及信号imf分量是需要获取的，即去噪后的接地网干扰信号。
[0087]
互相关系数小于第一分界点时对应的imf分量为噪声imf分量，互相关系数大于或者等于第一分界点、且欧氏距离大于或者等于第二分界点时对应的imf分量为混合imf分量，欧氏距离小于第二分界点时对应的imf分量为信号imf分量。
[0088]
本发明实施例通过步骤s4和s5，提出了一种信号分量分界点的确立方法，即利用互相关系数和欧氏距离准则将emd分解的imf分量分解为噪声imf、混合imf和信号imf，解决了相关技术中emd分解中模态分解点不明确的问题。
[0089]
s6，将每个混合imf分量进行归一化最小均方差自适应降噪处理。
[0090]
具体地，在利用互相关系数和欧氏距离将emd分解的imf分量分解为噪声imf分量、混合imf分量和信号imf分量之后，可将噪声imf分量移除，将每个混合imf分量进行归一化最小均方差(normalized least mean square，简称nlms)自适应降噪处理，以去除混合imf分量中的噪声。
[0091]
通过对混合有噪声和信号的imf分量进行自适应降噪处理，根据不同imf分量的含噪能量比重设置不同的滤波参数，实现定制化的降噪处理，且计算复杂度低，有效简化信号的处理过程
[0092]
s7，将处理后的各个混合imf分量和信号imf分量进行重构，以得到去除噪声后的接地网干扰信号。
[0093]
具体地，在去除混合imf分量中的噪声分量之后，将处理后的各个混合imf分量和信号imf分量进行重构，即可最终获取到去除噪声后的接地网干扰信号。
[0094]
总的来说，如图2所示，通过执行上述步骤，首先通过形态学滤波器对接地网干扰信号(也可称为原始信号)进行形态学滤波处理，以消除其中的脉冲干扰，之后，对滤波处理后的干扰信号进行emd分解得到多个imf分量和一个余项，利用互相关系数和欧氏距离准则将emd分解的imfs(ims分量)分解为噪声imfs、混合imfs和信号imfs，将噪声imfs移除，将混合imfs进行nlms自适应滤波处理，之后，将处理后的混合imfs和信号imfs进行重构，即可得到去噪后的接地网干扰信号
[0095]
由此，本发明实施例的接地网干扰信号的去噪方法，可以有效地去除瞬时脉冲干扰，降低随机噪声对接地网干扰信号的影响，提高接地网干扰匹配精度；且可以通过互相关系数和欧氏距离准则有效确定emd分解中模态分界点，从而可以提高去噪可靠性。
[0096]
在一个实施例中，上述步骤s2可包括以下步骤：
[0097]
s21，利用形态学滤波器对接地网干扰信号进行腐蚀运算及膨胀运算，运算公式分别为；
[0098]
(b
s
θg)(n)＝min{b
s
(n m)
‑
g(m):m∈d
g
}，
[0099][0100]
其中，θ为腐蚀运算，为膨胀运算，b
s
(n)为接地网干扰信号，b
s
(n)的定义域为d
f
＝{0,1,2,3,
……
，n}，g(n)为一维结构元素序列，g(n)的定义域为d
g
＝{0,1,2,
……
，p}，其中p和n均是整数，且n≥p，n为接地网干扰信号的长度，g(m)为形态学滤波器的结构元素，m为结构元素长度。
[0101]
s22，对接地网干扰信号进行形态学开运算，并对接地网干扰信号进行形态学闭运算，运算公式分别为：
[0102][0103][0104]
其中，为开运算，
·
为闭运算。
[0105]
s23，利用交替混合滤波器对接地网干扰信号进行交替混合滤波运算。
[0106]
具体地，为同时发挥各基本算子的优点，本发明实时采取交替混合滤波器的组合方式。其实现方式为：
[0107]
交替滤波器：
[0108][0109][0110]
交替混合滤波器：
[0111]
[(b
s
)altmix(g)](n)＝[(b
s
)oc(g) (b
s
)co(g)](n)/2
[0112]
其中，altmix为交替混合滤波运算，oc为先开后闭的交替滤波运算，co为先闭后开的交替滤波运算。
[0113]
由此，提出了一种改进的形态学滤波处理方式，可以有效去除瞬时脉冲干扰，降低随机噪声对接地网干扰信号的影响，提高接地网干扰匹配精度。
[0114]
在一个实施例实施例中，上述步骤s3可包括以下步骤：
[0115]
s31，令k＝1，i＝1，选取滤波信号b
s
(t)
′
中的所有极大值点和所有极小值点，并将所有极大值点连接起来形成上包络线、将所有极小值点连接起来形成下包络线，确定上包络线与下包络线间的均值信号m(t)。
[0116]
具体地，在得到滤波信号b
s
(t)
′
，根据该信号的所有极大值、极小值，对应画出上、下包络线b
s
(i)和b
s
‑
(j)，进而确定上、下包络线间的均值得到均值包络线，该均值包络线即对应均值信号m(t)＝(b
s
(i) b
s
‑
(j))/2。
[0117]
s32，将滤波信号b
s
(t)
′
减去均值信号m(t)，得到中间信号h
i
(t)，判断中间信号h
i
(t)是否满足emd分解的条件：
[0118][0119]
且中间信号h
i
(t)的极值点个数和过零点个数间的差值是否满足不大于1，若满足，则记第i个本征模态函数imf
i
(t)为：
[0120]
imf
i
(t)＝h
i
(t)，
[0121]
若不满足，则令k＝k 1，i＝1，返回步骤s31中。
[0122]
s33，将滤波信号b
s
(t)
′
减去第i个本征模态函数imf
i
(t)，得到一个余项信号r
i
(t)，然后判断r
i
(t)是否满足余项条件，若满足，令emd分解的最后一个余项为r
l
(t)；若不满足，则令i＝i 1，i＝1，b
s
(t)'＝r
i
(t)返回步骤s32继续分解；最终分解后的信号为：
[0123][0124]
其中，h
i
(t)代表第i个imf分量，l表示所有imf分量的个数，r
l
(t)为余项。
[0125]
余项条件可以是：不能分解，也就是说，在得到一个余项信号之后，可判断该余项信号是否还能进行分解，如果能，则令i＝i 1，i＝1，b
s
(t)'＝r
i
(t)返回步骤s32继续分解；如果不能，则记emd分解的最后一个余项为r
l
(t)。
[0126]
在一个实施例中，上述步骤s4可包括以下步骤：
[0127]
s41，计算接地网干扰信号b
s
(t)分别与各个imf分量间的互相关系数，分别记为ρ1,ρ2,
……
ρ
l
。
[0128]
s42，计算接地网干扰信号b
s
(t)及各个imf分量的概率密度函数，分别记为pdf(b
s
(t))和pdf(h
i
(t))，计算pdf(b
s
(t))分别和各个pdf(h
i
(t))之间的欧氏距离，分别记为s(1),s(2),
……
，s(l)。
[0129]
进一步地，上述步骤s5可包括以下步骤：
[0130]
s51，确定各个互相关系数中的极小值ρ
min
＝[ρ1,ρ2,
……
ρ
l
]，及确定各个欧氏距离中的极大值s
max
＝[s(1),s(2),
……
，s(l)]；
[0131]
s52，令第一分界点m1＝ρ
min
、第二分界点m2＝s
max
1，确定互相关系数小于m1时对应的imf分量为噪声imf分量，确定互相关系数大于或者等于m1、且欧氏距离大于或者等于m2时对应的imf分量为混合imf分量，确定欧氏距离小于m2时对应的imf分量为信号imf分量。
[0132]
也就是说，通过m1和m2将各imfs分解为噪声imfs，混合imfs和信号imfs。当imf分量对应的互相关系数小于m1时，确定该imf分量为噪声imf分量；当imf分量对应的互相关系数大于或者等于m1、且对应的欧氏距离大于或者等于m2时，确定该imf分量为混合imf分量；当imf分量对应的欧氏距离小于m2时，确定该imf分量为信号imf分量。
[0133]
在一个实施例中，上述步骤s6可包括以下步骤：
[0134]
通过变步长nlms滤波算法对每个混合imf分量进行自适应降噪处理，处理过程为：
[0135]
设任意一个混合imf分量为信号x
imfi
(n)，对应的滤波输出信号为：
[0136]
y
imfi
(n)＝w
t
(n)x
imfi
(n)，
[0137]
误差信号为：
[0138]
e(n)＝d(n)
‑
y
imfi
(n)，
[0139]
权系数更新信号为：
[0140][0141]
其中，α(0＜α＜1)为收敛因子，β(0＜β≤1)是正常数，||x(n)||2为抽头输入向量x
imf
(n)的欧氏范数平方。
[0142]
在一个示例中，将nlms降噪处理后的混合imfs与信号imfs进行重构，得到的去除
噪声后的接地网干扰信号为：
[0143][0144]
其中，m1为第一分界点，m2为第二分界点，y
imfi
(t)为将混合imf分量进行归一化最小均方差自适应降噪处理后得到的第i个信号，h
i
(t)代表第i个imf分量，为信号imf分量之和，l表示所有imf分量的个数，r
l
(t)为余项。
[0145]
综上所述，本发明实施例的接地网干扰信号的去噪方法，可以有效地去除瞬时脉冲干扰，降低随机噪声对接地网干扰信号的影响，提高接地网干扰匹配精度；针对经验模态分解中模态分界点不明确的问题，提出一种新的确立准则，利用互相关系数和欧氏距离准则将emd分解的本征模态函数分量分解为噪声imf、混合imf和信号imf；对混合有噪声和信号的imf进行自适应降噪处理，根据不同imf分量的含噪能量比重设置不同的滤波参数，实现定制化的降噪处理，且计算复杂度低，有效简化信号的处理过程。
[0146]
图3是本发明实施例的接地网干扰信号的去噪装置的方框示意图。
[0147]
如图3所示，该接地网干扰信号的去噪装置100包括：采集模块10、滤波模块20、分解模块30、确定模块40、选取模块50、降噪模块60及重构模块70。
[0148]
其中，采集模块10，用于采集接地网干扰信号b
s
(t)；滤波模块20，用于利用形态学滤波器对接地网干扰信号b
s
(t)进行交替混合形态滤波处理，以得到滤波信号b
s
(t)
′
；分解模块30，用于对滤波信号b
s
(t)
′
进行经验模态分解处理，以得到多个本征模态函数imf分量和一个余项；确定模块40，用于确定接地网干扰信号b
s
(t)分别与各个imf分量间的互相关系数，并确定接地网干扰信号b
s
(t)分别与各个imf分量间的欧氏距离；选取模块50，用于选取各个互相关系数中的极小值，作为第一分界点，并选取各个欧氏距离中的极大值点的下一个点，作为第二分界点，根据第一分界点及第二分界点，确定多个imf分量中的噪声imf分量、混合及信号imf分量，其中，第二分界点与第一分界点不同，噪声imf分量中仅包含噪声，混合imf分量中包含噪声和信号，信号imf分量中仅包含信号；降噪模块60，用于将每个混合imf分量进行归一化最小均方差自适应降噪处理；重构模块70，用于将处理后的各个混合imf分量和信号imf分量进行重构，以得到去除噪声后的接地网干扰信号。
[0149]
在一个实施例中，滤波模块20，可包括：
[0150]
第一运算单元，用于利用形态学滤波器对接地网干扰信号进行腐蚀运算及膨胀运算，运算公式分别为；
[0151]
(b
s
θg)(n)＝min{b
s
(n m)
‑
g(m):m∈d
g
}，
[0152][0153]
其中，θ为腐蚀运算，为膨胀运算，b
s
(n)为接地网干扰信号，b
s
(n)的定义域为d
f
＝{0,1,2,3,
……
，n}，g(n)为一维结构元素序列，g(n)的定义域为d
g
＝{0,1,2,
……
，p}，其中p和n均是整数，且n≥p，n为接地网干扰信号的长度，g(m)为形态学滤波器的结构元素，m为结构元素长度；
[0154]
第二运算单元，用于对接地网干扰信号进行形态学开运算，并对接地网干扰信号
b
s
(n)进行形态学闭运算，运算公式分别为：
[0155][0156][0157]
其中，为开运算，
·
为闭运算；
[0158]
第三运算单元，用于利用交替混合滤波器对接地网干扰信号进行交替混合滤波运算，运算公式为：
[0159]
[(b
s
)altmix(g)](n)＝[(b
s
)oc(g) (b
s
)co(g)](n)/2，
[0160]
其中，altmix为交替混合滤波运算，oc为先开后闭的交替滤波运算，co为先闭后开的交替滤波运算，
[0161][0162][0163]
在一个实施例中，分解模块30，可包括：
[0164]
第一确定单元，用于令k＝1，i＝1，选取滤波信号b
s
(t)
′
中的所有极大值点和所有极小值点，并将所有极大值点连接起来形成上包络线、将所有极小值点连接起来形成下包络线，确定上包络线与下包络线间的均值信号m(t)；
[0165]
第四运算单元，用于将滤波信号b
s
(t)
′
减去均值信号m(t)，得到中间信号h
i
(t)，判断中间信号h
i
(t)是否满足emd分解的条件：
[0166][0167]
且中间信号h
i
(t)的极值点个数和过零点个数间的差值是否满足不大于1，若满足，则记第i个本征模态函数imf
i
(t)为：
[0168]
imf
i
(t)＝h
i
(t)，
[0169]
若不满足，则令k＝k 1，i＝1，返回步骤s31中；
[0170]
第五运算单元，用于将滤波信号b
s
(t)
′
减去第i个本征模态函数imf
i
(t)，得到一个余项信号r
i
(t)，然后判断r
i
(t)是否满足余项条件，若满足，令emd分解的最后一个余项为r
l
(t)；若不满足，则令i＝i 1，i＝1，b
s
(t)'＝r
i
(t)返回步骤s32继续分解；最终分解后的信号为：
[0171][0172]
其中，h
i
(t)代表第i个imf分量，l表示所有imf分量的个数，r
l
(t)为余项。
[0173]
在一个实施例中，确定模块40，可包括：
[0174]
第六运算单元，用于计算接地网干扰信号b
s
(t)分别与各个imf分量间的互相关系数，分别记为ρ1,ρ2,
……
ρ
l
；
[0175]
第七运算单元，用于计算接地网干扰信号b
s
(t)及各个imf分量的概率密度函数，分别记为pdf(b
s
(t))和pdf(h
i
(t))，计算pdf(b
s
(t))分别和各个pdf(h
i
(t))之间的欧氏距离，分别记为s(1),s(2),
……
，s(l)。
[0176]
在一个实施例中，选取模块50，可包括：第二确定单元，用于确定各个互相关系数中的极小值ρ
min
＝[ρ1,ρ2,
……
ρ
l
]，及确定各个欧氏距离中的极大值s
max
＝[s(1),s(2),
……
，s(l)]；
[0177]
第三确定单元，用于令第一分界点m1＝ρ
min
、第二分界点m2＝s
max
1，确定互相关系数小于m1时对应的imf分量为噪声imf分量，确定互相关系数大于或者等于m1、且欧氏距离大于或者等于m2时对应的imf分量为混合imf分量，确定欧氏距离小于m2时对应的imf分量为信号imf分量。
[0178]
在一个实施例中，降噪模块60，具体可用于：通过变步长nlms滤波算法对每个混合imf分量进行自适应降噪处理，处理过程为：
[0179]
设任意一个混合imf分量为信号x
imfi
(n)，对应的滤波输出信号为：
[0180]
y
imfi
(n)＝w
t
(n)x
imfi
(n)，
[0181]
误差信号为：
[0182]
e(n)＝d(n)
‑
y
imfi
(n)，
[0183]
权系数更新信号为：
[0184][0185]
其中，α(0＜α＜1)为收敛因子，β(0＜β≤1)是正常数，||x(n)||2为抽头输入向量x
imf
(n)的欧氏范数平方。
[0186]
在一个实施例中，去除噪声后的接地网干扰信号为：
[0187][0188]
其中，m1为第一分界点，m2为第二分界点，y
imfi
(t)为将混合imf分量进行归一化最小均方差自适应降噪处理后得到的第i个信号，h
i
(t)代表第i个imf分量，为信号imf分量之和，l表示所有imf分量的个数，r
l
(t)为余项。
[0189]
需要说明的是，接地网干扰信号的去噪装置的其他具体实施方式可参见接地网干扰信号的去噪方法的具体实施方式，为避免冗余，此处不再赘述。
[0190]
本发明实施例的接地网干扰信号的去噪装置，可以有效地去除瞬时脉冲干扰，降低随机噪声对接地网干扰信号的影响，提高接地网干扰匹配精度；且可以通过互相关系数和欧氏距离准则有效确定emd分解中模态分界点，从而可以提高去噪可靠性。
[0191]
对应上述实施例，本发明还提出一种计算机设备。
[0192]
本发明实施例的计算机设备包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行该计算机程序时，可实现根据本发明上述实施例所述的接地网干扰信号的去噪方法。
[0193]
根据本发明实施例的计算机设备，处理器执行存储在存储器上的计算机程序时，首先通过形态学滤波器对接地网干扰信号(也可称为原始信号)进行形态学滤波处理，以消除其中的脉冲干扰，之后，对滤波处理后的干扰信号进行emd分解得到多个imf分量和一个
余项，利用互相关系数和欧氏距离准则将emd分解的imfs(ims分量)分解为噪声imfs、混合imfs和信号imfs，将混合imfs进行nlms自适应滤波处理，之后，将处理后的混合imfs和信号imfs进行重构，即可得到去噪后的接地网干扰信号。
[0194]
本发明实施例的计算机设备，在执行计算机程序时，可以有效地去除瞬时脉冲干扰，降低随机噪声对接地网干扰信号的影响，提高接地网干扰匹配精度；且可以通过互相关系数和欧氏距离准则有效确定emd分解中模态分界点，从而可以提高去噪可靠性。
[0195]
对应上述实施例，本发明还提出一种非临时性计算机可读存储介质。
[0196]
本发明实施例的非临时性计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时可实现根据本发明上述实施例所述的接地网干扰信号的去噪方法。
[0197]
根据本发明实施例的非临时性计算机可读存储介质，处理器执行存储在其上的计算机程序时，首先首先通过形态学滤波器对接地网干扰信号(也可称为原始信号)进行形态学滤波处理，以消除其中的脉冲干扰，之后，对滤波处理后的干扰信号进行emd分解得到多个imf分量和一个余项，利用互相关系数和欧氏距离准则将emd分解的imfs(ims分量)分解为噪声imfs、混合imfs和信号imfs，将混合imfs进行nlms自适应滤波处理，之后，将处理后的混合imfs和信号imfs进行重构，即可得到去噪后的接地网干扰信号。
[0198]
本发明实施例的非临时性计算机可读存储介质，在执行计算机程序时，可以有效地去除瞬时脉冲干扰，降低随机噪声对接地网干扰信号的影响，提高接地网干扰匹配精度；且可以通过互相关系数和欧氏距离准则有效确定emd分解中模态分界点，从而可以提高去噪可靠性。
[0199]
在本发明的描述中，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。“多个”的含义是两个或两个以上，除非另有明确具体的限定。
[0200]
在本发明中，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或成一体；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通或两个元件的相互作用关系。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本发明中的具体含义。
[0201]
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
[0202]
流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为，表示包括一个或更多个用于实现特定逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分，并且本发明的优选实施方式的范围包括另外的实现，其中可以不按所示出或讨论的顺序，包括根据所涉及的功能按基本同时的方式或按相反的顺序，来执行功能，这应被本发明的实施例所属技术领域的技术人员所理解。
[0203]
在流程图中表示或在此以其他方式描述的逻辑和/或步骤，例如，可以被认为是用
于实现逻辑功能的可执行指令的定序列表，可以具体实现在任何计算机可读介质中，以供指令执行系统、装置或设备(如基于计算机的系统、包括处理器的系统或其他可以从指令执行系统、装置或设备取指令并执行指令的系统)使用，或结合这些指令执行系统、装置或设备而使用。就本说明书而言，“计算机可读介质”可以是任何可以包含、存储、通信、传播或传输程序以供指令执行系统、装置或设备或结合这些指令执行系统、装置或设备而使用的装置。计算机可读介质的更具体的示例(非穷尽性列表)包括以下：具有一个或多个布线的电连接部(电子装置)，便携式计算机盘盒(磁装置)，随机存取存储器(ram)，只读存储器(rom)，可擦除可编辑只读存储器(eprom或闪速存储器)，光纤装置，以及便携式光盘只读存储器(cdrom)。另外，计算机可读介质甚至可以是可在其上打印所述程序的纸或其他合适的介质，因为可以例如通过对纸或其他介质进行光学扫描，接着进行编辑、解译或必要时以其他合适方式进行处理来以电子方式获得所述程序，然后将其存储在计算机存储器中。
[0204]
应当理解，本发明的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施方式中，多个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件来实现。例如，如果用硬件来实现，和在另一实施方式中一样，可用本领域公知的下列技术中的任一项或他们的组合来实现：具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路，具有合适的组合逻辑门电路的专用集成电路，可编程门阵列(pga)，现场可编程门阵列(fpga)等。
[0205]
本技术领域的普通技术人员可以理解实现上述实施例方法携带的全部或部分步骤是可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，该程序在执行时，包括方法实施例的步骤之一或其组合。
[0206]
此外，在本发明各个实施例中的各功能单元可以集成在一个处理模块中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个模块中。上述集成的模块既可以采用硬件的形式实现，也可以采用软件功能模块的形式实现。所述集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时，也可以存储在一个计算机可读取存储介质中。
[0207]
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。