原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
2011电工杯建模竞赛——用灰色理论预测风电功率
2011“电工杯”建模竞赛问题一
——运用灰色理论对风电功率进行实时预测
因为风具有间歇性的特点，而且短时间内风速具有一定的稳定性，我们取待测时点前短期的数据为历史数据，采用灰色预测GM（1，1）模型给出下一时段的风电功率。整个预测的输入、求解与输出全过程均可通过编程实现，具体程序代码见附录。
GM（1，1）是一阶的微分方程，设有风电功率的原始时间序列（待输入）：
，预测得到的风电功率记为。
对原始时间序列进行累加生成递推数列：
其中，由灰色模型理论，即可组建微分方程：
其中，a与u，并对其求解可得计算式：
由于a和u是通过最小二乘法求出的近似值，故是一个近似表达式，符号上将其区别于。
对函数表达式及进行离散，并将两者做差以便还原原序列，得到近似数据序列
通过以上的模型可以得到一个新的预测序列：
基于以上的预测原理，我们将待测时点前4小时（共16个时点）内的风电功率构成原始序列，用以上步骤得到实时预测所要求的，未来4小时内16时点的预测风电功率。
这是其中一次预测产生的图像。
图中横坐标为时间（简化为时点，函数有所平移），纵坐标为风电功率值。散点为原始数据，曲线为序列按时间描点后连线得到的平滑曲线。图中显示，在原始数据的时间区间内，预测曲线与实际相符。未来四小时内的风电功率由程序自动生成。
若想求出5月31日到6月6日之间所有时点的实时预测，工作量过于庞大。由于部分可以反映整体，我们同样抽样预测。
我们抽取的时间点为：
5月31日的第20个时点，6月2日的第40个时点，6月4日的第60个时点，6月6日的第80个时点。分别对四个对象进行灰色预测，四台机组的总功率由求和即可生成。再用灰色理论预测58台机组的总功率。
求解过程中相关MATLAB程序代码如下：
灰色预测的程序
clear
syms a b;
c=[a b];
A=[323.8125 338.625 411.9375 364.6875 405.375 496.2188 375.6563 380.8125 410.625 106.9688 12 -3.5625 10.78125 -1.875 -2.25 -3.5625]; 注：矩阵A代表原始序列
B=cumsum(A);
n=length(A);
for i=1:(n-1)
C(i)=(B(i)+B(i+1))/2;
end
%calculate the diameters
D=A;D(1)=[];
D=D;
E=[-C;ones(1,n-1)];
c=inv(E*E)*E*D;
c=c;
a=c(1);b=c(2);
%predict
F=[];F(1)=A(1);
for i=2:(n+16)
F(i)=(A(1)-b/a)/exp(a*(i-1))+b/a;
end
G=[];G(1)=A(1);
for i=2:(n+16)
G(i)=F(i)-F(i-1);
end
t1=-15:0;
t2=-15:16;
G 注：矩阵G代表预测序列
plot(t1,A,o,t2,G)
文档评论（0）