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

xinzhiyong

GM(1,1)灰色模型的程序实现function GM1=fungry1(x0) %输入原始数据x0 8 g# q' Z+ P7 R- x- }T=input('T=');%从键盘输入从最后一个历史数据算起的第T时点 x1=zeros(1,length(x0));B=zeros(length(x0)-1,2); yn=zeros(length(x0)-1,1);Hatx0=zeros(1,length(x0)+T); ; I' {! y' J( j# I$ V) nHatx00=zeros(1,length(x0));Hatx1=zeros(1,length(x0)+T); epsilon=zeros(length(x0),1);omega=zeros(length(x0),1); $ X/ W. H: i7 Y, d, @4 ufor i=1:length(x0)     for j=1:i ) w/ M. @! i) r" g' H( L. d        x1(i)=x1(i)+x0(j);     end 0 z7 ~- R( }# [" w4 |5 J0 ^end 2 A% H9 E% V& |% j8 Dfor i=1:length(x0)-1 # t1 }& e, {* h4 n7 ]8 ?    B(i,1)=(-1/2)*(x1(i)+x1(i+1));     B(i,2)=1;     yn(i)=x0(i+1); 5 @- g$ f; m& x9 {1 \: @' Rend HatA=(inv(B'*B))*B'*yn % GM(1,1)模型参数估计   h, |' `) h+ _% B5 i/ D$ mfor 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); 7 Y) @# {- ~+ _$ ]for k=2:length(x0)+T     Hatx0(k)=Hatx1(k)-Hatx1(k-1);%累计还原得到历史数据的模拟值 & o" G9 f# M# [! W3 cend ( y# I6 |6 z: N: O: i% Tfor i=1:length(x0) %开始模型检验 ' r: Y6 ?2 Z: `( m/ x& Q    epsilon(i)=x0(i)-Hatx0(i);     omega(i)=(epsilon(i)/x0(i))*100; end 3 b* |' q6 H* E* l6 Y+ b( J% x0;Hatx0;epsilon;omega;  %必要时去掉％得到各种数据 c=std(epsilon)/std(x0);p=0; 3 i$ S/ l4 f$ ~# z8 J) k2 a5 N0 [for i=1:length(x0)     if abs(epsilon(i)-mean(epsilon))<0.6745*std(x0) + e* M3 f* Q8 U, p7 Q        p=p+1; ! _5 G  ^. Z! j    end end ! q6 k/ t4 Y+ w+ ]% @; M: j! x0 bp=p/length(x0) if p>0.95 & c<0.35 9 T. j6 P! E* M3 g$ n    disp('The model is good,and the forecast is:'), * o9 y, A6 T! I; |9 F% B: l! @    disp(Hatx0(length(x0)+T)) ( d3 T" _- I' j" k9 K5 {elseif p>0.85 & c<0.5 ' n6 d  {  r/ \- }    disp('The model is eligibility,and the forecast is:'),     disp(Hatx0(length(x0)+T)) / f8 ~6 c& j1 e" ?* c$ J/ belseif p>0.7 & c>0.65     disp('The model is not good,and the forecast is:'), 8 @* o; H. _7 l9 ?  _- J: G    disp(Hatx0(length(x0)+T)) ; h: T+ Z& A  r9 Q) Pelse p<=0.7 & c>0.65   K3 Y" }4 M  ]! D" W    disp('The model is bad and try again') end ! Q. x8 t4 v, E. y/ `4 c# h) M. Mfor 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') endT=input('T=');%从键盘输入从最后一个历史数据算起的第T时点????是指什么啊，请大哥们，大姐们教一下，我急用，请快，谢谢我的初始值x0=[1.620938526 4 w) |$ r. J( M0.07925621 0.052318818 . H2 C. K6 H$ _0.041252502 0.021800479 7 o) l6 n1 b- J5 o3 h/ [- X0.053132975 0.089908836 0.109153219 0.079331832 0.342192598 ! U% i1 y, T5 ?7 w5 j6 g. V) w0.099718142 0 y5 Q+ V1 b9 P4 A, F! S0.135194823 * J7 ]6 V- I$ N3 {8 m* ~0.109274037 - K* q! F. v3 h- o4 N* `0.08152013 0.067876355 0.064706843 % b: w* N" q1 X+ W0.055562197 3 F3 N) }# g2 |5 |, x0.050848544 ]';

daiqiang5566

楼主莫急..

yysclshi

a= -0.0080
u= 0.0713
预测值
0 X) @: T$ u% `% c5 a" |4 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
4 s1 [, A3 @  Q) k2 \初始值
    1.6209    0.0793    0.0523    0.0413    0.0218    0.0531
    0.0899    0.1092    0.0793    0.3422    0.0997    0.1352
! K( t' d+ d; ]9 f; I& {8 ^    0.1093    0.0815    0.0679    0.0647    0.0556    0.0508
残差
         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
! _$ w" t9 l2 I0 W相对误差
         0    0.0672    0.6297    1.0835    2.9741    0.6437
    0.0209    0.1870    0.1276    0.7365    0.0885    0.3223
! G1 D" b% j, C; Y( z% u8 V    0.1548    0.1420    0.3826    0.4619    0.7162    0.8903
% N' S4 H, b- K/ S- y方差比
    0.1869
p =
0 o+ A1 {9 Y* `     1

yysclshi

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

杨晓敬

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

gxj820

很好啊！！顶顶！

一剑卡卡

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

jshzncd

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