CN107895203A - 一种基于信号稀疏表示的建筑分项冷负荷获取方法 - Google Patents

一种基于信号稀疏表示的建筑分项冷负荷获取方法 Download PDF

Info

Publication number
CN107895203A
CN107895203A CN201711027913.2A CN201711027913A CN107895203A CN 107895203 A CN107895203 A CN 107895203A CN 201711027913 A CN201711027913 A CN 201711027913A CN 107895203 A CN107895203 A CN 107895203A
Authority
CN
China
Prior art keywords
msub
msup
refrigeration duty
matrix
mrow
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
Application number
CN201711027913.2A
Other languages
English (en)
Other versions
CN107895203B (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.)
Tianjin University
Original Assignee
Tianjin 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 Tianjin University filed Critical Tianjin University
Priority to CN201711027913.2A priority Critical patent/CN107895203B/zh
Publication of CN107895203A publication Critical patent/CN107895203A/zh
Application granted granted Critical
Publication of CN107895203B publication Critical patent/CN107895203B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/213Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
    • G06F18/2136Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on sparsity criteria, e.g. with an overcomplete basis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315Needs-based resource requirements planning or analysis

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Theoretical Computer Science (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Marketing (AREA)
  • Development Economics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Game Theory and Decision Science (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Educational Administration (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Evolutionary Biology (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Air Conditioning Control Device (AREA)

Abstract

一种基于信号稀疏表示的建筑分项冷负荷获取方法:字典学习,即:将总冷负荷分为四项分项冷负荷,将四项分项冷负荷的逐时历史数据分别写作矩阵形式;对所述的矩阵形式分别采用非负K‑SVD字典学习算法,得到四个字典矩阵;将总冷负荷分解,是利用稀疏表示相关算法,由总冷负荷和字典矩阵得到四个稀疏系数矩阵;分项冷负荷信号重组,是利用矩阵相乘的方法,将求得的对应于温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷的四个字典矩阵及四个稀疏系数矩阵对应相乘,获取由总冷负荷分解所得的温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷。本发明可以高效快捷地实现建筑总冷负荷信号的分解及建筑分项冷负荷动态数据的获取。

Description

一种基于信号稀疏表示的建筑分项冷负荷获取方法
技术领域
本发明涉及一种建筑分项冷负荷获取方法。特别是涉及一种将建筑总的逐时冷负荷分解为具有不同物理意义的分项逐时冷负荷的基于信号稀疏表示的建筑分项冷负荷获取方法。
背景技术
建筑节能问题已受到社会的广泛关注,其中,由于夏季供冷能耗在建筑整体能耗中占很大比例,供冷节能问题也已成为当今社会的关注焦点。在探讨建筑供冷节能问题时,不仅应关注供冷系统的总能耗,更应该分析各供冷子系统的用能规律,从而掌握建筑各部件的作用与表现,以优化建筑设计与指导建筑运行。供冷系统总能耗可以用建筑总冷负荷来表示,建筑总冷负荷由四项分项冷负荷组成,包括由于围护结构内外温差传热形成的冷负荷(简称温差冷负荷),由于通过窗户的太阳辐射形成的冷负荷(简称辐射冷负荷),由于引入室外新风形成的冷负荷(简称新风冷负荷),由于人员、灯光、设备等内扰因素形成的冷负荷(简称内扰冷负荷),对分项冷负荷的获取与分析对建筑的节能设计与运行有着重要的现实意义。然而,建筑的分项冷负荷是无法直接监测与测量的,只有建筑总冷负荷可以通过供冷系统的实际运行数据测得,因此如果能通过建筑的动态运行数据反向解析得到建筑的分项冷负荷,则可以高效便利地获取建筑分项冷负荷的动态数据,从而对建筑节能实现重要指导意义。
发明内容
本发明所要解决的技术问题是,提供一种通过建筑总冷负荷的动态数据,反向分解得到建筑分项冷负荷动态数据的基于信号稀疏表示的建筑分项冷负荷获取方法。
本发明所采用的技术方案是:一种基于信号稀疏表示的建筑分项冷负荷获取方法,包括如下步骤:
1)字典学习,包括:
(1)将总冷负荷分为四项分项冷负荷,包括:温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷,将四项分项冷负荷的逐时历史数据分别写作矩阵形式;
(2)对所述的矩阵形式分别采用非负K-SVD字典学习算法,得到四个字典矩阵;
2)将总冷负荷分解,是利用稀疏表示相关算法,由总冷负荷和字典矩阵得到四个稀疏系数矩阵;
3)分项冷负荷信号重组,是利用矩阵相乘的方法,将步骤1)和步骤2)求得的对应于温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷的四个字典矩阵及四个稀疏系数矩阵对应相乘,获取由总冷负荷分解所得的温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷。
步骤1)第(1)步具体是:将温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷四项分项冷负荷的逐时历史数据分别写作矩阵X1、X2、X3和X4的形式,其中X1、X2、X3和X4分别对应表示温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷的历史数据,且X1、X2、X3和X4矩阵的行数为24,表示一天24小时,列数为历史数据所包含的天数。
步骤1)第(2)步具体是:分别对矩阵X1、X2、X3和X4采用非负K-SVD字典学习算法,得到四个字典矩阵B1、B2、B3和B4,其中,四个字典矩阵B1、B2、B3和B4分别对应表示温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷的字典,且四个字典矩阵B1、B2、B3和B4矩阵的行数为24,表示一天24小时,列数为大于24且小于矩阵X1、X2、X3和X4列数的数值。
步骤2)包括:
(1)将用于分解的总冷负荷的逐时数据写作矩阵的形式,总冷负荷矩阵的行数为24,表示一天24小时,列数分别为总冷负荷逐时数据所包含的天数;
(2)设置矩阵A'1、A'2、A'3和A'4为任意非负矩阵,其中矩阵A'1、A'2、A'3和A'4的行数等于字典矩阵B1、B2、B3和B4的列数,矩阵A'1、A'2、A'3和A'4的列数等于总冷负荷矩阵的列数,矩阵A'1、A'2、A'3和A'4分别为温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷的稀疏系数矩阵;
(3)设置稀疏度λ为0至1之间的任意数值;
(4)循环迭代以下各步直至收敛
将A'1、A'2、A'3和A'4矩阵中的所有负数项赋值为0;
(5)得到各分项冷负荷的稀疏系数矩阵A'1、A'2、A'3和A'4
步骤3)包括:
(1)分别将四个字典矩阵B1、B2、B3和B4与四个稀疏系数矩阵A'1、A'2、A'3和A'4对应相乘;
(2)乘积所得的四个矩阵写作其中分别表示分解所得的温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷的估计值,的行数和列数与用于分解的总冷负荷矩阵相同,即,行数为24,列数分别为总冷负荷逐时数据所包含的天数,且每一列代表一天中24小时的逐时冷负荷。
本发明的一种基于信号稀疏表示的建筑分项冷负荷获取方法,所需数据较少,且算法运行速度较快,可以高效快捷地实现建筑总冷负荷信号的分解及建筑分项冷负荷动态数据的获取。解决了建筑分项冷负荷无法直接测量的缺陷,有利于优化建筑的节能设计与运行,本发明具有以下优点:
1、创新性地提出了通过将建筑总冷负荷信号分解来获取建筑各分项冷负荷的动态数据的方法。
2、该方法在获取建筑分项冷负荷动态数据的过程中无需得知气象参数、建筑结构与热工参数等不易获知的数据,在应用时有着便捷与高效的优势。
3、使用该方法所分解得到的各个分项冷负荷各自具有特定的物理意义,而不是单纯基于不同频率的冷负荷分项分解。
附图说明
图1是本发明一种基于信号稀疏表示的建筑分项冷负荷获取方法的流程图;
图2a是用于分解的总冷负荷示意图;
图2b是分解所得温差冷负荷的估计值与实际值的对比示意图;
图2c是分解所得辐射冷负荷的估计值与实际值的对比示意图;
图2d是分解所得新风冷负荷的估计值与实际值的对比示意图;
图2e是分解所得内扰冷负荷的估计值与实际值的对比示意图。
具体实施方式
下面结合实施例和附图对本发明的一种基于信号稀疏表示的建筑分项冷负荷获取方法做出详细说明。
本发明的一种基于信号稀疏表示的建筑分项冷负荷获取方法,是由非负字典学习与非负稀疏表示相关算法所实现的一种将建筑总冷负荷信号分解为各个分项冷负荷信号的方法。该方法的核心思想是,首先采用字典学习方法分别为每项分项冷负荷找到一个可以特定表现其形态特征的冗余字典;而后应用已得到的字典,将总冷负荷信号写成四个字典矩阵与相应稀疏系数矩阵的线性组合,求得四个稀疏系数矩阵的最优解;最后分别将四个字典矩阵与四个稀疏系数矩阵相乘,求得的四个矩阵即为四项分项冷负荷的估计值。本发明的方法由字典学习,总冷负荷信号分解,分项冷负荷信号重组三个步骤组成。
如图1所示,本发明的一种基于信号稀疏表示的建筑分项冷负荷获取方法,包括如下步骤:
1)字典学习,包括:
(1)将总冷负荷分为四项分项冷负荷,包括:温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷,将四项分项冷负荷的逐时历史数据分别写作矩阵形式;具体是:
将温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷四项分项冷负荷的逐时历史数据分别写作矩阵X1、X2、X3和X4的形式,其中X1、X2、X3和X4分别对应表示温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷的历史数据,且X1、X2、X3和X4矩阵的行数为24,表示一天24小时,列数为历史数据所包含的天数。
(2)对所述的矩阵形式分别采用非负K-SVD字典学习算法,得到四个字典矩阵;具体是:
分别对矩阵X1、X2、X3和X4采用非负K-SVD字典学习算法(可详见[1]),得到四个字典矩阵B1、B2、B3和B4,其中,四个字典矩阵B1、B2、B3和B4分别对应表示温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷的字典,且四个字典矩阵B1、B2、B3和B4矩阵的行数为24,表示一天24小时,列数为大于24且小于矩阵X1、X2、X3和X4列数的数值。
2)将总冷负荷分解,是利用稀疏表示相关算法,由总冷负荷和字典矩阵得到四个稀疏系数矩阵;包括:
(1)将用于分解的总冷负荷的逐时数据写作矩阵的形式,总冷负荷矩阵的行数为24,表示一天24小时,列数分别为总冷负荷逐时数据所包含的天数;
(2)设置矩阵A'1、A'2、A'3和A'4为任意非负矩阵,其中矩阵A'1、A'2、A'3和A'4的行数等于字典矩阵B1、B2、B3和B4的列数,矩阵A'1、A'2、A'3和A'4的列数等于总冷负荷矩阵的列数,矩阵A'1、A'2、A'3和A'4分别为温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷的稀疏系数矩阵;
(3)设置稀疏度λ为0至1之间的任意数值;
(4)循环迭代以下各步直至收敛
将A'1、A'2、A'3和A'4矩阵中的所有负数项赋值为0;
(5)得到各分项冷负荷的稀疏系数矩阵A'1、A'2、A'3和A'4
3)分项冷负荷信号重组,是利用矩阵相乘的方法,将步骤1)和步骤2)求得的对应于温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷的四个字典矩阵及四个稀疏系数矩阵对应相乘,获取由总冷负荷分解所得的温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷。包括:
(1)分别将四个字典矩阵B1、B2、B3和B4与四个稀疏系数矩阵A'1、A'2、A'3和A'4对应相乘;
(2)乘积所得的四个矩阵写作其中分别表示分解所得的温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷的估计值,的行数和列数与用于分解的总冷负荷矩阵相同,即,行数为24,列数分别为总冷负荷逐时数据所包含的天数,且每一列代表一天中24小时的逐时冷负荷。
下面给出一使用本发明方法的实例:
(1)应用Trnsys模拟软件对一栋6层居住建筑在1990年至2010年天津的气象条件下的冷负荷进行逐时模拟,得到相应总冷负荷及各项分项冷负荷数据。
(2)将建筑从1990年至2009年共20年的四项分项冷负荷数据作为历史数据,并将四组逐时历史数据分别写作矩阵X1、X2、X3和X4的形式。其中X1、X2、X3和X4分别代表温差冷负荷、辐射冷负荷、新风冷负荷、内扰冷负荷的历史数据,且X1、X2、X3和X4矩阵的行数和列数分别为24和1000。
(3)分别对X1、X2、X3和X4采用非负K-SVD字典学习算法,得到四个字典矩阵B1、B2、B3和B4。其中,B1、B2、B3和B4分别代表温差冷负荷、辐射冷负荷、新风冷负荷、内扰冷负荷的字典,且B1、B2、B3和B4矩阵的行数和列数分别为24和150。
(4)将建筑在2010年夏季所有工作日的逐时总冷负荷数据作为用于分解的逐时数据。将该逐时总冷负荷数据写作矩阵的形式。且矩阵的行数和列数分别为24和40。
(5)设置矩阵A'1、A'2、A'3和A'4为任意非负矩阵,其中矩阵A'1、A'2、A'3和A'4分别代表温差冷负荷、辐射冷负荷、新风冷负荷、内扰冷负荷的稀疏系数矩阵,且矩阵A'1、A'2、A'3和A'4的行数和列数分别为150和40。设置稀疏度λ为0.5。
(6)循环迭代以下各步直至收敛,以获得最优的稀疏系数矩阵A'1、A'2、A'3和A'4
将A'1、A'2、A'3和A'4矩阵中的所有负数项赋值为0。
(7)分别将字典矩阵B1、B2、B3和B4与相对应的稀疏系数矩阵A'1、A'2、A'3和A'4相乘。将乘积所得的四个矩阵写作其中分别代表分解所得的温差冷负荷、辐射冷负荷、新风冷负荷、内扰冷负荷的估计值。的行数和列数分别为24和40,且每一列代表一天中24小时的逐时冷负荷,列数代表所得分解结果的天数为40天。
(8)将模拟所得建筑在2010年夏季所有工作日的四项分项冷负荷数据作为各分项冷负荷的实际值,并将其与所得分解结果中四项分项冷负荷数据的估计值相比较,对比结果如图2所示。
由对比结果图2可以看出,四项分项冷负荷数据的估计值与实际值在形态与幅值上非常相近,因此该分解结果也可以证明该发明所提供方法的可行性与准确性。

Claims (5)

1.一种基于信号稀疏表示的建筑分项冷负荷获取方法,其特征在于,包括如下步骤:
1)字典学习,包括:
(1)将总冷负荷分为四项分项冷负荷,包括:温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷,将四项分项冷负荷的逐时历史数据分别写作矩阵形式;
(2)对所述的矩阵形式分别采用非负K-SVD字典学习算法,得到四个字典矩阵;
2)将总冷负荷分解,是利用稀疏表示相关算法,由总冷负荷和字典矩阵得到四个稀疏系数矩阵;
3)分项冷负荷信号重组,是利用矩阵相乘的方法,将步骤1)和步骤2)求得的对应于温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷的四个字典矩阵及四个稀疏系数矩阵对应相乘,获取由总冷负荷分解所得的温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷。
2.根据权利要求1所述的一种基于信号稀疏表示的建筑分项冷负荷获取方法,其特征在于,步骤1)第(1)步具体是:将温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷四项分项冷负荷的逐时历史数据分别写作矩阵X1、X2、X3和X4的形式,其中X1、X2、X3和X4分别对应表示温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷的历史数据,且X1、X2、X3和X4矩阵的行数为24,表示一天24小时,列数为历史数据所包含的天数。
3.根据权利要求1所述的一种基于信号稀疏表示的建筑分项冷负荷获取方法,其特征在于,步骤1)第(2)步具体是:分别对矩阵X1、X2、X3和X4采用非负K-SVD字典学习算法,得到四个字典矩阵B1、B2、B3和B4,其中,四个字典矩阵B1、B2、B3和B4分别对应表示温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷的字典,且四个字典矩阵B1、B2、B3和B4矩阵的行数为24,表示一天24小时,列数为大于24且小于矩阵X1、X2、X3和X4列数的数值。
4.根据权利要求1所述的一种基于信号稀疏表示的建筑分项冷负荷获取方法,其特征在于,步骤2)包括:
(1)将用于分解的总冷负荷的逐时数据写作矩阵的形式,总冷负荷矩阵的行数为24,表示一天24小时,列数分别为总冷负荷逐时数据所包含的天数;
(2)设置矩阵A'1、A'2、A'3和A'4为任意非负矩阵,其中矩阵A'1、A'2、A'3和A'4的行数等于字典矩阵B1、B2、B3和B4的列数,矩阵A'1、A'2、A'3和A'4的列数等于总冷负荷矩阵的列数,矩阵A'1、A'2、A'3和A'4分别为温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷的稀疏系数矩阵;
(3)设置稀疏度λ为0至1之间的任意数值;
(4)循环迭代以下各步直至收敛
<mrow> <msup> <msub> <mi>A</mi> <mn>1</mn> </msub> <mo>&amp;prime;</mo> </msup> <mo>&amp;LeftArrow;</mo> <msup> <msub> <mi>A</mi> <mn>1</mn> </msub> <mo>&amp;prime;</mo> </msup> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <msup> <msub> <mi>B</mi> <mn>1</mn> </msub> <mi>T</mi> </msup> <mrow> <mo>(</mo> <msup> <mover> <mi>X</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mo>-</mo> <msub> <mi>B</mi> <mn>2</mn> </msub> <msup> <msub> <mi>A</mi> <mn>2</mn> </msub> <mo>&amp;prime;</mo> </msup> <mo>-</mo> <msub> <mi>B</mi> <mn>3</mn> </msub> <msup> <msub> <mi>A</mi> <mn>3</mn> </msub> <mo>&amp;prime;</mo> </msup> <mo>-</mo> <msub> <mi>B</mi> <mn>4</mn> </msub> <msup> <msub> <mi>A</mi> <mn>4</mn> </msub> <mo>&amp;prime;</mo> </msup> <mo>)</mo> </mrow> </mrow> <mrow> <msup> <msub> <mi>B</mi> <mn>1</mn> </msub> <mi>T</mi> </msup> <msub> <mi>B</mi> <mn>1</mn> </msub> <msup> <msub> <mi>A</mi> <mn>1</mn> </msub> <mo>&amp;prime;</mo> </msup> <mo>+</mo> <mi>&amp;lambda;</mi> </mrow> </mfrac> </mrow>
<mrow> <msup> <msub> <mi>A</mi> <mn>2</mn> </msub> <mo>&amp;prime;</mo> </msup> <mo>&amp;LeftArrow;</mo> <msup> <msub> <mi>A</mi> <mn>2</mn> </msub> <mo>&amp;prime;</mo> </msup> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <msup> <msub> <mi>B</mi> <mn>2</mn> </msub> <mi>T</mi> </msup> <mrow> <mo>(</mo> <msup> <mover> <mi>X</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mo>-</mo> <msub> <mi>B</mi> <mn>1</mn> </msub> <msup> <msub> <mi>A</mi> <mn>1</mn> </msub> <mo>&amp;prime;</mo> </msup> <mo>-</mo> <msub> <mi>B</mi> <mn>3</mn> </msub> <msup> <msub> <mi>A</mi> <mn>3</mn> </msub> <mo>&amp;prime;</mo> </msup> <mo>-</mo> <msub> <mi>B</mi> <mn>4</mn> </msub> <msup> <msub> <mi>A</mi> <mn>4</mn> </msub> <mo>&amp;prime;</mo> </msup> <mo>)</mo> </mrow> </mrow> <mrow> <msup> <msub> <mi>B</mi> <mn>2</mn> </msub> <mi>T</mi> </msup> <msub> <mi>B</mi> <mn>2</mn> </msub> <msup> <msub> <mi>A</mi> <mn>2</mn> </msub> <mo>&amp;prime;</mo> </msup> <mo>+</mo> <mi>&amp;lambda;</mi> </mrow> </mfrac> </mrow>
<mrow> <msup> <msub> <mi>A</mi> <mn>3</mn> </msub> <mo>&amp;prime;</mo> </msup> <mo>&amp;LeftArrow;</mo> <msup> <msub> <mi>A</mi> <mn>3</mn> </msub> <mo>&amp;prime;</mo> </msup> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <msup> <msub> <mi>B</mi> <mn>3</mn> </msub> <mi>T</mi> </msup> <mrow> <mo>(</mo> <msup> <mover> <mi>X</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mo>-</mo> <msub> <mi>B</mi> <mn>1</mn> </msub> <msup> <msub> <mi>A</mi> <mn>1</mn> </msub> <mo>&amp;prime;</mo> </msup> <mo>-</mo> <msub> <mi>B</mi> <mn>2</mn> </msub> <msup> <msub> <mi>A</mi> <mn>2</mn> </msub> <mo>&amp;prime;</mo> </msup> <mo>-</mo> <msub> <mi>B</mi> <mn>4</mn> </msub> <msup> <msub> <mi>A</mi> <mn>4</mn> </msub> <mo>&amp;prime;</mo> </msup> <mo>)</mo> </mrow> </mrow> <mrow> <msup> <msub> <mi>B</mi> <mn>3</mn> </msub> <mi>T</mi> </msup> <msub> <mi>B</mi> <mn>3</mn> </msub> <msup> <msub> <mi>A</mi> <mn>3</mn> </msub> <mo>&amp;prime;</mo> </msup> <mo>+</mo> <mi>&amp;lambda;</mi> </mrow> </mfrac> </mrow>
<mrow> <msup> <msub> <mi>A</mi> <mn>4</mn> </msub> <mo>&amp;prime;</mo> </msup> <mo>&amp;LeftArrow;</mo> <msup> <msub> <mi>A</mi> <mn>4</mn> </msub> <mo>&amp;prime;</mo> </msup> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <msup> <msub> <mi>B</mi> <mn>4</mn> </msub> <mi>T</mi> </msup> <mrow> <mo>(</mo> <msup> <mover> <mi>X</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;prime;</mo> </msup> <mo>-</mo> <msub> <mi>B</mi> <mn>2</mn> </msub> <msup> <msub> <mi>A</mi> <mn>2</mn> </msub> <mo>&amp;prime;</mo> </msup> <mo>-</mo> <msub> <mi>B</mi> <mn>3</mn> </msub> <msup> <msub> <mi>A</mi> <mn>3</mn> </msub> <mo>&amp;prime;</mo> </msup> <mo>-</mo> <msub> <mi>B</mi> <mn>4</mn> </msub> <msup> <msub> <mi>A</mi> <mn>4</mn> </msub> <mo>&amp;prime;</mo> </msup> <mo>)</mo> </mrow> </mrow> <mrow> <msup> <msub> <mi>B</mi> <mn>4</mn> </msub> <mi>T</mi> </msup> <msub> <mi>B</mi> <mn>4</mn> </msub> <msup> <msub> <mi>A</mi> <mn>4</mn> </msub> <mo>&amp;prime;</mo> </msup> <mo>+</mo> <mi>&amp;lambda;</mi> </mrow> </mfrac> </mrow>
将A'1、A'2、A'3和A'4矩阵中的所有负数项赋值为0;
(5)得到各分项冷负荷的稀疏系数矩阵A'1、A'2、A'3和A'4
5.根据权利要求1所述的一种基于信号稀疏表示的建筑分项冷负荷获取方法,其特征在于,步骤3)包括:
(1)分别将四个字典矩阵B1、B2、B3和B4与四个稀疏系数矩阵A'1、A'2、A'3和A'4对应相乘;
(2)乘积所得的四个矩阵写作其中分别表示分解所得的温差冷负荷、辐射冷负荷、新风冷负荷和内扰冷负荷的估计值,的行数和列数与用于分解的总冷负荷矩阵相同,即,行数为24,列数分别为总冷负荷逐时数据所包含的天数,且每一列代表一天中24小时的逐时冷负荷。
CN201711027913.2A 2017-10-28 2017-10-28 一种基于信号稀疏表示的建筑分项冷负荷获取方法 Active CN107895203B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711027913.2A CN107895203B (zh) 2017-10-28 2017-10-28 一种基于信号稀疏表示的建筑分项冷负荷获取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711027913.2A CN107895203B (zh) 2017-10-28 2017-10-28 一种基于信号稀疏表示的建筑分项冷负荷获取方法

Publications (2)

Publication Number Publication Date
CN107895203A true CN107895203A (zh) 2018-04-10
CN107895203B CN107895203B (zh) 2021-06-25

Family

ID=61803328

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711027913.2A Active CN107895203B (zh) 2017-10-28 2017-10-28 一种基于信号稀疏表示的建筑分项冷负荷获取方法

Country Status (1)

Country Link
CN (1) CN107895203B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113191418A (zh) * 2021-04-27 2021-07-30 华中科技大学 基于室外气象参数的非侵入式建筑分项冷负荷监测方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100256958A1 (en) * 2007-11-12 2010-10-07 The Industry & Academic Cooperation In Chungnam National University Method for predicting cooling load
CN102073999A (zh) * 2011-01-20 2011-05-25 西安电子科技大学 基于双冗余字典学习的自然图像去噪方法
CN104573738A (zh) * 2013-10-28 2015-04-29 北京大学 信号处理方法及其装置
CN107169606A (zh) * 2017-05-18 2017-09-15 天津大学 一种办公建筑冷负荷的预测方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100256958A1 (en) * 2007-11-12 2010-10-07 The Industry & Academic Cooperation In Chungnam National University Method for predicting cooling load
CN102073999A (zh) * 2011-01-20 2011-05-25 西安电子科技大学 基于双冗余字典学习的自然图像去噪方法
CN104573738A (zh) * 2013-10-28 2015-04-29 北京大学 信号处理方法及其装置
CN107169606A (zh) * 2017-05-18 2017-09-15 天津大学 一种办公建筑冷负荷的预测方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113191418A (zh) * 2021-04-27 2021-07-30 华中科技大学 基于室外气象参数的非侵入式建筑分项冷负荷监测方法

Also Published As

Publication number Publication date
CN107895203B (zh) 2021-06-25

Similar Documents

Publication Publication Date Title
CN108694467B (zh) 一种对配电网线损率进行预测的方法及系统
CN102622530B (zh) 基于改进遗传算法的抽汽供热机组热电负荷分配优化方法
Lu et al. Bi-population based discrete bat algorithm for the low-carbon job shop scheduling problem
Geng et al. Multi‐Objective Reentrant Hybrid Flowshop Scheduling with Machines Turning on and off Control Strategy Using Improved Multi‐Verse Optimizer Algorithm
CN102819641A (zh) 适于电磁暂态仿真的大规模配电网络整体模型化简方法
CN103049617B (zh) 保留无源性的大规模配电网络电磁暂态仿真模型化简方法
CN110264012A (zh) 基于经验模态分解的可再生能源功率组合预测方法及系统
Xiangkun et al. Synthesis of multi-stream heat exchanger network for multi-period operation with genetic/simulated annealing algorithms
Yuan et al. Irregular distribution of wind power prediction
Van Schijndel A review of the application of SimuLink S-functions to multi domain modelling and building simulation
CN110443404A (zh) 一种基于热量组分解析的建筑空调负荷预测方法
Pan et al. Probabilistic Short‐Term Wind Power Forecasting Using Sparse Bayesian Learning and NWP
CN107895203A (zh) 一种基于信号稀疏表示的建筑分项冷负荷获取方法
CN105024645A (zh) 一种基于矩阵进化的光伏阵列故障定位方法
Wei et al. The optimal reactive power flow model in mixed polar form based on transformer dummy nodes
CN105939026A (zh) 基于混合Laplace分布的风电功率波动量概率分布模型建立方法
CN116542385A (zh) 一种针对风、光功率预测的集成学习方法及系统
CN112784439A (zh) 一种基于离散化模型的能源互联网规划方法及装置
CN110020506B (zh) 基于电热型综合能源系统运行优化的差分格式选择方法
Lu et al. Wind farm layout design optimization through multi-scenario decomposition with complementarity constraints
Lee et al. On-policy learning-based deep reinforcement learning assessment for building control efficiency and stability
Xiangkun et al. Synthesis of flexible multi-stream heat exchanger networks based on stream pseudo-temperature with genetic/simulated annealing algorithms
Odonkor et al. Optimization of energy use strategies in building clusters using pareto bands
Ke et al. Optimizing and controlling building electric energy using cat boost under the energy internet of things
Ye et al. A combined prediction method of wind farm power

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