CN108776843A - 基于波尔兹曼机的光伏超短期功率预测方法 - Google Patents
基于波尔兹曼机的光伏超短期功率预测方法 Download PDFInfo
- Publication number
- CN108776843A CN108776843A CN201810325176.2A CN201810325176A CN108776843A CN 108776843 A CN108776843 A CN 108776843A CN 201810325176 A CN201810325176 A CN 201810325176A CN 108776843 A CN108776843 A CN 108776843A
- Authority
- CN
- China
- Prior art keywords
- rbm
- photovoltaic
- formula
- data
- training
- 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.)
- Pending
Links
- 238000013277 forecasting method Methods 0.000 title claims abstract description 8
- 238000000034 method Methods 0.000 claims abstract description 26
- 238000012549 training Methods 0.000 claims description 38
- 238000009826 distribution Methods 0.000 claims description 19
- 230000006870 function Effects 0.000 claims description 15
- 230000004913 activation Effects 0.000 claims description 9
- 230000005855 radiation Effects 0.000 claims description 9
- 239000006185 dispersion Substances 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 238000013480 data collection Methods 0.000 claims description 4
- 238000005192 partition Methods 0.000 claims description 4
- 238000001556 precipitation Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 239000007788 liquid Substances 0.000 claims description 3
- 238000010248 power generation Methods 0.000 abstract description 8
- 238000013135 deep learning Methods 0.000 abstract description 4
- 238000005516 engineering process Methods 0.000 abstract description 3
- 230000000694 effects Effects 0.000 description 9
- 238000012360 testing method Methods 0.000 description 7
- 241001269238 Data Species 0.000 description 4
- 241000208340 Araliaceae Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 210000004027 cell Anatomy 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012417 linear regression Methods 0.000 description 2
- 238000010801 machine learning Methods 0.000 description 2
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 210000002569 neuron Anatomy 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 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/06—Energy or water supply
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及光伏发电功率预测、深度学习技术领域，为提出一种基于深度学习中的限制波尔兹曼机模型对光伏功率进行预测的方法。为精确预测光伏功率提供有力支持。为此，本发明采用的技术方案是，基于波尔兹曼机的光伏超短期功率预测方法，利用影响光伏的因素，通过DBNS模型对数据之间进行训练以及拟合，获取其与光伏功率之间的相关性，从而对光伏超短期功率进行预测。本发明主要应用于光伏发电场合。
Description
技术领域
本发明涉及光伏发电功率预测、深度学习技术领域，具体讲,涉及基于波尔兹曼机的光伏超短期功率预测方法。
背景技术
太阳能已被公认为是未来最具竞争力能源之一，具有取之不尽、用之不竭绿色环保等特点。目前，光伏发电(Photovoltaic power generation)是利用太阳能的有效方式之一，分布式光伏并网发电系统与大规模集中式光伏并网系统是国内光伏利用的两种重要的方式。在环境和能源等压力下，分布式光伏发电是未来重点的发展方向。由于太阳能受各种气象因素的影响，使用光伏系统进行发电的功率是不稳定的，会造成对电力系统的不利影响。光伏功率短期预测的时间尺度一般为1～3天，超短期预测的时间一般为0～4小时，短期和超短期功率预测对电网实时调度和稳定等具有重要的决定作用。
针对于光伏发电功率的预测，国内外学者做了大量的研究。目前所有短期太阳能光伏发电预测方法都是基于相同的思路，即：首先利用数学和物理学理论及相关数据建立预测公式或模型，再通过求得的预测公式或模型对光伏电站发电量进行预测。
限制波尔兹曼机(Restricted Boltzmann Machine,RBM)是一种基于能量的模型，是Hinton 于1986年提出的波尔兹曼机(Boltzmann Machine,BM)的一种改进，当其仅有显元的分布为高斯分布时称为高斯波尔兹曼机(Gaussian-Bernoulli RBM,GBRBM)。同样RBM也是深度学习 (deep learning,DL)中的一个重要的模型DBN的基本单元。RBM可以看作一种自动编码器 (autoencoder)，可利用其强大的无监督学习的能力，将输入的数据进行重构，将数据从原始空间变换到其他维数的空间，易于发现其内部之间隐含的关系，从而便于构建后续的机器学习模型。
发明内容
为克服现有技术的不足，本发明旨在提出一种基于深度学习中的限制波尔兹曼机模型对光伏功率进行预测的方法。为精确预测光伏功率提供有力支持。为此，本发明采用的技术方案是，基于波尔兹曼机的光伏超短期功率预测方法，利用影响光伏的因素，通过DBNS模型对数据之间进行训练以及拟合，获取其与光伏功率之间的相关性，从而对光伏超短期功率进行预测。
具体步骤细化如下：
步骤1：原始数据的获取
首先获取原始的光伏数据，数据包括影响光伏的12个相关因素：具体有总柱液态水TCLW、总柱冰状水TCIW、表面压力SP、1000毫巴相对湿度RH、总云量TCC、10米高风的U分量10U、10米高风的V分量10V、2米处温度2T、表面太阳辐射累积量SSRD、表面热辐射累积量STRD、大气顶部太阳辐射累积量TSR、总降水TP，以及相对应的光伏发电数据；
步骤2：数据集的划分
剔除每天9时～20时夜间光伏发电为0的情况，分析其中剩余光伏发电不为0的情况，并将数据集分为A和B两个部分，B数据集包括某个月数据，A数据集包括其余数据，并且利用RBM将数据集A进行了重构得到数据集A'；
步骤3：模型的构建
训练模型是由波尔兹曼机BM为基础单元构建的,通过波尔兹曼机的组合，构成训练的网络模型，通过各个BM的显元与隐元之间参数的训练对其传递的权值进行调整，并通过选择合适的决策模型获取最优的预测。
在步骤(3)模型的构建过程中可以对波尔兹曼机进行优选，改进方式如下：
采用限制玻尔兹曼机RBM；或者，采用高斯玻尔兹曼机GBRBM，对于一个有n个显元，m个隐元的RBM，v代表显元的状态，h代表隐元的状态，那么RBM的概率分布P(v,h)满足Boltzmann分布；高斯RBM，其能量函数如式(1)所示：
其中：vi表示二值显元的取值；hj表示隐元的取值；ai,bj表示偏置(bias)；wij表示第i个显元和第j个隐元间的权重；σi表示显元的标准差；
全概率分布P(v,h)如式(2)所示：
其中：Z为配分函数，表达式为式(3)：
表示v和h的所有可能的组合；
根据公式(2)得到高斯RBM关于v的分布P(v)，即P(v,h)的边缘分布，定义P(v)为高斯RBM的似然函数，表达式为：
由于高斯RBM各隐元之间的激活条件是独立的，所以第j个隐元的激活概率如式(5) 所示：
同理，第i个显元的激活概率如式(6)所示：
其中：是罗杰斯特函数(logistics sigmoid function)。
训练RBM的任务就是确定RBM的参数的值，即θ＝{wij,ai,bj}，参数θ可以通过最大化 RBM在含有T个样本训练集上的对数似然函数学习得到，描述式如式(7)所示：
RBM参数计算方法有随机梯度上升法和对比散度法；
1)随机梯度上升法
利用随机梯度上升法可以获取最优的参数θ*，即求出：
对于模型参数的偏导数计算如式(9)～(11)所示：
其中：v0表示输入样本；<·>P表示关于分布P的数学期望；P(h|v0,θ)表示显元状态为已知输入样本v0时隐层的概率分布；采用对比散度法CD(Contrastive Divergence)算法获取其近似值；
2)对比散度法CD-k(Contrastive Divergence)
在CD-k算法开始时，设定该训练样本为显元的状态，利用式(5)计算隐元的状态；然后根据式(6)再重新计算显元的状态，即显元取值为1的概率，这样就产生了显层的一个重构 (reconstruction)，k的取值即为重构的次数；
式(9)～(11)的近似计算结果，如式(12)～(14)表示：
此时，RBM参数的增量为：
ε表示学习率；
3)假设在利用训练集样本训练RBM前的参数为θ，训练之后RBM的参数为θ*，RBM 参数的增量为Δθ，则对应有：
θ*＝θ+Δθ (17)
由步骤1)～3)计算，最终可得在已知训练集上所构建的RBM的参数。
本发明的特点及有益效果是：
本发明基于GEFCom2014所给出的光伏出力的相关数据，在构建的训练集与测试集的基础上，针对GBRBM重构前的数据集A构建了线性回归LR模型和SVR模型，并在重构后的数据集A'的基础上构建了LR模型，称其为GBRBM-LR模型，然后对三种模型在未参与训练的数据集B上进行了测试，对发电功率取得了较好的预测效果。
附图说明：
图1 RBM和BM结构。
图中，(a)BM结构 (b)RBM结构。
图2 CD-k算法示意。
图3 2014年5月该块光伏面板发电功率。
图4 GBRBM-LR模型预测结果。
图5 LR模型预测结果。
图6 SVR模型预测结果。
具体实施方式
基于波尔兹曼机对光伏进行超短期预测需要进行以下的步骤：
步骤1：原始数据的获取
首先获取原始的光伏数据，数据来源于The Global Energy ForecastingCompetition 2014 (GEFCom2014)，是从2012年4月1日至2014年6月1日对某块光伏面板的相关影响因素进行逐小时采集得到的。数据包括12个相关因素，包括总柱液态水(TCLW)、总柱冰状水(TCIW)、表面压力(SP)、1000毫巴相对湿度(RH)、总云量(TCC)、10米高风的U分量(10U)、 10米高风的V分量(10V)、2米处温度(2T)、表面太阳辐射累积量(SSRD)、表面热辐射累积量(STRD)、大气顶部太阳辐射累积量(TSR)、总降水(TP)，涉及风速、气温、气压、辐照、降水等，以及相对应的光伏发电数据。
步骤2：数据集的划分
影响光伏发电功率预测的影响因素有很多关系也相当复杂，剔除每天9时～20时(世界时)夜间光伏发电为0的情况，分析其中剩余光伏发电不为0的情况，并将数据集分为A和 B两个部分，B数据集包括2014年5月份数据，A数据集包括其余数据，并且利用RBM将数据集A进行了重构得到数据集A'。
为了实现模型的预测效果，在训练模型时，输入的影响因素的标签为未来1小时的光伏发电数据。在这种情况下，训练得到的模型输入的是当前12个相关影响因素，输出的是未来 1小时光伏的发电功率，从而实现了对光伏面板超短期(1小时)发电功率预测。
步骤3：模型的构建
训练模型是由波尔兹曼机BM为基础单元构建的,通过波尔兹曼机的组合，构成了训练的网络模型。通过各个BM的显元与隐元之间参数的训练对其传递的权值进行调整，并通过选择合适的决策模型获取最优的预测。
2014年5月该块光伏面板发电功共计372条数据，作为模型的测试集B。将2012年4月-2014年4月共计9120条数据作为训练集A，利用RBM将数据集A进行了重构得到数据集A'，然后分别训练了GBRBM-LR模型,LR模型和SVR模型。最后在测试集上对光伏面板发电功率进行了测试。
将模型预测输出得到的结果(generated value)用蓝色交叉点所示，三种模型在测试集上得到的结果如图4～图6所示，其中横坐标代表时间(年-月-日)，纵坐标代表光伏发电状况(标幺值)。MAS、MSE、EV、R2结果如表1所示利用平均绝对误差MAE(Mean AbsoluteError)、平均平方误差MSE(Mean Squared Error)、解释变异(Explained variance，EV)以及决定系数 (Coefficient of determination,又称R2分数)来分别评价三个光伏面板超短期发电功率预测模型的表现。
表1三种模型评价指标
在步骤(3)模型的构建过程中可以对波尔兹曼机进行优选。改进方式如下:
1.RBM(限制玻尔兹曼机)
RBM是BM的一种改进，BM具有强大的无监督学习能力，RBM在继承了BM强大的无监督学习能力的前提下对其进行了改进，限制了BM显层和隐层的连接。BM各层单元之间是相互连接的，而RBM各层单元之间并不存在相互连接的关系。
相对于BM，可以通过建立RBM来较好地重构原始数据，发现数据间隐含存在的关系，并且避免训练时间过长，从而可以便于后续机器学习模型的建立。
2.GBRBM(高斯玻尔兹曼机)
对于一个有n个显元，m个隐元的RBM，v代表显元(visible units)的状态，h代表隐元 (hidden units)的状态，那么RBM的概率分布P(v,h)满足Boltzmann分布。
实际的数据往往近似满足高斯分布，此时假设显元具有独立高斯噪音的线性单元，隐元为随机二值神经元，此时RBM即高斯-伯努利RBM，简称高斯RBM，其能量函数如式(1)所示：
其中：vi表示二值显元的取值；hj表示隐元的取值；ai,bj表示偏置(bias)；wij表示第i个显元和第j个隐元间的权重(weight)；σi表示显元的标准差。
全概率分布P(v,h)如式(2)所示：
其中：Z为配分函数(partition function)，表达式为式(3)：
表示v和h的所有可能的组合。
根据公式(2)可以得到高斯RBM关于v的分布P(v)，即P(v,h)的边缘分布。定义P(v)为高斯RBM的似然函数(likelihood)，表达式为：
由于高斯RBM各隐元之间的激活条件是独立的，所以第j个隐元的激活概率如式(5) 所示：
同理，第i个显元的激活概率如式(6)所示：
其中：是罗杰斯特函数(logistics sigmoid function)。
3.参数确定
训练RBM的任务就是确定RBM的参数的值，即θ＝{wij,ai,bj}。参数θ可以通过最大化 RBM在含有T个样本训练集上的对数似然函数学习得到，描述式如式(7)所示：
下面具体描述RBM参数计算方法。
4)随机梯度上升法(stochastic gradient ascent)
利用随机梯度上升法可以获取最优的参数θ*，即求出：
一般情况下，RBM的初始参数是随机初始化而得到的，此时只需要在训练集上计算出 RBM模型各个参数(即wij,ai,bj)的增量即能得到最终的RBM参数。因此计算logP(vt|θ)对模型的各个参数偏导数对于求得参数θ的增量至关重要。
对于模型参数(wij或ai或bj)的偏导数计算如式(9)～(11)所示：
其中：v0表示输入样本；<·>P表示关于分布P的数学期望；P(h|v0,θ)表示显元状态为已知输入样本v0时隐层的概率分布。
由于Z的存在(计算Z至少需要2n+m次计算)，P(v,h|θ)难以计算，因此本文采用对比散度法(Contrastive Divergence,CD)算法获取其近似值。
5)对比散度法(Contrastive Divergence,CD-k)
为了避免RBM训练时间过长，Hinton提出RBM的一个快速学习算法——CD-k算法。CD-k算法流程如图2所示。
在CD-k算法开始时，设定该训练样本为显元的状态，利用式(5)计算隐元的状态；然后根据式(6)再重新计算显元的状态，即显元取值为1的概率，这样就产生了显层的一个重构 (reconstruction)。k的取值即为重构的次数。
Hinton发现，k＝1时就能取得较好的效果。同时在利用CD-1算法训练RBM时，若用“data”来简记P(h|v0,θ)，“recon”来简记P(v,h|θ)，CD-1算法就可以得出式(9)～(11)的近似计算结果，如式(12)～(14)表示：
此时，RBM参数的增量为：
ε表示学习率。
6)假设在利用训练集样本训练RBM前的参数为θ，训练之后RBM的参数为θ*，RBM 参数的增量为Δθ，则对应有：
θ*＝θ+Δθ (17)
由步骤1)～3)计算，最终可得在已知训练集上所构建的RBM的参数。
4.GBRBM-LR训练模型
GBRBM有着强大的无监督学习的能力，当一个样本经过训练好的RBM，输出就是对应的重构后的数据。能够发现数据内在的结构，也即特征提取和选择的过程。本文在重构后的数据A′上构建LR模型，即为GBRBM-LR模型。这相当于GBRBM-LR模型输入的是当前 12个相关影响因素，输出的是未来1小时光伏的发电功率。
本发明基于GEFCom2014所给出的光伏出力的相关数据，在构建的训练集与测试集的基础上，针对GBRBM重构前的数据集A构建了线性回归LR模型和SVR模型，并在重构后的数据集A'的基础上构建了LR模型，称其为GBRBM-LR模型，然后对三种模型在未参与训练的数据集B上进行了测试。
由图4～6可以看出：GBRBM-LR模型的拟合效果远远优于LR模型的拟合效果，与SVR模型的拟合效果相当。说明GBRBM对LR模型的拟合能力有着极大的提高作用，达到甚至优于SVR模型的拟合能力，从而在测试集上，GBRBM-LR模型对光伏面板超短期(1小时) 发电功率有着较好的预测效果。
从表1所示模型评价指标比较可得，GBRBM-LR模型的MAE和MSE的值小于LR模型和SVR模型的值，说明GBRBM-LR输出的预测值与真实值更加接近；EV的值较大说明 GBRBM-LR输出的预测值的分布与真实值的分布更加相似。
从表1还可知：GBRBM-LR模型的R2的值更接近于1，说明其对数据的拟合程度更高，同样也说明了RBM重构之后的数据(自变量)，应用于LR模型上，比在原始数据上直接应用LR模型，对光伏发电数据(因变量)的解释程度更高，即RBM对数据所进行的重构，克服了数据在低维空间内难以构造良好的线性模型的问题，但是在相应的高维空间中构建线性模型就可以取得良好的拟合效果。
Claims (3)
1.一种基于波尔兹曼机的光伏超短期功率预测方法，其特征是，利用影响光伏的因素，通过DBNS模型对数据之间进行训练以及拟合，获取其与光伏功率之间的相关性，从而对光伏超短期功率进行预测。
2.如权利要求1所述的基于波尔兹曼机的光伏超短期功率预测方法，其特征是，具体步骤细化如下：
步骤1：原始数据的获取
首先获取原始的光伏数据，数据包括影响光伏的12个相关因素：具体有总柱液态水TCLW、总柱冰状水TCIW、表面压力SP、1000毫巴相对湿度RH、总云量TCC、10米高风的U分量10U、10米高风的V分量10V、2米处温度2T、表面太阳辐射累积量SSRD、表面热辐射累积量STRD、大气顶部太阳辐射累积量TSR、总降水TP，以及相对应的光伏发电数据；
步骤2：数据集的划分
剔除每天9时～20时夜间光伏发电为0的情况，分析其中剩余光伏发电不为0的情况，并将数据集分为A和B两个部分，B数据集包括某个月数据，A数据集包括其余数据，并且利用RBM将数据集A进行了重构得到数据集A'；
步骤3：模型的构建
训练模型是由波尔兹曼机BM为基础单元构建的,通过波尔兹曼机的组合，构成训练的网络模型，通过各个BM的显元与隐元之间参数的训练对其传递的权值进行调整，并通过选择合适的决策模型获取最优的预测。
3.如权利要求1所述的基于波尔兹曼机的光伏超短期功率预测方法，其特征是，在步骤(3)模型的构建过程中可以对波尔兹曼机进行优选，改进方式如下：采用限制玻尔兹曼机RBM；或者，采用高斯玻尔兹曼机GBRBM，对于一个有n个显元，m个隐元的RBM，v代表显元的状态，h代表隐元的状态，那么RBM的概率分布P(v,h)满足Boltzmann分布；高斯RBM，其能量函数如式(1)所示：
其中：vi表示二值显元的取值；hj表示隐元的取值；ai,bj表示偏置(bias)；wij表示第i个显元和第j个隐元间的权重；σi表示显元的标准差；
全概率分布P(v,h)如式(2)所示：
其中：Z为配分函数，表达式为式(3)：
表示v和h的所有可能的组合；
根据公式(2)得到高斯RBM关于v的分布P(v)，即P(v,h)的边缘分布，定义P(v)为高斯RBM的似然函数，表达式为：
由于高斯RBM各隐元之间的激活条件是独立的，所以第j个隐元的激活概率如式(5)所示：
同理，第i个显元的激活概率如式(6)所示：
其中：是罗杰斯特函数(logistics sigmoid function)。
训练RBM的任务就是确定RBM的参数的值，即θ＝{wij,ai,bj}，参数θ可以通过最大化RBM在含有T个样本训练集上的对数似然函数学习得到，描述式如式(7)所示：
RBM参数计算方法有随机梯度上升法和对比散度法；
1)随机梯度上升法
利用随机梯度上升法可以获取最优的参数θ*，即求出：
对于模型参数的偏导数计算如式(9)～(11)所示：
其中：v0表示输入样本；<·>P表示关于分布P的数学期望；P(h|v0,θ)表示显元状态为已知输入样本v0时隐层的概率分布；采用对比散度法CD(Contrastive Divergence)算法获取其近似值；
2)对比散度法CD-k(Contrastive Divergence)
在CD-k算法开始时，设定该训练样本为显元的状态，利用式(5)计算隐元的状态；然后根据式(6)再重新计算显元的状态，即显元取值为1的概率，这样就产生了显层的一个重构(reconstruction)，k的取值即为重构的次数；
式(9)～(11)的近似计算结果，如式(12)～(14)表示：
此时，RBM参数的增量为：
ε表示学习率；
3)假设在利用训练集样本训练RBM前的参数为θ，训练之后RBM的参数为θ*，RBM参数的增量为△θ，则对应有：
θ*＝θ+△θ (17)
由步骤1)～3)计算，最终可得在已知训练集上所构建的RBM的参数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810325176.2A CN108776843A (zh) | 2018-04-12 | 2018-04-12 | 基于波尔兹曼机的光伏超短期功率预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810325176.2A CN108776843A (zh) | 2018-04-12 | 2018-04-12 | 基于波尔兹曼机的光伏超短期功率预测方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108776843A true CN108776843A (zh) | 2018-11-09 |
Family
ID=64033733
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810325176.2A Pending CN108776843A (zh) | 2018-04-12 | 2018-04-12 | 基于波尔兹曼机的光伏超短期功率预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108776843A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113297791A (zh) * | 2021-05-18 | 2021-08-24 | 四川大川云能科技有限公司 | 一种基于改进dbn的风功率组合预测方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20160252405A1 (en) * | 2014-07-31 | 2016-09-01 | Huazhong University Of Science And Technology | Noninvasive measuring method for rapid temperature variation under dc excitation magnetic field |
| CN106503461A (zh) * | 2016-10-26 | 2017-03-15 | 广东产品质量监督检验研究院 | 一种基于深度学习法构建的光伏组件加速退化模型及光伏组件寿命预测方法 |
| CN106600041A (zh) * | 2016-12-02 | 2017-04-26 | 天津大学 | 基于受限玻尔兹曼机的光伏发电短期功率概率预测方法 |
| CN107256437A (zh) * | 2017-05-15 | 2017-10-17 | 内蒙古电力(集团)有限责任公司 | 一种光伏电站超短期辐照度预测方法及系统 |
-
2018
- 2018-04-12 CN CN201810325176.2A patent/CN108776843A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20160252405A1 (en) * | 2014-07-31 | 2016-09-01 | Huazhong University Of Science And Technology | Noninvasive measuring method for rapid temperature variation under dc excitation magnetic field |
| CN106503461A (zh) * | 2016-10-26 | 2017-03-15 | 广东产品质量监督检验研究院 | 一种基于深度学习法构建的光伏组件加速退化模型及光伏组件寿命预测方法 |
| CN106600041A (zh) * | 2016-12-02 | 2017-04-26 | 天津大学 | 基于受限玻尔兹曼机的光伏发电短期功率概率预测方法 |
| CN107256437A (zh) * | 2017-05-15 | 2017-10-17 | 内蒙古电力(集团)有限责任公司 | 一种光伏电站超短期辐照度预测方法及系统 |
Non-Patent Citations (1)
| Title |
|---|
| 路志英 等: "基于深度学习的短时强降水天气识别", 《天津大学学报（自然科学与工程技术版）》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113297791A (zh) * | 2021-05-18 | 2021-08-24 | 四川大川云能科技有限公司 | 一种基于改进dbn的风功率组合预测方法 |
| CN113297791B (zh) * | 2021-05-18 | 2024-02-06 | 四川大川云能科技有限公司 | 一种基于改进dbn的风功率组合预测方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Puri et al. | A hybrid artificial intelligence and internet of things model for generation of renewable resource of energy | |
| CN109214575B (zh) | 一种基于小波长短期记忆网络的超短期风电功率预测方法 | |
| CN111008728B (zh) | 一种用于分布式光伏发电系统短期出力的预测方法 | |
| Tang et al. | Entropy method combined with extreme learning machine method for the short-term photovoltaic power generation forecasting | |
| US20230035108A1 (en) | Photovoltaic cell parameter identification method based on improved equilibrium optimizer algorithm | |
| CN109165774A (zh) | 一种短期光伏功率预测方法 | |
| CN108734331A (zh) | 基于lstm的短期光伏发电功率预测方法及系统 | |
| Guo et al. | Study on short-term photovoltaic power prediction model based on the Stacking ensemble learning | |
| CN110942205B (zh) | 一种基于himvo-svm的短期光伏发电功率预测方法 | |
| CN106251001A (zh) | 一种基于改进模糊聚类算法的光伏功率预测方法 | |
| CN103489038A (zh) | 基于lm-bp神经网络的光伏超短期功率预测方法 | |
| Wang et al. | An improved model combining evolutionary algorithm and neural networks for PV maximum power point tracking | |
| CN107562992B (zh) | 一种基于svm和粒子群算法的光伏列阵最大功率跟踪方法 | |
| CN109086928A (zh) | 基于saga-fcm-lssvm模型的光伏电站实时功率预测方法 | |
| Liao | Genetic k-means algorithm based RBF network for photovoltaic MPP prediction | |
| CN112215428B (zh) | 基于误差修正和模糊逻辑的光伏发电功率预测方法及系统 | |
| CN109272134A (zh) | 一种考虑系统损耗的区域光伏电站并网功率预测方法 | |
| Chen et al. | Research on wind power prediction method based on convolutional neural network and genetic algorithm | |
| Natarajan et al. | Survey on renewable energy forecasting using different techniques | |
| CN116629416A (zh) | 光伏电站功率预测方法及装置 | |
| CN115759389A (zh) | 基于天气类型的相似日组合策略的日前光伏功率预测方法 | |
| CN106405683B (zh) | 基于g-l混合噪声特性核岭回归技术的风速预报方法及装置 | |
| Vohra et al. | End-to-end learning with multiple modalities for system-optimised renewables nowcasting | |
| CN116777039A (zh) | 基于训练集分段和误差修正的双层神经网络风速预测方法 | |
| CN108776843A (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 | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20181109 |