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

xinzhiyong

GM(1,1)灰色模型的程序实现function GM1=fungry1(x0) %输入原始数据x0 T=input('T=');%从键盘输入从最后一个历史数据算起的第T时点   H5 l, ]# h) m- y2 Ex1=zeros(1,length(x0));B=zeros(length(x0)-1,2); & T8 f8 ^) p. [' z0 Xyn=zeros(length(x0)-1,1);Hatx0=zeros(1,length(x0)+T); Hatx00=zeros(1,length(x0));Hatx1=zeros(1,length(x0)+T); & i+ i6 d+ S- R) K" Depsilon=zeros(length(x0),1);omega=zeros(length(x0),1); 0 v' e( ]! W; S: Q1 l% Xfor i=1:length(x0) 4 J( k( E4 Z1 Q8 N; P2 e    for j=1:i 1 X" V3 _/ J8 V# Q% O        x1(i)=x1(i)+x0(j);     end 0 e( y/ j9 @4 _* ]4 pend & \9 y9 _+ f+ z# Qfor i=1:length(x0)-1     B(i,1)=(-1/2)*(x1(i)+x1(i+1)); ' q' j- i% g( `  ]! `    B(i,2)=1;     yn(i)=x0(i+1); end $ k' G/ c: _! T9 {5 @" pHatA=(inv(B'*B))*B'*yn % GM(1,1)模型参数估计 for k=1:length(x0)+T 9 P- X2 b1 N2 U6 {" @    Hatx1(k)=(x0(1)-HatA(2)/HatA(1))*exp(-HatA(1)*(k-1))+HatA(2)/HatA(1); end 0 C% }7 p9 ~) aHatx0(1)=Hatx1(1); 4 {, y. a, o; ]- t7 G9 H$ yfor k=2:length(x0)+T     Hatx0(k)=Hatx1(k)-Hatx1(k-1);%累计还原得到历史数据的模拟值 end for i=1:length(x0) %开始模型检验 4 v! x, C/ Y2 M/ ^    epsilon(i)=x0(i)-Hatx0(i);     omega(i)=(epsilon(i)/x0(i))*100; # [1 @' i* C3 i3 _1 Yend % x0;Hatx0;epsilon;omega;  %必要时去掉％得到各种数据 1 }! N7 z0 a; W7 C6 o. |! ^" kc=std(epsilon)/std(x0);p=0; + A* M3 n5 S8 L  c: u+ \$ L. kfor i=1:length(x0) . V1 o& v3 F+ {7 x1 P" A" r    if abs(epsilon(i)-mean(epsilon))<0.6745*std(x0)         p=p+1; % R' k& C* g# c: `2 `, D6 |! O    end ! f: N5 |+ r1 L3 Z* _. z+ Kend 5 E  a5 W% m' M. n. Ap=p/length(x0) # `+ n) q5 e5 \+ |  @2 pif p>0.95 & c<0.35   a% z5 A) l$ {- E( w    disp('The model is good,and the forecast is:'),     disp(Hatx0(length(x0)+T)) elseif p>0.85 & c<0.5     disp('The model is eligibility,and the forecast is:'), % q5 @, @( ~/ u6 o, d' h- J    disp(Hatx0(length(x0)+T)) 8 c% p3 k3 S& p8 w1 F- F7 j3 N2 ?. Selseif p>0.7 & c>0.65     disp('The model is not good,and the forecast is:'),     disp(Hatx0(length(x0)+T)) else p<=0.7 & c>0.65     disp('The model is bad and try again') + T4 X' E7 j- }+ ^& U, aend * q2 z% r) O% c+ _for i=1:length(x0) 3 x! O, ^" {) r1 x" O    Hatx00(i)=Hatx0(i); end z=1:length(x0); 1 t( h9 O1 R1 [. q2 h0 kplot(z,x0,'-',z,Hatx00,':') %将原始数据和模拟值画在一个图上帮助观察 1 u; `# v' h& z  _9 ltext(2,x0(2),'History data: real line') $ w3 M+ y1 t* G( Ptext(length(x0)/2,Hatx00(length(x0))/2,'Simulation data:broken line') endT=input('T=');%从键盘输入从最后一个历史数据算起的第T时点????是指什么啊，请大哥们，大姐们教一下，我急用，请快，谢谢我的初始值x0=[1.620938526 0.07925621 0.052318818 0.041252502 0.021800479   `* r) ^$ d7 Q/ k! c% Y: c8 ^, Q0.053132975   ]# E0 V( s) N# C* e, W, h, F0.089908836 0.109153219 8 Y3 T1 e& q3 x- ~0.079331832 0.342192598 0.099718142 : P3 F# D( {" a7 m7 ^0.135194823 3 d. G  @" r, x& a( g  f. u: n0.109274037 9 |) _* H# {4 R& N0.08152013 6 D- ^0 _' M% u9 }0.067876355 0.064706843 ) {  v2 E4 h+ {0.055562197 & Q) y) c7 l( b- s+ ^0.050848544 . O" }" G" e) [0 C" i) ~* w, _! C]';

jshzncd

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

一剑卡卡

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

gxj820

很好啊！！顶顶！

杨晓敬

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

yysclshi

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

yysclshi

a= -0.0080
: Y. K9 T: j. R& gu= 0.0713
预测值
    1.6209    0.0846    0.0853    0.0859    0.0866    0.0873
9 E7 |. G0 Z; t: W8 W    0.0880    0.0887    0.0895    0.0902    0.0909    0.0916
    0.0924    0.0931    0.0938    0.0946    0.0954    0.0961
) M* |9 L* B' z/ ]初始值
    1.6209    0.0793    0.0523    0.0413    0.0218    0.0531
+ U" g6 f- W7 g1 q: ^+ U6 f5 }    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
9 |- K3 g1 N6 V5 u; q3 o7 S/ N    0.0019    0.0204   -0.0101    0.2520    0.0088    0.0436
. {1 y; g5 @* A. j    0.0169   -0.0116   -0.0260   -0.0299   -0.0398   -0.0453
相对误差
2 i* J, {4 L! ]2 w         0    0.0672    0.6297    1.0835    2.9741    0.6437
$ `* {& `/ `  {* D    0.0209    0.1870    0.1276    0.7365    0.0885    0.3223
2 }- C4 M+ e4 N# l: S  X5 E    0.1548    0.1420    0.3826    0.4619    0.7162    0.8903
方差比
    0.1869
) l, Q. K  ?. ?2 O* r9 ap =
3 t$ I+ k/ f" q- S5 ~9 u% Y$ b: P     1

daiqiang5566

楼主莫急..