CN111934716B - 一种电力线通信信号滤波方法及系统 - Google Patents

一种电力线通信信号滤波方法及系统 Download PDF

Info

Publication number
CN111934716B
CN111934716B CN202010971971.6A CN202010971971A CN111934716B CN 111934716 B CN111934716 B CN 111934716B CN 202010971971 A CN202010971971 A CN 202010971971A CN 111934716 B CN111934716 B CN 111934716B
Authority
CN
China
Prior art keywords
sequence
lagrangian
power line
line communication
matrix
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010971971.6A
Other languages
English (en)
Other versions
CN111934716A (zh
Inventor
翟明岳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangdong University of Petrochemical Technology
Original Assignee
Guangdong University of Petrochemical Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guangdong University of Petrochemical Technology filed Critical Guangdong University of Petrochemical Technology
Publication of CN111934716A publication Critical patent/CN111934716A/zh
Application granted granted Critical
Publication of CN111934716B publication Critical patent/CN111934716B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B3/00Line transmission systems
    • H04B3/02Details
    • H04B3/46Monitoring; Testing
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B3/00Line transmission systems
    • H04B3/54Systems for transmission via power distribution lines

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Power Engineering (AREA)
  • Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
  • Noise Elimination (AREA)

Abstract

本发明涉及一种电力线通信信号滤波方法及系统,涉及电力系统技术领域。该方法包括:利用实时的电力线通信信号序列中元素的差值确定信号差分序列;利用电力线通信信号序列和信号差分序列确定期望序列;利用信号差分序列和期望序列确定拉格朗日函数矩阵;利用拉格朗日函数矩阵确定拉格朗日因子矢量;利用拉格朗日函数矩阵和拉格朗日因子矢量确定拉格朗日最佳矩阵;利用拉格朗日最佳矩阵对电力线通信信号序列进行滤波,得到滤波后的信号序列。本发明通过引入拉格朗日因子矢量,利用拉格朗日函数矩阵和拉格朗日因子矢量确定拉格朗日最佳矩阵,利用拉格朗日最佳矩阵对电力线通信信号序列进行滤波,有效地滤除了PLC信号中的脉冲噪声。

Description

一种电力线通信信号滤波方法及系统
技术领域
本发明涉及电力系统技术领域,特别是涉及一种电力线通信信号滤波方法及系统。
背景技术
电力线通信(Power Line Communications,PLC),相比各种有线通信技术,具有无需重新布线和易于组网等优点,具有广阔的应用前景。电力线通信技术分为窄带电力线通信(Narrowband over power line,NPL)和宽带电力线通信(Broadband over power line,BPL)。窄带电力线通信是指带宽限定在3kHz~500kHz的电力线载波通信技术,包括欧洲电工标准化委员会(European Committee for Electrotechnical Standardization,CENELEC)的规定带宽(3~148.5kHz),美国联邦通讯委员会(FCC)的规定带宽(9~490kHz),日本无线工业及商贸联合会(Association of Radio Industries and Businesses,ARIB)的规定带宽(9~450kHz),和中国的规定带宽(3~500kHz)。窄带电力线通信技术多采用单载波调制技术,如相移键控(Phase-Shift Keying,PSK)技术,直接序列扩频(DirectSequence Spread Spectrum,DSSS)技术和线性调频信号(Chirp)等技术,通信速率小于1Mbits/s。宽带电力线通信技术指带宽限定在1.6MHz~30MHz之间、通信速率通常在1Mbps以上的电力线载波通信技术,采用以正交频分复用技术(Orthogonal Frequency DivisionMultiplexing,OFDM)为核心的多种扩频通信技术。
电力线通信(Power Line Communications,PLC)系统是利用中低压电力线网络传输高频通信信号的系统。虽然电力线通信系统有着广泛的应用,且技术相对成熟,但是电力线通信系统中大量的分支和电气设备,会在电力线信道中产生大量的噪声。其中随机脉冲噪声具有很大的随机性,噪声强度高,对电力线通信系统造成严重破坏,所以随机脉冲噪声的抑制技术,一直是国内外学者研究的重点。而且,随机脉冲噪声的噪声模型并不符合高斯分布,因此,传统的针对高斯噪声设计的通信系统不适用于电力线载波通信系统。在实际应用中,一些简单的非线性技术经常被应用于消除电力线信道噪声,如Clipping、Blanking和Clipping/Blanking技术。但是,这些非线性技术的研究方法都必须在一定的信噪比情况下才能良好工作,而且仅仅考虑了冲击噪声的消除。在电力线通信系统中,某些商用电力线发送器的特征是低发射功率,在一些特殊情况,发射功率甚至可能会低于18w。所以,在某些特殊情况时,PLC信号将会淹没在大量噪声中,导致电力线通信系统低信噪比情况。因此,在PLC通信系统中,对通信信号进行滤波是很重要的一步,常用的滤波方法是低通滤波器。
随着非线性电器的应用和普及(例如使用马达的果汁机、咖啡机和豆浆机等),PLC通信信号中出现了越来越多的脉冲噪声,此类噪声的瞬时功率很大,呈现出较为明显的非平稳性和非高斯特性,但常用的低通滤波器在滤除脉冲噪声方面表现不佳,严重影响了PLC通信系统的性能。
发明内容
本发明的目的是提供一种电力线通信信号滤波方法及系统,有效地滤除PLC信号中的脉冲噪声,保证电力线通信系统质量。
为实现上述目的,本发明提供了如下方案:
一种电力线通信信号滤波方法,包括:
获取实时的电力线通信信号序列;
利用所述电力线通信信号序列中相邻两个元素的差值确定信号差分序列;
利用所述电力线通信信号序列和所述信号差分序列确定期望序列;
利用所述信号差分序列和所述期望序列确定拉格朗日函数矩阵;
利用所述拉格朗日函数矩阵确定拉格朗日因子矢量;
利用所述拉格朗日函数矩阵和所述拉格朗日因子矢量确定拉格朗日最佳矩阵;
利用所述拉格朗日最佳矩阵对所述电力线通信信号序列进行滤波,得到滤波后的信号序列。
可选的,所述利用所述电力线通信信号序列中相邻两个元素的差值确定信号差分序列,具体包括:
计算所述电力线通信信号序列中相邻两个元素的差值;
通过所述差值确定信号差分序列,所述信号差分序列为:ΔSN=[0,s2-s1,...,sn-sn-1,...,sN-sN-1];
式中,所述电力线通信信号序列为S=[s1,s2,...,sn,...,sN],N表示所述电力线通信信号序列的长度;sn表示所述电力线通信信号序列的第n个元素,n=1,2,...,N。
可选的,所述利用所述电力线通信信号序列和所述信号差分序列确定期望序列,具体包括:
利用所述电力线通信信号序列和所述信号差分序列求解公式
Figure GDA0003158678640000031
得到期望序列;
式中,XEXP表示期望序列,ΔSN表示所述信号差分序列,σS表示所述电力线通信信号序列的均方差;σΔS为所述信号差分序列的均方差,S表示所述电力线通信信号序列。
可选的,所述利用所述信号差分序列和所述期望序列确定拉格朗日函数矩阵,具体包括:
利用所述信号差分序列和所述期望序列求解公式K={[ΔSN]TXEXP+1}3,得到拉格朗日函数矩阵;
式中,K表示拉格朗日函数矩阵,ΔSN表示所述信号差分序列,XEXP表示所述期望序列。
可选的,所述利用所述拉格朗日函数矩阵确定拉格朗日因子矢量,具体包括:
利用所述拉格朗日函数矩阵求解公式
Figure GDA0003158678640000032
得到拉格朗日因子矢量;
上式的约束条件为
Figure GDA0003158678640000033
式中,aOPT表示拉格朗日因子矢量;a表示中间矢量,a=[α12,...,αj,...,αN],αj表示中间矢量的第j个元素,N表示中间矢量a的元素数量;K表示拉格朗日函数矩阵;λMAX表示所述拉格朗日函数矩阵的最大特征值。
可选的,所述利用所述拉格朗日函数矩阵和所述拉格朗日因子矢量确定拉格朗日最佳矩阵,具体包括:
利用所述拉格朗日函数矩阵和所述拉格朗日因子矢量求解公式
Figure GDA0003158678640000041
得到拉格朗日最佳矩阵;
式中,WOPT表示拉格朗日最佳矩阵;Wn表示所述拉格朗日因子矢量aOPT的第n个元素;Zn表示所述拉格朗日函数矩阵的第n个特征矢量,N1表示所述拉格朗日函数矩阵的特征矢量的数量,所述拉格朗日因子矢量的元素的数量与所述拉格朗日函数矩阵的特征矢量的数量相等。
可选的,所述利用所述拉格朗日最佳矩阵对所述电力线通信信号序列进行滤波,得到滤波后的信号序列,具体包括:
利用所述拉格朗日最佳矩阵求解公式SNEW=WOPTS,对所述电力线通信信号序列进行滤波,得到滤波后的信号序列;
式中,SNEW表示滤波后的信号序列,WOPT表示所述拉格朗日最佳矩阵,S表示所述电力线通信信号序列。
一种电力线通信信号滤波系统,包括:
电力线通信信号序列获取模块,用于获取实时的电力线通信信号序列;
信号差分序列确定模块,用于利用所述电力线通信信号序列中相邻两个元素的差值确定信号差分序列;
期望序列确定模块,用于利用所述电力线通信信号序列和所述信号差分序列确定期望序列;
拉格朗日函数矩阵确定模块,用于利用所述信号差分序列和所述期望序列确定拉格朗日函数矩阵;
拉格朗日因子矢量确定模块,用于利用所述拉格朗日函数矩阵确定拉格朗日因子矢量;
最佳矩阵确定模块,用于利用所述拉格朗日函数矩阵和所述拉格朗日因子矢量确定拉格朗日最佳矩阵;
滤波模块,用于利用所述拉格朗日最佳矩阵对所述电力线通信信号序列进行滤波,得到滤波后的信号序列。
可选的,所述拉格朗日因子矢量确定模块,具体包括:
拉格朗日因子矢量确定单元,用于利用所述拉格朗日函数矩阵求解公式
Figure GDA0003158678640000051
得到拉格朗日因子矢量;
上式的约束条件为
Figure GDA0003158678640000052
式中,aOPT表示拉格朗日因子矢量;a表示中间矢量,a=[α12,...,αj,...,αN],αj表示中间矢量的第j个元素,N表示中间矢量a的元素数量;K表示拉格朗日函数矩阵;λMAX表示所述拉格朗日函数矩阵的最大特征值。
可选的,所述滤波模块,具体包括:
滤波单元,用于利用所述拉格朗日最佳矩阵求解公式SNEW=WOPTS,对所述电力线通信信号序列进行滤波,得到滤波后的信号序列;
式中,SNEW表示滤波后的信号序列,WOPT表示所述拉格朗日最佳矩阵,S表示所述电力线通信信号序列。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明提供了一种电力线通信信号滤波方法及系统。该方法包括:获取实时的电力线通信信号序列;利用电力线通信信号序列中相邻两个元素的差值确定信号差分序列;利用电力线通信信号序列和信号差分序列确定期望序列;利用信号差分序列和期望序列确定拉格朗日函数矩阵;利用拉格朗日函数矩阵确定拉格朗日因子矢量;利用拉格朗日函数矩阵和拉格朗日因子矢量确定拉格朗日最佳矩阵;利用拉格朗日最佳矩阵对电力线通信信号序列进行滤波,得到滤波后的信号序列。本发明利用拉格朗日函数矩阵确定拉格朗日因子矢量,通过引入拉格朗日因子矢量,利用拉格朗日函数矩阵和拉格朗日因子矢量确定拉格朗日最佳矩阵,利用拉格朗日最佳矩阵对电力线通信信号序列进行滤波,可以有效地滤除PLC信号中的脉冲噪声,保证电力线通信系统质量。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例所提供的电力线通信信号滤波方法的流程图;
图2为本发明实施例所提供的电力线通信信号滤波系统的结构图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种电力线通信信号滤波方法及系统,有效地滤除PLC信号中的脉冲噪声,保证电力线通信系统质量。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
本实施例提供一种电力线通信信号滤波方法,图1为本发明实施例所提供的电力线通信信号滤波方法的流程图,参见图1,电力线通信信号滤波方法包括:
步骤101,获取实时的电力线通信信号序列(PLC信号序列)。输入实测的PLC信号序列S=[s1,s2,…,sN],N为PLC信号序列的长度。
步骤102,利用电力线通信信号序列中相邻两个元素的差值确定信号差分序列。
步骤102具体包括:
计算电力线通信信号序列中相邻两个元素的差值。
通过差值确定信号差分序列ΔSN,信号差分序列为:ΔSN=[0,s2-s1,s3-s2,…,sN-sN-1]。
式中,电力线通信信号序列为S=[s1,s2,...,sn,...,sN],N表示电力线通信信号序列的长度;sn表示电力线通信信号序列的第n个元素,n=1,2,...,N。
步骤103,利用电力线通信信号序列和信号差分序列确定期望序列。
步骤103具体包括:
利用电力线通信信号序列和信号差分序列求解公式(1),得到期望序列。
Figure GDA0003158678640000071
式中,XEXP表示期望序列,ΔSN表示信号差分序列,σS表示电力线通信信号序列的均方差;σΔS为信号差分序列的均方差,S表示电力线通信信号序列。
步骤104,利用信号差分序列和期望序列确定拉格朗日函数矩阵。
步骤104具体包括:
利用信号差分序列和期望序列求解公式(2),得到拉格朗日函数矩阵。
K={[ΔSN]TXEXP+1}3 (2)
式中,K表示拉格朗日函数矩阵,ΔSN表示信号差分序列,XEXP表示期望序列。
步骤105,利用拉格朗日函数矩阵确定拉格朗日因子矢量。
步骤105具体包括:
利用拉格朗日函数矩阵求解公式(3),得到拉格朗日因子矢量aOPT
Figure GDA0003158678640000072
subject to
Figure GDA0003158678640000073
Figure GDA0003158678640000074
式中,aOPT表示拉格朗日因子矢量;a表示中间矢量,a=[α12,...,αj,...,αN],αj表示中间矢量的第j个元素,N表示中间矢量a的元素数量,j表示αj的序号,j=1,2,...,N;K表示拉格朗日函数矩阵;λMAX表示拉格朗日函数矩阵的最大特征值;subjctto表示使得αj满足
Figure GDA0003158678640000075
条件。中间矢量a的元素数量与PLC信号序列的元素数量相等。
通过正交三角分解(QR decomposition)法、奇异值分解(Singular ValueDecomposition,SVD)法和雅可比迭代法等方法求解拉格朗日函数矩阵的特征值,然后比较所有特征值确定最大特征值λMAX
步骤106,利用拉格朗日函数矩阵和拉格朗日因子矢量确定拉格朗日最佳矩阵。
步骤106具体包括:
利用拉格朗日函数矩阵和拉格朗日因子矢量求解公式(4),得到拉格朗日最佳矩阵(Lagrange最佳矩阵)。
Figure GDA0003158678640000081
式中,WOPT表示Lagrange最佳矩阵;Wn表示拉格朗日因子矢量aOPT的第n个元素;Zn表示拉格朗日函数矩阵的第n个特征矢量,N1表示拉格朗日函数矩阵的特征矢量的数量,拉格朗日因子矢量的元素的数量与拉格朗日函数矩阵的特征矢量的数量相等,拉格朗日函数矩阵的特征矢量的数量与PLC信号序列的元素数量相等。Zn为拉格朗日函数矩阵的特征矢量,可以通过对拉格朗日函数矩阵进行SVD分解得到,Zn是SVD分解得到的左酋矩阵中的第n列元素,是一个矢量。Lagrange因子矢量aOPT为一个N维矢量,里面有N个元素,第n个即为Wn。因为步骤105求解得到了拉格朗日因子矢量的最小值,即最佳拉格朗日因子矢量,所以将拉格朗日函数矩阵和拉格朗日因子矢量代入公式(4)能够得到Lagrange最佳矩阵。
步骤107,利用拉格朗日最佳矩阵对电力线通信信号序列进行滤波,得到滤波后的信号序列。
步骤107具体包括:
利用拉格朗日最佳矩阵求解公式(5),对电力线通信信号序列进行滤波,得到滤波后的信号序列。
SNEW=WOPTS (5)
式中,SNEW表示滤波后的信号序列,WOPT表示Lagrange最佳矩阵,S表示电力线通信信号序列。
步骤102-106的目的是使滤除噪声前后的信号之间的功率差值最小化,进而得到最佳逼近矩阵WOPT,通过将WOPT与S相乘可以滤除电力线通信信号序列中的噪声。
本实施例还提供一种电力线通信信号滤波系统,图2为本发明实施例所提供的电力线通信信号滤波系统的结构图。参见图2,电力线通信信号滤波系统包括:
电力线通信信号序列获取模块201,用于获取实时的电力线通信信号序列(PLC信号序列)。输入实测的PLC信号序列S=[s1,s2,…,sN],N为PLC信号序列的长度。
信号差分序列确定模块202,用于利用电力线通信信号序列中相邻两个元素的差值确定信号差分序列。
信号差分序列确定模块202具体包括:
差值计算单元,用于计算电力线通信信号序列中相邻两个元素的差值。
信号差分序列确定单元,用于通过差值确定信号差分序列ΔSN,信号差分序列为:ΔSN=[0,s2-s1,s3-s2,…,sN-sN-1]。
式中,电力线通信信号序列为S=[s1,s2,...,sn,...,sN],N表示电力线通信信号序列的长度;sn表示电力线通信信号序列的第n个元素,n=1,2,...,N。
期望序列确定模块203,用于利用电力线通信信号序列和信号差分序列确定期望序列。
期望序列确定模块203具体包括:
期望序列确定单元,用于利用电力线通信信号序列和信号差分序列求解公式(1),得到期望序列。
Figure GDA0003158678640000091
式中,XEXP表示期望序列,ΔSN表示信号差分序列,σS表示电力线通信信号序列的均方差;σΔS为信号差分序列的均方差,S表示电力线通信信号序列。
拉格朗日函数矩阵确定模块204,用于利用信号差分序列和期望序列确定拉格朗日函数矩阵。
拉格朗日函数矩阵确定模块204具体包括:
拉格朗日函数矩阵确定单元,用于利用信号差分序列和期望序列求解公式(2),得到拉格朗日函数矩阵。
K={[ΔSN]TXEXP+1}3 (2)
式中,K表示拉格朗日函数矩阵,ΔSN表示信号差分序列,XEXP表示期望序列。
拉格朗日因子矢量确定模块205,用于利用拉格朗日函数矩阵确定拉格朗日因子矢量。
拉格朗日因子矢量确定模块205具体包括:
拉格朗日因子矢量确定单元,用于利用拉格朗日函数矩阵求解公式(3),得到拉格朗日因子矢量aOPT
Figure GDA0003158678640000101
subject to
Figure GDA0003158678640000102
Figure GDA0003158678640000103
式中,aOPT表示拉格朗日因子矢量;a表示中间矢量,a=[α12,...,αj,...,αN],αj表示中间矢量的第j个元素,N表示中间矢量a的元素数量,j表示αj的序号,j=1,2,...,N;K表示拉格朗日函数矩阵;λMAX表示拉格朗日函数矩阵的最大特征值;subjctto表示使得αj满足
Figure GDA0003158678640000104
条件。中间矢量a的元素数量与PLC信号序列的元素数量相等。
最佳矩阵确定模块206,用于利用拉格朗日函数矩阵和拉格朗日因子矢量确定拉格朗日最佳矩阵。
最佳矩阵确定模块206具体包括:
最佳矩阵确定单元,用于利用拉格朗日函数矩阵和拉格朗日因子矢量求解公式(4),得到拉格朗日最佳矩阵。
Figure GDA0003158678640000111
式中,WOPT表示Lagrange最佳矩阵;Wn表示拉格朗日因子矢量aOPT的第n个元素;Zn表示拉格朗日函数矩阵的第n个特征矢量,N1表示拉格朗日函数矩阵的特征矢量的数量,拉格朗日因子矢量的元素的数量与拉格朗日函数矩阵的特征矢量的数量相等,拉格朗日函数矩阵的特征矢量的数量与PLC信号序列的元素数量相等。
滤波模块207,用于利用拉格朗日最佳矩阵对电力线通信信号序列进行滤波,得到滤波后的信号序列。
滤波模块207具体包括:
滤波单元,用于利用拉格朗日最佳矩阵求解公式(5),对电力线通信信号序列进行滤波,得到滤波后的信号序列。
SNEW=WOPTS (5)
式中,SNEW表示滤波后的信号序列,WOPT表示Lagrange最佳矩阵,S表示电力线通信信号序列。
本发明利用拉格朗日函数矩阵确定拉格朗日因子矢量,通过引入拉格朗日因子矢量,利用拉格朗日函数矩阵和拉格朗日因子矢量确定拉格朗日最佳矩阵,利用拉格朗日最佳矩阵对电力线通信信号序列进行滤波,可以有效地滤除PLC信号中的脉冲噪声,保证电力线通信系统质量;同时对有用信号具有高保真作用,不会造成有用信号的失真,所提出的方法较为简单,计算速度较快;滤除脉冲噪声之后,PLC信号的信噪比可以提高3dB左右。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。

Claims (3)

1.一种电力线通信信号滤波方法,其特征在于,包括:
获取实时的电力线通信信号序列;
利用所述电力线通信信号序列中相邻两个元素的差值确定信号差分序列;
利用所述电力线通信信号序列和所述信号差分序列确定期望序列,具体包括:利用所述电力线通信信号序列和所述信号差分序列求解公式
Figure FDA0003158678630000011
得到期望序列;
利用所述信号差分序列和所述期望序列确定拉格朗日函数矩阵;具体包括:利用所述信号差分序列和所述期望序列求解公式K={[ΔSN]TXEXP+1}3,得到拉格朗日函数矩阵;
利用所述拉格朗日函数矩阵确定拉格朗日因子矢量,具体包括:利用所述拉格朗日函数矩阵求解公式
Figure FDA0003158678630000012
得到拉格朗日因子矢量;所述拉格朗日函数矩阵求解公式的约束条件为
Figure FDA0003158678630000013
利用所述拉格朗日函数矩阵和所述拉格朗日因子矢量确定拉格朗日最佳矩阵,具体包括:
利用所述拉格朗日函数矩阵和所述拉格朗日因子矢量求解公式
Figure FDA0003158678630000014
得到拉格朗日最佳矩阵;
利用所述拉格朗日最佳矩阵对所述电力线通信信号序列进行滤波,得到滤波后的信号序列,具体包括:利用所述拉格朗日最佳矩阵求解公式SNEW=WOPTS,对所述电力线通信信号序列进行滤波,得到滤波后的信号序列;
式中,XEXP表示期望序列,ΔSN表示所述信号差分序列,σS表示所述电力线通信信号序列的均方差;σΔS为所述信号差分序列的均方差,S表示所述电力线通信信号序列,K表示拉格朗日函数矩阵,ΔSN表示所述信号差分序列,XEXP表示所述期望序列,aOPT表示拉格朗日因子矢量;a表示中间矢量,a=[α12,...,αj,...,αN],αj表示中间矢量的第j个元素,N表示中间矢量a的元素数量;K表示拉格朗日函数矩阵;λMAX表示所述拉格朗日函数矩阵的最大特征值,WOPT表示拉格朗日最佳矩阵;Wn表示所述拉格朗日因子矢量aOPT的第n个元素;Zn表示所述拉格朗日函数矩阵的第n个特征矢量,N1表示所述拉格朗日函数矩阵的特征矢量的数量,所述拉格朗日因子矢量的元素的数量与所述拉格朗日函数矩阵的特征矢量的数量相等,SNEW表示滤波后的信号序列。
2.根据权利要求1所述的电力线通信信号滤波方法,其特征在于,所述利用所述电力线通信信号序列中相邻两个元素的差值确定信号差分序列,具体包括:
计算所述电力线通信信号序列中相邻两个元素的差值;
通过所述差值确定信号差分序列,所述信号差分序列为:ΔSN=[0,s2-s1,...,sn-sn-1,...,sN-sN-1];
式中,所述电力线通信信号序列为S=[s1,s2,...,sn,...,sN],N表示所述电力线通信信号序列的长度;sn表示所述电力线通信信号序列的第n个元素,n=1,2,...,N。
3.一种电力线通信信号滤波系统,其特征在于,包括:
电力线通信信号序列获取模块,用于获取实时的电力线通信信号序列;
信号差分序列确定模块,用于利用所述电力线通信信号序列中相邻两个元素的差值确定信号差分序列;
期望序列确定模块,用于利用所述电力线通信信号序列和所述信号差分序列确定期望序列;
期望序列确定模块具体包括:
期望序列确定单元,用于利用电力线通信信号序列和信号差分序列求解公式
Figure FDA0003158678630000021
得到期望序列;
式中,XEXP表示期望序列,ΔSN表示信号差分序列,σS表示电力线通信信号序列的均方差;σΔS为信号差分序列的均方差,S表示电力线通信信号序列;
拉格朗日函数矩阵确定模块,用于利用所述信号差分序列和所述期望序列确定拉格朗日函数矩阵;
拉格朗日函数矩阵确定模块具体包括:
拉格朗日函数矩阵确定单元,用于利用信号差分序列和期望序列求解公式K={[ΔSN]TXEXP+1}3,得到拉格朗日函数矩阵;
式中,K表示拉格朗日函数矩阵,ΔSN表示信号差分序列,XEXP表示期望序列;
拉格朗日因子矢量确定模块,用于利用所述拉格朗日函数矩阵确定拉格朗日因子矢量;
所述拉格朗日因子矢量确定模块,具体包括:
拉格朗日因子矢量确定单元,用于利用所述拉格朗日函数矩阵求解公式
Figure FDA0003158678630000031
得到拉格朗日因子矢量;
上式的约束条件为
Figure FDA0003158678630000032
式中,aOPT表示拉格朗日因子矢量;a表示中间矢量,a=[α12,...,αj,...,αN],αj表示中间矢量的第j个元素,N表示中间矢量a的元素数量;K表示拉格朗日函数矩阵;λMAX表示所述拉格朗日函数矩阵的最大特征值;
最佳矩阵确定模块,用于利用所述拉格朗日函数矩阵和所述拉格朗日因子矢量确定拉格朗日最佳矩阵;
最佳矩阵确定模块具体包括:
最佳矩阵确定单元,用于利用拉格朗日函数矩阵和拉格朗日因子矢量求解公式
Figure FDA0003158678630000033
得到拉格朗日最佳矩阵;
式中,WOPT表示Lagrange最佳矩阵;Wn表示拉格朗日因子矢量aOPT的第n个元素;Zn表示拉格朗日函数矩阵的第n个特征矢量,N1表示拉格朗日函数矩阵的特征矢量的数量,拉格朗日因子矢量的元素的数量与拉格朗日函数矩阵的特征矢量的数量相等,拉格朗日函数矩阵的特征矢量的数量与PLC信号序列的元素数量相等;
滤波模块,用于利用所述拉格朗日最佳矩阵对所述电力线通信信号序列进行滤波,得到滤波后的信号序列;
所述滤波模块,具体包括:
滤波单元,用于利用所述拉格朗日最佳矩阵求解公式SNEW=WOPTS,对所述电力线通信信号序列进行滤波,得到滤波后的信号序列;
式中,SNEW表示滤波后的信号序列。
CN202010971971.6A 2019-11-02 2020-09-16 一种电力线通信信号滤波方法及系统 Active CN111934716B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201911062084.0A CN110784246A (zh) 2019-11-02 2019-11-02 一种利用Lagrande因子的PLC信号滤波方法和系统
CN2019110620840 2019-11-02

Publications (2)

Publication Number Publication Date
CN111934716A CN111934716A (zh) 2020-11-13
CN111934716B true CN111934716B (zh) 2021-09-17

Family

ID=69388539

Family Applications (2)

Application Number Title Priority Date Filing Date
CN201911062084.0A Withdrawn CN110784246A (zh) 2019-11-02 2019-11-02 一种利用Lagrande因子的PLC信号滤波方法和系统
CN202010971971.6A Active CN111934716B (zh) 2019-11-02 2020-09-16 一种电力线通信信号滤波方法及系统

Family Applications Before (1)

Application Number Title Priority Date Filing Date
CN201911062084.0A Withdrawn CN110784246A (zh) 2019-11-02 2019-11-02 一种利用Lagrande因子的PLC信号滤波方法和系统

Country Status (1)

Country Link
CN (2) CN110784246A (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111082837B (zh) * 2019-11-25 2021-04-23 宁波大学 基于递归压缩感知的电力线通信系统脉冲噪声抑制方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101140467A (zh) * 2000-12-22 2008-03-12 美国标准国际公司 制造和测试通信系统
CN107729845A (zh) * 2017-10-20 2018-02-23 开沃新能源汽车集团有限公司 一种基于子空间特征值分解的实测频响函数降噪方法
CN108551357A (zh) * 2018-04-19 2018-09-18 南方电网科学研究院有限责任公司 一种电力线通信系统的编码解码方法及装置
CN110017894A (zh) * 2019-05-14 2019-07-16 广东石油化工学院 变压器运行状态振声检测中随机噪声的滤除方法和装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9537697B2 (en) * 2014-08-12 2017-01-03 Bae Systems Information And Electronic System Integration Inc. Constrained optimization approach to compander design for OFDM PAPR reduction
CN108918932B (zh) * 2018-09-11 2021-01-15 广东石油化工学院 负荷分解中功率信号自适应滤波方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101140467A (zh) * 2000-12-22 2008-03-12 美国标准国际公司 制造和测试通信系统
CN107729845A (zh) * 2017-10-20 2018-02-23 开沃新能源汽车集团有限公司 一种基于子空间特征值分解的实测频响函数降噪方法
CN108551357A (zh) * 2018-04-19 2018-09-18 南方电网科学研究院有限责任公司 一种电力线通信系统的编码解码方法及装置
CN110017894A (zh) * 2019-05-14 2019-07-16 广东石油化工学院 变压器运行状态振声检测中随机噪声的滤除方法和装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于EMD_TFPF算法的电力线通信噪声消除技术研究;翟明岳;《电力系统保护与控制》;20150401;51-56 *

Also Published As

Publication number Publication date
CN111934716A (zh) 2020-11-13
CN110784246A (zh) 2020-02-11

Similar Documents

Publication Publication Date Title
CN108880621B (zh) 一种电力线通信信号自适应滤波方法
CN109257068B (zh) 一种电力线通信信号自适应滤波方法
CN111934716B (zh) 一种电力线通信信号滤波方法及系统
CN111756405A (zh) 一种利用f平方模的plc信道脉冲噪声检测方法和系统
CN112118027A (zh) 一种plc信道脉冲噪声检测方法及系统
CN111641435A (zh) 一种利用Fenchel共轭的PLC信号滤波方法和系统
CN111628804A (zh) 一种利用Gilbert优化的PLC信号滤波方法和系统
CN111756456A (zh) 一种利用偏移量的plc信道脉冲噪声检测方法和系统
CN110336591B (zh) 一种利用信号分离的plc信号滤波方法和系统
CN109039379B (zh) 电力线通信信号滤波方法
CN109117807B (zh) 一种plc通信信号自适应时频峰值滤波方法及系统
CN110635824B (zh) 一种利用分类回归树的plc信道脉冲噪声检测方法和系统
CN111641434A (zh) 一种利用完全矢量的plc信号滤波方法和系统
CN111800165A (zh) 一种利用奇异值矩阵的plc信号滤波方法和系统
CN111641436A (zh) 一种利用lp优化的plc信号滤波方法和系统
CN110739986B (zh) 一种利用投影累积量的plc信道脉冲噪声检测方法和系统
CN111541635A (zh) 一种利用t分布的PLC信号滤波方法和系统
CN112165342B (zh) 一种利用模式特征矢量的噪声检测方法和系统
CN110572189A (zh) 一种利用Fermal理论的PLC信号滤波方法和系统
CN112104392B (zh) 一种利用状态矩阵的plc信道脉冲噪声检测方法和系统
CN111610747A (zh) 一种利用残差相关的plc信号重构方法和系统
CN111669201A (zh) 一种利用Dantzig优化的PLC信号滤波方法和系统
CN112187319A (zh) 一种利用Dantzig选择器的PLC信号重构方法和系统
CN111786700A (zh) 一种利用克氏能量尺度的plc信号滤波方法和系统
CN112350747B (zh) 一种利用状态估计的plc信道脉冲噪声检测方法和系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant