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

xinzhiyong

GM(1,1)灰色模型的程序实现function GM1=fungry1(x0) %输入原始数据x0 T=input('T=');%从键盘输入从最后一个历史数据算起的第T时点 ) Z) B* M+ n7 V* b; A# c3 G! d% k! U8 Cx1=zeros(1,length(x0));B=zeros(length(x0)-1,2); 3 U2 M8 s! Q6 Wyn=zeros(length(x0)-1,1);Hatx0=zeros(1,length(x0)+T); Hatx00=zeros(1,length(x0));Hatx1=zeros(1,length(x0)+T); : d# [) X3 [" Q+ Wepsilon=zeros(length(x0),1);omega=zeros(length(x0),1); for i=1:length(x0)     for j=1:i         x1(i)=x1(i)+x0(j); % G8 G$ C# |( q0 X* @& n    end end for i=1:length(x0)-1     B(i,1)=(-1/2)*(x1(i)+x1(i+1));     B(i,2)=1;     yn(i)=x0(i+1); : A# m+ `  h( I7 hend HatA=(inv(B'*B))*B'*yn % GM(1,1)模型参数估计 6 n5 Z' i  B5 P( Sfor k=1:length(x0)+T , `; N" S8 M4 T: S    Hatx1(k)=(x0(1)-HatA(2)/HatA(1))*exp(-HatA(1)*(k-1))+HatA(2)/HatA(1); end 2 y! j5 t4 s5 K7 \7 C1 nHatx0(1)=Hatx1(1); for k=2:length(x0)+T 5 j3 h) B" ]* _+ a$ X2 m    Hatx0(k)=Hatx1(k)-Hatx1(k-1);%累计还原得到历史数据的模拟值 , j9 i" a4 f; R' E/ ^end 5 N6 Y+ q* w+ ?$ a' M# {0 s2 hfor i=1:length(x0) %开始模型检验     epsilon(i)=x0(i)-Hatx0(i);     omega(i)=(epsilon(i)/x0(i))*100;   ~8 w) M7 n' }, Hend % x0;Hatx0;epsilon;omega;  %必要时去掉％得到各种数据 2 I. P; H+ t/ e7 Kc=std(epsilon)/std(x0);p=0; $ ~% Q  C1 d0 r' ?/ Mfor i=1:length(x0) # t' v" R+ {1 E/ v! M  B    if abs(epsilon(i)-mean(epsilon))<0.6745*std(x0)         p=p+1;     end , k8 J' X! [/ Y7 D2 I. l" p& Q( lend p=p/length(x0) # |  F+ H$ t  l# J; y  a7 `$ P  ~if p>0.95 & c<0.35 ) S* {8 q! M# I2 z. u% n2 {    disp('The model is good,and the forecast is:'), 1 V; l# y( [1 F) t. N% F+ y* W    disp(Hatx0(length(x0)+T)) , D6 A; h4 \' E4 N7 felseif p>0.85 & c<0.5     disp('The model is eligibility,and the forecast is:'),     disp(Hatx0(length(x0)+T)) elseif p>0.7 & c>0.65     disp('The model is not good,and the forecast is:'), $ e7 K( m2 t/ T4 c. {* F    disp(Hatx0(length(x0)+T)) else p<=0.7 & c>0.65     disp('The model is bad and try again') end & f. N) h# ~: ?2 q8 Zfor i=1:length(x0)     Hatx00(i)=Hatx0(i); end 9 }. S0 [2 Y8 u& J5 T4 Fz=1:length(x0); plot(z,x0,'-',z,Hatx00,':') %将原始数据和模拟值画在一个图上帮助观察 text(2,x0(2),'History data: real line') : Z2 f9 K7 C7 a) R4 l! Atext(length(x0)/2,Hatx00(length(x0))/2,'Simulation data:broken line') 5 [+ A- Q/ z  U. k4 QendT=input('T=');%从键盘输入从最后一个历史数据算起的第T时点????是指什么啊，请大哥们，大姐们教一下，我急用，请快，谢谢我的初始值x0=[1.620938526 $ i3 e# J5 ^. V/ ?$ B4 e0.07925621 0.052318818 0.041252502 6 g, S1 p* _. Q1 H# \) [0.021800479 + _# ?7 X7 c. o) ~: ^& R) u0.053132975 , w# u7 y. W8 m7 Q0.089908836 6 v; x$ h# \: b0.109153219   F6 z0 B1 g% s0.079331832 0.342192598 0.099718142 - O/ U2 \6 I& \: N7 N7 z0 ~0.135194823 & d2 w# o1 N5 b# Q1 A0 g' a5 ^7 P0.109274037 1 a) C% _8 S  K6 @. A0.08152013 0.067876355 0.064706843 0.055562197 0.050848544 0 i6 a. x+ K: D2 K' ^! D" g]';

jshzncd

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

一剑卡卡

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

gxj820

很好啊！！顶顶！

杨晓敬

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

yysclshi

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

yysclshi

a= -0.0080
. z/ n# t& y! _% Nu= 0.0713
预测值
) x4 ?2 T9 X2 L: L3 d( Y    1.6209    0.0846    0.0853    0.0859    0.0866    0.0873
+ P7 b6 q0 Z) E5 N4 Q    0.0880    0.0887    0.0895    0.0902    0.0909    0.0916
" r7 H' b2 r0 G  [    0.0924    0.0931    0.0938    0.0946    0.0954    0.0961
9 @: z9 O5 J+ b2 Z初始值
: c! ]* U8 }9 A( n; f/ D    1.6209    0.0793    0.0523    0.0413    0.0218    0.0531
    0.0899    0.1092    0.0793    0.3422    0.0997    0.1352
6 ?/ q- ^/ l: ^4 r, D& P, D    0.1093    0.0815    0.0679    0.0647    0.0556    0.0508
; _9 q& U5 N0 f5 n残差
         0   -0.0053   -0.0329   -0.0447   -0.0648   -0.0342
    0.0019    0.0204   -0.0101    0.2520    0.0088    0.0436
* R0 o- d' {% G* B# M; s; [    0.0169   -0.0116   -0.0260   -0.0299   -0.0398   -0.0453
相对误差
1 s4 P+ s# n: i8 n9 K0 F# _         0    0.0672    0.6297    1.0835    2.9741    0.6437
; u) M8 ^  q: Y. l: U    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
* |+ ~8 v- Z: f& D& {p =
2 J0 i% P/ [6 e" ~: h" Z7 Z     1

daiqiang5566

楼主莫急..