GM(1,1)预测模型的MATLAB程序求助，急！！！

xinzhiyong

GM(1,1)灰色模型的程序实现function GM1=fungry1(x0) %输入原始数据x0 T=input('T=');%从键盘输入从最后一个历史数据算起的第T时点 x1=zeros(1,length(x0));B=zeros(length(x0)-1,2); # W! y# ?& g% A. k' Q( \. kyn=zeros(length(x0)-1,1);Hatx0=zeros(1,length(x0)+T); - f9 @' x, L) O! aHatx00=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) ; F5 p" k  m! d' Y- _    for j=1:i ' @# g3 G# a( c0 `2 n% s        x1(i)=x1(i)+x0(j);     end end   V0 z5 c1 l, U6 c/ _7 G/ x1 mfor i=1:length(x0)-1     B(i,1)=(-1/2)*(x1(i)+x1(i+1));     B(i,2)=1;     yn(i)=x0(i+1); 8 e8 h  _/ \/ w6 r5 uend 5 L3 F9 c: g# U' s! |- sHatA=(inv(B'*B))*B'*yn % GM(1,1)模型参数估计 . j# f0 V3 B. k) xfor k=1:length(x0)+T ; i6 r. f& V# Q. w6 ]  |0 t    Hatx1(k)=(x0(1)-HatA(2)/HatA(1))*exp(-HatA(1)*(k-1))+HatA(2)/HatA(1); end 0 e; {& D/ B( b& P6 f  r: }) tHatx0(1)=Hatx1(1); for k=2:length(x0)+T - I6 T. r+ q% @6 \7 D    Hatx0(k)=Hatx1(k)-Hatx1(k-1);%累计还原得到历史数据的模拟值 % t2 g6 Y9 A( n! @" e/ v0 n8 Xend 3 o8 u& H$ t. i9 ~& x0 W6 Cfor i=1:length(x0) %开始模型检验     epsilon(i)=x0(i)-Hatx0(i);     omega(i)=(epsilon(i)/x0(i))*100; : [: |* f- {' I: C+ Q: Kend 0 Z3 b2 @% \$ X  o$ t% 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) 1 N; j5 L3 ~9 |        p=p+1;     end end ' h: E5 B7 B! O' c: Lp=p/length(x0) if p>0.95 & c<0.35     disp('The model is good,and the forecast is:'),     disp(Hatx0(length(x0)+T)) . f' P5 q6 D- U5 \- N, t/ @0 b3 Felseif p>0.85 & c<0.5 / {; I' D0 u  q$ l  k    disp('The model is eligibility,and the forecast is:'),     disp(Hatx0(length(x0)+T)) 4 X) C. [& y" i. ^. @/ W8 a8 lelseif p>0.7 & c>0.65     disp('The model is not good,and the forecast is:'), 7 s  u. _& E( I8 K8 h+ `    disp(Hatx0(length(x0)+T)) else p<=0.7 & c>0.65 ) L" d3 u# F7 W% [5 t4 b. C    disp('The model is bad and try again') * h* A+ i2 g" y9 m% S& C6 Kend for i=1:length(x0) - I! ]* d- s5 H5 D# B, I    Hatx00(i)=Hatx0(i); end # G8 U* s# W8 U$ K- a0 L. _z=1:length(x0); . W9 L8 _. ^1 |  m) D* _' Aplot(z,x0,'-',z,Hatx00,':') %将原始数据和模拟值画在一个图上帮助观察 , K0 j, S- o" S6 V7 x" |text(2,x0(2),'History data: real line') 4 j/ w# z, z& Y1 E2 Ktext(length(x0)/2,Hatx00(length(x0))/2,'Simulation data:broken line') endT=input('T=');%从键盘输入从最后一个历史数据算起的第T时点????是指什么啊，请大哥们，大姐们教一下，我急用，请快，谢谢我的初始值x0=[1.620938526 0.07925621 + [2 @' G) @( r1 w; Y2 y, G0.052318818 0.041252502 , [8 y# B' A* ~+ r9 ^6 V3 I- O: B8 g  g0.021800479 0.053132975 5 k% A/ C* B" N# d' M0.089908836 0.109153219 ) Z& _- a3 @* j, B% e8 [. w' a) \0.079331832 0.342192598 0.099718142 0.135194823 ) q% r$ v/ n! q6 O. ?, g3 i6 H0.109274037 ; K7 h- _$ B% w2 q0 O0.08152013 0.067876355 # M- q3 ]  M( w; \& S9 m$ U( z0.064706843 1 E8 m: A5 I0 |" F$ V0.055562197 1 j' o( R8 j* L1 L- h  M0.050848544 ]';

daiqiang5566

楼主莫急..

yysclshi

a= -0.0080
u= 0.0713
1 J1 k. V, x* G7 @8 Z/ \- Y预测值
2 y2 u: B# V# O    1.6209    0.0846    0.0853    0.0859    0.0866    0.0873
    0.0880    0.0887    0.0895    0.0902    0.0909    0.0916
    0.0924    0.0931    0.0938    0.0946    0.0954    0.0961
初始值
# V: R' P. b" E& S' r0 F- m    1.6209    0.0793    0.0523    0.0413    0.0218    0.0531
    0.0899    0.1092    0.0793    0.3422    0.0997    0.1352
5 o* Q2 U; e+ ^" {, i) S; [4 o    0.1093    0.0815    0.0679    0.0647    0.0556    0.0508
) i0 v- V/ \8 F" ~& J  d# }残差
         0   -0.0053   -0.0329   -0.0447   -0.0648   -0.0342
    0.0019    0.0204   -0.0101    0.2520    0.0088    0.0436
! M& J- g  B& B  {    0.0169   -0.0116   -0.0260   -0.0299   -0.0398   -0.0453
2 _: s( Q8 u8 F相对误差
# h( t2 T9 B3 X' m- P+ m, g, l9 N         0    0.0672    0.6297    1.0835    2.9741    0.6437
6 Z* Q' M" O0 a  M, [, m    0.0209    0.1870    0.1276    0.7365    0.0885    0.3223
* |1 ^9 I- d% l2 p5 Z% \7 k) Y" e    0.1548    0.1420    0.3826    0.4619    0.7162    0.8903
方差比
3 n- e" x0 j; B9 x( ]4 A5 [    0.1869
p =
     1

yysclshi

这个图片我不会发，真难为情

杨晓敬

不好意思，不会啊！帮不上忙啊

gxj820

很好啊！！顶顶！

一剑卡卡

不懂。。。。。。。。。。。。。。

jshzncd

不懂啊，楼主~不好意思是哈