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); $ m( }% I/ }2 E. fyn=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); 8 G3 E/ X3 F0 z6 r* @for i=1:length(x0)     for j=1:i 2 L, C* X+ r8 _) \+ y        x1(i)=x1(i)+x0(j); 9 L% f, h: E) C! i& j& j" M    end end ; T4 w- ^( o6 O+ Y6 Efor i=1:length(x0)-1 ( u9 s* c' H0 c# V9 B3 |: E- N- z    B(i,1)=(-1/2)*(x1(i)+x1(i+1)); : g! M. t2 h# N0 A/ j, T5 E    B(i,2)=1; 2 M6 |( x5 C9 n( s& B    yn(i)=x0(i+1); - _7 y0 x$ ]3 B1 k  Zend HatA=(inv(B'*B))*B'*yn % GM(1,1)模型参数估计 ; a8 }. g+ V/ a7 J4 K! l, a: Nfor k=1:length(x0)+T % r; a% D6 x- q    Hatx1(k)=(x0(1)-HatA(2)/HatA(1))*exp(-HatA(1)*(k-1))+HatA(2)/HatA(1); end Hatx0(1)=Hatx1(1); for k=2:length(x0)+T     Hatx0(k)=Hatx1(k)-Hatx1(k-1);%累计还原得到历史数据的模拟值 4 u9 e  {2 \7 \3 A& Bend for i=1:length(x0) %开始模型检验     epsilon(i)=x0(i)-Hatx0(i);     omega(i)=(epsilon(i)/x0(i))*100; end ( m) C8 ?& N& f% x0;Hatx0;epsilon;omega;  %必要时去掉％得到各种数据 c=std(epsilon)/std(x0);p=0; for i=1:length(x0) 1 O3 Q$ v1 D' u9 D5 P8 x1 H9 H& ~    if abs(epsilon(i)-mean(epsilon))<0.6745*std(x0)         p=p+1;     end 8 ?4 i" |: f8 vend ; |/ O* _# m* Z! wp=p/length(x0) if p>0.95 & c<0.35   F2 ]/ U% L; Z9 [9 i( J: j  [    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:'), % _+ f3 X  U8 q; K9 }    disp(Hatx0(length(x0)+T)) elseif p>0.7 & c>0.65 * C1 n& D9 B  |- Z    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') , o- q' G$ d! G# aend for i=1:length(x0)     Hatx00(i)=Hatx0(i); end % O; `. N' o( S4 G) }# lz=1:length(x0); plot(z,x0,'-',z,Hatx00,':') %将原始数据和模拟值画在一个图上帮助观察 text(2,x0(2),'History data: real line') / W/ t! E6 }$ P% o9 l0 jtext(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 & K+ n+ K  E9 I5 ^0.053132975 0.089908836   s, b1 L: _& k( w+ P+ m" P7 J0.109153219 5 p; h; |, s4 C% @: q6 I0.079331832 . p# @2 n- s& s# c0.342192598 0.099718142 0.135194823 ! t$ `! J) i4 g6 A9 c0.109274037 / q# J3 l0 ^) _/ k/ w' h0.08152013 ' Z% X9 v; E; r' G# ]7 C; N0 A0.067876355 8 _  C; a: f+ Z7 V& a# J0.064706843 8 I7 }% e! E& r) p0.055562197 0.050848544 ]';

daiqiang5566

楼主莫急..

yysclshi

a= -0.0080
u= 0.0713
5 U8 ^) b2 |7 ~4 [& O7 K预测值
  ^5 h) v5 N  X) }  b    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
9 e" t2 \$ a) v初始值
& R0 v% F/ _- |1 a3 D    1.6209    0.0793    0.0523    0.0413    0.0218    0.0531
7 p5 g1 j  Y; ?    0.0899    0.1092    0.0793    0.3422    0.0997    0.1352
    0.1093    0.0815    0.0679    0.0647    0.0556    0.0508
6 h& E& s$ [2 o! x  h残差
         0   -0.0053   -0.0329   -0.0447   -0.0648   -0.0342
# Z7 R/ |; c9 @% e    0.0019    0.0204   -0.0101    0.2520    0.0088    0.0436
# `1 `$ J& h. D  u4 x. m2 k    0.0169   -0.0116   -0.0260   -0.0299   -0.0398   -0.0453
相对误差
         0    0.0672    0.6297    1.0835    2.9741    0.6437
+ H* A  V( t: m: d4 d$ ?! M& G    0.0209    0.1870    0.1276    0.7365    0.0885    0.3223
) A" t1 L* G0 x  L1 H    0.1548    0.1420    0.3826    0.4619    0.7162    0.8903
$ D7 X7 u9 r) Q+ k, P方差比
8 W7 ?& n& d& {) d/ f    0.1869
p =
     1

yysclshi

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

杨晓敬

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

gxj820

很好啊！！顶顶！

一剑卡卡

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

jshzncd

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