CN108090323A - 基于无限混合高斯和样条回归的风电功率曲线拟合方法 - Google Patents
基于无限混合高斯和样条回归的风电功率曲线拟合方法 Download PDFInfo
- Publication number
- CN108090323A CN108090323A CN201810023362.0A CN201810023362A CN108090323A CN 108090323 A CN108090323 A CN 108090323A CN 201810023362 A CN201810023362 A CN 201810023362A CN 108090323 A CN108090323 A CN 108090323A
- Authority
- CN
- China
- Prior art keywords
- power curve
- distribution
- power
- spline
- model
- 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
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000012549 training Methods 0.000 claims abstract description 9
- 239000000203 mixture Substances 0.000 claims description 7
- 238000007781 pre-processing Methods 0.000 claims description 3
- 230000005611 electricity Effects 0.000 abstract description 2
- 238000010276 construction Methods 0.000 abstract 1
- 230000000977 initiatory effect Effects 0.000 abstract 1
- 230000002159 abnormal effect Effects 0.000 description 3
- 230000014509 gene expression Effects 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000013528 artificial neural network Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 238000012706 support-vector machine Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000007418 data mining Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000007637 random forest analysis Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000006467 substitution reaction 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
- 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
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Wind Motors (AREA)
Abstract
本发明公开了一种基于无限混合高斯和样条回归的风电功率曲线拟合方法，包括进行数据预处理、构造鲁棒的样条回归模型、利用变分贝叶斯优化鲁棒的样条回归模型、得出功率曲线以及概率功率曲线的步骤，本发明的风电功率曲线拟合方法只需要设置一些初始化参数，方法简单，能够容忍训练数据中存在一些非一致样本，精度高、误差小，可以进一步提高风电预报的精度。
Description
技术领域
本发明涉及新能源领域和机器学习领域，特别涉及一种基于无限混合高斯和样条回归的风电功率曲线拟合方法。
背景技术
当前，随着全球的能源危机的出现以及使用传统能源如煤，石油等带来的日趋严峻的环境问题，找到替代原有的传统能源的任务越来越紧迫。风电作为一种清洁、可再生的能源受到了越来越多的关注。大规模的风电并网将会在一定程度上缓解能源危机，并且能带来经济效益和减少缓解污染。然而，由于风电本身的随机性和间歇性，从而导致在大规模风电并网后对整个电力系统的完全性和稳定性产生巨大影响。因此，从风电并网的角度来说，准确的功率预报是很有必要的，此外还可以减少电力系统的运行成本。
由于风速和功率之间的关系可以用功率曲线来表示。因此风电预报常用两步法，第一步先得到风速的预报值，第二步是利用功率曲线得到功率的预报。通常情况下，功率曲线是由风机的生产厂商提供的。然而，该功率曲线是一条理论的功率曲线，没有考虑环境等因素(如气温和湿度等)的影响。在实际中，随着地理环境和气候的不同，实际功率曲线会产生变化。因此，直接利用理论功率曲线来实现功率预报会带来额外的预报误差。因此，很多学者也研究如何去得到一个准确的实际功率曲线来进一步提高风电预报的精度。除此以外，准确的功率曲线还能够在线监测风机的运行状态和减少风机的运行和维护成本等。
目前，可以将风电功率曲线建模技术分成两大类：参数模型和非参数模型。一般地，参数模型主要由带有几个参数的数学表达式构成，主要包括线性分割法、多项式功率曲线、理想功率曲线、概率模型、动态功率曲线、4参数、5参数logistics模型以及修正的双曲正切模型等。对于多项式模型而言，常用三次功率曲线、二次功率曲线、6阶以及9阶多项式模型来拟合功率曲线。参数模型的缺点是在描述功率曲线的动态特性上其性能有限。与参数模型不同，非参数模型不需要利用一些数学表达式也不需要对功率曲线的形状有先验知识，仅仅利用历史的功率以及风速数据就可以拟合出各种各样的功率曲线。非参数模型主要包括样条回归，人工神经网络法，模糊法和一些数据挖掘的方法如支持向量机，随机森林和K近邻等。尽管非参数模型比参数模型灵活，但是其计算量也相对较大。
除了功率曲线的建模方法外，另一个影响我们获得一个精确的功率曲线的因素是数据的质量。然而，在实际中，获得的风速与功率的数据中经常会存在很多不一致的样本。造成以上现象的原因包括传感器误差、停机维护、弃风限电以及环境因素如结冰等。不一致样本的特性是在给定一个风速，真实的功率数据远离功率曲线。为了提高数据的质量，通常用提前去除这些不一致样本，然后再利用处理后的数据构建功率曲线模型。然而，此类方法的缺陷在于我们无法保证所有的不一致样本都能被检测到。
在数据中存在一些非一致样本时，由于实际功率与功率曲线上的功率距离较远，因此此处的误差较大，误差分布呈现出非高斯特性，具有长尾现象。此种情况下高斯分布很难描述以上特性的误差分布。然而，当前的一些模型如样条回归、多项式模型等都假设误差服从高斯分布，在训练数据中存在一些非一致样本时，真实的误差分布特性与假设的误差分布之间并不一致。因次，在存在非一致样本的情况下，假设误差服从高斯分布的功率曲线模型并不合适。
发明内容
本发明的目的是解决现有风电功率曲线的精度较低误差较大的技术问题，本发明提供一种基于无限混合高斯和样条回归的风电功率曲线拟合方法。
本发明解决技术问题采用如下技术方案：
一种基于无限混合高斯和样条回归的风电功率曲线拟合方法，包括以下步骤：
1)数据预处理：
根据实际的风速和功率数据画出经验功率曲线，然后将明显的异常数据去除，处理后的样本表示为所述的xi，yi分别表示风速和功率值，N表示训练样本的长度；
2)构造鲁棒的样条回归模型：
构造鲁棒样条回归模型yi＝Zi(xi)β+ei，所述的Zi(xi)表示根据样条基计算出的输入向量，β表示回归系数，ei表示服从无限混合高斯模型的回归误差：
所述的表示第k个高斯分布的方差；所述的回归系数β服从一个高斯分布；所述的πk表示第k个高斯分布的权值，且πk是关于的一个函数，表示为
所述的参数是一个变量，满足一个参数为v的Beta分布超参数v为满足参数为e0、f0的Gamma分布，v＝Gamma(v|e0，f0)，所述的e0、f0的初始值设置为0.0001并根据变分贝叶斯得到的参数的后验分布来更新e0、f0的值；
3)利用变分贝叶斯优化鲁棒的样条回归模型：
根据步骤2)中对鲁棒样条回归模型中各个参数的先验分布，构造最终的似然函数：p(R，β，w，ρ，τ，v|Y，Z)＝p(Y|R，Z，β，τ)p(τ)p(R|w)p(β|ρ)p(ρ)p(w|v)p(v)，所述p(.)是变量的概率分布；
令Θ＝{R，β，w，ρ，τ，v}，根据变分贝叶斯的原理，求出鲁棒样条回归模型中的所有参数的后验分布,
所述的<·>是期望运算；
4)得出功率曲线以及概率功率曲线：
根据步骤3)得到的参数的后验分布，给定一个新样本(x*，y*)，所述x*、y*分别表示未知风速值和待预测的真实功率，通过下式来推断y*的概率分布：
所述μβ,Σβ分别表示参数β的后验分布中的期望和方差，Z*表示的是根据样条基计算的与风速x*相对应的输入向量，y*的预测分布是混合高斯模型，所述的Z*μβ是y*的点预测值；
根据y＝Z(x)μβ得到预测的功率曲线；
根据每个预测值所相应的概率分布，求出相应的概率功率曲线。
本发明具有如下有益效果：本发明的风电功率曲线拟合方法由于采用了无限混合高斯分布，理论上能拟合任意复杂的分布，因此能够容忍训练数据中存在一些非一致样本，而无需采用一些其他的算法将非一致样本全部去除；采用变分贝叶斯的优化方法，不仅能得到确定的功率曲线，也可以得到概率功率曲线，进而可以利用概率功率曲线去识别出数据中的非一致样本；本发明的风电功率曲线拟合方法只需要设置一些初始化参数即可，没有任何待优化的参数，方法简单、精度高、误差小，可进一步提高风电预报的精度。
附图说明
图1为本发明基于无限混合高斯和样条回归的风电功率曲线拟合方法的流程图；
图2为本发明实施例中的原始的风速与功率数据；
图3为本发明实施例中不同模型得到的数据集C的功率曲线图；
图4为本发明实施例中不同数据集的概率功率曲线。
具体实施方式
下面结合实施方式和实施例对本发明的技术方案作进一步阐述。
具体实施方式：本实施方式是一种基于无限混合高斯和样条回归的风电功率曲线拟合方法，如图1所示，具体步骤如下：
1)数据预处理：
根据实际的风速和功率数据画出经验功率曲线，然后将明显的异常数据去除，处理后的样本表示为所述的xi，yi分别表示风速和功率值，N表示训练样本的长度；
2)构造鲁棒的样条回归模型：
构造鲁棒样条回归模型yi＝Zi(xi)β+ei，所述的Zi(xi)表示根据样条基计算出的输入向量，β表示回归系数，ei表示服从无限混合高斯模型的回归误差：
所述的表示第k个高斯分布的方差；所述的回归系数β服从一个高斯分布；所述的πk表示第k个高斯分布的权值，且πk是关于的一个函数，表示为
所述的参数是一个变量，满足一个参数为v的Beta分布超参数v为满足参数为e0、f0的Gamma分布，v＝Gamma(v|e0，f0)，所述的e0、f0的初始值设置为0.0001并根据变分贝叶斯得到的参数的后验分布来更新e0、f0的值；
3)利用变分贝叶斯优化鲁棒的样条回归模型：
根据步骤2)中对鲁棒样条回归模型中各个参数的先验分布，构造最终的似然函数：p(R，β，w，ρ，τ，v|Y，Z)＝p(Y|R，Z，β，τ)p(τ)p(R|w)p(β|ρ)p(ρ)p(w|v)p(v)，所述p(.)是变量的概率分布；
令Θ＝{R，β，w，ρ，τ，v}，根据变分贝叶斯的原理，求出鲁棒样条回归模型中的所有参数的后验分布,
所述的<·>是期望运算；
4)得出功率曲线以及概率功率曲线：
根据步骤3)得到的参数的后验分布，给定一个新样本(x*，y*)，所述x*、y*分别表示未知风速值和待预测的真实功率，通过下式来推断y*的概率分布：
所述μβ,Σβ分别表示参数β的后验分布中的期望和方差，Z*表示的是根据样条基计算的与风速x*相对应的输入向量，y*的预测分布是混合高斯模型，所述的Z*μβ是y*的点预测值；
根据y＝Z(x)μβ得到预测的功率曲线；
根据每个预测值所相应的概率分布，求出相应的概率功率曲线。
下面通过具体实施例，对本发明做进一步详细说明，应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
实施例采用以下三个数据集：数据集A来自于宁夏风电场，其中包含10000对样本(风速S和相应的风功率P)。在进行实验之前，我们将数据集分成2部分，前8000个样本作为训练样本，剩余的2000个样本作为测试样本。另外两个数据集(B和C)是来自同一风电场的两个不同风机上的数据。数据集B中包含6000对样本，其中前5000作为训练样本，剩余的1000样本作为测试样本。数据集C中包含7500对样本，前6000作为训练样本，剩余的1500样本作为测试样本。以上三个数据集的的样本采样频率为10min，原始风速和功率数据如图2所示。
在实际中，为了降低非一致样本对功率曲线建模的影响，通常情况下需要将这些非一致样本点剔除。然而，由于当前用于识别非一致样本的算法能力限制，我们并不能确保所有的非一致样本点都能剔除。为了模拟此种情况，我们首先采用简单的方法剔除原始数据中的非一致样本点。
在数据集A中，风速大于10m/s且功率小于14000kW时的样本被认为是非一致样本，且直接去除。在数据集B中，当风速大于6m/s且功率小于50kW的样本，风速大于13m/s功率小于700kW的样本，都被视为非一致样本且被去除。在数据集C中，风速大于6m/s功率小于50，风速大于12m/s功率小于600kW中的样本被去除。上述简单的方法只能去除较为明显的非一致样本，因此处理后的数据中仍然存在一些异常点。
为了说明本发明所提出的方法的有效性，我们对比以下功率曲线拟合方法常用的模型：四种参数模型(4参数logistic模型、5参数logistic模型、6阶和9阶多项式回归模型)，一个非参数模型(样条回归模型)和两个基于学习的模型(神经网络法和支持向量机)。为了表述方便，上述对比模型分别简写为4-PLM、5-PLM、6-PRM、9-PRM、SRM、ANN和SVM。本发明采用的模型简写为RSRM。为了定量的刻画各个模型的性能，一般采用以下两个指标，即平均绝对误差(MAE)和均方根误差(RMSE)，计算公式为
其中yi，分别表示为真实的功率和预测的功率。
本发明所提出的基于鲁棒回归模型的功率曲线拟合方法能容忍一些非一致样本点在样本中。不同模型的功率曲线拟合方法的拟合结果如表1所示，利用不同模型的功率曲线方法得到的数据集C的功率曲线如图3所示。
表1基于各种模型的功率曲线拟合方法的性能对比
从表1可以看出，四个参数模型的性能较为有限，两个基于学习的模型(MLP和SVM)的性能大体上要优于参数模型。本发明所提出的方法RSRM在数据中存在非一致样本时得到的功率曲线要好于参数模型和基于学习的模型。从图3也可以看到RSRM的功率曲线能够更好的拟合原始的功率数据。
本发明提出的方法不仅能提供确定性的功率曲线，也能够提供概率功率曲线。我们可以利用概率功率曲线去识别原始数据中的非一致样本点。用前面所述的简单方法从数据集A、B和C中识别出的非一致样本点分别为144、309和931个。在95％的置信度下，用本发明提出的方法能从数据集A、B和C中识别出445、401和1102个非一致样本，如图4中不同数据集的概率功率曲线所示。从上述结果可以看出，本发明提出的方法能够识别出更多的非一致样本。
最后应说明的是：以上实施方式和实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施方式和实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述实施方式和实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明实施方式和实施例技术方案的精神和范围。
Claims (1)
1.一种基于无限混合高斯和样条回归的风电功率曲线拟合方法，其特征在于，包括以下步骤：
1)数据预处理：
根据实际的风速和功率数据画出经验功率曲线，然后将明显的异常数据去除，处理后的样本表示为所述的xi，yi分别表示风速和功率值，N表示训练样本的长度；
2)构造鲁棒的样条回归模型：
构造鲁棒样条回归模型yi＝Zi(xi)β+ei，所述的Zi(xi)表示根据样条基计算出的输入向量，β表示回归系数，ei表示服从无限混合高斯模型的回归误差：
<mrow>
<msub>
<mi>e</mi>
<mi>t</mi>
</msub>
<mo>~</mo>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</msubsup>
<msub>
<mi>&pi;</mi>
<mi>k</mi>
</msub>
<mi>N</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>e</mi>
<mi>t</mi>
</msub>
<mo>|</mo>
<mn>0</mn>
<mo>,</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>k</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
所述的表示第k个高斯分布的方差；所述的回归系数β服从一个高斯分布；所述的πk表示第k个高斯分布的权值，且πk是关于的一个函数，表示为
所述的参数是一个变量，满足一个参数为v的Beta分布超参数v为满足参数为e0、f0的Gamma分布，v＝Gamma{v|e0，f0)，所述的e0、f0的初始值设置为0.0001并根据变分贝叶斯得到的参数的后验分布来更新e0、f0的值；
3)利用变分贝叶斯优化鲁棒的样条回归模型：
根据步骤2)中对鲁棒样条回归模型中各个参数的先验分布，构造最终的似然函数：p(R，β，w，ρ，τ，ν|Y，Z)＝p(Y|R，Z，β，τ)p(τ)p(R|w)p(β|ρ)p(ρ)p(w|v)p(v)，所述p(.)是变量的概率分布；
令Θ＝{R，β，w，ρ，τ，ν}，根据变分贝叶斯的原理，求出鲁棒样条回
<mrow>
<mi>Q</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&Theta;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>exp</mi>
<mo>{</mo>
<mo><</mo>
<mi>ln</mi>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>&Theta;</mi>
<mo>|</mo>
<mi>Y</mi>
<mo>,</mo>
<mi>Z</mi>
<mo>)</mo>
</mrow>
<msub>
<mo>></mo>
<mrow>
<mo>\</mo>
<msub>
<mi>&Theta;</mi>
<mrow>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</msub>
<mo>}</mo>
</mrow>
<mrow>
<mo>&Integral;</mo>
<mi>exp</mi>
<mo>{</mo>
<mo><</mo>
<mi>ln</mi>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>&Theta;</mi>
<mo>|</mo>
<mi>Y</mi>
<mo>,</mo>
<mi>Z</mi>
<mo>)</mo>
</mrow>
<msub>
<mo>></mo>
<mrow>
<mo>\</mo>
<msub>
<mi>&Theta;</mi>
<mrow>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</msub>
<mo>}</mo>
<msub>
<mi>d&Theta;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
</mrow>
归模型中的所有参数的后验分布，所述的<·>是期望运算；
4)得出功率曲线以及概率功率曲线：
根据步骤3)得到的参数的后验分布，给定一个新样本(x*，y*)，所述x*、y*分别表示未知风速值和待预测的真实功率，通过下式来
<mrow>
<mo>&Integral;</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<mo><</mo>
<msub>
<mi>&pi;</mi>
<mi>k</mi>
</msub>
<mo>></mo>
<mi>N</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>y</mi>
<mo>*</mo>
</msup>
<mo>|</mo>
<msup>
<mi>Z</mi>
<mo>*</mo>
</msup>
<mi>&beta;</mi>
<mo>,</mo>
<mo><</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>k</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>></mo>
<mo>)</mo>
</mrow>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>&beta;</mi>
<mo>|</mo>
<msub>
<mi>&mu;</mi>
<mi>&beta;</mi>
</msub>
<mo>,</mo>
<msub>
<mo>&Sigma;</mo>
<mi>&beta;</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>&beta;</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<mo><</mo>
<msub>
<mi>&pi;</mi>
<mi>k</mi>
</msub>
<mo>></mo>
<mi>N</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>y</mi>
<mo>*</mo>
</msup>
<mo>|</mo>
<msup>
<mi>Z</mi>
<mo>*</mo>
</msup>
<msub>
<mi>&mu;</mi>
<mi>&beta;</mi>
</msub>
<mo>,</mo>
<mo><</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>k</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>></mo>
<mo>+</mo>
<msup>
<mi>Z</mi>
<mo>*</mo>
</msup>
<msub>
<mo>&Sigma;</mo>
<mi>&beta;</mi>
</msub>
<msup>
<mi>Z</mi>
<mrow>
<mo>*</mo>
<mi>T</mi>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
推断y*的概率分布：
所述μβ,Σβ分别表示参数β的后验分布中的期望和方差，Z*表示的是根据样条基计算的与风速x*相对应的输入向量，y*的预测分布是混合高斯模型，所述的Z*μβ是y*的点预测值；
根据y＝Z(x)μβ得到预测的功率曲线；
根据每个预测值所相应的概率分布，求出相应的概率功率曲线。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810023362.0A CN108090323B (zh) | 2018-01-10 | 2018-01-10 | 基于无限混合高斯和样条回归的风电功率曲线拟合方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810023362.0A CN108090323B (zh) | 2018-01-10 | 2018-01-10 | 基于无限混合高斯和样条回归的风电功率曲线拟合方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108090323A true CN108090323A (zh) | 2018-05-29 |
| CN108090323B CN108090323B (zh) | 2020-10-09 |
Family
ID=62181932
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810023362.0A Active CN108090323B (zh) | 2018-01-10 | 2018-01-10 | 基于无限混合高斯和样条回归的风电功率曲线拟合方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108090323B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111027229A (zh) * | 2019-12-26 | 2020-04-17 | 中南大学 | 基于稀疏异方差多样条回归的风功率曲线拟合的方法 |
| CN117688793A (zh) * | 2024-02-04 | 2024-03-12 | 中国地质大学（武汉） | 一种分布鲁棒机组组合建模与求解方法、设备及存储设备 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103345577A (zh) * | 2013-06-27 | 2013-10-09 | 江南大学 | 变分贝叶斯概率假设密度多目标跟踪方法 |
| KR20140100783A (ko) * | 2013-02-07 | 2014-08-18 | 경희대학교 산학협력단 | 비모수적 베이지언 모션 인식 방법 및 그 장치 |
| CN104008574A (zh) * | 2014-06-16 | 2014-08-27 | 浙江大学 | 一种基于无限高斯混合模型的高光图图像解混方法 |
| WO2015054295A1 (en) * | 2013-10-07 | 2015-04-16 | Arizona Board Of Regents On Behalf Of Arizona State University | Kidney glomeruli measurement systems and methods |
| CN106683122A (zh) * | 2016-12-16 | 2017-05-17 | 华南理工大学 | 一种基于高斯混合模型和变分贝叶斯的粒子滤波方法 |
-
2018
- 2018-01-10 CN CN201810023362.0A patent/CN108090323B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20140100783A (ko) * | 2013-02-07 | 2014-08-18 | 경희대학교 산학협력단 | 비모수적 베이지언 모션 인식 방법 및 그 장치 |
| CN103345577A (zh) * | 2013-06-27 | 2013-10-09 | 江南大学 | 变分贝叶斯概率假设密度多目标跟踪方法 |
| WO2015054295A1 (en) * | 2013-10-07 | 2015-04-16 | Arizona Board Of Regents On Behalf Of Arizona State University | Kidney glomeruli measurement systems and methods |
| CN104008574A (zh) * | 2014-06-16 | 2014-08-27 | 浙江大学 | 一种基于无限高斯混合模型的高光图图像解混方法 |
| CN106683122A (zh) * | 2016-12-16 | 2017-05-17 | 华南理工大学 | 一种基于高斯混合模型和变分贝叶斯的粒子滤波方法 |
Non-Patent Citations (10)
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111027229A (zh) * | 2019-12-26 | 2020-04-17 | 中南大学 | 基于稀疏异方差多样条回归的风功率曲线拟合的方法 |
| CN111027229B (zh) * | 2019-12-26 | 2021-12-07 | 中南大学 | 基于稀疏异方差多样条回归的风功率曲线拟合的方法 |
| CN117688793A (zh) * | 2024-02-04 | 2024-03-12 | 中国地质大学（武汉） | 一种分布鲁棒机组组合建模与求解方法、设备及存储设备 |
| CN117688793B (zh) * | 2024-02-04 | 2024-05-10 | 中国地质大学（武汉） | 一种分布鲁棒机组组合建模与求解方法、设备及存储设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108090323B (zh) | 2020-10-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108022019B (zh) | 基于风电机组分类的风电场调度方法和系统 | |
| CN111881617A (zh) | 数据处理方法、风力发电机组的性能评估方法及系统 | |
| CN111190349A (zh) | 船舶机舱设备状态监测及故障诊断方法、系统及介质 | |
| CN108304350B (zh) | 基于大数据集近邻策略的风机指标预测与故障预警方法 | |
| CN108090323B (zh) | 基于无限混合高斯和样条回归的风电功率曲线拟合方法 | |
| CN112669168B (zh) | 一种风电功率短期预测方法 | |
| CN111091236A (zh) | 一种按桨距角分类的多分类深度学习短期风功率预测方法 | |
| CN111695724A (zh) | 一种基于混合神经网络模型的风速预测方法 | |
| CN103984986B (zh) | 实时校正的自学习arma模型风电功率超短期预测方法 | |
| CN111340307B (zh) | 预测风机风力发电功率的方法以及相关装置 | |
| CN110991701A (zh) | 一种基于数据融合的风电场风机风速预测方法及系统 | |
| CN111027229B (zh) | 基于稀疏异方差多样条回归的风功率曲线拟合的方法 | |
| Pan et al. | Probabilistic Short‐Term Wind Power Forecasting Using Sparse Bayesian Learning and NWP | |
| CN114119289A (zh) | 综合能源监控数据的处理方法及装置 | |
| CN107527093B (zh) | 一种风电机组运行状态诊断方法及装置 | |
| CN109190845A (zh) | 一种两阶段动态优化短期风电功率预测方法 | |
| CN116885711A (zh) | 一种风电功率预测方法、装置、设备和可读存储介质 | |
| CN109872248B (zh) | 一种风电场集群出力计算方法及系统 | |
| CN113761692B (zh) | 一种基于迁移成分分析的多风电机组运行状态辨识方法 | |
| CN113048012B (zh) | 基于混合高斯模型的风电机组偏航角度识别方法和装置 | |
| CN116050072A (zh) | 一种基于随机采样一致性的风电机组理论功率曲线识别方法和装置 | |
| CN115510961A (zh) | 一种基于主动学习的社区综合能源系统运行安全评估方法 | |
| CN114239920A (zh) | 一种风电功率短期预测方法 | |
| CN112070320A (zh) | 一种基于动态谐波回归的超短期风电功率预测方法及系统 | |
| CN117117858B (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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |