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

xinzhiyong

GM(1,1)灰色模型的程序实现function GM1=fungry1(x0) %输入原始数据x0 9 q! c# `) ?" BT=input('T=');%从键盘输入从最后一个历史数据算起的第T时点 x1=zeros(1,length(x0));B=zeros(length(x0)-1,2); $ s5 f* p) N6 T: m% X( Syn=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); ' i1 o: g. Q# p$ q4 _for i=1:length(x0)     for j=1:i         x1(i)=x1(i)+x0(j); ( A0 v& {: G4 K- _* _! ~    end ! x6 O7 h! J' jend + l8 e! e: v9 U% n9 j: A- Q2 s3 pfor i=1:length(x0)-1     B(i,1)=(-1/2)*(x1(i)+x1(i+1)); 4 x- v' X9 f8 H  `0 F6 M    B(i,2)=1;     yn(i)=x0(i+1); ' r& ]! i8 |8 ]) A7 ?end ( ]6 m" D- M: \$ N* X' @  K9 b* qHatA=(inv(B'*B))*B'*yn % GM(1,1)模型参数估计 ) R' T5 E- Q: ]+ r4 Bfor 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); 3 d/ g6 r7 g/ [1 V0 }/ A5 dfor k=2:length(x0)+T     Hatx0(k)=Hatx1(k)-Hatx1(k-1);%累计还原得到历史数据的模拟值 . }) Q) X8 M7 f' r6 t9 y8 O; @end 6 P" R% ~# d, B- g  v5 Y; hfor i=1:length(x0) %开始模型检验 ( y' H4 ]  b' O    epsilon(i)=x0(i)-Hatx0(i); % ?0 P) ]* {7 A, d. N9 g    omega(i)=(epsilon(i)/x0(i))*100; end % x0;Hatx0;epsilon;omega;  %必要时去掉％得到各种数据 c=std(epsilon)/std(x0);p=0; 9 O/ ?* ~/ k+ P. x* ffor i=1:length(x0) - \$ M! x2 G, h3 {    if abs(epsilon(i)-mean(epsilon))<0.6745*std(x0)         p=p+1;     end end 5 Y+ g4 d/ N- d4 L7 Sp=p/length(x0) * y4 }, h! O( xif 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 - D7 O9 T3 O* M+ t% o( N    disp('The model is eligibility,and the forecast is:'), ( c0 ]' w; i, h9 e9 K    disp(Hatx0(length(x0)+T)) * A. f4 O1 @+ b) T! c$ }' e- lelseif p>0.7 & c>0.65 ) Y& ^" H9 _" ?    disp('The model is not good,and the forecast is:'), * v7 J7 P  J" s' N6 t    disp(Hatx0(length(x0)+T))   H; G; I% G6 y8 m4 m4 @else p<=0.7 & c>0.65 ( V, U, D9 p! K0 e# e! w    disp('The model is bad and try again') end ' w3 _; P+ v  M, Afor i=1:length(x0) # N" y. j# a0 t: d& P- Z    Hatx00(i)=Hatx0(i); 1 B2 C, A- R* z% v/ c9 fend 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 ; V8 }# m0 O" k5 v# i$ w% |0.021800479 0.053132975 0.089908836 0.109153219 9 {1 Y# v4 b6 r# o7 r0.079331832 0 T% F& q1 m5 H3 K" S0.342192598 0.099718142 0.135194823   J! \6 \1 F# P! h0.109274037 7 u8 P. _  e  ~; U! \2 g! h+ r0.08152013 ' W" Q& @( I* z# s" k$ Y0.067876355 - \5 i3 y/ p6 e2 L& Z) _0.064706843 3 S/ c' E+ G. U' B: q' a0.055562197 , o6 p2 q% G; }% C' j0.050848544 + D4 Y( z6 x7 l# S1 n% z6 {]';

daiqiang5566

楼主莫急..

yysclshi

a= -0.0080
7 S8 ^+ {- a  j; iu= 0.0713
预测值
  ?5 `2 ]' q& T% F; M2 C4 b    1.6209    0.0846    0.0853    0.0859    0.0866    0.0873
7 N& n, v2 i: h5 g! h5 q    0.0880    0.0887    0.0895    0.0902    0.0909    0.0916
    0.0924    0.0931    0.0938    0.0946    0.0954    0.0961
初始值
3 Q7 {& m! g; t% c, o$ J4 v    1.6209    0.0793    0.0523    0.0413    0.0218    0.0531
    0.0899    0.1092    0.0793    0.3422    0.0997    0.1352
) T& e) p/ x! t% @    0.1093    0.0815    0.0679    0.0647    0.0556    0.0508
残差
9 V+ |9 v. Q4 }         0   -0.0053   -0.0329   -0.0447   -0.0648   -0.0342
    0.0019    0.0204   -0.0101    0.2520    0.0088    0.0436
3 ?' A2 q9 @9 s+ i. N    0.0169   -0.0116   -0.0260   -0.0299   -0.0398   -0.0453
! u  f* X0 `! Z相对误差
         0    0.0672    0.6297    1.0835    2.9741    0.6437
    0.0209    0.1870    0.1276    0.7365    0.0885    0.3223
    0.1548    0.1420    0.3826    0.4619    0.7162    0.8903
方差比
- G9 U7 Q$ D* g( R; Z    0.1869
p =
& J5 P! d% H5 ]: K  w; u+ o2 y     1

yysclshi

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

杨晓敬

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

gxj820

很好啊！！顶顶！

一剑卡卡

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

jshzncd

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