一种利用奇异值矩阵的plc信号滤波方法和系统
阅读说明：本技术 一种利用奇异值矩阵的plc信号滤波方法和系统 (PLC signal filtering method and system using singular value matrix ) 是由 翟明岳 于 2020-07-08 设计创作，主要内容包括：本发明的实施例公开一种利用奇异值矩阵的PLC信号滤波方法和系统,所述方法包括：步骤101获取按时间顺序采集的信号序列S；步骤102求取信号延迟矩阵；步骤103对信号延迟矩阵进行奇异值分解；步骤104求取信号阈值；步骤105求取有效特征值个数；步骤106求取恢复矩阵；步骤107求取滤除噪声后的信号序列。(The embodiment of the invention discloses a PLC signal filtering method and a system by using a singular value matrix, wherein the method comprises the following steps: step 101, acquiring a signal sequence S acquired according to a time sequence; step 102, obtaining a signal delay matrix; 103, performing singular value decomposition on the signal delay matrix; step 104, calculating a signal threshold; step 105, calculating the number of effective characteristic values; step 106, obtaining a recovery matrix; step 107 is to find the signal sequence after noise filtering.)
技术领域
本发明涉及通信领域，尤其涉及一种PLC信号滤波方法和系统。
背景技术
电力线通信，相比各种有线通信技术，无需重新布线，易于组网等优点，具有广阔的应用前景。电力线通信技术分为窄带电力线通信(Narrowband over power line，NPL)和宽带电力线通信(Broadband over power line，BPL)；窄带电力线通信是指带宽限定在3k至500kHz的电力线载波通信技术；电力线通信技术包括欧洲CENELEC的规定带宽(3148.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为核心的多种扩频通信技术。
虽然电力线通信系统有着广泛的应用，且技术相对成熟，但是电力线通信系统中大量的分支和电气设备，会在电力线信道中产生大量的噪声；而其中随机脉冲噪声具有很大的随机性，噪声强度高，对电力线通信系统造成严重破坏，因此，针对随机脉冲噪声的抑制技术，一直是国内外学者研究的重点；而且噪声模型并不符合高斯分布。因此，传统的针对高斯噪声设计的通信系统不再适用于电力线载波通信系统，必须研究相应的噪声抑制技术，以提高电力线通信系统信噪比，降低误码率，保证电力线通信系统质量。
在实际应用中，一些简单的非线性技术经常被应用于消除电力线信道噪声，如Clip-ping、Blanking和Clipping/Blanking技术，但是这些研究方法都必须在一定的信噪比情况下才能良好工作，仅仅考虑了冲击噪声的消除，在电力线通信系统中，某些商用电力线发送器的特征是低发射功率，在一些特殊情况，发射功率甚至可能会低于18w，因此，在某些特殊情况，信号将会淹没在大量噪声中，导致电力线通信系统低信噪比情况。
发明内容
随着非线性电器的应用和普及，中低压输配电网络中背景噪声呈现出较为明显的非平稳性和非高斯特性，常用的低通滤波器在非平稳和非高斯噪声环境中难以达到理想的滤波效果，很难滤除非平稳非高斯噪声，严重影响了PLC通信系统的性能。
本发明的目的是提供一种利用奇异值矩阵的PLC信号滤波方法和系统，所提出的方法利用了PLC调制信号、脉冲噪声和背景噪声在信号混合表示领域中的差异，通过奇异值矩阵性质区分PLC调制信号、脉冲噪声和背景噪声。所提出的方法具有较好的噪声滤除性能，计算也非常简单。
为实现上述目的，本发明提供了如下方案：
一种利用奇异值矩阵的PLC信号滤波方法，包括：
步骤101获取按时间顺序采集的信号序列S；
步骤102求取信号延迟矩阵，具体为：所述信号延迟矩阵D的第i行第j列元素为di,j，其计算公式为di,j＝s|i+j+κ|N；其中，s|i+j+κ|N为所述信号序列S的第|i+j+κ|N个元素；| |N表示以N为模取余数运算；κ为延迟因子，其计算公式为
步骤103对信号延迟矩阵进行奇异值分解，具体为：所述信号延迟矩阵D的奇异值分解为D＝PVQT；其中，P为所述信号延迟矩阵D的左奇异矢量矩阵；Q为所述信号延迟矩阵D的右奇异矢量矩阵；V为所述信号延迟矩阵D的奇异值矩阵；
步骤104求取信号阈值，所述信号阈值ε0的求取公式为其中，m0为所述信号序列S的均值；σ0为所述信号序列S的均方差；λmax为平均相关矩阵B的最大特征值；λmin为所述平均相关矩阵B的非零最
小特征值；所述平均相关矩阵B的计算公式为B＝[S-m0]T[S-m0]；
步骤105求取有效特征值个数，所述有效特征值个数M的求取过程为：所述信号延迟矩阵D的特征值大于或者等于所述信号阈值ε0的个数即为所述有效特征值M的值；
步骤106求取恢复矩阵，具体为：所述恢复矩阵ω的求取公式为：其中，vm为所述信号延迟矩阵D的第m个奇异值；pm为所述左奇异值矩阵P的第m列矢量；qm为所述右奇异值矩阵Q的第m列矢量；
步骤107求取滤除噪声后的信号序列，具体为：所述滤除噪声后的信号序列Snew的计算公式为Snew＝ωS。
一种利用奇异值矩阵的PLC信号滤波系统，包括：
模块201获取按时间顺序采集的信号序列S；
模块202求取信号延迟矩阵，具体为：所述信号延迟矩阵D的第i行第j列元素为di,j，其计算公式为di,j＝s|i+j+κ|N；其中，s|i+j+κ|N为所述信号序列S的第|i+j+κ|N个元素；| |N表示以N为模取余数运算；κ为延迟因子，其计算公式为snr为所述信号序列S的信噪比；表示下取整运算；
模块203对信号延迟矩阵进行奇异值分解，具体为：所述信号延迟矩阵D的奇异值分解为D＝PVQT；其中，P为所述信号延迟矩阵D的左奇异矢量矩阵；Q为所述信号延迟矩阵D的右奇异矢量矩阵；V为所述信号延迟矩阵D的奇异值矩阵；
模块204求取信号阈值，所述信号阈值ε0的求取公式为
模块205求取有效特征值个数，所述有效特征值个数M的求取过程为：所述信号延迟矩阵D的特征值大于或者等于所述信号阈值ε0的个数即为所述有效特征值M的值；
模块206求取恢复矩阵，具体为：所述恢复矩阵ω的求取公式为：其中，vm为所述信号延迟矩阵D的第m个奇异值；pm为所述左奇异值矩阵P的第m列矢量；qm为所述右奇异值矩阵Q的第m列矢量；
模块207求取滤除噪声后的信号序列，具体为：所述滤除噪声后的信号序列Snew的计算公式为Snew＝ωS。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
随着非线性电器的应用和普及，中低压输配电网络中背景噪声呈现出较为明显的非平稳性和非高斯特性，常用的低通滤波器在非平稳和非高斯噪声环境中难以达到理想的滤波效果，很难滤除非平稳非高斯噪声，严重影响了PLC通信系统的性能。
本发明的目的是提供一种利用奇异值矩阵的PLC信号滤波方法和系统，所提出的方法利用了PLC调制信号、脉冲噪声和背景噪声在信号混合表示领域中的差异，通过奇异值矩阵性质区分PLC调制信号、脉冲噪声和背景噪声。所提出的方法具有较好的噪声滤除性能，计算也非常简单。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍。显而易见，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明的方法流程示意图；
图2为本发明的系统流程示意图；
图3为本发明的具体实施案例流程示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1一种利用奇异值矩阵的PLC信号滤波方法的流程示意图
图1为本发明一种利用奇异值矩阵的PLC信号滤波方法的流程示意图。如图1所示，所述的一种利用奇异值矩阵的PLC信号滤波方法具体包括以下步骤：
步骤101获取按时间顺序采集的信号序列S；
步骤102求取信号延迟矩阵，具体为：所述信号延迟矩阵D的第i行第j列元素为di,j，其计算公式为di,j＝s|i+j+κ|N；其中，s|i+j+κ|N为所述信号序列S的第|i+j+κ|N个元素；| |N表示以N为模取余数运算；κ为延迟因子，其计算公式为snr为所述信号序列S的信噪比；表示下取整运算；
步骤103对信号延迟矩阵进行奇异值分解，具体为：所述信号延迟矩阵D的奇异值分解为D＝PVQT；其中，P为所述信号延迟矩阵D的左奇异矢量矩阵；Q为所述信号延迟矩阵D的右奇异矢量矩阵；V为所述信号延迟矩阵D的奇异值矩阵；
步骤104求取信号阈值，所述信号阈值ε0的求取公式为
步骤105求取有效特征值个数，所述有效特征值个数M的求取过程为：所述信号延迟矩阵D的特征值大于或者等于所述信号阈值ε0的个数即为所述有效特征值M的值；
步骤106求取恢复矩阵，具体为：所述恢复矩阵ω的求取公式为：
步骤107求取滤除噪声后的信号序列，具体为：所述滤除噪声后的信号序列Snew的计算公式为Snew＝ωS。
图2一种利用奇异值矩阵的PLC信号滤波系统的结构意图
图2为本发明一种利用奇异值矩阵的PLC信号滤波系统的结构示意图。如图2所示，所述一种利用奇异值矩阵的PLC信号滤波系统包括以下结构：
模块201获取按时间顺序采集的信号序列S；
模块202求取信号延迟矩阵，具体为：所述信号延迟矩阵D的第i行第j列元素为di,j，其计算公式为di,j＝s|i+j+κ|N；其中，s|i+j+κ|N为所述信号序列S的第|i+j+κ|N个元素；| |N表示以N为模取余数运算；κ为延迟因子，其计算公式为snr为所述信号序列S的信噪比；表示下取整运算；
模块203对信号延迟矩阵进行奇异值分解，具体为：所述信号延迟矩阵D的奇异值分解为D＝PVQT；其中，P为所述信号延迟矩阵D的左奇异矢量矩阵；Q为所述信号延迟矩阵D的右奇异矢量矩阵；V为所述信号延迟矩阵D的奇异值矩阵；
模块204求取信号阈值，所述信号阈值ε0的求取公式为其中，m0为所述信号序列S的均值；σ0为所述信号序列S的均方差；λmax为平均相关矩阵B的最大特征值；λmin为所述平均相关矩阵B的非零最小特征值；所述平均相关矩阵B的计算公式为B＝[S-m0]T[S-m0]；
模块205求取有效特征值个数，所述有效特征值个数M的求取过程为：所述信号延迟矩阵D的特征值大于或者等于所述信号阈值ε0的个数即为所述有效特征值M的值；
模块206求取恢复矩阵，具体为：所述恢复矩阵ω的求取公式为：
模块207求取滤除噪声后的信号序列，具体为：所述滤除噪声后的信号序列Snew的计算公式为Snew＝ωS。
下面提供一个具体实施案例，进一步说明本发明的方案
图3为本发明具体实施案例的流程示意图。如图3所示，具体包括以下步骤：
步骤301获取按时间顺序采集的信号序列S；
步骤302求取信号延迟矩阵，具体为：所述信号延迟矩阵D的第i行第j列元素为di,j，其计算公式为di,j＝s|i+j+κ|N；其中，s|i+j+κ|N为所述信号序列S的第|i+j+κ|N个元素；| |N表示以N为模取余数运算；κ为延迟因子，其计算公式为snr为所述信号序列S的信噪比；表示下取整运算；
步骤303对信号延迟矩阵进行奇异值分解，具体为：所述信号延迟矩阵D的奇异值分解为D＝PVQT；其中，P为所述信号延迟矩阵D的左奇异矢量矩阵；Q为所述信号延迟矩阵D的右奇异矢量矩阵；V为所述信号延迟矩阵D的奇异值矩阵；
步骤304求取信号阈值，所述信号阈值ε0的求取公式为
步骤305求取有效特征值个数，所述有效特征值个数M的求取过程为：所述信号延迟矩阵D的特征值大于或者等于所述信号阈值ε0的个数即为所述有效特征值M的值；
步骤306求取恢复矩阵，具体为：所述恢复矩阵ω的求取公式为：
步骤307求取滤除噪声后的信号序列，具体为：所述滤除噪声后的信号序列Snew的计算公式为Snew＝ωS。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述较为简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。
- 上一篇：一种医用注射器针头装配设备
- 下一篇：一种HPLC互联互通测试装置及测试方法