CN104463381B - 基于kpca与wlssvm的建筑能耗预测方法 - Google Patents
基于kpca与wlssvm的建筑能耗预测方法 Download PDFInfo
- Publication number
- CN104463381B CN104463381B CN201510000613.XA CN201510000613A CN104463381B CN 104463381 B CN104463381 B CN 104463381B CN 201510000613 A CN201510000613 A CN 201510000613A CN 104463381 B CN104463381 B CN 104463381B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msup
- mtd
- munderover
- 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
Links
- 238000005265 energy consumption Methods 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000005457 optimization Methods 0.000 claims abstract description 17
- 238000004458 analytical method Methods 0.000 claims abstract description 6
- 239000002245 particle Substances 0.000 claims description 25
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000013507 mapping Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000012549 training Methods 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 238000010998 test method Methods 0.000 claims description 3
- 230000007786 learning performance Effects 0.000 abstract 1
- 238000004134 energy conservation Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 201000004569 Blindness Diseases 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005183 dynamical system Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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/08—Construction
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Tourism & Hospitality (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Medicines That Contain Protein Lipid Enzymes And Other Medicines (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种基于KPCA与WLSSVM的建筑能耗预测方法,该方法利用核主元分析(KPCA)消除样本共线性,降低维数,进而建立加权最小二乘支持向量(WLSSVM)模型,并结合粒子群(PSO)算法优化模型参数,提高模型的学习性能和泛化能力。本发明可有效的应用于建筑能耗预测,具有良好的预测精度和鲁棒性能。
Description
技术领域
本发明涉及一种基于KPCA与WLSSVM的建筑能耗预测方法。
背景技术
建筑能耗在我国能源消耗中的比重日益增大,建筑节能已经迫在眉睫。建筑能耗的准确预测不仅是实现建筑能耗分析和节能评估的前提条件,也是进行建筑节能优化设计的重要依据之一。
建筑能耗***作为一个典型的具有多变量、强耦合和不确定特性的复杂动态***,建立其精确的预测模型仍然面临很大的困难。一方面,由于建筑室内***或设备启停等情况的存在,使得用于建模的能耗数据样本中不可避免包含异常值,这势必影响模型预测结果的可靠性。另一方面,由于建筑能耗和环境温度、气象条件、建筑围护结构以及建筑节能参数等多种影响因素之间存在着复杂的非线性关系,而且各因素之间存在很强的相关性,在建模过程中会降低建模精度,增加计算复杂度。
发明内容
本发明的目的在于提供一种基于KPCA与WLSSVM的建筑能耗预测方法,能够有效消除变量之前的相关性,降低输入变量的维数,简化模型结构,提高模型训练速度,具有较高预测精度和泛化能力。
为实现上述目的,本发明采用如下技术方案:一种基于KPCA与WLSSVM的建筑能耗预测模型方法,其特征在于,包括以下步骤:
步骤S1:通过DEST-C动态计算采用正交试验方法,获取建筑物全年能耗数据样本集,讲所述建筑物全年能耗数据样本集进行归一化处理后作为建筑能耗预测模型数据;
步骤S2:通过核主元分析消除变量之间的相关性,并通过提取主成分降低样本维数;
步骤S3:采用粒子群优化算法对所述建筑能耗预测模型的参数进行优化;
步骤S4:求解b、α可建立LSSVM的建筑能耗预测模型,再根据模型重新计算各个样本的拟合误差ξi,各个样本的权值vi以及各个建筑能耗样本的权值vi,求解b、α*建立WLSSVM的建筑能耗预测模型;
步骤S5:输入新的建筑能耗样本,得到模型输出将数据进行反归一化。
在本发明一实施例中,所述步骤S2的具体方法如下:
将原空间xi(i=1,2,…,N)映射到某个高维特征空间,并在所述高维特征空间实现PCA,所述高维特征空间的协方差矩阵R为:
求所述协方差矩阵R的特征值λi和特征向量ui,有:
Rui=λiui (2)
所述特征向量ui可以表示为φ(x)的线性组合:
以上三式两边左乘φ(xk),整理得:
式中,k=1,2,…,N,定义一个N×N维的核函数矩阵K:
Kij=K(xi,xj)=φ(xi)T·φ(xj) (5)
式(4)进一步化简为:
Nλia=Ka (6)
式中,a=(a1,a2,…,aN)T,提取样本数据x映射后的第t非线性主元pt:
在本发明一实施例中,所述步骤S3的具体方法如下:
定义LSSVM的核函数为高斯径向基函数:
式中,σ为核宽参数;
再采用粒子群优化算法来优化LSSVM的正则化参数C和核宽参数σ的值,优化关键步骤如下:
步骤S31:建立待优化目标函数:
式中,yi为第i个样本实际值,为第i个样本预测值,设定约束集C∈(Cmin,Cmax)和σ∈(σmin,σmax);
步骤S32:初始化粒子种群,将参数(C,σ)的值表示为粒子的位置;
步骤S33:计算所述粒子的适应度值,以更新所述粒子的位置和速度;
步骤S34:评价所述粒子种群是否满足最优解条件,如果满足则输出最优参数(Cbest,σbest),否则返回步骤S33再次计算粒子适应度。
在本发明一实施例中,所述步骤S4的具体方法如下:
通过WLSSVM将样本误差的二范数定义为损失函数,并将不等式约束改为等式约束,同时,依据样本训练的重要性分别赋予其不同权重,设所述二范数的权值为vi,则其优化问题为:
s.t. yi=ωT·φ(xi)+b+ξi (11)
式中,ω是权系数向量,φ(·)是输入空间到高维空间的映射,C是惩罚因子,b是阈值;
引入Lagrange函数:
式中,为Lagrange乘子,根据优化条件,引入核函数:
K(xi,xj)=φ(xi)T·φ(xj) (13)
则式(10)的优化问题变为以下线性方程组的求解问题:
式中,l1×N是1×N的单位行向量,lN×1是N×1的单位列向量,
y=[y1,y2,…,yN]T;由此可得回归函数的形式:
本发明与现有技术相比具有以下有益效果:
1、通过核主元分析能够消除建筑能耗影响因素间的冗余性和共线性,正确提取输入样本的特征信息,简化模型的输入样本,提高建模效率;
2、WLSSVM模型具有良好的非线性处理能力,使得基于WLSSVM的建筑能耗预测模型具有更好的学习精度和泛化能力;
3、采用PSO优化算法对模型的参数进行寻优,避免了模型参数选择的盲目性,进一步提高建筑能耗预测模型的预测精度。
附图说明
图1是本发明方法流程图。
具体实施方式
下面结合附图及实施例对本发明做进一步说明。
请参照图1,本发明提供一种基于KPCA与WLSSVM的建筑能耗预测方法,其特征在于,包括以下步骤:
步骤S1:通过DEST-C动态计算采用正交试验方法,获取建筑物全年能耗数据样本集,讲所述建筑物全年能耗数据样本集进行归一化处理后作为建筑能耗预测模型数据;
步骤S2:通过核主元分析消除变量之间的相关性,并通过提取主成分降低样本维数;具体方法如下:
将原空间xi(i=1,2,…,N)映射到某个高维特征空间,并在所述高维特征空间实现PCA,所述高维特征空间的协方差矩阵R为:
求所述协方差矩阵R的特征值λi和特征向量ui,有:
Rui=λiui (2)
所述特征向量ui可以表示为φ(x)的线性组合:
以上三式两边左乘φ(xk),整理得:
式中,k=1,2,…,N,定义一个N×N维的核函数矩阵K:
Kij=K(xi,xj)=φ(xi)T·φ(xj) (5)
式(4)进一步化简为:
Nλia=Ka (6)
式中,a=(a1,a2,…,aN)T,提取样本数据x映射后的第t非线性主元pt:
步骤S3:采用粒子群优化算法对所述建筑能耗预测模型的参数进行优化;具体方法如下:
在核函数选择方面,高斯径向基函数具有良好的处理样本输入与输出之间复杂非线性关系的能力,而且需要确定的参数少,计算效率高,因此定义LSSVM的核函数为高斯径向基函数:
式中,σ为核宽参数;
对于C和σ的选择,一般是依靠试算和经验的方法,不仅耗时且不准确,粒子群优化算法是一种通用的启发式搜索技术,该算法具有快速和全局优化的特点,因此,本发明采用粒子群优化算法来优化LSSVM的正则化参数C和核宽参数σ的值,优化关键步骤如下:
步骤S31:建立待优化目标函数:
式中,yi为第i个样本实际值,为第i个样本预测值,设定约束集C∈(Cmin,Cmax)和σ∈(σmin,σmax);
步骤S32:初始化粒子种群,将参数(C,σ)的值表示为粒子的位置;
步骤S33:计算所述粒子的适应度值,以更新所述粒子的位置和速度;
步骤S34:评价所述粒子种群是否满足最优解条件,如果满足则输出最优参数(Cbest,σbest),否则返回步骤S33再次计算粒子适应度。
步骤S4:求解b、α可建立LSSVM的建筑能耗预测模型,再根据模型重新计算各个样本的拟合误差ξi,各个样本的权值vi以及各个建筑能耗样本的权值vi,求解b、α*建立WLSSVM的建筑能耗预测模型;具体方法如下:
通过WLSSVM将样本误差的二范数定义为损失函数,并将不等式约束改为等式约束,同时,依据样本训练的重要性分别赋予其不同权重,设所述二范数的权值为vi,则其优化问题为:
s.t. yi=ωT·φ(xi)+b+ξi (11)
式中,ω是权系数向量,φ(·)是输入空间到高维空间的映射,C是惩罚因子,b是阈值;
引入Lagrange函数:
式中,为Lagrange乘子,根据优化条件,引入核函数:
K(xi,xj)=φ(xi)T·φ(xj) (13)
则式(10)的优化问题变为以下线性方程组的求解问题:
式中,l1×N是1×N的单位行向量,lN×1是N×1的单位列向量,
y=[y1,y2,…,yN]T;由此可得回归函数的形式:
步骤S5:输入新的建筑能耗样本,得到模型输出将数据进行反归一化。
以上所述仅为本发明的较佳实施例,凡依本发明申请专利范围所做的均等变化与修饰,皆应属本发明的涵盖范围。
Claims (1)
1.一种基于KPCA与WLSSVM的建筑能耗预测方法,其特征在于,包括以下步骤:
步骤S1:通过DEST-C动态计算采用正交试验方法,获取建筑物全年能耗数据样本集,讲所述建筑物全年能耗数据样本集进行归一化处理后作为建筑能耗预测模型数据;
步骤S2:通过核主元分析消除变量之间的相关性,并通过提取主成分降低样本维数;
步骤S3:采用粒子群优化算法对所述建筑能耗预测模型的参数进行优化;
步骤S4:求解b、α可建立LSSVM的建筑能耗预测模型,再根据模型重新计算各个样本的拟合误差ξi,各个样本的权值vi以及各个建筑能耗样本的权值vi,求解b、α*建立WLSSVM的建筑能耗预测模型;
步骤S5:输入新的建筑能耗样本,得到模型输出将数据进行反归一化;
所述步骤S2的具体方法如下:
将原空间xi(i=1,2,…,N)映射到某个高维特征空间,并在所述高维特征空间实现PCA,所述高维特征空间的协方差矩阵R为:
<mrow>
<mi>R</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>&phi;</mi>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
求所述协方差矩阵R的特征值λi和特征向量ui,有:
Rui=λiui (2)
所述特征向量ui可以表示为φ(x)的线性组合:
<mrow>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>a</mi>
<mi>j</mi>
</msub>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
以上三式两边左乘φ(xk),整理得:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mo>{</mo>
<mo>&lsqb;</mo>
<mi>&phi;</mi>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>&CenterDot;</mo>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>&CenterDot;</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>a</mi>
<mi>j</mi>
</msub>
<mo>&lsqb;</mo>
<mi>&phi;</mi>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>&CenterDot;</mo>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>a</mi>
<mi>j</mi>
</msub>
<mo>&lsqb;</mo>
<mi>&phi;</mi>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>&CenterDot;</mo>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,k=1,2,…,N,定义一个N×N维的核函数矩阵K:
Kij=K(xi,xj)=φ(xi)T·φ(xj) (5)
式(4)进一步化简为:
Nλia=Ka (6)
式中,a=(a1,a2,…,aN)T,提取样本数据x映射后的第t非线性主元pt:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>p</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<msup>
<msub>
<mi>u</mi>
<mi>t</mi>
</msub>
<mi>T</mi>
</msup>
<mo>&CenterDot;</mo>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>a</mi>
<mi>j</mi>
</msub>
<mi>&phi;</mi>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>&CenterDot;</mo>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>a</mi>
<mi>j</mi>
</msub>
<mi>K</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
所述步骤S3的具体方法如下:
定义LSSVM的核函数为高斯径向基函数:
<mrow>
<mi>K</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<mi>&sigma;</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,σ为核宽参数;
再采用粒子群优化算法来优化LSSVM的正则化参数C和核宽参数σ的值,优化关键步骤如下:
步骤S31:建立待优化目标函数:
<mrow>
<mi>min</mi>
<mi> </mi>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mo>,</mo>
<mi>&sigma;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msubsup>
<mi>y</mi>
<mi>i</mi>
<mo>*</mo>
</msubsup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,yi为第i个样本实际值,为第i个样本预测值,设定约束集C∈(Cmin,Cmax)和σ∈(σmin,σmax);
步骤S32:初始化粒子种群,将参数(C,σ)的值表示为粒子的位置;
步骤S33:计算所述粒子的适应度值,以更新所述粒子的位置和速度;
步骤S34:评价所述粒子种群是否满足最优解条件,如果满足则输出最优参数(Cbest,σbest),否则返回步骤S33再次计算粒子适应度;
所述步骤S4的具体方法如下:
通过WLSSVM将样本误差的二范数定义为损失函数,并将不等式约束改为等式约束,同时,依据样本训练的重要性分别赋予其不同权重,设所述二范数ξi 2的权值为vi,则其优化问题为:
<mrow>
<mi>min</mi>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mi>&omega;</mi>
<mi>T</mi>
</msup>
<mi>&omega;</mi>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mi>C</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>v</mi>
<mi>i</mi>
</msub>
<msubsup>
<mi>&xi;</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
s.t.yi=ωT·φ(xi)+b+ξi (11)
式中,ω是权系数向量,φ(·)是输入空间到高维空间的映射,C是惩罚因子,b是阈值;
引入Lagrange函数:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mo>,</mo>
<mi>b</mi>
<mo>,</mo>
<mi>&xi;</mi>
<mo>,</mo>
<msup>
<mi>&alpha;</mi>
<mo>*</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mi>&omega;</mi>
<mi>T</mi>
</msup>
<mi>&omega;</mi>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mi>C</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>v</mi>
<mi>i</mi>
</msub>
<msubsup>
<mi>&xi;</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>&alpha;</mi>
<mi>i</mi>
<mo>*</mo>
</msubsup>
<mo>&lsqb;</mo>
<msup>
<mi>&omega;</mi>
<mi>T</mi>
</msup>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>b</mi>
<mo>+</mo>
<msub>
<mi>&xi;</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,为Lagrange乘子,根据优化条件,引入核函数:
K(xi,xj)=φ(xi)T·φ(xj) (13)
则式(10)的优化问题变为以下线性方程组的求解问题:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<msub>
<mi>l</mi>
<mrow>
<mn>1</mn>
<mo>&times;</mo>
<mi>N</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>l</mi>
<mrow>
<mi>N</mi>
<mo>&times;</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mrow>
<mi>R</mi>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mi>C</mi>
</mfrac>
<mi>V</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>b</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<msup>
<mi>&alpha;</mi>
<mo>*</mo>
</msup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>y</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,l1×N是1×N的单位行向量,lN×1是N×1的单位列向量,R={K(xi,xj)|i,j=1,2,…,N},y=[y1,y2,…,yN]T;由此可得回归函数的形式:
<mrow>
<mi>y</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>&alpha;</mi>
<mi>i</mi>
<mo>*</mo>
</msubsup>
<mi>K</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>b</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
2
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510000613.XA CN104463381B (zh) | 2015-01-04 | 2015-01-04 | 基于kpca与wlssvm的建筑能耗预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510000613.XA CN104463381B (zh) | 2015-01-04 | 2015-01-04 | 基于kpca与wlssvm的建筑能耗预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104463381A CN104463381A (zh) | 2015-03-25 |
CN104463381B true CN104463381B (zh) | 2017-12-08 |
Family
ID=52909385
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510000613.XA Expired - Fee Related CN104463381B (zh) | 2015-01-04 | 2015-01-04 | 基于kpca与wlssvm的建筑能耗预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104463381B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104807039B (zh) * | 2015-04-23 | 2017-01-11 | 华北电力大学 | 一种锅炉燃烧优化的变量降维建模方法 |
CN106250905B (zh) * | 2016-07-08 | 2020-05-12 | 复旦大学 | 一种结合高校建筑结构特征的实时能耗异常检测方法 |
CN106779137A (zh) * | 2016-11-14 | 2017-05-31 | 上海海事大学 | 一种根据海况与操纵条件预测船舶油耗的方法 |
CN106910144B (zh) * | 2017-01-11 | 2020-04-28 | 华南理工大学 | 基于分时刻实际用能系数的大型建筑逐时能耗在线预测方法 |
CN107730006B (zh) * | 2017-09-13 | 2021-01-05 | 重庆电子工程职业学院 | 基于可再生能源大数据深度学习的建筑近零能耗控制方法 |
CN108204944A (zh) * | 2018-01-13 | 2018-06-26 | 福州大学 | 基于apso优化的lssvm的埋地管道腐蚀速率预测方法 |
CN108764305A (zh) * | 2018-05-14 | 2018-11-06 | 浙江大学 | 一种改进的群智能机器学习故障诊断*** |
CN109284877A (zh) * | 2018-11-19 | 2019-01-29 | 福州大学 | 基于aiga-wlssvm埋地管道腐蚀速率预测方法 |
CN110610274A (zh) * | 2019-09-18 | 2019-12-24 | 福州大学 | 一种建筑能耗预测方法及*** |
CN111523710A (zh) * | 2020-04-10 | 2020-08-11 | 三峡大学 | 基于pso-lssvm在线学习的电力设备温度预测方法 |
CN111832839B (zh) * | 2020-07-24 | 2021-04-30 | 河北工业大学 | 基于充分增量学习的能耗预测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102194018A (zh) * | 2010-03-05 | 2011-09-21 | 朗德华信(北京)自控技术有限公司 | 基于云计算的建筑物能耗预测分析***及方法 |
CN102938092A (zh) * | 2012-10-08 | 2013-02-20 | 珠海派诺科技股份有限公司 | 一种基于神经网络的建筑节假日能耗预测方法 |
-
2015
- 2015-01-04 CN CN201510000613.XA patent/CN104463381B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102194018A (zh) * | 2010-03-05 | 2011-09-21 | 朗德华信(北京)自控技术有限公司 | 基于云计算的建筑物能耗预测分析***及方法 |
CN102938092A (zh) * | 2012-10-08 | 2013-02-20 | 珠海派诺科技股份有限公司 | 一种基于神经网络的建筑节假日能耗预测方法 |
Non-Patent Citations (2)
Title |
---|
Predicting hourly cooling load in the building: A comparison of support vector machine and different artificial neural networks;Qiong Li 等;《Energy Conversion and Management》;20090228;第50卷(第1期);第90-91页 * |
飞灰含碳量自适应校正WLSSVM软测量模型;王伟 等;《热力发电》;20130831;第42卷(第8期);第75-76页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104463381A (zh) | 2015-03-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104463381B (zh) | 基于kpca与wlssvm的建筑能耗预测方法 | |
CN106600059B (zh) | 基于改进rbf神经网络的智能电网短期负荷预测方法 | |
Andrés et al. | Efficient aerodynamic design through evolutionary programming and support vector regression algorithms | |
Han et al. | Surrogate-based aerodynamic shape optimization with application to wind turbine airfoils | |
CN110866592B (zh) | 模型训练方法、装置、能效预测方法、装置和存储介质 | |
Zhou et al. | Comparison of different deep neural network architectures for isothermal indoor airflow prediction | |
CN106650920A (zh) | 一种基于优化极限学习机的预测模型 | |
CN105260786A (zh) | 一种电力推进***仿真可信度评估模型综合优化方法 | |
CN109978253B (zh) | 一种基于增量学习的电力***短期负荷预测方法 | |
CN110781595B (zh) | 能源使用效率pue的预测方法、装置、终端及介质 | |
Kyriacou et al. | Efficient PCA-driven EAs and metamodel-assisted EAs, with applications in turbomachinery | |
CN115310727B (zh) | 一种基于迁移学习的建筑冷热电负荷预测方法及*** | |
CN112381282A (zh) | 基于宽度学习***的光伏发电功率预测方法 | |
CN116796141A (zh) | 一种基于gbdt回归模型的办公建筑能耗预测方法 | |
You et al. | A novel mixed uncertainty support vector machine method for structural reliability analysis | |
Gilan et al. | Sustainable building design: A challenge at the intersection of machine learning and design optimization | |
CN201716727U (zh) | 基于遥感与gis的地理模拟*** | |
Du et al. | Super resolution generative adversarial networks for multi-fidelity pressure distribution prediction | |
Wang et al. | Mutual information-weighted principle components identified from the depth features of stacked autoencoders and original variables for oil dry point soft sensor | |
CN117034808A (zh) | 一种基于图注意力网络的天然气管网压力估计方法 | |
Zhang et al. | RSVRs based on feature extraction: a novel method for prediction of construction projects’ costs | |
Selvan | On the effect of shape parameterization on aerofoil shape optimization | |
Hailong et al. | Data fusion method for temperature monitoring of bio-oxidation with wireless sensor networks | |
CN115034432A (zh) | 一种风电场风力发电机组风速预测方法 | |
Wang et al. | Neural network and PSO-based structural approximation analysis for blade of wind turbine |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171208 Termination date: 20210104 |