CN112347922B - 一种利用Hankerl矩阵的功率信号滤波方法和*** - Google Patents

一种利用Hankerl矩阵的功率信号滤波方法和*** Download PDF

Info

Publication number
CN112347922B
CN112347922B CN202011229098.XA CN202011229098A CN112347922B CN 112347922 B CN112347922 B CN 112347922B CN 202011229098 A CN202011229098 A CN 202011229098A CN 112347922 B CN112347922 B CN 112347922B
Authority
CN
China
Prior art keywords
matrix
hankerl
signal
signal sequence
specifically
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.)
Expired - Fee Related
Application number
CN202011229098.XA
Other languages
English (en)
Other versions
CN112347922A (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.)
North China Electric Power University
Original Assignee
North China Electric Power University
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 North China Electric Power University filed Critical North China Electric Power University
Priority to CN202011229098.XA priority Critical patent/CN112347922B/zh
Publication of CN112347922A publication Critical patent/CN112347922A/zh
Application granted granted Critical
Publication of CN112347922B publication Critical patent/CN112347922B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/02Preprocessing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/02Preprocessing
    • G06F2218/04Denoising

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Complex Calculations (AREA)

Abstract

本发明的实施例公开一种利用Hankerl矩阵的功率信号滤波方法和***,所述方法包括:步骤101获取按时间顺序采集的信号序列S;步骤102生成Hankel矩阵;步骤103求取最佳秩判断阈值;步骤104求取最佳秩;步骤105求取收缩化特征值;步骤106求取低秩逼近信号矩阵;步骤107求取滤除噪声后的信号序列。

Description

一种利用Hankerl矩阵的功率信号滤波方法和***
技术领域
本发明涉及电力领域,尤其涉及一种功率信号的滤波方法和***。
背景技术
随着智能电网的发展,家庭用电负荷的分析变得越来越重要。通过用电负荷的分析,家庭用户可以及时获得每个电器的用电信息,以及电费的精细化清单;电力部门可以获得更详尽的用户用电信息,并可以提高用电负荷预测的准确度,为电力部门提供统筹规划的依据。同时,利用每个电器的用电信息,可获知用户的用电行为,这对于家庭能耗评估和节能策略的研究具有指导意义。
当前用电负荷分解主要分为侵入式负荷分解和非侵入式负荷分解两种方法。非侵入式负荷分解方法不需要在负荷的内部用电设备上安装监测设备,只需要根据用电负荷总信息即可获得每个用电设备的负荷信息。非侵入式负荷分解方法具有投入少、方便使用等特点,因此,该方法适用于家庭负荷用电的分解。
非侵入式负荷分解算法中,电气设备的开关事件检测是其中最重要的环节。最初的开关事件检测以有功功率P的变化值作为开关事件检测的判断依据,方便且直观。这是因为任何一个用电设备的运行状态发生变化,其所消耗的功率值也必然发生改变,并且该改变也将会在所有电器所消耗的总功率中体现出来。这种方法除了需要设置功率变化值的合理阈值,还需要解决事件检测方法在实际应用中存在的问题,例如某些电器启动时刻的瞬时功率值会出现较大的尖峰(马达启动电流远大于额定电流),会造成电器稳态功率变化值不准确,从而影响对开关事件检测的判断;而且不同家用电器的暂态过程或长或短(脉冲噪声的持续时间和发生频率相差较大),因此功率变化值的确定变得较为困难;由于电能质量的变化(如电压突降)有功功率会出现突变的情况,这样很可能会出现误判。
因此,开关事件检测过程中,所使用的实测功率信号常常受到噪声的影响,利用这些不完善的功率信号是不能正确地进行开关事件检测的。因此如何有效地重构不完整的功率信号,滤除噪声的影响,是此方法能否成功的关键。现在常用的方法,对此问题重视不够,还未采取有效的措施解决此问题。
发明内容
开关事件检测过程中,所使用的实测功率信号常常受到噪声的影响,利用这些不完善的功率信号是不能正确地进行开关事件检测的。因此如何有效地重构不完整的功率信号,滤除噪声的影响,是此方法能否成功的关键。现在常用的方法,对此问题重视不够,还未采取有效的措施解决此问题。
本发明的目的是提供一种利用Hankerl矩阵的功率信号滤波方法和***,所提出的方法利用了功率信号与噪声在发生机制方面的差异,根据Hankerl矩阵性质,实现功率信号的滤波。所提出的方法具有较好的鲁棒性,计算也较为简单。
为实现上述目的,本发明提供了如下方案:
一种利用Hankerl矩阵的功率信号滤波方法,包括:
步骤101获取按时间顺序采集的信号序列S;
步骤102生成Hankel矩阵,具体为:Hankel矩阵记为H,其第i行第j列元素记为hij,所用求取公式为:
Figure BDA0002764572510000025
其中:
i=1,2,···,N为行序号,
j=1,2,···,N为列序号,
Figure BDA0002764572510000026
表示所述信号序列S的第|i-1+j|N个元素,
|i-1+j|N表示以N为模对i-1+j取余数,
N为所述信号序列S的长度;
步骤103求取最佳秩判断阈值,具体为:最佳秩判断阈值记为ε,所用求取公式为:
ε=ln(SNR+1)+2lnσ0
其中:
σ0为所述信号序列S的均方差,
SNR为所述信号序列S的信噪比;
步骤104求取最佳秩,具体为:最佳秩记为R,所用求取公式为:
Figure BDA0002764572510000021
其中:
count{*}表示集合*中元素的个数,
*表示任一集合,
σj(H)为Hankel矩阵H的第j个特征值(特征值按降序排列),
σmed(H)为Hankel矩阵H的中值特征值;
步骤105求取收缩化特征值,具体为:第r个收缩化特征值记为
Figure BDA0002764572510000022
所用求取公式为:
Figure BDA0002764572510000023
Figure BDA0002764572510000024
为幂收缩因子,
r=1,2,···,R为收缩序号,
σR+1(H)为Hankel矩阵H的第R+1个特征值(特征值按降序排列),
σr(H)为Hankel矩阵H的第r个特征值(特征值按降序排列);
步骤106求取低秩逼近信号矩阵,具体为:低秩逼近信号矩阵记为
Figure BDA0002764572510000031
所用求取公式为:
Figure BDA0002764572510000032
其中:
Figure BDA0002764572510000033
为收缩化的特征值矩阵,
U为所述Hankel矩阵H的左特征矢量矩阵,
V为所述Hankel矩阵H的右特征矢量矩阵;
步骤107求取滤除噪声后的信号序列,具体为:滤除噪声后的信号序列为Snew,其第l个元素记为
Figure BDA0002764572510000034
所用求取公式为:
Figure BDA0002764572510000035
其中:
m=1,2,···,N为第一求和参数,
n=1,2,···,N为第二求和参数,
Figure BDA0002764572510000036
表示所述低秩逼近信号矩阵
Figure BDA0002764572510000037
的第m行第n列元素,
l=1,2,···,N为元素序号。
一种利用Hankerl矩阵的功率信号滤波***,包括:
模块201获取按时间顺序采集的信号序列S;
模块202生成Hankel矩阵,具体为:Hankel矩阵记为H,其第i行第j列元素记为hij,所用求取公式为:
Figure BDA0002764572510000038
其中:
i=1,2,···,N为行序号,
j=1,2,···,N为列序号,
Figure BDA0002764572510000039
表示所述信号序列S的第|i-1+j|N个元素,
|i-1+j|N表示以N为模对i-1+j取余数,
N为所述信号序列S的长度;
模块203求取最佳秩判断阈值,具体为:最佳秩判断阈值记为ε,所用求取公式为:
ε=ln(SNR+1)+2lnσ0
其中:
σ0为所述信号序列S的均方差,
SNR为所述信号序列S的信噪比;
模块204求取最佳秩,具体为:最佳秩记为R,所用求取公式为:
Figure BDA0002764572510000041
其中:
count{*}表示集合*中元素的个数,
*表示任一集合,
σj(H)为Hankel矩阵H的第j个特征值(特征值按降序排列),
σmed(H)为Hankel矩阵H的中值特征值;
模块205求取收缩化特征值,具体为:第r个收缩化特征值记为
Figure BDA0002764572510000042
所用求取公式为:
Figure BDA0002764572510000043
Figure BDA0002764572510000044
为幂收缩因子,
r=1,2,···,R为收缩序号,
σR+1(H)为Hankel矩阵H的第R+1个特征值(特征值按降序排列),
σr(H)为Hankel矩阵H的第r个特征值(特征值按降序排列);
模块206求取低秩逼近信号矩阵,具体为:低秩逼近信号矩阵记为
Figure BDA0002764572510000045
所用求取公式为:
Figure BDA0002764572510000046
其中:
Figure BDA0002764572510000047
为收缩化的特征值矩阵,
U为所述Hankel矩阵H的左特征矢量矩阵,
V为所述Hankel矩阵H的右特征矢量矩阵;
模块207求取滤除噪声后的信号序列,具体为:滤除噪声后的信号序列为Snew,其第l个元素记为
Figure BDA0002764572510000048
所用求取公式为:
Figure BDA0002764572510000049
其中:
m=1,2,···,N为第一求和参数,
n=1,2,···,N为第二求和参数,
Figure BDA00027645725100000410
表示所述低秩逼近信号矩阵
Figure BDA00027645725100000411
的第m行第n列元素,
l=1,2,···,N为元素序号。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
开关事件检测过程中,所使用的实测功率信号常常受到噪声的影响,利用这些不完善的功率信号是不能正确地进行开关事件检测的。因此如何有效地重构不完整的功率信号,滤除噪声的影响,是此方法能否成功的关键。现在常用的方法,对此问题重视不够,还未采取有效的措施解决此问题。
本发明的目的是提供一种利用Hankerl矩阵的功率信号滤波方法和***,所提出的方法利用了功率信号与噪声在发生机制方面的差异,根据Hankerl矩阵性质,实现功率信号的滤波。所提出的方法具有较好的鲁棒性,计算也较为简单。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍。显而易见,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明的方法流程示意图;
图2为本发明的***流程示意图;
图3为本发明的具体实施案例流程示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述。显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1一种利用Hankerl矩阵的功率信号滤波方法的流程示意图
图1为本发明一种利用Hankerl矩阵的功率信号滤波方法的流程示意图。如图1所示,所述的一种利用Hankerl矩阵的功率信号滤波方法具体包括以下步骤:
步骤101获取按时间顺序采集的信号序列S;
步骤102生成Hankel矩阵,具体为:Hankel矩阵记为H,其第i行第j列元素记为hij,所用求取公式为:
Figure BDA0002764572510000051
其中:
i=1,2,···,N为行序号,
j=1,2,···,N为列序号,
Figure BDA0002764572510000052
表示所述信号序列S的第|i-1+j|N个元素,
|i-1+j|N表示以N为模对i-1+j取余数,
N为所述信号序列S的长度;
步骤103求取最佳秩判断阈值,具体为:最佳秩判断阈值记为ε,所用求取公式为:
ε=ln(SNR+1)+2lnσ0
其中:
σ0为所述信号序列S的均方差,
SNR为所述信号序列S的信噪比;
步骤104求取最佳秩,具体为:最佳秩记为R,所用求取公式为:
Figure BDA0002764572510000061
其中:
count{*}表示集合*中元素的个数,
*表示任一集合,
σj(H)为Hankel矩阵H的第j个特征值(特征值按降序排列),
σmed(H)为Hankel矩阵H的中值特征值;
步骤105求取收缩化特征值,具体为:第r个收缩化特征值记为
Figure BDA0002764572510000062
所用求取公式为:
Figure BDA0002764572510000063
Figure BDA0002764572510000064
为幂收缩因子,
r=1,2,···,R为收缩序号,
σR+1(H)为Hankel矩阵H的第R+1个特征值(特征值按降序排列),
σr(H)为Hankel矩阵H的第r个特征值(特征值按降序排列);
步骤106求取低秩逼近信号矩阵,具体为:低秩逼近信号矩阵记为
Figure BDA0002764572510000065
所用求取公式为:
Figure BDA0002764572510000066
其中:
Figure BDA0002764572510000067
为收缩化的特征值矩阵,
U为所述Hankel矩阵H的左特征矢量矩阵,
V为所述Hankel矩阵H的右特征矢量矩阵;
步骤107求取滤除噪声后的信号序列,具体为:滤除噪声后的信号序列为Snew,其第l个元素记为
Figure BDA0002764572510000068
所用求取公式为:
Figure BDA0002764572510000069
其中:
m=1,2,···,N为第一求和参数,
n=1,2,···,N为第二求和参数,
Figure BDA0002764572510000071
表示所述低秩逼近信号矩阵
Figure BDA0002764572510000072
的第m行第n列元素,
l=1,2,···,N为元素序号。
图2一种利用Hankerl矩阵的功率信号滤波***的结构意图
图2为本发明一种利用Hankerl矩阵的功率信号滤波***的结构示意图。如图2所示,所述一种利用Hankerl矩阵的功率信号滤波***包括以下结构:
模块201获取按时间顺序采集的信号序列S;
模块202生成Hankel矩阵,具体为:Hankel矩阵记为H,其第i行第j列元素记为hij,所用求取公式为:
Figure BDA0002764572510000074
其中:
i=1,2,···,N为行序号,
j=1,2,···,N为列序号,
Figure BDA0002764572510000075
表示所述信号序列S的第|i-1+j|N个元素,
|i-1+j|N表示以N为模对i-1+j取余数,
N为所述信号序列S的长度;
模块203求取最佳秩判断阈值,具体为:最佳秩判断阈值记为ε,所用求取公式为:
ε=ln(SNR+1)+2lnσ0
其中:
σ0为所述信号序列S的均方差,
SNR为所述信号序列S的信噪比;
模块204求取最佳秩,具体为:最佳秩记为R,所用求取公式为:
Figure BDA0002764572510000073
其中:
count{*}表示集合*中元素的个数,
*表示任一集合,
σj(H)为Hankel矩阵H的第j个特征值(特征值按降序排列),
σmed(H)为Hankel矩阵H的中值特征值;
模块205求取收缩化特征值,具体为:第r个收缩化特征值记为
Figure BDA0002764572510000081
所用求取公式为:
Figure BDA0002764572510000082
Figure BDA0002764572510000083
为幂收缩因子,
r=1,2,···,R为收缩序号,
σR+1(H)为Hankel矩阵H的第R+1个特征值(特征值按降序排列),
σr(H)为Hankel矩阵H的第r个特征值(特征值按降序排列);
模块206求取低秩逼近信号矩阵,具体为:低秩逼近信号矩阵记为
Figure BDA0002764572510000084
所用求取公式为:
Figure BDA0002764572510000085
其中:
Figure BDA0002764572510000086
为收缩化的特征值矩阵,
U为所述Hankel矩阵H的左特征矢量矩阵,
V为所述Hankel矩阵H的右特征矢量矩阵;
模块207求取滤除噪声后的信号序列,具体为:滤除噪声后的信号序列为Snew,其第l个元素记为
Figure BDA0002764572510000087
所用求取公式为:
Figure BDA0002764572510000088
其中:
m=1,2,···,N为第一求和参数,
n=1,2,···,N为第二求和参数,
Figure BDA0002764572510000089
表示所述低秩逼近信号矩阵
Figure BDA00027645725100000810
的第m行第n列元素,
l=1,2,···,N为元素序号。
下面提供一个具体实施案例,进一步说明本发明的方案
图3为本发明具体实施案例的流程示意图。如图3所示,具体包括以下步骤:
步骤301获取按时间顺序采集的信号序列S;
步骤302生成Hankel矩阵,具体为:Hankel矩阵记为H,其第i行第j列元素记为hij,所用求取公式为:
Figure BDA00027645725100000811
其中:
i=1,2,···,N为行序号,
j=1,2,···,N为列序号,
Figure BDA0002764572510000098
表示所述信号序列S的第|i-1+j|N个元素,
|i-1+j|N表示以N为模对i-1+j取余数,
N为所述信号序列S的长度;
步骤303求取最佳秩判断阈值,具体为:最佳秩判断阈值记为ε,所用求取公式为:
ε=ln(SNR+1)+2lnσ0
其中:
σ0为所述信号序列S的均方差,
SNR为所述信号序列S的信噪比;
步骤304求取最佳秩,具体为:最佳秩记为R,所用求取公式为:
Figure BDA0002764572510000091
其中:
count{*}表示集合*中元素的个数,
*表示任一集合,
σj(H)为Hankel矩阵H的第j个特征值(特征值按降序排列),
σmed(H)为Hankel矩阵H的中值特征值;
步骤305求取收缩化特征值,具体为:第r个收缩化特征值记为
Figure BDA0002764572510000092
所用求取公式为:
Figure BDA0002764572510000093
Figure BDA0002764572510000094
为幂收缩因子,
r=1,2,···,R为收缩序号,
σR+1(H)为Hankel矩阵H的第R+1个特征值(特征值按降序排列),
σr(H)为Hankel矩阵H的第r个特征值(特征值按降序排列);
步骤306求取低秩逼近信号矩阵,具体为:低秩逼近信号矩阵记为
Figure BDA0002764572510000095
所用求取公式为:
Figure BDA0002764572510000096
其中:
Figure BDA0002764572510000097
为收缩化的特征值矩阵,
U为所述Hankel矩阵H的左特征矢量矩阵,
V为所述Hankel矩阵H的右特征矢量矩阵;
步骤307求取滤除噪声后的信号序列,具体为:滤除噪声后的信号序列为Snew,其第l个元素记为
Figure BDA0002764572510000101
所用求取公式为:
Figure BDA0002764572510000102
其中:
m=1,2,···,N为第一求和参数,
n=1,2,···,N为第二求和参数,
Figure BDA0002764572510000103
表示所述低秩逼近信号矩阵
Figure BDA0002764572510000104
的第m行第n列元素,
l=1,2,···,N为元素序号。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的***而言,由于其与实施例公开的方法相对应,所以描述较为简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。

Claims (2)

1.一种利用Hankerl矩阵的功率信号滤波方法,其特征在于,包括:
步骤101获取按时间顺序采集的信号序列S;
步骤102生成Hankerl矩阵,具体为:Hankerl矩阵记为H,其第i行第j列元素记为hij,所用求取公式为:
Figure FDA0003353832180000011
其中:
i=1,2,…,N为行序号,
j=1,2,…,N为列序号,
Figure FDA0003353832180000012
表示所述信号序列S的第|i-1+j|N个元素,
|i-1+j|N表示以N为模对i-1+j取余数,
N为所述信号序列S的长度;
步骤103求取最佳秩判断阈值,具体为:最佳秩判断阈值记为ε,所用求取公式为:
ε=ln(SNR+1)+2lnσ0
其中:
σ0为所述信号序列S的均方差,
SNR为所述信号序列S的信噪比;
步骤104求取最佳秩,具体为:最佳秩记为R,所用求取公式为:
Figure FDA0003353832180000013
其中:
count{*}表示集合*中元素的个数,
*表示任一集合,
σj(H)为Hankerl矩阵H的第j个特征值,特征值按降序排列;
σmed(H)为Hankerl矩阵H的中值特征值;
步骤105求取收缩化特征值,具体为:第r个收缩化特征值记为
Figure FDA0003353832180000014
所用求取公式为:
Figure FDA0003353832180000015
Figure FDA0003353832180000016
为幂收缩因子,
r=1,2,…,R为收缩序号,
σR+1(H)为Hankerl矩阵H的第R+1个特征值,特征值按降序排列;
σr(H)为Hankerl矩阵H的第r个特征值,特征值按降序排列;
步骤106求取低秩逼近信号矩阵,具体为:低秩逼近信号矩阵记为
Figure FDA0003353832180000021
所用求取公式为:
Figure FDA0003353832180000022
其中:
Figure FDA0003353832180000023
为收缩化的特征值矩阵,
U为所述Hankerl矩阵H的左特征矢量矩阵,
V为所述Hankerl矩阵H的右特征矢量矩阵;
步骤107求取滤除噪声后的信号序列,具体为:滤除噪声后的信号序列为Snew,其第l个元素记为
Figure FDA0003353832180000024
所用求取公式为:
Figure FDA0003353832180000025
其中:
m=1,2,…,N为第一求和参数,
n=1,2,…,N为第二求和参数,
Figure FDA0003353832180000026
表示所述低秩逼近信号矩阵
Figure FDA0003353832180000027
的第m行第n列元素,
l=1,2,…,N为元素序号。
2.一种利用Hankerl矩阵的功率信号滤波***,其特征在于,包括:
模块201获取按时间顺序采集的信号序列S;
模块202生成Hankerl矩阵,具体为:Hankerl矩阵记为H,其第i行第j列元素记为hij,所用求取公式为:
Figure FDA0003353832180000028
其中:
i=1,2,…,N为行序号,
j=1,2,…,N为列序号,
Figure FDA0003353832180000029
表示所述信号序列S的第|i-1+j|N个元素,
|i-1+j|N表示以N为模对i-1+j取余数,
N为所述信号序列S的长度;
模块203求取最佳秩判断阈值,具体为:最佳秩判断阈值记为ε,所用求取公式为:
ε=ln(SNR+1)+2lnσ0
其中:
σ0为所述信号序列S的均方差,
SNR为所述信号序列S的信噪比;
模块204求取最佳秩,具体为:最佳秩记为R,所用求取公式为:
Figure FDA0003353832180000031
其中:
count{*}表示集合*中元素的个数,
*表示任一集合,
σj(H)为Hankerl矩阵H的第j个特征值,特征值按降序排列;
σmed(H)为Hankerl矩阵H的中值特征值;
模块205求取收缩化特征值,具体为:第r个收缩化特征值记为
Figure FDA0003353832180000032
所用求取公式为:
Figure FDA0003353832180000033
Figure FDA0003353832180000034
为幂收缩因子,
r=1,2,…,R为收缩序号,
σR+1(H)为Hankerl矩阵H的第R+1个特征值,特征值按降序排列;
σr(H)为Hankerl矩阵H的第r个特征值,特征值按降序排列;
模块206求取低秩逼近信号矩阵,具体为:低秩逼近信号矩阵记为
Figure FDA0003353832180000035
所用求取公式为:
Figure FDA0003353832180000036
其中:
Figure FDA0003353832180000037
为收缩化的特征值矩阵,
U为所述Hankerl矩阵H的左特征矢量矩阵,
V为所述Hankerl矩阵H的右特征矢量矩阵;
模块207求取滤除噪声后的信号序列,具体为:滤除噪声后的信号序列为Snew,其第l个元素记为
Figure FDA0003353832180000038
所用求取公式为:
Figure FDA0003353832180000039
其中:
m=1,2,…,N为第一求和参数,
n=1,2,…,N为第二求和参数,
Figure FDA00033538321800000310
表示所述低秩逼近信号矩阵
Figure FDA00033538321800000311
的第m行第n列元素,
l=1,2,…,N为元素序号。
CN202011229098.XA 2020-11-06 2020-11-06 一种利用Hankerl矩阵的功率信号滤波方法和*** Expired - Fee Related CN112347922B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011229098.XA CN112347922B (zh) 2020-11-06 2020-11-06 一种利用Hankerl矩阵的功率信号滤波方法和***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011229098.XA CN112347922B (zh) 2020-11-06 2020-11-06 一种利用Hankerl矩阵的功率信号滤波方法和***

Publications (2)

Publication Number Publication Date
CN112347922A CN112347922A (zh) 2021-02-09
CN112347922B true CN112347922B (zh) 2022-02-08

Family

ID=74428378

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011229098.XA Expired - Fee Related CN112347922B (zh) 2020-11-06 2020-11-06 一种利用Hankerl矩阵的功率信号滤波方法和***

Country Status (1)

Country Link
CN (1) CN112347922B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111666870A (zh) * 2020-06-04 2020-09-15 广东石油化工学院 一种利用二次约束的功率信号重构方法和***
CN111680590A (zh) * 2020-05-26 2020-09-18 广东石油化工学院 一种利用收缩梯度的功率信号滤波方法和***
CN111737645A (zh) * 2020-06-28 2020-10-02 广东石油化工学院 一种利用预测矩阵的功率信号重构方法和***

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101667141B1 (ko) * 2015-03-11 2016-10-25 한국과학기술원 소멸필터를 이용한 고속 mr 영상 복원 알고리듬 개발
CN105205461B (zh) * 2015-09-18 2018-05-11 中国石油大学(华东) 一种用于模态参数识别的信号降噪方法
CN106446868A (zh) * 2016-10-13 2017-02-22 成都芯安尤里卡信息科技有限公司 一种基于emd与奇异值差分谱的侧信道信号特征提取方法
US11138513B2 (en) * 2017-06-13 2021-10-05 Princeton University Dynamic learning system
CN110826017A (zh) * 2019-09-25 2020-02-21 中国地质大学(武汉) 一种基于参数优化Hankel矩阵和奇异值分解的信号去噪方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111680590A (zh) * 2020-05-26 2020-09-18 广东石油化工学院 一种利用收缩梯度的功率信号滤波方法和***
CN111666870A (zh) * 2020-06-04 2020-09-15 广东石油化工学院 一种利用二次约束的功率信号重构方法和***
CN111737645A (zh) * 2020-06-28 2020-10-02 广东石油化工学院 一种利用预测矩阵的功率信号重构方法和***

Also Published As

Publication number Publication date
CN112347922A (zh) 2021-02-09

Similar Documents

Publication Publication Date Title
CN111680590A (zh) 一种利用收缩梯度的功率信号滤波方法和***
CN111666870A (zh) 一种利用二次约束的功率信号重构方法和***
CN112434567B (zh) 一种利用噪声抖动性质的功率信号滤波方法和***
CN112347922B (zh) 一种利用Hankerl矩阵的功率信号滤波方法和***
CN111832474A (zh) 一种利用能量尺度的功率信号滤波方法和***
CN110244115B (zh) 一种基于信号连接性的负载开关事件检测方法及***
CN111639606A (zh) 一种利用Dantzig总梯度最小化的功率信号滤波方法和***
CN110221119B (zh) 基于功率与akie融合信息的负载开关事件检测方法和***
CN111830405A (zh) 一种利用频率差异量的负荷开关事件检测方法和***
CN111929608A (zh) 一种负载开关事件检测方法及***
CN110542855B (zh) 基于离散余弦变换的负荷开关事件检测方法和***
CN111737645A (zh) 一种利用预测矩阵的功率信号重构方法和***
CN112257576B (zh) 一种利用Maha距离测度的负荷开关事件检测方法和***
CN112307986B (zh) 一种利用高斯梯度的负荷开关事件检测方法和***
CN110702981A (zh) 一种利用分类树的负荷开关事件检测方法和***
CN112183410A (zh) 一种利用ds理论的功率信号滤波方法和***
CN110514884B (zh) 一种基于延迟矢量的功率信号滤波方法和***
CN111679120A (zh) 一种利用解耦合的功率信号滤波方法和***
CN112329637B (zh) 一种利用模式特征的负荷开关事件检测方法和***
CN112270282B (zh) 一种利用矩阵谱模的功率信号滤波方法和***
CN112307997B (zh) 一种利用主模式分解的功率信号重构方法和***
CN111665389A (zh) 一种利用随机投影的功率信号滤波方法和***
CN112180153A (zh) 一种利用KUllback-Leibler距离的负荷开关事件检测方法和***
CN110784019A (zh) 一种利用Mercer投影的功率信号滤波方法和***
CN112180155A (zh) 一种利用紧支集的负荷开关事件检测方法和***

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20220208