数学建模社区-数学中国

标题: GM(1,1)预测模型的MATLAB程序求助，急！！！ [打印本页]

作者: xinzhiyong 时间: 2009-8-27 23:08
标题: GM(1,1)预测模型的MATLAB程序求助，急！！！
GM(1,1)灰色模型的程序实现function GM1=fungry1(x0) %输入原始数据x0
T=input('T=');%从键盘输入从最后一个历史数据算起的第T时点
x1=zeros(1,length(x0));B=zeros(length(x0)-1,2);
yn=zeros(length(x0)-1,1);Hatx0=zeros(1,length(x0)+T);
Hatx00=zeros(1,length(x0));Hatx1=zeros(1,length(x0)+T);
epsilon=zeros(length(x0),1);omega=zeros(length(x0),1);
for i=1:length(x0)
for j=1:i
      x1(i)=x1(i)+x0(j);
end
end
for i=1:length(x0)-1
B(i,1)=(-1/2)*(x1(i)+x1(i+1));
B(i,2)=1;
yn(i)=x0(i+1);
end
HatA=(inv(B'*B))*B'*yn % GM(1,1)模型参数估计
for k=1:length(x0)+T
Hatx1(k)=(x0(1)-HatA(2)/HatA(1))*exp(-HatA(1)*(k-1))+HatA(2)/HatA(1);
end
Hatx0(1)=Hatx1(1);
for k=2:length(x0)+T
Hatx0(k)=Hatx1(k)-Hatx1(k-1);%累计还原得到历史数据的模拟值
end
for i=1:length(x0) %开始模型检验
epsilon(i)=x0(i)-Hatx0(i);
omega(i)=(epsilon(i)/x0(i))*100;
end
% x0;Hatx0;epsilon;omega;  %必要时去掉％得到各种数据
c=std(epsilon)/std(x0);p=0;
for i=1:length(x0)
if abs(epsilon(i)-mean(epsilon))<0.6745*std(x0)
      p=p+1;
end
end
p=p/length(x0)
if p>0.95 & c<0.35
disp('The model is good,and the forecast is:'),
disp(Hatx0(length(x0)+T))
elseif p>0.85 & c<0.5
disp('The model is eligibility,and the forecast is:'),
disp(Hatx0(length(x0)+T))
elseif p>0.7 & c>0.65
disp('The model is not good,and the forecast is:'),
disp(Hatx0(length(x0)+T))
else p<=0.7 & c>0.65
disp('The model is bad and try again')
end
for i=1:length(x0)
Hatx00(i)=Hatx0(i);
end
z=1:length(x0);
plot(z,x0,'-',z,Hatx00,':') %将原始数据和模拟值画在一个图上帮助观察
text(2,x0(2),'History data: real line')
text(length(x0)/2,Hatx00(length(x0))/2,'Simulation data:broken line')
endT=input('T=');%从键盘输入从最后一个历史数据算起的第T时点????是指什么啊，请大哥们，大姐们教一下，我急用，请快，谢谢我的初始值x0=[1.620938526
0.07925621
0.052318818
0.041252502
0.021800479
0.053132975
0.089908836
0.109153219
0.079331832
0.342192598
0.099718142
0.135194823
0.109274037
0.08152013
0.067876355
0.064706843
0.055562197
0.050848544
]';

作者: starbinbin 时间: 2010-3-10 20:57
楼主什么意思啊？？？？？？？？？

作者: 数学者 时间: 2010-3-11 23:03
就是预测期的长度的意思嘛，比如历史数据是由2000年到2007年的，T=3。则预测2008，2009，2010这三年的数据~

作者: xuchangshou0604 时间: 2010-4-11 10:18
好东西，多谢楼主分享。。。。。

作者: alair009 时间: 2012-1-26 09:32
(*^__^*) 指点系词……激扬文字……5129454739650846259430356733215641124313647551654634064881650094230240059449553

作者: 风过也 时间: 2012-5-10 22:57
程序我就先收下了

作者: Netbird_China 时间: 2012-5-14 23:20

作者: lalala263071405 时间: 2016-8-31 08:30
非常好的代码，帮了我大忙，感谢

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)