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

xinzhiyong

GM(1,1)灰色模型的程序实现function GM1=fungry1(x0) %输入原始数据x0   c" ]- x' r2 n0 E# l  I# ?T=input('T=');%从键盘输入从最后一个历史数据算起的第T时点 0 V% N' h1 H2 Lx1=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 6 d# K7 H1 ^7 T5 H2 T5 @* N        x1(i)=x1(i)+x0(j); + ^' D- |" k- O* T% b# E6 ^    end $ ~& ~9 \8 G- U8 e  E, D4 rend ; K4 ?- ]  f8 x4 M. X8 d  gfor i=1:length(x0)-1     B(i,1)=(-1/2)*(x1(i)+x1(i+1)); / w7 I! A- r* l" e; R    B(i,2)=1;     yn(i)=x0(i+1); . n7 m+ I$ n# _% ?% T8 [, Aend " F3 B5 v+ c( M4 eHatA=(inv(B'*B))*B'*yn % GM(1,1)模型参数估计 9 X3 @" q# v2 U( s. }8 Wfor k=1:length(x0)+T     Hatx1(k)=(x0(1)-HatA(2)/HatA(1))*exp(-HatA(1)*(k-1))+HatA(2)/HatA(1); ' \8 i4 z/ B3 ?8 a" t2 Eend Hatx0(1)=Hatx1(1); for k=2:length(x0)+T     Hatx0(k)=Hatx1(k)-Hatx1(k-1);%累计还原得到历史数据的模拟值 end 5 r0 \+ O' D: @/ Z4 s+ \. }* Zfor i=1:length(x0) %开始模型检验 7 R. c- \  ]4 t0 I5 {    epsilon(i)=x0(i)-Hatx0(i); / o+ r+ Z/ V4 ^8 r, ]9 G9 y    omega(i)=(epsilon(i)/x0(i))*100; end . k+ l: ?5 N# D  A3 G1 x% x0;Hatx0;epsilon;omega;  %必要时去掉％得到各种数据 c=std(epsilon)/std(x0);p=0; 0 j5 Y& B6 m7 m: J8 S0 ?for i=1:length(x0)     if abs(epsilon(i)-mean(epsilon))<0.6745*std(x0) 8 X* S3 N& \- G6 z2 Y4 m& E        p=p+1;     end end ; x  v' p' p# E9 w3 W) @p=p/length(x0) & d! ]4 A/ m# }+ i$ g4 ?if p>0.95 & c<0.35     disp('The model is good,and the forecast is:'),     disp(Hatx0(length(x0)+T)) / [! f. }9 i. t6 U/ K1 qelseif p>0.85 & c<0.5 3 m; G0 @# j+ W- K! d    disp('The model is eligibility,and the forecast is:'), # p2 p+ l; _# H0 J    disp(Hatx0(length(x0)+T))   m8 i- P& S  u1 q1 O+ Zelseif p>0.7 & c>0.65 ( J* o7 k8 H1 a$ k    disp('The model is not good,and the forecast is:'), ; F3 E& N- D6 f* A# i8 N: w8 {7 r    disp(Hatx0(length(x0)+T)) 3 w) W& S, _2 X) L  Uelse p<=0.7 & c>0.65 * R: q( M  Z0 E9 O    disp('The model is bad and try again') 0 g% c, {  g4 g) E2 jend / ~: a$ ^2 g2 e6 l: m3 L! w+ Kfor 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') . I* c: a5 z& q2 b0 S, A5 |/ GendT=input('T=');%从键盘输入从最后一个历史数据算起的第T时点????是指什么啊，请大哥们，大姐们教一下，我急用，请快，谢谢我的初始值x0=[1.620938526 0.07925621 0 O% X! e/ w; \% A, U0.052318818 . }/ S  x+ ~1 }8 }0 T0 s* C+ y0.041252502 8 U- t" u: ?8 _. t4 M0.021800479 0.053132975 + Z2 U) t: s. N6 \  C0.089908836 * L3 ^8 B6 x7 n, ^& }& R8 k1 @8 ]2 t* g0.109153219 6 F: k- Q0 j: T! p0.079331832 ' Z# X* R7 v5 Z0.342192598 0.099718142 0.135194823 % k4 Y- v  c+ |0 n4 Y/ r/ A& ], v0.109274037 1 K& M. S. u* K3 t2 m% `0.08152013 0.067876355 0.064706843 ; t, ?& P1 P% X; ?2 K0.055562197 + ?! G: u- p  V! b9 A0.050848544 + z* P2 G$ t( h8 y]';

jshzncd

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

一剑卡卡

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

gxj820

很好啊！！顶顶！

杨晓敬

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

yysclshi

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

yysclshi

a= -0.0080
5 @/ F6 r4 |3 B' {5 ^$ s0 qu= 0.0713
, E: c9 R+ @7 \" t预测值
$ u# \1 Q" W( `# w# f) l    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
3 t+ c1 @8 k+ U# T$ @' a! N初始值
    1.6209    0.0793    0.0523    0.0413    0.0218    0.0531
1 A, }5 t( w% h1 q0 R    0.0899    0.1092    0.0793    0.3422    0.0997    0.1352
    0.1093    0.0815    0.0679    0.0647    0.0556    0.0508
残差
         0   -0.0053   -0.0329   -0.0447   -0.0648   -0.0342
2 X- m  o0 G: W' d' @$ u  n7 h    0.0019    0.0204   -0.0101    0.2520    0.0088    0.0436
    0.0169   -0.0116   -0.0260   -0.0299   -0.0398   -0.0453
相对误差
         0    0.0672    0.6297    1.0835    2.9741    0.6437
    0.0209    0.1870    0.1276    0.7365    0.0885    0.3223
0 p1 E4 [# |) M$ T5 Y    0.1548    0.1420    0.3826    0.4619    0.7162    0.8903
& J# m4 v+ E! V2 ?9 z' m* R% V方差比
    0.1869
! u+ |& G9 O- b  S$ Ep =
     1

daiqiang5566

楼主莫急..