CN106600059B - 基于改进rbf神经网络的智能电网短期负荷预测方法 - Google Patents

基于改进rbf神经网络的智能电网短期负荷预测方法 Download PDF

Info

Publication number
CN106600059B
CN106600059B CN201611148874.7A CN201611148874A CN106600059B CN 106600059 B CN106600059 B CN 106600059B CN 201611148874 A CN201611148874 A CN 201611148874A CN 106600059 B CN106600059 B CN 106600059B
Authority
CN
China
Prior art keywords
center
output
neural network
basis function
hidden layer
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.)
Active
Application number
CN201611148874.7A
Other languages
English (en)
Other versions
CN106600059A (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.)
Nanchang University
Beijing University of Posts and Telecommunications
Original Assignee
Nanchang University
Beijing University of Posts and Telecommunications
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 Nanchang University, Beijing University of Posts and Telecommunications filed Critical Nanchang University
Priority to CN201611148874.7A priority Critical patent/CN106600059B/zh
Publication of CN106600059A publication Critical patent/CN106600059A/zh
Application granted granted Critical
Publication of CN106600059B publication Critical patent/CN106600059B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/048Activation functions
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • 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
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Strategic Management (AREA)
  • Human Resources & Organizations (AREA)
  • Software Systems (AREA)
  • Tourism & Hospitality (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Linguistics (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • General Business, Economics & Management (AREA)
  • Molecular Biology (AREA)
  • Artificial Intelligence (AREA)
  • Marketing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • Public Health (AREA)
  • Development Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开一种基于改进RBF神经网络的智能电网短期负荷预测方法,涉及智能电网技术领域,用于确定基函数中心,提升智能电网负荷预测精度。该预测方法包括:网络初始化;S2、计算基函数中心ci;S3、根据基函数中心ci,计算方差ζi;S4、根据基函数中心ci以及方差ζi,计算隐含层的输出Ri;S5、根据隐含层的输出Ri,计算输出层的输出;S6、根据均方误差和函数计算预测误差E;S7、对神经网络中隐含层神经元与输出层神经元的连接权重进行更新;S8、对预测误差E进行判断,若预测误差E在预期内,则迭代计算结束;反之,返回步骤S4,重新迭代计算预测误差E。本发明用于对电网负荷进行预测。

Description

基于改进RBF神经网络的智能电网短期负荷预测方法
技术领域
本发明涉及智能电网技术领域,尤其涉及一种基于改进RBF神经网络的智能电网短期负荷预测方法。
背景技术
智能电网的快速发展产生了大量的用电数据(亦称为样本数据),对这些样本数据进行分析具有重大意义。利用预测方法将样本数据应用到短期负荷预测中,从而提升负荷预测精度,这对电力***的安全调度和经济运行起着重要的作用。径向基函数(RadialBasis Function,以下简称RBF)神经网络是应用于负荷预测中最为广泛的一种预测方法,因为它是一种局部逼近网络,能以任意精度逼近任意连续函数,具有唯一最佳逼近特性且无局部极小问题,而且拓扑结构简单、学习速率快。RBF神经网络预测方法中主要有三个参数影响预测精度,分别为基函数中心、基函数半径以及网络隐含层与输出层的连接权值。其中连接权值常采用梯度下降法求得。基函数中心以及基函数半径对预测精度的影响非常大,因此现有研究主要集中在如何确定RBF神经网络的基函数中心以及基函数半径上。现有技术中主要采用如下方法计算基函数中心以及基函数半径:
第一种采用聚类分析法(例如,K-means方法和FCM方法)计算基函数中心以及基函数半径;第二种采用启发式方法(例如,如遗传方法和粒子群方法)计算基函数中心以及基函数半径。其中,启发式方法复杂度高,在智能电网大规模负荷数据下预测时间较长,因此在智能电网大规模负荷预测下,聚类法更适合确定RBF神经网络基函数中心以及基函数半径。
另外,对于RBF神经网络基函数中心,主要采用FCM方法确定,但是在智能电网负荷预测中,FCM方法负荷数据规模大、维数多,方法较复杂,最终使得智能电网负荷预测精确度较低。
发明内容
本发明的目的在于提供一种基于改进RBF神经网络的智能电网短期负荷预测方法,用于确定基函数中心,提升智能电网负荷预测精度。
为达到上述目的,本发明采用如下技术方案:
本发明提供一种基于改进RBF神经网络的智能电网短期负荷预测方法,基于改进RBF神经网络的智能电网短期负荷预测方法包括:
S1、网络初始化;
S2、计算基函数中心ci
S3、根据基函数中心ci,计算方差ζi
S4、根据基函数中心ci以及方差ζi,计算隐含层的输出Ri
S5、根据隐含层的输出Ri,计算输出层的输出;
S6、根据均方误差和函数计算预测误差E;
S7、对神经网络中隐含层神经元与输出层神经元的连接权重进行更新;
S8、对预测误差E进行判断,若预测误差E在预期内,则迭代计算结束;反之,返回步骤S4,重新迭代计算预测误差E。
步骤S1包括:确定输入层神经元个数N,隐含层神经元个数N,输出层神经元个数N,以及初始化学习速率η和基函数重叠系数η;其中,输出层神经元个数N=1,隐含层神经元个数N即为基函数中心个数。
步骤S2包括:
S21、S21、输入模糊指数m,迭代停止阈值ε以及PCA累积贡献率因子δ,基函数中心个数N,聚类的属性权值ωn,原始数据X;;其中,样本数据X={x1,x2,…,xN},N为样本点的个数,xj为每一个样本点,xj={xj1,xj2,…,xjs},j=1,2,…N,s表示每一个样本点包含的属性的个数,即每一个样本点的维度,在PCA-WFCM方法中将样本数据X划分为K类,聚类中心V={v1,v2,…,vk},共有K个,K的范围为2≤K≤N。
S22、对样本数据X进行PCA属性降维处理:根据如下公式计算出每一个样本点降维后的维度L
Figure BDA0001179393020000031
保留维度L之前的所有维度作为聚类属性,式中,S表示每一个样本点的原始维度,λn表示PCA中协方差矩阵的特征值,δ表示PCA累积贡献率因子;根据如下公式计算每一个聚类的属性权值:
Figure BDA0001179393020000032
n=1,2,…,L;根据聚类的属性权值ωn得到聚类属性的权重向量W={ωn}以及降维后的样本数据Xnew;。
S23、初始化隶属度矩阵U:令0≤uij≤1,
Figure BDA0001179393020000033
U={uij},式中,uij表示第j个样本点属于第i类的隶属度,K表示聚类中心V的个数。
S24、根据隶属度矩阵U以及降维后的样本数据Xnew,计算聚类中心V:
Figure BDA0001179393020000034
V={vi},式中,m为模糊指数,xj表示第j个样本点。
S25、根据降维后的样本数据Xnew,迭代计算隶属度矩阵U:
Figure BDA0001179393020000035
U={uij},式中,m为模糊指数,表示隶属度矩阵U的模糊度,m的值越大,则隶属度矩阵U的模糊度越高,令m=2,K表示聚类中心V的个数,dij表示每一个样本点xj到聚类中心vi的加权欧式距离,由如下公式计算得出:,L表示每一个样本点降维后的维度,ωn表示聚类的属性权值。
S26、根据隶属度矩阵U和聚类中心V计算目标函数J:
Figure BDA0001179393020000036
式中,m为模糊指数,dij表示每一个样本点xj到聚类中心vi的加权欧式距离,K表示聚类中心V的个数,N表示样本点个数。
S27、对目标函数J进行判断:若|J(t)-J(t-1)|<ε,则输出聚类中心V,也就是基函数中心ci;反之,则返回步骤S24,直到满足公式|J(t)-J(t-1)|<ε,则停止迭代计算,计算出聚类中心V;式中,ε表示迭代停止阈值,t表示迭代次数。
步骤S3包括:根据如下公式计算隐含层神经元的半径ζi:ζi=λminj‖ci-cj‖,i,j=1,2,…N;式中,ci表示基函数中心,N表示隐含层神经元个数,η表示基函数重叠系数。
步骤S4包括:根据降维后的样本数据Xnew、基函数中心ci以及方差ζi,得到隐含层的输出:
Figure BDA0001179393020000041
i=1,2,…N,式中,N表示隐含层神经元个数。
步骤S5包括:根据降维后的样本数据Xnew、得到输出层的输出:
Figure BDA0001179393020000042
式中,N表示隐含层神经元个数;W表示第i个隐含层神经元与输出层神经元的连接权重,Ri表示隐含层的输出。
步骤S6包括:利用均方误差和函数计算预测误差E:令将一组输入向量{xj,j=1,2…O}和对应输出值{yj,j=1,2,…O}作为训练样本,
式中,O为样本数,预测误差
Figure BDA0001179393020000043
步骤S7包括:根据如下公式对连接权重W进行更新:
Figure BDA0001179393020000044
式中,η表示学习速率,E为预测误差,q为更新次数。
聚类加权值Z的计算方法如下:
计算出协方差矩阵C:令一个S维空间数据映射到L维子空间中,其中,L<<S,假设X={xn}为零均值数据,即
Figure BDA0001179393020000045
其中,n=1,2,…,N,
Figure BDA0001179393020000046
T为转秩符号;对协方差矩阵C进行特征值分解:
Figure BDA0001179393020000047
将一个S维的数据xi向L维主成分方向投影,即,Y=XQL;设协方差矩阵C的特征根λ1≥λ2≥…≥λS,且
Figure BDA0001179393020000048
为第l个主成分的贡献率,
Figure BDA0001179393020000051
为前L个主成分的累积贡献率;降维后第k个属性的聚类属性权值:
Figure BDA0001179393020000052
l=1,2…,L;特征向量合集Q=[q1,q2,…,qS],特征值Л=diag(λ12,…,λS),令累积贡献率大于95%。
采用PCA-WFCM确定RBF基函数中心,负荷数据先通过PCA降维处理,获得较少的不相关的预测输入,从而减小RBF基函数的重叠,更好的确定基函数中心并降低预测算法复杂度;另外,采用加权FCM基函数聚类,将PCA处理后获得的不同属性的方差贡献率给降维后属性加权,从而提升聚类准确度,获得更准确的基函数中心,从而提高负荷预测精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例中RBF神经网络基本结构;
图2为本实施例中基于改进RBF神经网络的智能电网短期负荷预测方法的流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明是基于PCA降维的加权FCM聚类算法,合称PCA-FCM,得到智能电网短期负荷预测结果的。
另外,本发明中涉及到的迭代计算,指的是是逐次逼近,先取一个粗糙的近似值,然后用同一个或几个公式,反复校正此初值,直至达到预定精度要求为止。例如,本发明中涉及的迭代建立预测误差E的训练模型,其中涉及隐含层的输出、输出层的输出、加权值的更新等。
本发明中隐含层中心ci亦或称为基函数中心ci;隐含层神经元的半径ζi亦或称为方差ζi
实施例一
本实施例提出一种基于改进RBF神经网络的智能电网短期负荷预测方法。该方法采用PCA-WFCM聚类算法确定RBF基函数中心ci,采用梯度下降法确定RBF神经网络隐含层与输出层连接权值。下面详细的描述基于改进RBF神经网络的智能电网短期负荷预测方法:
首先建立样本数据X={x1,x2,…,xN},xj为每一个样本点,xj={xj1,xj2,…,xjs},j=1,2,…N,共有N个样本点,xj表示每一个样本点包含s个属性,在聚类分析中将样本数据X划分为K类,K的范围是2≤K≤N;聚类中心V={v1,v2,…,vk},共K个聚类中心。
径向基函数(RBF)神经网络是一种基于函数逼近理论的前馈神经网络,与BP(BackPropagation)神经网络相比,它具有函数逼近特性好、结构简单和训练速度快等优点。如图1所示,RBF神经网络是由输入层1、隐含层2、输出层3构成的三层神经网络,它的结构如下图所示。
RBF神经网络的核心思想是将径向基函数作为隐含层单元的基,构成隐含层空间,在隐含层中对输入矢量进行变换,使得在低维空间内线性不可分的数据在高维空间中线性可分。输入层到隐含层的变换是非线性变换,而隐含层到输出层是线性变换。隐含层变换函数即为径向基函数,它是局部分布中心径向对称的非负非线性函数。RBF神经网络输入层和隐含层之间的连接权值为1,隐含层完成对激活函数的参数调整,输出层对连接权值进行调整。在RBF神经网络中,需要求解的参数有3个:基函数中心ci、隐含层的宽度、隐含层到输出层的连接权重。高斯函数是RBF神经网络中常用的基函数,因此隐含层神经元的输出为:
Figure BDA0001179393020000071
其中,ci为隐含层中心,即第i个高斯函数的中心,N为隐含层神经元数,σi是RBF隐含层神经元的半径,可以表示为:
ζi=τminj‖ci-cj‖,i,j=1,2,…N; (2)
其中,ci表示隐含层中心,N表示隐含层神经元个数,η表示基函数重叠系数。
RBF神经网络的输出为所有隐含层神经元输出的线性组合,可以表示为:
Figure BDA0001179393020000072
其中,W是第j个隐含层神经元与输出神经元的连接权重,Ri表示隐含层的输出。
RBF神经网络非线性逼近过程中,给定训练样例后,算法需要解决以下两个关键问题:1)确定网络结构,即确定RBF神经网络的基函数中心ci;2)调整隐含层与输出层的连接权重ωn。这些参数的选取会影响RBF神经网络的预测性能,因此,做预测之前,我们要选取最优ωn和ci来提升RBF神经网络的预测性能。
隐含层与输出层的连接权重ωn一般采用梯度下降法训练。将一组输入向量{xj,j=1,2…O}和对应输出值{yj,j=1,2,…O}作为训练样本,其中K为样本数。则均方误差和函数为:
Figure BDA0001179393020000073
为了最小化误差函数,连接权重W为:
Figure BDA0001179393020000074
基函数中心ci采用PCA-WFCM聚类方法确定。下面详细介绍PCA-WFCM聚类算法。
PCA-WFCM聚类算法是基于传统FCM聚类算法,通过对样本数据做属性降维,降低算法复杂度,同时将降维后的每个属性方差贡献率作为属性权值做加权FCM聚类,进一步提升聚类性能。该算法包括两个步骤,先做PCA降维,然后做加权FCM聚类,下面分别介绍各个步骤。
主成分分析(PCA)是线性降维算法,可以将一个S维空间数据映射到L维子空间中,其中L<<S。数学计算上要求计算出原始数据的协方差矩阵的特征向量。假设X={xn},n=1,2,…,N为零均值数据,即
Figure BDA0001179393020000081
定义协方差矩阵C为:
Figure BDA0001179393020000082
做特征值分解得:
Figure BDA0001179393020000083
其中Q=[q1,q2,…,qS]为特征向量集合,Λ=diag(λ12,…,λS)为特征值。可以利用前L个特征向量UL=[u1,u2,…,uL]将一个S维的数据xi向L维主成分方向投影为Y=XQL。设协方差矩阵C的特征根λ1≥λ2≥…≥λS≥0,定义
Figure BDA0001179393020000084
为第l个主成分的贡献率,
Figure BDA0001179393020000085
为前L个主成分的累积贡献率。为了达到降维的目的,并希望信息损失小,一般取累积贡献率大于95%。
PCA-WFCM算法是基于FCM算法的,它是一种模糊聚类算法,即软划分方法。每一个样本点不能严格地被划分到某一类,而是以一定的隶属度属于某一类。令uij表示第j个样本点属于第i类的隶属度,则隶属度矩阵和聚类中心分别为U={uij}和V={vi}。PCA-WFCM算法是在FCM算法基础上通过给降维后属性赋予不同权值来将属性重要度考虑进去。属性的权值采用PCA对原始数据属性降维后每一属性的方差贡献率,贡献率越大的属性表示其在数据集中发挥的作用也越大。降维后第k个属性的权值:
Figure BDA0001179393020000091
l=1,2…,L。
聚类算法的目标是使类内相似性最大、类间相似性最小,且相似性采用欧式距离度量。因此算法通过最小化目标函数来确定聚类中心V和模糊矩阵U,该目标函数为
Figure BDA0001179393020000092
其中
Figure BDA0001179393020000093
其中dij是样本xj到聚类中心vi的加权欧式距离
Figure BDA0001179393020000094
式(8)中m≥1是模糊加权指数,表示隶属度矩阵U的模糊度,m越大,分类的模糊度越高,通常取m=2;L表示每一个样本点降维后的维度,Z表示聚类加权值。
通过对(8),(9)微分计算,我们可以得到uij和vi的计算公式为
Figure BDA0001179393020000095
Figure BDA0001179393020000096
其中xj表示第j个样本。
实施例二
根据实施例一中的思想,本实施例中神经网络预测需要事先确定好神经网络输入输出以及隐含层节点数。网络输入由影响预测值的一系列参数决定。因为智能电网用户负荷曲线具有良好的周期特性,因此对某一时刻的负荷值影响可考虑日周期特性及周周期特性,即,选择预测时刻前一天的同一时刻负荷值以及前一周同一时刻负荷值。具体预测输入层神经元个数N=9,包括8个负荷点:预测点前一时刻负荷值L(t-1)、预测点前两时刻负荷值L(t-2)、前一天同一预测点负荷值L(t-48)、前一天同一预测点前一时刻负荷值L(t-49)、前一天同一预测点后一时刻负荷值L(t-47)、前一周同一预测点负荷值L(t-48×7)、前一周同一预测点前一时刻负荷值L(t-48×7-1)、前一周同一预测点后一时刻负荷值L(t-48×7+1)。还包括一个日类型参数,即是否为周末。预测输出神经元数为N=1,即预测某一时刻的负荷。隐含层神经元数目N根据预测误差最小确定。所有网络输入都采用最大最小归一化处理。
如图2所示,基于改进RBF神经网络的智能电网短期负荷预测方法包括:
(1)网络初始化。根据***输入输出序列确定网络输入层神经元个数N,隐含层神经元个数N,输出层神经元个数N,以及初始化学习速率η和基函数重叠系数η。
(2)求RBF基函数中心ci。采用PCA-WFCM聚类算法确定基函数中心,具体流程为:
S21、输入模糊指数m,迭代停止阈值ε以及主成分累积贡献率因子δ,聚类数目即基函数中心个数N,连接权重ωn,样本数据X;
S22、对样本数据做PCA属性降维处理。根据式
Figure BDA0001179393020000101
求得L,保留维度L之前的所有维度作为聚类属性,并根据式
Figure BDA0001179393020000102
n=1,2,…,L;初始化每一个聚类属性,得到聚类属性的连接权重向量W={ωn}以及降维后的样本数据Xnew;其中,Xnew={x1,x2,…,xg},g为样本点是个数,xg为每一个样本点,xg={xg1,xg2,…,xgs},s表示每一个样本点包含的属性的个数,即每一个样本点的维度。
S23、根据式(9)初始化隶属度矩阵U;
S24、根据式(12)计算聚类中心V;
S25、根据式(11)计算隶属度矩阵U;
S26、根据式(8)计算目标函数J;
S27、如果|J(t)-J(t-1)|<ε,则得到聚类中心V,也就是基函数中心ci;反之,返回步骤S24;
(3)根据式(2)求解方差ζi
(4)计算隐含层的输出。根据降维后的样本数据Xnew,隐含层中心ci以及方差ζi,根据式(1)计算隐含层的输出。
(5)根据式(3)计算输出层的输出。
(6)根据式(4)计算预测误差E。
(7)根据式(5)进行聚类加权值的更新。
(8)判断算法迭代是否结束,若没结束,返回步骤(4)。
采用PCA-WFCM确定RBF基函数中心,负荷数据先通过PCA降维处理,获得较少的不相关的预测输入,从而减小RBF基函数的重叠,更好的确定基函数中心并降低预测算法复杂度;另外,采用加权FCM基函数聚类,将PCA处理后获得的不同属性的方差贡献率给降维后属性加权,从而提升聚类准确度,获得更准确的基函数中心,从而提高负荷预测精度。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。

Claims (7)

1.一种基于改进RBF神经网络的智能电网短期负荷预测方法,其特征在于,包括:
S1、网络初始化;
将智能电网样本数据进行降维处理;将PCA降维处理后的样本数据进PCA行加权FCM基函数聚类,其中,在对PCA降维处理后的智能电网样本数据进行加权FCM基函数聚类时,将PCA降维处理中获得的不同属性的方差贡献率给降维后属性聚类加权,得到基函数中心ci,具体如下:
S2、计算基函数中心ci
S3、根据基函数中心ci,计算方差ζi
聚类加权值的计算方法如下:
计算出协方差矩阵C:令一个S维空间数据映射到L维子空间中,其中,L<<S,假设X={xn}为零均值数据,即
Figure FDA0002385717240000011
其中,n=1,2,…,N,
Figure FDA0002385717240000012
T为转秩符号;
对协方差矩阵C进行特征值分解:
Figure FDA0002385717240000013
将一个S维的数据xi向L维主成分方向投影,即,Y=XQL
设协方差矩阵C的特征根λ1≥λ2≥…≥λS,且
Figure FDA0002385717240000014
为第l个主成分的贡献率,
Figure FDA0002385717240000015
为前L个主成分的累积贡献率;
降维后第k个属性的聚类属性权值:
Figure FDA0002385717240000016
特征向量合集Q=[q1,q2,…,qS],特征值Л=diag(λ12,…,λS),令累积贡献率大于95%;
依据得到的基函数中心ci和方差ζi,对加权FCM基函数聚类后的样本数据进行RBF神经网络预测得到输出层输出,具体如下:
S4、根据基函数中心ci以及方差ζi,计算隐含层的输出Ri
步骤S4包括:
根据降维后的样本数据Xnew、基函数中心ci以及方差ζi,得到隐含层的输出:
Figure FDA0002385717240000021
式中,N表示隐含层神经元个数;
S5、根据隐含层的输出Ri,计算输出层的输出,从而得到智能电网短期负荷预测结果;
S6、根据均方误差和函数计算预测误差E;
S7、对神经网络中隐含层神经元与输出层神经元的连接权重进行更新;
S8、对预测误差E进行判断,若预测误差E在预期内,则迭代计算结束;反之,返回步骤S4,重新迭代计算预测误差E。
2.根据权利要求1所述的基于改进RBF神经网络的智能电网短期负荷预测方法,其特征在于,步骤S1包括:
确定输入层神经元个数N,隐含层神经元个数N,输出层神经元个数N,以及初始化学习速率η和基函数重叠系数τ;其中,输出层神经元个数N=1,隐含层神经元个数N即为基函数中心个数。
3.根据权利要求1所述的基于改进RBF神经网络的智能电网短期负荷预测方法,其特征在于,步骤S2包括:
S21、输入模糊指数m,迭代停止阈值ε以及PCA累积贡献率因子δ,基函数中心个数N,聚类的属性权值ωn,原始数据X;
其中,样本数据X={x1,x2,…,xN},N为样本点的个数,xj为每一个样本点,xj={xj1,xj2,…,xjs},j=1,2,…N,s表示每一个样本点包含的属性的个数,即每一个样本点的维度,在PCA-WFCM方法中将样本数据X划分为K类,聚类中心V={v1,v2,…,vk},共有K个,K的范围为2≤K≤N;
S22、对样本数据X进行PCA属性降维处理:根据如下公式计算出每一个样本点降维后的维度L:
Figure FDA0002385717240000031
保留维度L之前的所有维度作为聚类属性;
式中,S表示每一个样本点的原始维度,λn为特征值,表示PCA中方差贡献率,δ表示PCA累积贡献率因子;
根据如下公式计算每一个聚类的属性权值:
Figure FDA0002385717240000032
根据聚类的属性权值ωn得到聚类属性的权重向量W={ωn}以及降维后的样本数据Xnew
S23、初始化隶属度矩阵U:
令0≤uij≤1,
Figure FDA0002385717240000033
U={uij};
式中,uij表示第j个样本点属于第i类的隶属度,K表示聚类中心V的个数;
S24、根据隶属度矩阵U以及降维后的样本数据Xnew,计算聚类中心V:
Figure FDA0002385717240000034
V={vi};
式中,m为模糊指数,xj表示第j个样本点;
S25、根据降维后的样本数据Xnew,迭代计算隶属度矩阵U:
Figure FDA0002385717240000035
U={uij};
式中,m为模糊指数,表示隶属度矩阵U的模糊度,m的值越大,则隶属度矩阵U的模糊度越高,令m=2,K表示聚类中心V的个数,dij表示每一个样本点xj到聚类中心vi的加权欧式距离,由如下公式计算得出:
Figure FDA0002385717240000036
L表示每一个样本点降维后的维度,ωn表示聚类的属性权值;
S26、根据隶属度矩阵U和聚类中心V计算目标函数J:
Figure FDA0002385717240000037
式中,m为模糊指数,dij表示每一个样本点xj到聚类中心vi的加权欧式距离,K表示聚类中心V的个数,N表示样本点个数;
S27、对目标函数J进行判断:若|J(t)-J(t-1)|<ε,则输出聚类中心V,也就是基函数中心ci;反之,则返回步骤S24,直到满足公式|J(t)-J(t-1)|<ε,则停止迭代计算,计算出聚类中心V,式中,ε表示迭代停止阈值,t表示迭代次数。
4.根据权利要求1所述的基于改进RBF神经网络的智能电网短期负荷预测方法,其特征在于,步骤S3包括:
根据如下公式计算隐含层神经元的半径σi
σi=τminj||ci-cj||,i,j=1,2,…N
式中,ci表示基函数中心,N表示隐含层神经元个数,τ表示基函数重叠系数。
5.根据权利要求1所述的基于改进RBF神经网络的智能电网短期负荷预测方法,其特征在于,步骤S5包括:
根据降维后的样本数据Xnew、得到输出层的输出:
Figure FDA0002385717240000041
式中,N表示隐含层神经元个数;W表示第i个隐含层神经元与输出层神经元的连接权重,Ri表示隐含层的输出。
6.根据权利要求1所述的基于改进RBF神经网络的智能电网短期负荷预测方法,其特征在于,步骤S6包括:
利用均方误差和函数计算预测误差E:
令将一组输入向量{xj,j=1,2…O}和对应输出值{yj,j=1,2,…O}作为训练样本,
式中,O为样本数,预测误差
Figure FDA0002385717240000042
7.根据权利要求1所述的基于改进RBF神经网络的智能电网短期负荷预测方法,其特征在于,步骤S7包括:
根据如下公式对连接权重W进行更新:
Figure FDA0002385717240000043
式中,η表示学习速率,E为预测误差,q为更新次数。
CN201611148874.7A 2016-12-13 2016-12-13 基于改进rbf神经网络的智能电网短期负荷预测方法 Active CN106600059B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611148874.7A CN106600059B (zh) 2016-12-13 2016-12-13 基于改进rbf神经网络的智能电网短期负荷预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611148874.7A CN106600059B (zh) 2016-12-13 2016-12-13 基于改进rbf神经网络的智能电网短期负荷预测方法

Publications (2)

Publication Number Publication Date
CN106600059A CN106600059A (zh) 2017-04-26
CN106600059B true CN106600059B (zh) 2020-07-24

Family

ID=58802260

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611148874.7A Active CN106600059B (zh) 2016-12-13 2016-12-13 基于改进rbf神经网络的智能电网短期负荷预测方法

Country Status (1)

Country Link
CN (1) CN106600059B (zh)

Families Citing this family (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107330518A (zh) * 2017-06-21 2017-11-07 国家电网公司 基于温度调节负荷预测的能源管理控制方法和***
CN107403188A (zh) * 2017-06-28 2017-11-28 中国农业大学 一种水质评价方法及装置
CN107194524B (zh) * 2017-07-28 2020-05-22 合肥工业大学 一种基于rbf神经网络的煤与瓦斯突出预测方法
CN107657351A (zh) * 2017-10-26 2018-02-02 广州泰阳能源科技有限公司 一种基于plc与主元分析‑rbf神经网络的负荷预测***
CN108021935B (zh) * 2017-11-27 2024-01-23 中国电力科学研究院有限公司 一种基于大数据技术的维度约简方法及装置
CN108229754B (zh) * 2018-01-31 2021-12-10 杭州电子科技大学 基于相似日分段和lm-bp网络的短期负荷预测方法
CN108230121B (zh) * 2018-02-09 2022-06-10 艾凯克斯(嘉兴)信息科技有限公司 一种基于循环神经网络的产品设计方法
CN108680358A (zh) * 2018-03-23 2018-10-19 河海大学 一种基于轴承温度模型的风电机组故障预测方法
CN108631817B (zh) * 2018-05-10 2020-05-19 东北大学 一种基于时频分析和径向神经网络进行跳频信号频段预测的方法
CN108734355B (zh) * 2018-05-24 2022-03-08 国网福建省电力有限公司 一种应用于电能质量综合治理场景的短期电力负荷并行预测方法及***
CN109309382B (zh) * 2018-09-13 2022-03-04 广东工业大学 一种短期电力负荷预测方法
CN109179133A (zh) * 2018-11-05 2019-01-11 常熟理工学院 用于预判故障的电梯智能维保预测方法及***
CN109583044B (zh) * 2018-11-09 2022-07-15 中国直升机设计研究所 基于rbf神经网络的直升机旋翼飞行载荷预测方法
CN109284876A (zh) * 2018-11-19 2019-01-29 福州大学 基于pca-rbf埋地管道腐蚀速率预测方法
CN110097126B (zh) * 2019-05-07 2023-04-21 江苏优聚思信息技术有限公司 基于dbscan聚类算法的核查重点人员、房屋漏登记的方法
CN110059824A (zh) * 2019-05-22 2019-07-26 杭州电子科技大学 一种基于主成分分析的神经网络预测方法
CN110365647B (zh) * 2019-06-13 2021-09-14 广东工业大学 一种基于pca和bp神经网络的虚假数据注入攻击检测方法
CN110443318B (zh) * 2019-08-09 2023-12-08 武汉烽火普天信息技术有限公司 一种基于主成分分析和聚类分析的深度神经网络方法
CN110796158A (zh) * 2019-09-10 2020-02-14 国网浙江省电力有限公司杭州供电公司 基于rbf径向基神经网络的电网公司分类方法
CN110824142B (zh) * 2019-11-13 2022-06-24 杭州鲁尔物联科技有限公司 一种地质灾害的预测方法、装置及设备
CN110852522B (zh) * 2019-11-19 2024-03-29 南京工程学院 一种短期电力负荷预测方法及***
JP7295431B2 (ja) * 2019-11-27 2023-06-21 富士通株式会社 学習プログラム、学習方法および学習装置
CN110991638B (zh) * 2019-11-29 2024-01-05 国网山东省电力公司聊城供电公司 基于聚类及神经网络的广义负荷建模方法
CN111259943A (zh) * 2020-01-10 2020-06-09 天津大学 基于机器学习的温跃层预测方法
CN112580853A (zh) * 2020-11-20 2021-03-30 国网浙江省电力有限公司台州供电公司 基于径向基神经网络的母线短期负荷预测方法
CN112766076B (zh) * 2020-12-31 2023-05-12 上海电机学院 一种电力负荷超短期预测方法、***、设备及存储介质
CN113408699A (zh) * 2021-06-16 2021-09-17 中国地质科学院 基于改进的径向基函数神经网络的岩性识别方法及***
CN114875196B (zh) * 2022-07-01 2022-09-30 北京科技大学 一种转炉出钢量的确定方法和***
CN115796057A (zh) * 2023-02-06 2023-03-14 广东电网有限责任公司中山供电局 基于bas-lstm的电缆接头温度预测方法和***

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102982393A (zh) * 2012-11-09 2013-03-20 山东电力集团公司聊城供电公司 一种输电线路动态容量的在线预测方法
CN103095494A (zh) * 2012-12-31 2013-05-08 北京邮电大学 一种电力通信网风险评估方法
CN103136598A (zh) * 2013-02-26 2013-06-05 福建省电力有限公司 基于小波分析的月度电力负荷计算机预测方法
CN103646354A (zh) * 2013-11-28 2014-03-19 国家电网公司 基于有效指标fcm和rbf神经网络的变电站负荷特性分类方法
CN105678404A (zh) * 2015-12-30 2016-06-15 东北大学 基于网购电和动态关联因子的微电网负荷预测***及方法
CN105787584A (zh) * 2016-01-28 2016-07-20 华北电力大学(保定) 一种基于云平台的风电机群故障预警方法
EP3098762A1 (en) * 2015-05-29 2016-11-30 Samsung Electronics Co., Ltd. Data-optimized neural network traversal

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102982393A (zh) * 2012-11-09 2013-03-20 山东电力集团公司聊城供电公司 一种输电线路动态容量的在线预测方法
CN103095494A (zh) * 2012-12-31 2013-05-08 北京邮电大学 一种电力通信网风险评估方法
CN103136598A (zh) * 2013-02-26 2013-06-05 福建省电力有限公司 基于小波分析的月度电力负荷计算机预测方法
CN103646354A (zh) * 2013-11-28 2014-03-19 国家电网公司 基于有效指标fcm和rbf神经网络的变电站负荷特性分类方法
EP3098762A1 (en) * 2015-05-29 2016-11-30 Samsung Electronics Co., Ltd. Data-optimized neural network traversal
CN105678404A (zh) * 2015-12-30 2016-06-15 东北大学 基于网购电和动态关联因子的微电网负荷预测***及方法
CN105787584A (zh) * 2016-01-28 2016-07-20 华北电力大学(保定) 一种基于云平台的风电机群故障预警方法

Also Published As

Publication number Publication date
CN106600059A (zh) 2017-04-26

Similar Documents

Publication Publication Date Title
CN106600059B (zh) 基于改进rbf神经网络的智能电网短期负荷预测方法
CN109063911B (zh) 一种基于门控循环单元网络的负荷聚合体分组预测方法
CN110175386B (zh) 变电站电气设备温度预测方法
CN111027772B (zh) 基于pca-dbilstm的多因素短期负荷预测方法
CN108985515B (zh) 一种基于独立循环神经网络的新能源出力预测方法及***
CN106067034B (zh) 一种基于高维矩阵特征根的配电网负荷曲线聚类方法
Li et al. A novel double incremental learning algorithm for time series prediction
CN109492748B (zh) 一种基于卷积神经网络的电力***的中长期负荷预测模型建立方法
CN109255726A (zh) 一种混合智能技术的超短期风功率预测方法
Chen et al. Construction of tunable radial basis function networks using orthogonal forward selection
CN112434848A (zh) 基于深度信念网络的非线性加权组合风电功率预测方法
CN110895772A (zh) 基于灰色关联分析与SA-PSO-Elman算法结合的售电量预测方法
CN113641722A (zh) 一种基于变体lstm的长期时间序列数据预测方法
CN111355633A (zh) 一种基于pso-delm算法的比赛场馆内手机上网流量预测方法
CN112766603A (zh) 一种交通流量预测方法、***、计算机设备及存储介质
CN110471768A (zh) 一种基于fastPCA-ARIMA的负载预测方法
CN116526450A (zh) 计及误差补偿的两阶段短期电力负荷组合预测方法
CN115186803A (zh) 一种考虑pue的数据中心算力负荷需求组合预测方法和***
CN113722980A (zh) 海洋浪高预测方法、***、计算机设备、存储介质、终端
CN110738363B (zh) 一种光伏发电功率预测方法
CN113836823A (zh) 一种基于负荷分解和优化双向长短期记忆网络的负荷组合预测方法
Han et al. A Hybrid BPNN-GARF-SVR PredictionModel Based on EEMD for Ship Motion.
CN111079902A (zh) 基于神经网络的分解模糊***优化方法及装置
CN114330119B (zh) 一种基于深度学习的抽蓄机组调节***辨识方法
Schiavo et al. Powerful and flexible fuzzy algorithm for nonlinear dynamic system identification

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