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

xinzhiyong

GM(1,1)灰色模型的程序实现function GM1=fungry1(x0) %输入原始数据x0 T=input('T=');%从键盘输入从最后一个历史数据算起的第T时点 9 w+ }4 {+ K6 y6 F3 H; m( Rx1=zeros(1,length(x0));B=zeros(length(x0)-1,2); : K9 n; L5 F. q5 xyn=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)); ' w7 W0 [. X1 G5 Z1 @  a    B(i,2)=1; + Q6 D% H8 I4 H% H    yn(i)=x0(i+1); ; [; @0 a5 P9 eend HatA=(inv(B'*B))*B'*yn % GM(1,1)模型参数估计 3 E: j& T% D/ t- c2 tfor k=1:length(x0)+T     Hatx1(k)=(x0(1)-HatA(2)/HatA(1))*exp(-HatA(1)*(k-1))+HatA(2)/HatA(1); 1 R9 @' z) r9 A/ V5 Qend 8 N& D% Q" l) j) w: V4 f) o* mHatx0(1)=Hatx1(1); 6 I" p% l6 y; Z& a- L+ ofor k=2:length(x0)+T - O, w: j+ K* O, S' P7 [0 }( g    Hatx0(k)=Hatx1(k)-Hatx1(k-1);%累计还原得到历史数据的模拟值 1 S9 J* V- S1 l- G! lend 8 z: i2 ~/ i, H: q. ?8 W' g% Ufor i=1:length(x0) %开始模型检验 9 {  e, |# ]$ J( a9 p6 \8 h    epsilon(i)=x0(i)-Hatx0(i); ; {" f# x4 o5 P: U/ [0 ~    omega(i)=(epsilon(i)/x0(i))*100; end + w3 X- x8 z1 P  m, G& L! T% x0;Hatx0;epsilon;omega;  %必要时去掉％得到各种数据 : U9 ?+ |3 g: Q( w8 G- k$ d/ hc=std(epsilon)/std(x0);p=0; for i=1:length(x0) " ^3 F2 f. Y) b  H4 B# u0 {% E    if abs(epsilon(i)-mean(epsilon))<0.6745*std(x0)         p=p+1; , A% K4 _6 \8 ^7 I    end + T( G, a2 ]1 ]) Q0 zend & E9 Y! d5 Z$ M, d' v# 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))   [- b! h: @8 ?$ b3 r7 L3 uelseif p>0.85 & c<0.5     disp('The model is eligibility,and the forecast is:'), : \. B+ i6 j1 Z0 R    disp(Hatx0(length(x0)+T)) 0 q) K+ {2 s! P: s$ I  belseif p>0.7 & c>0.65 ! z& Y1 K6 s& H# ?# D+ l; J" [    disp('The model is not good,and the forecast is:'), 3 `  l. a/ v' o; s+ Y2 d    disp(Hatx0(length(x0)+T)) 8 m2 E& H( H( d5 t( Velse p<=0.7 & c>0.65     disp('The model is bad and try again') 5 T2 ^. g  ^0 c9 b1 C2 ]. fend for i=1:length(x0)   }- H; z6 J+ P! [* \" [    Hatx00(i)=Hatx0(i); end % G- L, _& ~+ J5 q5 z' sz=1:length(x0); plot(z,x0,'-',z,Hatx00,':') %将原始数据和模拟值画在一个图上帮助观察   p8 s/ _3 I. N" ?, U' O) j9 qtext(2,x0(2),'History data: real line') text(length(x0)/2,Hatx00(length(x0))/2,'Simulation data:broken line') 6 H/ |+ x; \& \' _endT=input('T=');%从键盘输入从最后一个历史数据算起的第T时点????是指什么啊，请大哥们，大姐们教一下，我急用，请快，谢谢我的初始值x0=[1.620938526 ( h4 u' `1 C" z: Q0.07925621 ' i. P" _* o8 j5 A  ^  u, W  ]' k0.052318818 0.041252502 0.021800479 2 l. c; q# ^) n4 O& Y0.053132975 0.089908836 ( y) O+ U0 D+ M7 g+ @& d0.109153219 1 C6 V1 H5 O2 U1 D7 X0.079331832 5 c# c  {8 t" }. `8 U0.342192598 0.099718142 + }; ^  V  g  L; n; |0.135194823 0.109274037 0 ~0 u. C9 t$ Q3 Q1 W0.08152013 ' q( X- e8 i) e. k  i5 v( g& h0.067876355 0.064706843 0.055562197 ' \  }' n- M# c; q, I1 q9 \0.050848544 ]';

daiqiang5566

楼主莫急..

yysclshi

a= -0.0080
u= 0.0713
( B9 F' D) }$ Y/ G预测值
7 D) }  L6 }# a+ g6 e& B( s    1.6209    0.0846    0.0853    0.0859    0.0866    0.0873
9 O3 M- J5 d3 N. R0 B# S    0.0880    0.0887    0.0895    0.0902    0.0909    0.0916
7 y) e9 r  K4 d6 L+ g" h: y    0.0924    0.0931    0.0938    0.0946    0.0954    0.0961
! E1 W8 ]( d/ [8 M4 D: y初始值
    1.6209    0.0793    0.0523    0.0413    0.0218    0.0531
. O7 D1 _9 K; e    0.0899    0.1092    0.0793    0.3422    0.0997    0.1352
    0.1093    0.0815    0.0679    0.0647    0.0556    0.0508
残差
3 j. e! B) m% t5 Q$ C" u         0   -0.0053   -0.0329   -0.0447   -0.0648   -0.0342
    0.0019    0.0204   -0.0101    0.2520    0.0088    0.0436
    0.0169   -0.0116   -0.0260   -0.0299   -0.0398   -0.0453
$ a4 r; ?; B/ }& J相对误差
         0    0.0672    0.6297    1.0835    2.9741    0.6437
8 q2 W* P* T  I    0.0209    0.1870    0.1276    0.7365    0.0885    0.3223
: D# T8 W# H- z, ]' c9 ~    0.1548    0.1420    0.3826    0.4619    0.7162    0.8903
方差比
    0.1869
' p2 |& ]/ n- S' `$ r( K9 jp =
     1

yysclshi

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

杨晓敬

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

gxj820

很好啊！！顶顶！

一剑卡卡

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

jshzncd

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