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

xinzhiyong

GM(1,1)灰色模型的程序实现function GM1=fungry1(x0) %输入原始数据x0 0 m2 ?( T* E. v+ @( l) N1 B& IT=input('T=');%从键盘输入从最后一个历史数据算起的第T时点 4 u) J/ o" a& Y+ R9 \x1=zeros(1,length(x0));B=zeros(length(x0)-1,2); & n( L& M3 k9 T7 m8 Fyn=zeros(length(x0)-1,1);Hatx0=zeros(1,length(x0)+T); / a1 Y, O0 [" _$ zHatx00=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   E  p: W5 L0 d        x1(i)=x1(i)+x0(j);     end ; F6 |/ E/ Z' _$ ]1 W2 Zend - l; e" P  w7 l0 f  kfor i=1:length(x0)-1 % C4 M& Q  X- W6 {! R$ }    B(i,1)=(-1/2)*(x1(i)+x1(i+1));     B(i,2)=1;     yn(i)=x0(i+1); end ; a/ T8 f0 W/ ]* U' KHatA=(inv(B'*B))*B'*yn % GM(1,1)模型参数估计 4 z7 z8 m$ f8 u. ?5 @0 Dfor k=1:length(x0)+T     Hatx1(k)=(x0(1)-HatA(2)/HatA(1))*exp(-HatA(1)*(k-1))+HatA(2)/HatA(1); end 1 m. L7 n/ v0 f9 V- Z: h5 e; M; uHatx0(1)=Hatx1(1); for k=2:length(x0)+T     Hatx0(k)=Hatx1(k)-Hatx1(k-1);%累计还原得到历史数据的模拟值 end 2 ~2 h) J/ F6 j  X3 P& X& l% ofor i=1:length(x0) %开始模型检验 % ]+ ?7 S9 `4 o" s    epsilon(i)=x0(i)-Hatx0(i);     omega(i)=(epsilon(i)/x0(i))*100; * m. l6 w; N  _2 a# E1 v8 Jend 9 P; }% B) |! C! u- T% x0;Hatx0;epsilon;omega;  %必要时去掉％得到各种数据 ( `$ A5 I( X( G' lc=std(epsilon)/std(x0);p=0; 8 ?% G! w0 _3 d  j- Rfor i=1:length(x0)     if abs(epsilon(i)-mean(epsilon))<0.6745*std(x0) 1 I6 s1 m- `3 U6 z- b        p=p+1;     end 5 _  {( j" ]* [* E* Aend p=p/length(x0) if p>0.95 & c<0.35 ' \9 ]( V+ i" _( ^! Y    disp('The model is good,and the forecast is:'), : @$ c, q+ U2 ^% `! k0 n3 I) ]    disp(Hatx0(length(x0)+T)) * W& q# [# `7 ~8 O/ r; D& kelseif p>0.85 & c<0.5     disp('The model is eligibility,and the forecast is:'), ' |" V' u0 y  |    disp(Hatx0(length(x0)+T)) elseif p>0.7 & c>0.65     disp('The model is not good,and the forecast is:'), $ l9 ~: l3 W! h' o5 A    disp(Hatx0(length(x0)+T)) else p<=0.7 & c>0.65 ! m; D( T# e8 b4 ~  M7 |# j    disp('The model is bad and try again') end # V9 d! l: A. m; a8 Pfor i=1:length(x0)     Hatx00(i)=Hatx0(i); # r* u" l! q/ c2 wend z=1:length(x0); 8 X! p8 q! Z7 `5 c( R" gplot(z,x0,'-',z,Hatx00,':') %将原始数据和模拟值画在一个图上帮助观察 3 o+ n8 G: K6 h" Q9 Mtext(2,x0(2),'History data: real line') 1 P& y2 U2 i# z8 ptext(length(x0)/2,Hatx00(length(x0))/2,'Simulation data:broken line') endT=input('T=');%从键盘输入从最后一个历史数据算起的第T时点????是指什么啊，请大哥们，大姐们教一下，我急用，请快，谢谢我的初始值x0=[1.620938526 0.07925621 , O2 I9 W% L$ l, ^% m* G4 S0.052318818 & a  ?' @" S& x7 z! O5 I$ W0.041252502 0.021800479 " v7 {! y1 O, I9 p7 p0.053132975 7 x# t$ s/ f( }: a; u0.089908836 . X6 F2 X" }: p) |1 \% \' w0.109153219 0.079331832 0.342192598 2 ]: }; \- y& D, x% @3 t) H0.099718142 ( Q8 X$ n. E1 g7 F0.135194823 0.109274037 . a8 P8 h8 L6 ^' ?0.08152013 0.067876355 0.064706843 3 W+ Y) z) s" V0 m; ?0.055562197 0.050848544 # A' z. a* |5 a' E/ e]';

daiqiang5566

楼主莫急..

yysclshi

a= -0.0080
u= 0.0713
4 T8 O$ h0 m( c/ w  c, H预测值
    1.6209    0.0846    0.0853    0.0859    0.0866    0.0873
    0.0880    0.0887    0.0895    0.0902    0.0909    0.0916
& P( A8 @! V) h/ g4 s7 ?    0.0924    0.0931    0.0938    0.0946    0.0954    0.0961
$ F( j; R+ r6 w; a! e1 W初始值
    1.6209    0.0793    0.0523    0.0413    0.0218    0.0531
    0.0899    0.1092    0.0793    0.3422    0.0997    0.1352
: d: E% c( B$ p5 Y' y) m- x    0.1093    0.0815    0.0679    0.0647    0.0556    0.0508
残差
         0   -0.0053   -0.0329   -0.0447   -0.0648   -0.0342
, _- X" z9 \3 p. N    0.0019    0.0204   -0.0101    0.2520    0.0088    0.0436
    0.0169   -0.0116   -0.0260   -0.0299   -0.0398   -0.0453
相对误差
; w% X* _! d& Q4 T3 t3 a         0    0.0672    0.6297    1.0835    2.9741    0.6437
    0.0209    0.1870    0.1276    0.7365    0.0885    0.3223
/ Z2 [- F8 }/ c) D( y5 A* Z2 U    0.1548    0.1420    0.3826    0.4619    0.7162    0.8903
方差比
, n' `/ ^" d  z9 n: @8 `- t' d" o    0.1869
8 a: B; e/ Z3 m6 \6 r2 `p =
# v. X% ~: p/ E) z     1

yysclshi

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

杨晓敬

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

gxj820

很好啊！！顶顶！

一剑卡卡

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

jshzncd

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