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

xinzhiyong

GM(1,1)灰色模型的程序实现function GM1=fungry1(x0) %输入原始数据x0 T=input('T=');%从键盘输入从最后一个历史数据算起的第T时点 x1=zeros(1,length(x0));B=zeros(length(x0)-1,2); 1 t$ z( V5 W) T* hyn=zeros(length(x0)-1,1);Hatx0=zeros(1,length(x0)+T);   Q% z1 U; L4 H1 H' [% y$ b1 SHatx00=zeros(1,length(x0));Hatx1=zeros(1,length(x0)+T); 6 M. I; x- P1 nepsilon=zeros(length(x0),1);omega=zeros(length(x0),1); for i=1:length(x0)     for j=1:i ; {( P6 E) i- W- v, m% H        x1(i)=x1(i)+x0(j); * s/ X9 r3 X. }7 Q% n1 A9 j8 x    end end ; c9 }, Z/ g1 X2 w6 C8 F* j) m  [- Hfor i=1:length(x0)-1 0 U  U: J2 F5 z    B(i,1)=(-1/2)*(x1(i)+x1(i+1)); , D& g+ l, d* c4 V    B(i,2)=1;     yn(i)=x0(i+1); end ( V, f' a. c+ ?7 K5 T% a7 ]HatA=(inv(B'*B))*B'*yn % GM(1,1)模型参数估计 ( H- L2 C* Q: S% @2 z% t- `6 ]for k=1:length(x0)+T     Hatx1(k)=(x0(1)-HatA(2)/HatA(1))*exp(-HatA(1)*(k-1))+HatA(2)/HatA(1);   K$ j; |+ ]4 r( t  x% @8 Bend . A0 T/ W" h3 Q; @( ZHatx0(1)=Hatx1(1); - v3 S' G& |, l9 U' M& H5 `for k=2:length(x0)+T , f% ]; _' x- }9 A; o! p6 S    Hatx0(k)=Hatx1(k)-Hatx1(k-1);%累计还原得到历史数据的模拟值 + k3 v* ?- I  o% E" u- y3 b8 Jend for i=1:length(x0) %开始模型检验     epsilon(i)=x0(i)-Hatx0(i); 2 o9 a9 i$ m* l5 v+ |( _0 Y0 r4 j    omega(i)=(epsilon(i)/x0(i))*100; end % x0;Hatx0;epsilon;omega;  %必要时去掉％得到各种数据 c=std(epsilon)/std(x0);p=0; for i=1:length(x0) 8 f- q- p( a$ ]' W9 V    if abs(epsilon(i)-mean(epsilon))<0.6745*std(x0) 1 H( z7 p% F$ a2 |) b1 P        p=p+1; : o  B7 l- \, Y  k& D6 E2 b    end , J1 J3 M4 `* n' ^3 O# cend % y. x2 l/ @& Sp=p/length(x0) 1 I/ [( |8 R6 O1 }7 h* N3 nif p>0.95 & c<0.35     disp('The model is good,and the forecast is:'), - f0 j, \0 ^1 L    disp(Hatx0(length(x0)+T)) elseif p>0.85 & c<0.5     disp('The model is eligibility,and the forecast is:'),   m* w4 _  G' {& W4 Z* H    disp(Hatx0(length(x0)+T)) elseif 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') end , m) }6 H/ |' R7 M9 ^% e- w- _for i=1:length(x0) 4 n7 Y1 q4 \  Y4 ~2 \2 ]$ A    Hatx00(i)=Hatx0(i); end z=1:length(x0); * Q) Q. l. D% y% |0 K4 `plot(z,x0,'-',z,Hatx00,':') %将原始数据和模拟值画在一个图上帮助观察 text(2,x0(2),'History data: real line') + p1 f$ x, i/ K) r0 ^# dtext(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 0.053132975 0.089908836 7 p2 p8 f+ A7 O0.109153219 0.079331832 7 W1 l3 `: Y- N2 T4 x: [5 r0.342192598 , ]0 @6 K9 L: h+ l0.099718142 , k  F% c7 y9 L+ d1 d. F3 k9 P0.135194823 # f% i0 ~' X, Y9 n+ d3 n+ @0.109274037 0.08152013 0.067876355 0.064706843 " Q7 K& B/ ^8 n0.055562197 1 V3 ~9 g# b9 |; M0.050848544 ]';

daiqiang5566

楼主莫急..

yysclshi

a= -0.0080
u= 0.0713
/ [! r, t. Z- W* y7 E0 u( s预测值
    1.6209    0.0846    0.0853    0.0859    0.0866    0.0873
8 N) U$ z7 s8 q, h/ _! i* G    0.0880    0.0887    0.0895    0.0902    0.0909    0.0916
    0.0924    0.0931    0.0938    0.0946    0.0954    0.0961
初始值
    1.6209    0.0793    0.0523    0.0413    0.0218    0.0531
    0.0899    0.1092    0.0793    0.3422    0.0997    0.1352
    0.1093    0.0815    0.0679    0.0647    0.0556    0.0508
+ D/ I0 W' Q4 `+ B残差
" U  M' J2 B6 k2 R7 O6 v         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
相对误差
( I* o) j' _  ?+ b/ @         0    0.0672    0.6297    1.0835    2.9741    0.6437
  t$ J4 @  W# |4 n    0.0209    0.1870    0.1276    0.7365    0.0885    0.3223
    0.1548    0.1420    0.3826    0.4619    0.7162    0.8903
方差比
    0.1869
p =
5 a- J4 e. x( E" e7 ]5 r$ _     1

yysclshi

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

杨晓敬

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

gxj820

很好啊！！顶顶！

一剑卡卡

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

jshzncd

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