CN113469502A - 一种基于字典学习的区别分组非侵入式负荷分解方法 - Google Patents
一种基于字典学习的区别分组非侵入式负荷分解方法 Download PDFInfo
- Publication number
- CN113469502A CN113469502A CN202110637292.XA CN202110637292A CN113469502A CN 113469502 A CN113469502 A CN 113469502A CN 202110637292 A CN202110637292 A CN 202110637292A CN 113469502 A CN113469502 A CN 113469502A
- Authority
- CN
- China
- Prior art keywords
- dictionary
- grouping
- decomposition
- learning
- representing
- 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.)
- Granted
Links
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 78
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000012549 training Methods 0.000 claims abstract description 46
- 238000012360 testing method Methods 0.000 claims description 14
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000010606 normalization Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000007781 pre-processing Methods 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 4
- 230000015556 catabolic process Effects 0.000 claims description 3
- 238000006731 degradation reaction Methods 0.000 claims description 3
- 230000004069 differentiation Effects 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 5
- 238000012544 monitoring process Methods 0.000 abstract description 3
- 238000005070 sampling Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 6
- 238000005406 washing Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 238000005265 energy consumption Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000006399 behavior Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000011478 gradient descent method Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000013101 initial test Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 108010076504 Protein Sorting Signals Proteins 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000007787 long-term memory Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000006403 short-term memory Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Tourism & Hospitality (AREA)
- Evolutionary Biology (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Health & Medical Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Development Economics (AREA)
- Primary Health Care (AREA)
- Operations Research (AREA)
- Public Health (AREA)
- Educational Administration (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- General Health & Medical Sciences (AREA)
- Water Supply & Treatment (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Compression, Expansion, Code Conversion, And Decoders (AREA)
Abstract
本发明涉及电力负荷分解技术领域,公开了一种基于字典学习的区别分组非侵入式负荷分解方法,采用区别字典学习算法以及分组思想训练出每个设备的字典,再根据字典进行稀疏编码,实现负荷的非侵入式分解。利用字典学习方法能够较快的学习到各类用电设备的用电模式,使得该算法迁移到不同地方时具有较强的适应能力。同时使所训练的字典能够较好的针对特定设备,这样训练的字典保留了设备各自的独特性,能够降低因为特征相近所带来的误判,使得负荷分解结果准确性更好。对训练样本数据进行设备集分组,能够将非侵入式负荷监测问题由传统的整体分解任务转化为递归分解任务,提高负荷分解结果的准确性和稳定性,避免出现功率相近电器分解效果差的现象。
Description
技术领域
本发明涉及电力负荷分解技术领域,特别涉及一种基于字典学习的区别分组非侵入式负荷分解方法。
背景技术
随着现代社会的发展,人们对能源的需求也同步激增。能源危机的到来,使得节能问题引起了全社会的高度重视,提高能源效率是减少碳排放、缓解能源危机的重要途径。电能作为一种经济、实用、清洁的可再生二次能源,能够实现灵活调节和精确控制,使终端能源利用效率提高。智能电网是一个被广泛讨论的可以提高能源效率的概念。负荷分解技术作为智能电网实现节能的关键技术之一,主要是通过将负荷总用电信息分解为各时段每个负荷的运行状态和对应能耗信息,使得用户获取详细用电信息,以引导用户调整用电习惯达到节能降耗的目的。
目前负荷分解技术主要包括侵入式和非侵入式两类。侵入式负荷分解需要为每个用电设备安装传感器以获得详细用电数据,其优点是数据准确可靠,缺点是可操作性低、成本高、用户接受程度低。而非侵入式负荷分解只需采集总能耗信息就可以利用各用电设备的工作特征识别出单个设备能耗信息。非侵入式负荷分解技术弥补了传统侵入式的缺陷,是未来负荷分解的主要发展方向。
目前国内外都对非侵入式负荷分解方法进行了大量的研究,如组合优化算法,支持向量机,稀疏编码,隐马尔可夫模型,深度长短时记忆网络等算法。相较于其他方法深度神经网络和稀疏编码方法能够取得更好结果。但是深度神经网络随着网络层数增加,会导致参数量的增加,模型的训练难度也随之加大。稀疏编码算法在这方面具有一定优势能提升分解性能,但在识别的准确度上仍然有较大的提升空间。因此,如何提高现有非侵入式负荷分解方法的效率与准确性,是现有技术亟待解决的问题之一。
发明内容
本发明的主要目的为提供一种基于字典学习的区别分组非侵入式负荷分解方法,旨在解决现有技术中非侵入式负荷分解方法的效率与准确性难以提高的技术问题。
为实现上述发明目的,本发明采用了如下技术方案:
一种基于字典学习的区别分组非侵入式负荷分解方法,包括以下步骤:
S1:获取设备的数据集,对所述数据集进行数据预处理,并将预处理后的数据分为训练样本数据和测试样本数据;所述数据集包括原始测量的功率信号;
S2:对训练样本数据进行设备集分组,所述设备集分组包括第一分组与第二分组,所述第一分组包括一个设备,所述第二分组包括多个设备;
S3:根据训练样本数据,对第一分组和第二分组进行区别字典训练,分别得到对应的第一字典和第二字典;
S4:基于区别指标,对所述第一字典与所述第二字典区别衡量,并记录区别指标值;
S5:在第一分组和第二分组中选择出字典区别指标最小的分组,将其确定为最优分组;
S6:进行功率分解,利用最优分组训练好的第一字典和第二字典对预处理好的聚合功率信号数据进行稀疏编码,得到稀疏表示,并将稀疏表示与区别字典相乘得到分解结果;
S7:通过剔除最优分组里面被分出来的设备及其功率预测值来更新设备集和总功率;
S8:重复步骤S2-S7,得到每次迭代的分解结果,直到所述第二分组中只包括一个设备,进入S9;
S9:将每次递归迭代结果作为一个设备的预测结果,并对该预测结果进行反归一化处理及输出。
优选的,步骤S1中,所述数据预处理包括以下步骤:
(1-1)对所述数据集进行选取;
(1-2)将所述数据集的进行除燥;
(1-3)将除燥后的数据集进行归一化;
(1-4)生成样本数据。
优选的,步骤S2中,所述设备集分组包括第一分组与第二分组,所述第一分组包括一个设备,命名为one;所述第二分组包括多个设备,命名为rest。
优选的,步骤S3中,所述区别字典训练包括以下步骤:
(2-1)对每个设备的功率信号Xi∈Rm×n学习一个字典Di∈Rm×q和一个稀疏表示Zi∈Rq×n使得每个设备的功率信号等于一个字典和一个稀疏表示的乘积,其计算公式为:
Xi=DiZi,i=1...N (1)
其中,i表示设备编号,N表示总设备数,Xi为设备的功率信号,Di为字典,Zi为稀疏表示;
(2-2)利用预处理好的训练样本数据建立学习字典,所述学习字典的基本模型如下:
其中,λ1是权衡重构误差与稀疏性的参数,添加约束Di≥0,Zi≥0的原因是设备的功率信号总是非负的,dj表示字典Di的第j列,即Di=[d1,d2,...,dj,...,dq],亦称为字典原子,令||dj||2=1是为了防止出现退化解;||·||F表示F范数运算,||·||1表示1范数运算;
(2-3)加入第一正则化项,提高字典的辨别能力,所述第一正则化项的计算公式为:
(2-4)加入第二正则化项,提高字典间的区别,所述第二正则化项的计算公式为:
(2-5)联立方程(1)~(4),得到第一分组one和第二分组rest最终的区别字典学习的模型。
(2-6)根据第一分组one和第二分组rest最终的区别字典学习的模型,分别求解分组第一分组one和第二分组rest的字典:
优选的,步骤(2-5)中,所述第一分组one的区别字典学习的模型为:
其中,X1、Z1、分别表示第一分组one的设备功耗和第一稀疏表示,onedj表示第一字典oneD的原子;λ1是权衡重构误差与稀疏性的参数,α是第一正则化参数,β是第二正则化参数,||·||F表示F范数运算,||·||1表示1范数运算。
优选的,步骤(2-5)中,所述第二分组rest的区别字典学习的模型为:
其中,X2、Z2分别表示第二分组rest的设备功耗和第二稀疏表示,restdj分别表示第二字典restD的原子;λ1是权衡重构误差与稀疏性的参数,α是第一正则化参数,β是第二正则化参数,||·||F表示F范数运算,||·||1表示1范数运算。
优选的,步骤S4中,对所述第一字典与所述第二字典区别衡量包括定义区别矩阵A=oneDTrestD,且A∈Rq×q,并将训练所得的第一字典oneD和第二字典restD作为输入,计算得到两者的区别指标。
优选的,所述计算区别指标的计算公式为:
其中,DI表示区别指标,aij表示区别矩阵A第i列第j行的元素。
优选的,步骤S6中,所述功率分解包括利用聚合信号X、最优分组结构和训练所得第一字典oneD、第二字典restD建立稀疏分解模式。
优选的,所述的稀疏分解模式计算公式为:
其中,Z1、Z2分别表示第一分组与第二分组的中各设备的稀疏表示,oneD表示所述第一字典的原子,restD表示所述第二字典的原子;λ2是权衡重构误差与稀疏性的参数。
本发明的有益效果为:本发明提供了一种基于字典学习的区别分组非侵入式负荷分解方法,采用区别字典学习算法以及分组思想训练出每个设备的字典,再根据字典进行稀疏编码,实现负荷的非侵入式分解。利用字典学习方法能够较快的学习到各类用电设备的用电模式,使得该算法迁移到不同地方时具有较强的适应能力。同时,在训练阶段考虑了字典内部的非连贯性和字典间的区别性,使所训练的字典能够较好的针对特定设备。这样训练的字典保留了设备各自的独特性,能够降低因为特征相近所带来的误判,使得负荷分解结果准确性更好。对训练样本数据进行设备集分组,能够将非侵入式负荷监测问题由传统的整体分解任务转化为递归分解任务,提高负荷分解结果的准确性和稳定性,避免出现功率相近电器分解效果差的现象。本发明能够应用于负荷分解、智能电网、用电行为分析、负荷节能等多个技术领域,对居民、电力公司以及全社会都具有重要意义。
附图说明
图1是本发明实施例提供的一种基于字典学习的区别分组非侵入式负荷分解方法的整体流程图;
图2是是本发明实施例提供的一种基于字典学习的区别分组非侵入式负荷分解方法的具体流程图;
图3是是本发明实施例提供的一种基于字典学习的区别分组非侵入式负荷分解方法的分组结构示例图;
本发明目的的实现、功能特点及优点将结合实施例,参照附图做进一步说明。
具体实施方式
应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
请参阅图1、图2及图3,本发明实施例提供了一种基于字典学习的区别分组非侵入式负荷分解方法,包括以下步骤:
S1:获取设备的数据集,对所述数据集进行数据预处理,并将预处理后的数据分为训练样本数据和测试样本数据;所述数据集包括原始测量的功率信号;
在本发明实施例步骤S1中,所述测试样本数据计算设备中的设备总功率。
在本发明实施例步骤S1中,所述数据预处理的方法包括
(1-1)对所述数据集进行选取;选用REDD公开数据集进行实验,此数据集包含六个家庭的用电数据,电源级采样频率为1Hz,设备级的电器采样频率为1/3Hz。针对每个家庭的聚合功率信号和单设备功率信号进行数据对齐和数据填充操作。
(1-2)将所述数据集的进行除燥;测量获得的功率原始数据往往存在噪声,所以选择中值滤波方法进行去噪,保证在滤去噪声的同时保护信号边缘。
(1-3)将除燥后的数据集进行归一化;对信号进行归一化处理,通过采用最大值最小值归一化方法,使结果映射到区间[0,1]内,得到所需要的序列。
(1-4)生成样本数据。将归一化处理后的负荷序列分解为初始训练序列和初始测试序列;再对所述初始训练序列和初始测试序列通过设定窗口长度和步长进行滑动处理得到训练样本数据和测试样本数据。
具体的,在本发明实施例中,所述数据预处理是利用REDD数据集中house_1的数据来获取训练样本数据和测试样本数据,选用其中的洗碗机、冰箱、微波炉和洗衣机四种用电设备用做目标设备进行负荷分解,即设备集B={洗碗机,冰箱,微波炉,洗衣机},设备数N=4。
首先由于设备级的采样频率(1/3Hz)和电源级的采样频率(1Hz)不一致,需要先将设备级的数据和电源级的数据根据时间戳进行对齐,对于缺失的设备级数据进行后填充;
针对设备级的数据通过中值滤波去除噪声和毛刺,并保留信号边缘。例如对功率信号序列X进行滤波处理时,定义窗口长度为L=2×5+1=11。对窗口内的信号样本值按从小到大的顺序排列后得X(j-5),…,X(j),…,X(j+5),其中值X(j)被定义为中值滤波的输出值;
对经过中值滤波的数据进行下采样。例如存在功率信号X={X1,X2,…,Xn},下采样时间为1min,即是60个间隔采样一次,得到下采样之后的功率信号为X’={X1,X61,X121,…};
然后对数据进行归一化处理以消除指标间的纲量影响,采用最大值最小值归一化方法,使其映射到区间[0,1]。归一化函数为y=(x-xmin)/(xmax-xmin),其中,x为样本数据,xmax为样本数据最大值,xmin为样本数据的最小值,y为归一化后得到的样本数据;
为了增加样本数量,通过重叠滑动窗口获得训练样本。滑动窗口大小和步长分别为180和30,代表的时间分别是3小时和半小时。对于长为z的训练序列来说需要滑动(z-180)/30+1次来获取训练样本;
最后获取测试样本,测试样本采用非重叠滑动方式获取,长度为h的测试序列需要滑动h/180次。其中训练样本和测试样本没有重合,并且训练样本和测试样本的数据量比例为4:1;
对于设备集B中的四个目标设备均执行上述操作就可以得到四组一一对应的用于负荷分解的训练样本数据和测试样本数据,分别为Xa、Xb、Xc、Xd和Xa’、Xb’、Xc’、Xd’。
S2:对训练样本数据进行设备集分组,所述设备集分组包括第一分组与第二分组,所述第一分组包括一个设备,命名为one,所述第二分组包括多个设备,命名为rest;
如上述S2步骤所述的,针对设备集B选取其中一个设备作为第一分组,并命名为one。剩下的设备组成第二分组,命名为rest。例如one={洗碗机},rest={冰箱,微波炉,洗衣机}。由于用电设备之间的功耗具有相似性,这种相似性会对分解结果产生极大的影响,降低分解精度,通过将设备集进行分组能够消减这种相似性,以达到更好的分解效果。
S3:根据训练样本数据,对第一分组和第二分组进行区别字典训练,分别得到对应的第一字典和第二字典;
在本发明实施例步骤S3中,所述区别字典训练包括以下步骤:
(2-1)对每个设备的功率信号Xi∈Rm×n学习一个字典Di∈Rm×q和一个稀疏表示Zi∈Rq×n使得每个设备的功率信号等于一个字典和一个稀疏表示的乘积,其计算公式为:
Xi=DiZi,i=1...N (1)
其中,i表示设备编号,N表示总设备数,Xi为设备的功率信号,Di为字典,Zi为稀疏表示;
(2-2)利用预处理好的训练样本数据建立学习字典,所述学习字典的基本模型如下:
其中,λ1是权衡重构误差与稀疏性的参数,添加约束Di≥0,Zi≥0的原因是设备的功率信号总是非负的,dj表示字典Di的第j列,即Di=[d1,d2,...,dj,...,dq],亦称为字典原子,令||dj||2=1是为了防止出现退化解;||·||F表示F范数运算,||·||1表示1范数运算;
(2-3)加入第一正则化项,提高字典的辨别能力,所述第一正则化项的计算公式为:
(2-4)加入第二正则化项,提高字典间的区别,所述第二正则化项的计算公式为:
(2-5)联立方程(1)~(4),得到最终区别字典学习的模型:
所述第一分组one的区别字典学习的模型为:
其中,X1、Z1、分别表示第一分组one的设备功耗和第一稀疏表示,onedj表示第一字典oneD的原子;λ1是权衡重构误差与稀疏性的参数,α是第一正则化参数,β是第二正则化参数,||·||F表示F范数运算,||·||1表示1范数运算。
所述第二分组rest的区别字典学习的模型为:
其中,X2、Z2分别表示第二分组rest的设备功耗和第二稀疏表示,restdj分别表示第二字典restD的原子;λ1是权衡重构误差与稀疏性的参数,α是第一正则化参数,β是第二正则化参数,||·||F表示F范数运算,||·||1表示1范数运算。
(2-6)利用式(5)和式(6)分别求解第一分组one和第二分组rest的字典。
在本发明实施例中,利用交替最小化求解上述问题,当近似误差变为收敛或迭代次数达到指定数值时,迭代结束。将上述问题转化为如下优化问题:
首先,固定第一稀疏表示Z1求解第一字典oneD,式(5)转换为:
将式(7)采用梯度下降法解第一字典oneD,求解形式如下:
其中,μ表示梯度下降的步长,[·]+表示负数置零,以确保oneD≥0。
然后再对oneD的每列onedj进行单位化||onedj||2=1。
同理,固定oneD求解Z1,式(5)转换为:
上式属于l1最小化问题,且带有非负约束,采用非负矩阵分解求解,求解形式如下:
Z1=Z1.*(2oneDTX1)./(2oneDToneDZ1+λ1) (10)
同理,可求得第二字典restD和第二稀疏表示Z2,求解形式如下:
对第二字典restD的每列restdj进行单位化||restdj||2=1
Z2=Z2.*(2restDTX2)./(2restDTrestDZ2+λ1) (12)
通过交替最小化求解上述最小化问题可以得到第一字典oneD和第二字典restD。求解字典时使用梯度下降法,求解稀疏表示时使用非负矩阵分解,以保证非负约束。交替最小化迭代的结束条件是目标函数的变化率小于10-5,即上一次迭代的目标函数值减本次迭代的目标函数值的绝对值除以上一次迭代的目标函数值的值小于10-5。
S4:基于区别指标,对所述第一字典与所述第二字典区别衡量,并记录区别指标值;
在本发明实施例步骤S4中,对所述第一字典与所述第二字典区别衡量包括定义区别矩阵A=oneDTrestD,且A∈Rq×q,并将训练所得的第一字典oneD和第二字典restD作为输入,计算得到两者的区别指标。
在本发明实施例中,所述计算区别指标的计算公式为:
其中,DI表示区别指标,aij表示区别矩阵A第i列第j行的元素。DI越小说明字典oneD和restD之间的区别越大。
S5:在第一分组和第二分组中选择出字典区别指标最小的分组,将其确定为最优分组;
如上述S5步骤所述的,重复S2-S4,递归的从设备集B中选择一个设备作为one分组,分别进行区别字典学习,记录字典区别评估指标。从记录的评估指标中选择出最小值,将此分组作为最优分组,即是在此分组情况下设备集B能够训练出最具有区别性的两个字典,能够带来更高的分解精度。
S6:进行功率分解,利用最优分组训练好的第一字典和第二字典对预处理好的聚合功率信号数据进行稀疏编码,得到稀疏表示,并将稀疏表示与区别字典相乘得到分解结果;
在本发明实施例步骤S6中,所述功率分解包括利用聚合信号X、最优分组结构和训练所得第一字典oneD、第二字典restD建立稀疏分解模式。
所述的稀疏分解模式计算公式为:
其中,Z1、Z2分别表示第一分组与第二分组的中各设备的稀疏表示,oneD表示所述第一字典的原子,restD表示所述第二字典的原子;λ2是权衡重构误差与稀疏性的参数。
具体的,输入待分解的总功率信号X,利用最优分组结构和和训练好的字典oneD和restD对其进行拟合,分别求得稀疏编码系数,求解方法选择非负矩阵分解。然后将对应的字典和稀疏编码系数相乘得到本次迭代的分解结果。例如one分组的分解结果为rest分组的分解结果为其中Z′1和Z′2利用下式求解:
X表示待分解总功率,针对第一分组one={洗碗机},第二分组rest={冰箱,微波炉,洗衣机}有X=Xa’+Xb’+Xc’+Xd’。其中,λ2用于权衡重构误差和稀疏度,因为最终分解目的是使重构误差较小,所以将λ2设置为一个较小值0.001。
S7:通过剔除最优分组里面被分出来的设备及其功率预测值来更新设备集和总功率;
如上述S7所述的,针对选出的最优分组可以得到分组情况和该分组的分解效果。设备集B的更新是从B中剔除最优第一分组one所包含的设备,例如目前设备集={洗碗机,冰箱,微波炉,洗衣机},最优分组的情况是第一分组one={微波炉},第二分组rest={洗碗机,冰箱,洗衣机},那么更新后的设备集={洗碗机,冰箱,洗衣机}。总功率X的更新是利用原来的总功率X1减去最优第一分组one的分解结果oneDZ'1最优,即X=X1-oneDZ'1最优。
S8:重复步骤S2-S7,得到每次迭代的分解结果,直到所述第二分组中只包括一个设备,进入S9;
如上述S8所述的,将更新后的设备集和总功率以及训练样本数据作为输入,并重复步骤S2-S7,每次迭代都针对选择出来的最优分组结构进行功率分解,直到设备集中只剩下一个设备时停止迭代。
S9:将每次递归迭代结果作为一个设备的预测结果,并对该预测结果进行反归一化处理及输出。
如图3所示是一个分组递归分解的分组结构示例图,具体的,每次迭代都能得到一个设备的预测结果,对该预测结果进行反归一化处理,即是将得到的预测数据从[0,1]区间还原到设备实际功率消耗值。反归一化的公式为:y’=x’(xmax-xmin)+xmin,其中,x’为预测数据,xmax为样本数据最大值,xmin为样本数据的最小值,y’为反归一化后得到的预测数据,即是最终预测结果。
本发明提供了一种基于字典学习的区别分组非侵入式负荷分解方法,采用区别字典学习算法以及分组思想训练出每个设备的字典,再根据字典进行稀疏编码,实现负荷的非侵入式分解。利用字典学习方法能够较快的学习到各类用电设备的用电模式,使得该算法迁移到不同地方时具有较强的适应能力。同时,在训练阶段考虑了字典内部的非连贯性和字典间的区别性,使所训练的字典能够较好的针对特定设备。这样训练的字典保留了设备各自的独特性,能够降低因为特征相近所带来的误判,使得负荷分解结果准确性更好。对训练样本数据进行设备集分组,能够将非侵入式负荷监测问题由传统的整体分解任务转化为递归分解任务,提高负荷分解结果的准确性和稳定性,避免出现功率相近电器分解效果差的现象。本发明能够应用于负荷分解、智能电网、用电行为分析、负荷节能等多个技术领域,对居民、电力公司以及全社会都具有重要意义。
上列详细说明是针对本发明可行实施例的具体说明,以上实施例并非用以限制本发明的专利范围,凡未脱离本发明所为的等效实施或变更,均应包含于本案的专利范围中。
Claims (10)
1.一种基于字典学习的区别分组非侵入式负荷分解方法,其特征在于,包括以下步骤:
S1:获取设备的数据集,对所述数据集进行数据预处理,并将预处理后的数据分为训练样本数据和测试样本数据;所述数据集包括原始测量的功率信号;
S2:对训练样本数据进行设备集分组,所述设备集分组包括第一分组与第二分组,所述第一分组包括一个设备,所述第二分组包括多个设备;
S3:根据训练样本数据,对第一分组和第二分组进行区别字典训练,分别得到对应的第一字典和第二字典;
S4:基于区别指标,对所述第一字典与所述第二字典区别衡量,并记录区别指标值;
S5:在第一分组和第二分组中选择出字典区别指标最小的分组,将其确定为最优分组;
S6:进行功率分解,利用最优分组训练好的第一字典和第二字典对预处理好的聚合功率信号数据进行稀疏编码,得到稀疏表示,并将稀疏表示与区别字典相乘得到分解结果;
S7:通过剔除最优分组里面被分出来的设备及其功率预测值来更新设备集和总功率;
S8:重复步骤S2-S7,得到每次迭代的分解结果,直到所述第二分组中只包括一个设备,进入S9;
S9:将每次递归迭代结果作为一个设备的预测结果,并对该预测结果进行反归一化处理及输出。
2.根据权利要求1所述的基于字典学习的区别分组非侵入式负荷分解方法,其特征在于,步骤S1中,所述数据预处理包括以下步骤:
(1-1)对所述数据集进行选取;
(1-2)将所述数据集的进行除燥;
(1-3)将除燥后的数据集进行归一化;
(1-4)生成样本数据。
3.根据权利要求1所述的基于字典学习的区别分组非侵入式负荷分解方法,其特征在于,步骤S2中,所述设备集分组包括第一分组与第二分组,所述第一分组包括一个设备,命名为one;所述第二分组包括多个设备,命名为rest。
4.根据权利要求1所述的基于字典学习的区别分组非侵入式负荷分解方法,其特征在于,步骤S3中,所述区别字典训练包括以下步骤:
(2-1)对每个设备的功率信号Xi∈Rm×n学习一个字典Di∈Rm×q和一个稀疏表示Zi∈Rq×n使得每个设备的功率信号等于一个字典和一个稀疏表示的乘积,其计算公式为:
Xi=DiZi,i=1...N (1)
其中,i表示设备编号,N表示总设备数,Xi为设备的功率信号,Di为字典,Zi为稀疏表示;
(2-2)利用预处理好的训练样本数据建立学习字典,所述学习字典的基本模型如下:
s.t.Di≥0,Zi≥0,||dj||2=1 i=1...N j=1...q (2)
其中,λ1是权衡重构误差与稀疏性的参数,添加约束Di≥0,Zi≥0的原因是设备的功率信号总是非负的,dj表示字典Di的第j列,即Di=[d1,d2,...,dj,...,dq],亦称为字典原子,令||dj||2=1是为了防止出现退化解;||·||F表示F范数运算,||·||1表示1范数运算;
(2-3)加入第一正则化项,提高字典的辨别能力,所述第一正则化项的计算公式为:
(2-4)加入第二正则化项,提高字典间的区别,所述第二正则化项的计算公式为:
(2-5)联立方程(1)~(4),得到第一分组one和第二分组rest最终的区别字典学习的模型。
(2-6)根据第一分组one和第二分组rest最终的区别字典学习的模型,分别求解分组第一分组one和第二分组rest的字典。
7.根据权利要求1所述的基于字典学习的区别分组非侵入式负荷分解方法,其特征在于,步骤S4中,对所述第一字典与所述第二字典区别衡量包括定义区别矩阵A=oneDTrestD,且A∈Rq×q,并将训练所得的第一字典oneD和第二字典restD作为输入,计算得到两者的区别指标。
9.根据权利要求1所述的基于字典学习的区别分组非侵入式负荷分解方法,其特征在于,步骤S6中,所述功率分解包括利用聚合信号X、最优分组结构和训练所得第一字典oneD、第二字典restD建立稀疏分解模式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110637292.XA CN113469502B (zh) | 2021-06-08 | 2021-06-08 | 一种基于字典学习的区别分组非侵入式负荷分解方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110637292.XA CN113469502B (zh) | 2021-06-08 | 2021-06-08 | 一种基于字典学习的区别分组非侵入式负荷分解方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113469502A true CN113469502A (zh) | 2021-10-01 |
CN113469502B CN113469502B (zh) | 2023-08-15 |
Family
ID=77869181
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110637292.XA Active CN113469502B (zh) | 2021-06-08 | 2021-06-08 | 一种基于字典学习的区别分组非侵入式负荷分解方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113469502B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116859169A (zh) * | 2023-09-01 | 2023-10-10 | 国网浙江省电力有限公司余姚市供电公司 | 一种非侵入式电力负荷分解方法及系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130289788A1 (en) * | 2012-04-25 | 2013-10-31 | Abhay Gupta | Energy Disaggregation Techniques for Low Resolution Whole-House Energy Consumption Data |
WO2018149133A1 (zh) * | 2017-02-17 | 2018-08-23 | 深圳大学 | 基于核非负矩阵分解的字典学习和稀疏特征表示的人脸识别方法及系统 |
CN110443404A (zh) * | 2019-06-28 | 2019-11-12 | 天津大学 | 一种基于热量组分解析的建筑空调负荷预测方法 |
CN111311434A (zh) * | 2020-02-14 | 2020-06-19 | 广州水沐青华科技有限公司 | 用电设备负荷分离方法、装置、计算机设备和存储介质 |
CN111461189A (zh) * | 2020-03-23 | 2020-07-28 | 东南大学 | 基于非侵入式量测技术的居民节律检测方法 |
CN112348096A (zh) * | 2020-11-11 | 2021-02-09 | 合肥工业大学 | 非侵入式负荷分解方法及系统 |
-
2021
- 2021-06-08 CN CN202110637292.XA patent/CN113469502B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130289788A1 (en) * | 2012-04-25 | 2013-10-31 | Abhay Gupta | Energy Disaggregation Techniques for Low Resolution Whole-House Energy Consumption Data |
WO2018149133A1 (zh) * | 2017-02-17 | 2018-08-23 | 深圳大学 | 基于核非负矩阵分解的字典学习和稀疏特征表示的人脸识别方法及系统 |
CN110443404A (zh) * | 2019-06-28 | 2019-11-12 | 天津大学 | 一种基于热量组分解析的建筑空调负荷预测方法 |
CN111311434A (zh) * | 2020-02-14 | 2020-06-19 | 广州水沐青华科技有限公司 | 用电设备负荷分离方法、装置、计算机设备和存储介质 |
CN111461189A (zh) * | 2020-03-23 | 2020-07-28 | 东南大学 | 基于非侵入式量测技术的居民节律检测方法 |
CN112348096A (zh) * | 2020-11-11 | 2021-02-09 | 合肥工业大学 | 非侵入式负荷分解方法及系统 |
Non-Patent Citations (1)
Title |
---|
ZIYUE JIA等: "Sequence to point learning based on bidirectional dilated residual network for non intrusive load monitoring", ARXIV, pages 1 - 10 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116859169A (zh) * | 2023-09-01 | 2023-10-10 | 国网浙江省电力有限公司余姚市供电公司 | 一种非侵入式电力负荷分解方法及系统 |
CN116859169B (zh) * | 2023-09-01 | 2023-11-17 | 国网浙江省电力有限公司余姚市供电公司 | 一种非侵入式电力负荷分解方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN113469502B (zh) | 2023-08-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110610280B (zh) | 一种电力负荷短期预测方法、模型、装置及系统 | |
CN110336534B (zh) | 一种基于光伏阵列电气参数时间序列特征提取的故障诊断方法 | |
CN109635928B (zh) | 基于深度学习模型融合的电压暂降原因识别方法 | |
CN109598451B (zh) | 一种基于pca与lstm神经网络的非侵入式负荷辨识方法 | |
CN110309603B (zh) | 一种基于风速特性的短期风速预测方法及系统 | |
CN111193256A (zh) | 基于变分模态分解和门控循环单元的电力负荷预测方法 | |
CN108376262B (zh) | 一种风电出力典型特性的分析模型构建方法 | |
CN109975673B (zh) | 一种光伏微电网直流侧故障电弧识别方法 | |
CN115130830B (zh) | 基于级联宽度学习和麻雀算法的非侵入式负荷分解方法 | |
Zhu et al. | New feature extraction method for photovoltaic array output time series and its application in fault diagnosis | |
CN113469502B (zh) | 一种基于字典学习的区别分组非侵入式负荷分解方法 | |
CN115100466A (zh) | 一种非侵入式负荷监测方法、装置及介质 | |
CN114139820A (zh) | 一种用于非侵入式电能负荷预测的改进模态分解方法 | |
CN117272230A (zh) | 一种基于多任务学习模型的非侵入式负荷监测方法和系统 | |
CN116561569A (zh) | 一种基于EO特征选择结合AdaBoost算法的工业电力负荷辨识方法 | |
CN116008747A (zh) | 基于小波阈值去噪的Yogi-mLSTM电缆局放识别方法与诊断系统 | |
CN114462306A (zh) | 基于变权重时域卷积网络的非侵入式用电负荷分解方法 | |
CN111008584B (zh) | 一种模糊自组织神经网络的电能质量量测缺失修复方法 | |
Liu et al. | Nonintrusive load disaggregation combining with external attention mechanism and Seq2Piont | |
Yu et al. | Nonintrusive load disaggregation method based on graph signal processing | |
El Bakali et al. | Data-Based Solar Radiation Forecasting with Pre-Processing Using Variational Mode Decomposition | |
CN112668770A (zh) | 一种基于信息与波形总体相似度的电力负荷预测方法 | |
CN113420436B (zh) | 基于gru和arima-t的非侵入式事件检测方法 | |
Xiong et al. | Research on Industrial Load Decomposition Based on Correlation Analysis and DNN | |
Tian et al. | An Efficient Short-Term Electricity Forecasting Approach Based on EEMD-LSTM Model with Feature Factors |
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 |