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); ; d7 y/ T/ {: p" R- Q4 G0 ^9 R! v0 nyn=zeros(length(x0)-1,1);Hatx0=zeros(1,length(x0)+T); ' `$ L) Q- D2 E! B* N+ }5 bHatx00=zeros(1,length(x0));Hatx1=zeros(1,length(x0)+T);   d# F, p$ u0 P8 L& X. M; O+ kepsilon=zeros(length(x0),1);omega=zeros(length(x0),1); . t7 o8 J# K- h4 k/ U& w  n8 Dfor i=1:length(x0) / I9 j1 x6 z; @: f6 k& u9 y3 `    for j=1:i / H4 I" V2 l6 y6 b- y- s2 B* ^        x1(i)=x1(i)+x0(j);     end ; @% Q7 {5 I+ `' `$ _- mend 2 V# u/ T/ K% c. `" l/ ~3 ?for i=1:length(x0)-1     B(i,1)=(-1/2)*(x1(i)+x1(i+1)); 9 y2 @; V1 n. R# H4 P$ O0 |    B(i,2)=1;     yn(i)=x0(i+1); end HatA=(inv(B'*B))*B'*yn % GM(1,1)模型参数估计 6 `0 h( Z; B; k1 ^for k=1:length(x0)+T ; i$ K5 K, W. T! ^3 S# k: C, ^    Hatx1(k)=(x0(1)-HatA(2)/HatA(1))*exp(-HatA(1)*(k-1))+HatA(2)/HatA(1); end 6 u" \* V; N. Y' R0 pHatx0(1)=Hatx1(1); ) V$ _" C* O1 _* f; _for k=2:length(x0)+T 0 }4 q5 }4 n' d! m9 \  ^    Hatx0(k)=Hatx1(k)-Hatx1(k-1);%累计还原得到历史数据的模拟值 ; `' ^  ^3 \- F' ]end for i=1:length(x0) %开始模型检验     epsilon(i)=x0(i)-Hatx0(i); , O/ W+ _# P6 M( l    omega(i)=(epsilon(i)/x0(i))*100; " Z) [: }5 m1 u7 Yend : o+ z7 U3 K+ u$ P% x0;Hatx0;epsilon;omega;  %必要时去掉％得到各种数据 " @4 e  D$ L" L( Z, C# e  J- Ec=std(epsilon)/std(x0);p=0; for i=1:length(x0)     if abs(epsilon(i)-mean(epsilon))<0.6745*std(x0)         p=p+1; ; I- @! }  R$ \1 ~! ?# Y    end end p=p/length(x0)   C, C/ @3 r0 ^0 @, c4 Uif p>0.95 & c<0.35     disp('The model is good,and the forecast is:'),     disp(Hatx0(length(x0)+T)) ) M- \1 `1 D0 `2 a3 r4 W$ A  }elseif p>0.85 & c<0.5 + J( H  t. c* t9 i- g    disp('The model is eligibility,and the forecast is:'), * Z  W' c4 q/ V' J    disp(Hatx0(length(x0)+T)) " ?9 |' r, F( Y, M  Kelseif p>0.7 & c>0.65     disp('The model is not good,and the forecast is:'), 9 `" \+ d8 a% i) {4 \    disp(Hatx0(length(x0)+T)) else p<=0.7 & c>0.65 / v0 d, h8 `! q1 L. a7 A    disp('The model is bad and try again') 3 @! G2 T0 R" Fend for i=1:length(x0) 8 \' @- x5 o7 V& C: l    Hatx00(i)=Hatx0(i); 4 j! V3 u" O8 x, o* [4 W* t7 gend 0 l; Z, X# L$ M4 L% f5 }- v7 Yz=1:length(x0); plot(z,x0,'-',z,Hatx00,':') %将原始数据和模拟值画在一个图上帮助观察 text(2,x0(2),'History data: real line') $ T8 }6 B' `9 K. d/ K- H4 utext(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 1 p# f7 I: Y6 j0.053132975 0.089908836 0 H  C+ D( h7 S; V% @0.109153219 0.079331832 9 S4 P0 `0 j8 i) ~1 n/ u5 K0.342192598 0.099718142 6 J4 K( y6 D- u6 g" A# p  A0.135194823 0.109274037 0.08152013 0.067876355 0.064706843 0.055562197 6 F+ R3 C9 w* U, B9 |0.050848544 6 m3 Y4 d8 C+ \+ t]';

daiqiang5566

楼主莫急..

yysclshi

a= -0.0080
u= 0.0713
预测值
    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
初始值
    1.6209    0.0793    0.0523    0.0413    0.0218    0.0531
( c) @; c+ l7 B/ ]    0.0899    0.1092    0.0793    0.3422    0.0997    0.1352
; \  S" ]; |& G9 s    0.1093    0.0815    0.0679    0.0647    0.0556    0.0508
2 i+ L/ I1 x' m7 J5 Q9 I残差
         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
相对误差
- b1 s) o. M' M         0    0.0672    0.6297    1.0835    2.9741    0.6437
    0.0209    0.1870    0.1276    0.7365    0.0885    0.3223
    0.1548    0.1420    0.3826    0.4619    0.7162    0.8903
% F% R. @% q* T$ O" N+ l; u0 Z方差比
3 a7 G  u; H' Z7 L6 A  F3 ^  |    0.1869
p =
     1

yysclshi

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

杨晓敬

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

gxj820

很好啊！！顶顶！

一剑卡卡

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

jshzncd

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