基本灰色负荷预测程序（Matlab）

qiuwei_002 · 发表于 2008-5-29 16:51:34

马上加入，结交更多好友，共享更多资料，让你轻松玩转电力研学社区！

您需要登录才可以下载或查看，没有账号？立即加入

×

function SGrey

X0 = input('请输入原始负荷数据：'); %输入原始数据
n = length(X0);  %原始n年数据

%累加生成
X1 = zeros(1,n);
for i = 1:n
if i == 1
      X1(1,i) = X0(1,i);
else
      X1(1,i) = X0(1,i) + X1(1,i-1);
end
end
X1

%计算数据矩阵B和数据向量Y
B = zeros(n-1,2);
Y = zeros(n-1,1);
for i = 1:n-1
B(i,1) = -0.5*(X1(1,i) + X1(1,i+1));
B(i,2) = 1;
Y(i,1) = X0(1,i+1);
end
B,Y

%计算GM(1,1)微分方程的参数a和u
A = zeros(2,1);
A = inv(B'*B)*B'*Y;
a = A(1,1);
u = A(2,1);
a,u

%建立灰色预测模型
XX0(1,1) = X0(1,1);
for i = 2:n
XX0(1,i) = (X0(1,1) - u/a)*(1-exp(a))*exp(-a*(i-1));
end
XX0
%模型精度的后验差检验
e = 0;       %求残差平均值
for i =1:n
e = e + (X0(1,i) - XX0(1,i));
end
e = e/n;
e
aver = 0;    %求历史数据平均值
for i = 1:n
aver = aver + X0(1,i);
end
aver = aver / n;
aver
s12 = 0;    %求历史数据方差
for i = 1:n
s12 = s12 + (X0(1,i)-aver)^2;
end
s12 = s12 / n;
s12
s22 = 0;    %求残差方差
for i = 1:n
s22 = s22 + ((X0(1,i) - XX0(1,i)) - e)^2;
end
s22 = s22 / n;
s22
C = s22 / s12; %求后验差比值
C
cout = 0; %求小误差概率
for i = 1:n
if abs((X0(1,i) - XX0(1,i)) - e) < 0.6754*sqrt(s12)
      cout = cout+1;
else
      cout = cout;
end
end
P = cout / n;
P
if (C < 0.35 & P > 0.95)
disp('预测精度为一级');
m = input('请输入需要预测的年数: m = '); %预测往后各年的负荷
disp('往后m各年负荷为：');
f = zeros(1,m);
for i = 1:m
      f(1,i) = (X0(1,1) - u/a)*(1-exp(a))*exp(-a*(i+n-1));
end
f
else
disp('灰色预测法不适用');
end

zhao06102642 · 发表于 2016-12-19 16:43:08

多谢楼主，好人啊！

qiuwei_002 · 发表于 2008-7-15 09:35:05

大家有没有改进的灰色预测程序呢？

乔丹 · 发表于 2009-4-3 15:53:11

毕业设计正在做这个，谢谢楼主分享，呵呵！

weilingyou · 发表于 2009-5-12 14:45:00

很好这个

no123456 · 发表于 2009-5-12 16:01:01

谢谢楼主分享！很好这个

snowbird · 发表于 2009-9-10 22:23:43

不错，谢谢楼主

daisyjiangd · 发表于 2010-3-29 11:58:26

谢谢，我正需要。

wangyuan01 · 发表于 2010-4-21 21:51:21

正在研究这个东西，谢谢！

family1420 · 发表于 2010-4-26 15:22:28

谢谢楼主，不知是否和神经网络结合在一起了。

zhpy · 发表于 2010-4-26 15:28:59

灰色预测以前数学建模时用过

账号		自动登录	找回密码
密码			立即加入

基本灰色负荷预测程序（Matlab）

马上加入，结交更多好友，共享更多资料，让你轻松玩转电力研学社区！

相关帖子