一种利用Gibuss效应的PLC信号滤波方法和系统与流程

本发明涉及电力领域，尤其涉及一种plc信号滤波方法和系统。

背景技术：

电力线通信，相比各种有线通信技术，无需重新布线，易于组网等优点，具有广阔的应用前景。电力线通信技术分为窄带电力线通信(narrowbandoverpowerline，npl)和宽带电力线通信(broadbandoverpowerline，bpl)；窄带电力线通信是指带宽限定在3k～500khz的电力线载波通信技术；电力线通信技术包括欧洲cenelec的规定带宽(3～148.5khz)，美国联邦通讯委员会(fcc)的规定带宽(9～490khz)，日本无线工业及商贸联合会(associationofradioindustriesandbusinesses，arib)的规定带宽(9～450khz)，和中国的规定带宽(3～500khz)。窄带电力线通信技术多采用单载波调制技术，如psk技术，dsss技术和线性调频chirp等技术，通信速率小于1mbits/s；宽带电力线通信技术指带宽限定在1.6～30mhz之间、通信速率通常在1mbps以上的电力线载波通信技术，采用以ofdm为核心的多种扩频通信技术。

虽然电力线通信系统有着广泛的应用，且技术相对成熟，但是电力线通信系统中大量的分支和电气设备，会在电力线信道中产生大量的噪声；而其中随机脉冲噪声具有很大的随机性，噪声强度高，对电力线通信系统造成严重破坏，因此，针对随机脉冲噪声的抑制技术，一直是国内外学者研究的重点；而且噪声模型并不符合高斯分布。因此，传统的针对高斯噪声设计的通信系统不再适用于电力线载波通信系统，必须研究相应的噪声抑制技术，以提高电力线通信系统信噪比，降低误码率，保证电力线通信系统质量。在实际应用中，一些简单的非线性技术经常被应用于消除电力线信道噪声，如clipping、blanking和clipping/blanking技术，但是这些研究方法都必须在一定的信噪比情况下才能良好工作，仅仅考虑了冲击噪声的消除，在电力线通信系统中，某些商用电力线发送器的特征是低发射功率，在一些特殊情况，发射功率甚至可能会低于18w，因此，在某些特殊情况，信号将会淹没在大量噪声中，导致电力线通信系统低信噪比情况。

随着非线性电器的应用和普及，中低压输配电网络中背景噪声呈现出较为明显的非平稳性和非高斯特性，常用的低通滤波器在非平稳和非高斯噪声环境中难以达到理想的滤波效果，很难滤除非平稳非高斯噪声，严重影响了plc通信系统的性能。

技术实现要素：

本发明的目的是提供一种利用gibuss效应的plc信号滤波方法和系统，所提出的方法利用了plc调制信号和背景噪声在gibuss变换领域中的差异，通过gibuss系数区分plc调制信号和背景噪声。所提出的方法具有较好的噪声滤除性能，计算也非常简单。

为实现上述目的，本发明提供了如下方案：

一种利用gibuss效应的plc信号滤波方法，包括：

步骤1，输入实测的plc信号序列s；

步骤2，根据gibuss效应对所述plc信号序列s进行滤除噪声处理，滤除噪声后的信号序列为snew；具体为，snew＝sopt+[s^ts]^-1(s-sopt)；其中，sopt为最佳预测矢量。

一种利用gibuss效应的plc信号滤波系统，包括：

获取模块，输入实测的plc信号序列s；

滤波模块，根据gibuss效应对所述plc信号序列s进行滤除噪声处理，滤除噪声后的信号序列为snew；具体为，snew＝sopt+[s^ts]^-1(s-sopt)；其中，sopt为最佳预测矢量。

根据本发明提供的具体实施例，本发明公开了以下技术效果：

虽然电力线通信系统有着广泛的应用，且技术相对成熟，但是电力线通信系统中大量的分支和电气设备，会在电力线信道中产生大量的噪声；随着非线性电器的应用和普及，中低压输配电网络中背景噪声呈现出较为明显的非平稳性和非高斯特性，常用的低通滤波器在非平稳和非高斯噪声环境中难以达到理想的滤波效果，很难滤除非平稳非高斯噪声，严重影响了plc通信系统的性能。

本发明的目的是提供一种利用gibuss效应的plc信号滤波方法和系统，所提出的方法利用了plc调制信号和背景噪声在gibuss变换领域中的差异，通过gibuss系数区分plc调制信号和背景噪声。所提出的方法具有较好的噪声滤除性能，计算也非常简单。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍。显而易见，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明的方法流程示意图；

图2为本发明的系统结构示意图；

图3为本发明具体实施案例的流程示意图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。

图1一种利用gibuss效应的plc信号滤波方法的流程示意图

图1为本发明一种利用gibuss效应的plc信号滤波方法的流程示意图。如图1所示，所述的一种利用gibuss效应的plc信号滤波方法具体包括以下步骤：

步骤1，输入实测的plc信号序列s；

步骤2，根据gibuss效应对所述plc信号序列s进行滤除噪声处理，滤除噪声后的信号序列为snew；具体为，snew＝sopt+[s^ts]^-1(s-sopt)；其中，sopt为最佳预测矢量。

所述步骤2之前，所述方法还包括：

步骤3，求取所述最佳预测矢量sopt。

所述步骤3包括：

步骤301，求取混叠矩阵p，具体为：

p＝diag[p1,p2,…,pn]

其中：

n：所述信号序列s的长度

i:虚数单位

tn：混叠因子，[1，n]范围内随机分布的整数

步骤302，求取延迟矢量，具体为：

p＝s^td

其中：

d＝ifft[p]

ifft[p]:矩阵p的逆傅里叶变换

步骤303，迭代初始化，具体为：

y1＝s：初始化的迭代矢量

n＝1:初始化的迭代控制参量

步骤304，迭代更新，具体为：

其中：

第n+1步的迭代矢量

第n步的迭代矢量

x:中间变量

gibuss因子

反因子

s＝[s1，s2，…，sn]：所述信号序列s

n：所述信号序列s的长度

ms：所述信号序列s的均值

σs：所述信号序列s的均方差

步骤305，迭代终止，具体为：

所述迭代控制参量n加1，并返回步骤304持续迭代，直至相邻两次的迭代矢量之差小于千分之一为止。此时所述迭代控制参量n＝k，并得到最佳预测矢量sopt＝yn+1。

图2一种利用gibuss效应的plc信号滤波系统的结构意图

图2为本发明一种利用gibuss效应的plc信号滤波系统的结构示意图。如图2所示，所述一种利用gibuss效应的plc信号滤波系统包括以下结构：

获取模块401，输入实测的plc信号序列s；

滤波模块402，根据gibuss效应对所述plc信号序列s进行滤除噪声处理，滤除噪声后的信号序列为snew；具体为，snew＝sopt+[s^ts]^-1(s-sopt)；其中，sopt为最佳预测矢量。

所述的系统，还包括：

计算模块403，求取所述最佳预测矢量sopt。

所述计算模块403，还包括下列单元，具体为：

第一计算单元4031，求取混叠矩阵p，具体为：

p＝diag[p1,p2,…,pn]

其中：

n：所述信号序列s的长度

i:虚数单位

tn：混叠因子，[1，n]范围内随机分布的整数

第二计算单元4032，求取延迟矢量，具体为：

p＝s^td

其中：

d＝ifft[p]

ifft[p]:矩阵p的逆傅里叶变换

迭代初始单元4033，迭代初始化，具体为：

y1＝s：初始化的迭代矢量

n＝1:初始化的迭代控制参量

迭代更新单元4034，迭代更新，具体为：

其中：

第n+1步的迭代矢量

第n步的迭代矢量

x:中间变量

gibuss因子

反因子

s＝[s1，s2，…，sn]：所述信号序列s

n：所述信号序列s的长度

ms：所述信号序列s的均值

σs：所述信号序列s的均方差

迭代终止单元4035，迭代终止，具体为：

所述迭代控制参量n加1，并返回迭代更新单元4034持续迭代，直至相邻两次的迭代矢量之差小于千分之一为止。此时所述迭代控制参量n＝k，并得到最佳预测矢量sopt＝yn+1。

下面提供一个具体实施案例，进一步说明本发明的方案

图3为本发明具体实施案例的流程示意图。如图3所示，具体包括以下步骤：

1.输入实测的plc信号序列

s＝[s1,s2,…,sn-1,sn]

其中：

s：实测的plc信号数据序列，长度为nsi,i＝1,2,…,n:序号为i的实测plc信号

2.求取混叠矩阵

p＝diag[p1,p2,…,pn]

其中：

n：所述信号序列s的长度

i:虚数单位

tn：混叠因子，[1，n]范围内随机分布的整数

3.求取延迟矢量

p＝s^td

其中：

d＝ifft[p]

ifft[p]:矩阵p的逆傅里叶变换

4.迭代初始化

y1＝s：初始化的迭代矢量

n＝1:初始化的迭代控制参量

5.迭代更新

其中：

第n+1步的迭代矢量

第n步的迭代矢量

x:中间变量

gibuss因子

反因子

s＝[s1，s2，…，sn]：所述信号序列s

n：所述信号序列s的长度

ms：所述信号序列s的均值

σs：所述信号序列s的均方差

6.迭代终止

所述迭代控制参量n加1，并返回迭代更新持续迭代，直至相邻两次的迭代矢量之差小于千分之一为止。此时所述迭代控制参量n＝k，并得到最佳预测矢量sopt＝yn+1。

7.滤波

根据gibuss效应对所述plc信号序列s进行滤除噪声处理，滤除噪声后的信号序列为snew；具体为，snew＝sopt+[s^ts]^-1(s-sopt)；其中，sopt为最佳预测矢量。

本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述较为简单，相关之处参见方法部分说明即可。

本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。