一种利用一阶统计量的plc信号滤波方法和系统
阅读说明:本技术 一种利用一阶统计量的plc信号滤波方法和系统 (PLC signal filtering method and system using first-order statistics ) 是由 翟明岳 于 2019-09-06 设计创作,主要内容包括:本发明的实施例公开一种利用一阶统计量的PLC信号滤波方法和系统,所述方法包括:步骤1,输入实测的PLC信号序列S;步骤2,根据统计量性质对所述信号序列S进行滤除噪声处理,滤除噪声后的信号序列为S<Sub>NEW</Sub>。具体为,<Image he="121" wi="646" file="DDA0002194784610000011.GIF" imgContent="drawing" imgFormat="GIF" orientation="portrait" inline="no"></Image>其中,V<Sub>OPT</Sub>为最佳预测矢量;γ为修正因子;<Image he="76" wi="186" file="DDA0002194784610000012.GIF" imgContent="drawing" imgFormat="GIF" orientation="portrait" inline="no"></Image>为修正矢量。(the embodiment of the invention discloses a PLC signal filtering method and a system utilizing first-order statistics, wherein the method comprises the following steps: step 1, inputting an actually measured PLC signal sequence S; and 2, carrying out noise filtering processing on the signal sequence S according to the statistical quality property, wherein the signal sequence after noise filtering is SNEW. Specifically, the VOPT is an optimal prediction vector; gamma is a correction factor; is a correction vector.)
技术领域
本发明涉及电力领域,尤其涉及一种PLC信号滤波方法和系统。
背景技术
电力线通信,相比各种有线通信技术,无需重新布线,易于组网等优点,具有广阔的应用前景。电力线通信技术分为窄带电力线通信(Narrowband over power line,NPL)和宽带电力线通信(Broadband over power line,BPL);窄带电力线通信是指带宽限定在3k~500kHz的电力线载波通信技术;电力线通信技术包括欧洲CENELEC的规定带宽(3~148.5kHz),美国联邦通讯委员会(FCC)的规定带宽(9~490kHz),日本无线工业及商贸联合会(Association of Radio Industries and Businesses,ARIB)的规定带宽(9~450kHz),和中国的规定带宽(3~500kHz)。窄带电力线通信技术多采用单载波调制技术,如PSK技术,DSSS技术和线性调频Chirp等技术,通信速率小于1Mbits/s;宽带电力线通信技术指带宽限定在1.6~30MHz之间、通信速率通常在1Mbps以上的电力线载波通信技术,采用以OFDM为核心的多种扩频通信技术。
虽然电力线通信系统有着广泛的应用,且技术相对成熟,但是电力线通信系统中大量的分支和电气设备,会在电力线信道中产生大量的噪声;而其中随机脉冲噪声具有很大的随机性,噪声强度高,对电力线通信系统造成严重破坏,因此,针对随机脉冲噪声的抑制技术,一直是国内外学者研究的重点;而且噪声模型并不符合高斯分布。因此,传统的针对高斯噪声设计的通信系统不再适用于电力线载波通信系统,必须研究相应的噪声抑制技术,以提高电力线通信系统信噪比,降低误码率,保证电力线通信系统质量。在实际应用中,一些简单的非线性技术经常被应用于消除电力线信道噪声,如Clipping、Blanking和Clipping/Blanking技术,但是这些研究方法都必须在一定的信噪比情况下才能良好工作,仅仅考虑了冲击噪声的消除,在电力线通信系统中,某些商用电力线发送器的特征是低发射功率,在一些特殊情况,发射功率甚至可能会低于18w,因此,在某些特殊情况,信号将会淹没在大量噪声中,导致电力线通信系统低信噪比情况。
随着非线性电器的应用和普及,中低压输配电网络中背景噪声呈现出较为明显的非平稳性和非高斯特性,常用的低通滤波器在非平稳和非高斯噪声环境中难以达到理想的滤波效果,很难滤除非平稳非高斯噪声,严重影响了PLC通信系统的性能。
发明内容
本发明的目的是提供一种利用一阶统计量的PLC信号滤波方法和系统,所提出的方法利用了PLC调制信号和背景噪声在信号统计量方面的差异,通过一阶统计量的性质分离PLC调制信号和背景噪声。所提出的方法具有较好的噪声滤除性能,计算也非常简单。
为实现上述目的,本发明提供了如下方案:
一种利用一阶统计量的PLC信号滤波方法,包括:
步骤1,输入实测的PLC信号序列S;
步骤2,根据统计量性质对所述信号序列S进行滤除噪声处理,滤除噪声后的信号序列为SNEW。具体为,其中,VOPT为最佳预测矢量;γ为修正因子;为修正矢量。
一种利用一阶统计量的PLC信号滤波系统,包括:
获取模块,输入实测的PLC信号序列S;
滤波模块,根据统计量性质对所述信号序列S进行滤除噪声处理,滤除噪声后的信号序列为SNEW。具体为,其中,VOPT为最佳预测矢量;γ为修正因子;为修正矢量。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
虽然电力线通信系统有着广泛的应用,且技术相对成熟,但是电力线通信系统中大量的分支和电气设备,会在电力线信道中产生大量的噪声;随着非线性电器的应用和普及,中低压输配电网络中背景噪声呈现出较为明显的非平稳性和非高斯特性,常用的低通滤波器在非平稳和非高斯噪声环境中难以达到理想的滤波效果,很难滤除非平稳非高斯噪声,严重影响了PLC通信系统的性能。
本发明的目的是提供一种利用一阶统计量的PLC信号滤波方法和系统,所提出的方法利用了PLC调制信号和背景噪声在信号统计量方面的差异,通过一阶统计量的性质分离PLC调制信号和背景噪声。所提出的方法具有较好的噪声滤除性能,计算也非常简单。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍。显而易见,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明的方法流程示意图;
图2为本发明的系统结构示意图;
图3为本发明具体实施案例的流程示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述。显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1一种利用一阶统计量的PLC信号滤波方法的流程示意图
图1为本发明一种利用一阶统计量的PLC信号滤波方法的流程示意图。如图1所示,所述的一种利用一阶统计量的PLC信号滤波方法具体包括以下步骤:
步骤1,输入实测的PLC信号序列S;
步骤2,根据统计量性质对所述信号序列S进行滤除噪声处理,滤除噪声后的信号序列为SNEW。具体为,其中,VOPT为最佳预测矢量;γ为修正因子;为修正矢量。
所述步骤2之前,所述方法还包括:
步骤3,求取所述最佳预测矢量VOPT、修正因子γ和修正矢量
所述步骤3包括:
步骤301,初始化重叠矩阵具体为:
k=0
其中:
重叠矩阵P0中的第i行第j列元素
行参数
列参数
k:迭代控制参数
SNR:所述信号序列S的信噪比
下取整运算
si:所述信号序列S中的第i个元素
步骤302,更新所述重叠矩阵,具体为:
Pk+1=[I+γkD2(Pk)]-1P0
其中:
I:单位矩阵
对角矩阵D(Pk)的第m行第m列元素
第k步迭代的修正因子
步骤303,迭代结束,具体为:
所述迭代控制参量k加1,返回到步骤302继续迭代更新,直到相邻两次迭代结果之差小于千分之一为止,此时所述迭代控制参量k=K,得到所述最佳预测矢量VOPT=VK+1。
步骤304,求取所述修正因子γ,具体为:
γ=max[γ1,γ2,…,γK+1]
步骤305,求取所述修正矢量具体为:
其中:
Ωi:所述修正矢量中的第i个元素[i=1,2,…,N]
所述最佳预测矢量VOPT中的第i个元素[i=1,2,…,N]
mS:所述信号序列S的均值
σS:所述信号序列S的均方差
图2一种利用一阶统计量的PLC信号滤波系统的结构意图
图2为本发明一种利用一阶统计量的PLC信号滤波系统的结构示意图。如图2所示,所述一种利用一阶统计量的PLC信号滤波系统包括以下结构:
获取模块401,输入实测的PLC信号序列S;
滤波模块402,根据统计量性质对所述信号序列S进行滤除噪声处理,滤除噪声后的信号序列为SNEW。具体为,其中,VOPT为最佳预测矢量;γ为修正因子;为修正矢量。
所述的系统,还包括:
计算模块403,求取所述最佳预测矢量VOPT、修正因子γ和修正矢量
所述计算模块403,还包括下列单元,具体为:
初始化单元4031,初始化重叠矩阵具体为:
k=0
其中:
重叠矩阵P0中的第i行第j列元素
行参数
列参数
k:迭代控制参数
SNR:所述信号序列S的信噪比
下取整运算
si:所述信号序列S中的第i个元素
迭代单元4032,更新所述重叠矩阵,具体为:
Pk+1=[I+γkD2(Pk)]-1P0
其中:
I:单位矩阵
对角矩阵D(Pk)的第m行第m列元素
第k步迭代的修正因子
判断单元4033,迭代结束,具体为:
所述迭代控制参量k加1,返回到迭代模块4032继续迭代更新,直到相邻两次迭代结果之差小于千分之一为止,此时所述迭代控制参量k=K,得到所述最佳预测矢量VOPT=VK +1。
第一计算单元4034,求取所述修正因子γ,具体为:
γ=max[γ1,γ2,…,γK+1]
第二计算单元4035,求取所述修正矢量具体为:
其中:
Ωi:所述修正矢量中的第i个元素[i=1,2,…,N]
所述最佳预测矢量VOPT中的第i个元素[i=1,2,...,N]
mS:所述信号序列S的均值
σS:所述信号序列S的均方差
下面提供一个具体实施案例,进一步说明本发明的方案
图3为本发明具体实施案例的流程示意图。如图3所示,具体包括以下步骤:
1.输入实测的PLC信号序列
S=[s1,s2,…,sN-1,sN]
其中:
S:实测的PLC信号数据序列,长度为N
si,i=1,2,…,N:序号为i的实测PLC信号
2.初始化重叠矩阵
k=0
其中:
重叠矩阵P0中的第i行第j列元素
行参数
列参数
3.更新重叠矩阵
Pk+1=[I+γkD2(Pk)]-1P0
其中:
I:单位矩阵
对角矩阵D(Pk)的第m行第m列元素
第k步迭代的修正因子
4.求取修正因子
γ=max[γ1,γ2,…,γK+1]
5.求取修正矢量
其中:
Ωi:所述修正矢量中的第i个元素[i=1,2,…,N]
所述最佳预测矢量VOPT中的第i个元素[i=1,2,…,N]
mS:所述信号序列S的均值
σS:所述信号序列S的均方差
6.滤波
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述较为简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。