ruonangu
2011-11-27, 13:07
实时预测是滚动地预测每个时点未来4小时内的16个时点（每15分钟一个时点）的风电功率数值现在要实时预测一天的数据。历史数据是之前一天的16个时点。
>> syms a b;
c=[a,b]';
A=[347.0625 303.6563 351.8438 405.9375 303 250.875 270.2813 206.25 170.625 290.625 267.75 203.25 307.6875 279.0938 425.1563 195.1875
];
B=cumsum(A);
n=length(A);
for i=1:(n-1);
C(i)=(B(i)+B(i+1))/2;
end
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);
F=[];
for j=1:112;
F(j)=A(j);
for i=(j+1):(j+n+15);
F(i)=(A(1)-b/a)/exp(a*(i-1)+b/a);
G[];G(j)=A(j);
G(i)=F(i)-F(i-1);
end
end
t1=1:16;t2=1:128)
plot(t1,A,'o',t2,G)
这是一个灰色预测模型。还有更好的方法没有的吗？我们的这个程序是错误的，我自己编的。请指教
>> syms a b;
c=[a,b]';
A=[347.0625 303.6563 351.8438 405.9375 303 250.875 270.2813 206.25 170.625 290.625 267.75 203.25 307.6875 279.0938 425.1563 195.1875
];
B=cumsum(A);
n=length(A);
for i=1:(n-1);
C(i)=(B(i)+B(i+1))/2;
end
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);
F=[];
for j=1:112;
F(j)=A(j);
for i=(j+1):(j+n+15);
F(i)=(A(1)-b/a)/exp(a*(i-1)+b/a);
G[];G(j)=A(j);
G(i)=F(i)-F(i-1);
end
end
t1=1:16;t2=1:128)
plot(t1,A,'o',t2,G)
这是一个灰色预测模型。还有更好的方法没有的吗？我们的这个程序是错误的，我自己编的。请指教