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

xinzhiyong

GM(1,1)灰色模型的程序实现function GM1=fungry1(x0) %输入原始数据x0 T=input('T=');%从键盘输入从最后一个历史数据算起的第T时点 6 U0 q, A3 J, nx1=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); ; v1 E% E3 r4 E& B! \for i=1:length(x0)     for j=1:i         x1(i)=x1(i)+x0(j);     end 8 H. u1 S% D# @. r3 @2 P) Eend for i=1:length(x0)-1     B(i,1)=(-1/2)*(x1(i)+x1(i+1));     B(i,2)=1; & }" ^, E; i7 a9 Q/ R    yn(i)=x0(i+1); " |1 M7 }& t' ^) E5 rend 0 o4 p( \0 y" A1 RHatA=(inv(B'*B))*B'*yn % GM(1,1)模型参数估计 * Q! \2 z! k% l( Ffor k=1:length(x0)+T 1 q2 Y1 H" o  U1 b$ r- [; z5 A    Hatx1(k)=(x0(1)-HatA(2)/HatA(1))*exp(-HatA(1)*(k-1))+HatA(2)/HatA(1); ) p9 I' x7 h( T: zend Hatx0(1)=Hatx1(1); for k=2:length(x0)+T . p5 g5 ^! e; T+ u" A    Hatx0(k)=Hatx1(k)-Hatx1(k-1);%累计还原得到历史数据的模拟值 - Y% A4 j# ^% l' f  b( l# Z; yend , J1 A1 _( i- m* l  Lfor i=1:length(x0) %开始模型检验 ) v* U" M# K/ d( \% E0 G    epsilon(i)=x0(i)-Hatx0(i);     omega(i)=(epsilon(i)/x0(i))*100; end % x0;Hatx0;epsilon;omega;  %必要时去掉％得到各种数据 : W* u& J& q0 i( B+ kc=std(epsilon)/std(x0);p=0; 3 K. Y5 g  S' k' j/ X+ gfor i=1:length(x0)     if abs(epsilon(i)-mean(epsilon))<0.6745*std(x0)         p=p+1; 3 }' E8 [" Z5 {6 ?: q. S" J    end end p=p/length(x0) % z. e2 X5 b$ ~+ mif p>0.95 & c<0.35     disp('The model is good,and the forecast is:'), : H0 s0 J5 C+ N/ x/ ]  D    disp(Hatx0(length(x0)+T)) 5 Y: L) g' l/ S. x3 Jelseif p>0.85 & c<0.5 5 O2 N! f0 P: B+ ^    disp('The model is eligibility,and the forecast is:'),     disp(Hatx0(length(x0)+T)) elseif p>0.7 & c>0.65 # c/ V1 a! T. P( k2 o    disp('The model is not good,and the forecast is:'),     disp(Hatx0(length(x0)+T)) $ D9 p2 P8 }- j! |& ?9 J  xelse p<=0.7 & c>0.65 7 v5 S) D) \! D: `    disp('The model is bad and try again') , R* ]4 Q; V5 ?- k0 P  K9 Pend for i=1:length(x0)     Hatx00(i)=Hatx0(i); end z=1:length(x0); 0 k4 w/ @6 B+ G  X5 N2 ~# l% v5 e  z2 vplot(z,x0,'-',z,Hatx00,':') %将原始数据和模拟值画在一个图上帮助观察 text(2,x0(2),'History data: real line') . R! K0 V  E1 Ntext(length(x0)/2,Hatx00(length(x0))/2,'Simulation data:broken line') endT=input('T=');%从键盘输入从最后一个历史数据算起的第T时点????是指什么啊，请大哥们，大姐们教一下，我急用，请快，谢谢我的初始值x0=[1.620938526 ) m! [" C% C2 o  j7 h" D! F' O0.07925621 ( y  i. O  J5 U3 w* o0 a. ?/ N8 Z0.052318818 ! S% b) R7 n8 B5 Z2 s$ W& j0.041252502 # C1 {. Z: P3 I. ~- P3 I7 J/ U: C0.021800479 0.053132975 $ c/ I! |/ U$ I& x" [1 W6 \0.089908836 0.109153219 ; r- z: p. s" e0.079331832 - T0 l* w8 }, B! _- c, h0 Q' p0.342192598 0.099718142 0.135194823 9 V- n: w2 U# X/ V) S  I0 s1 z0.109274037 0.08152013 0.067876355 0.064706843 0.055562197 0.050848544 ]';

daiqiang5566

楼主莫急..

yysclshi

a= -0.0080
% C( q: O; W- wu= 0.0713
预测值
    1.6209    0.0846    0.0853    0.0859    0.0866    0.0873
, e* n- U" r4 T  J7 F) J3 P    0.0880    0.0887    0.0895    0.0902    0.0909    0.0916
    0.0924    0.0931    0.0938    0.0946    0.0954    0.0961
* S/ K' H- O" ^- ?! E4 R初始值
    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
残差
         0   -0.0053   -0.0329   -0.0447   -0.0648   -0.0342
' A$ D5 s5 }# S3 ^) z9 t' B- M' ?. m    0.0019    0.0204   -0.0101    0.2520    0.0088    0.0436
! s: m0 s; K0 I6 W8 X8 k& N    0.0169   -0.0116   -0.0260   -0.0299   -0.0398   -0.0453
相对误差
; d5 N) }4 u4 t9 W         0    0.0672    0.6297    1.0835    2.9741    0.6437
  k' j- c5 m# D( q8 h    0.0209    0.1870    0.1276    0.7365    0.0885    0.3223
+ X0 t. M8 Y$ W3 _% g, @    0.1548    0.1420    0.3826    0.4619    0.7162    0.8903
- F: A& ]. R+ N( B- _方差比
$ C7 k7 C/ R8 s4 e% T% a    0.1869
+ G; i$ [. [& L/ r6 [p =
     1

yysclshi

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

杨晓敬

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

gxj820

很好啊！！顶顶！

一剑卡卡

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

jshzncd

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