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

xinzhiyong

GM(1,1)灰色模型的程序实现function GM1=fungry1(x0) %输入原始数据x0 / ~$ ?2 u  o% K0 pT=input('T=');%从键盘输入从最后一个历史数据算起的第T时点 / H: e; j* W# n5 }& h, x7 ax1=zeros(1,length(x0));B=zeros(length(x0)-1,2); yn=zeros(length(x0)-1,1);Hatx0=zeros(1,length(x0)+T); 4 N* Y, D5 A4 [# ?( _3 U. pHatx00=zeros(1,length(x0));Hatx1=zeros(1,length(x0)+T); epsilon=zeros(length(x0),1);omega=zeros(length(x0),1); . D' u- z6 t* N) Y/ Hfor i=1:length(x0)     for j=1:i         x1(i)=x1(i)+x0(j); 2 A1 v% C% e8 [3 {6 Q9 T    end end ; z% U$ O% Z% p( {' ffor i=1:length(x0)-1     B(i,1)=(-1/2)*(x1(i)+x1(i+1)); & n  }( K4 T9 \& \9 y! e% N2 }    B(i,2)=1;     yn(i)=x0(i+1); end 3 d5 N) H7 P1 a* n: w! @, p7 m4 bHatA=(inv(B'*B))*B'*yn % GM(1,1)模型参数估计 for k=1:length(x0)+T     Hatx1(k)=(x0(1)-HatA(2)/HatA(1))*exp(-HatA(1)*(k-1))+HatA(2)/HatA(1); end ) m3 T- U4 L+ s9 W* Q) i) E( Q' ~) YHatx0(1)=Hatx1(1); for k=2:length(x0)+T     Hatx0(k)=Hatx1(k)-Hatx1(k-1);%累计还原得到历史数据的模拟值 end % s1 L: R) Z9 Bfor i=1:length(x0) %开始模型检验 + V2 u9 P4 `; M, I    epsilon(i)=x0(i)-Hatx0(i); / k4 v( ~; q$ a    omega(i)=(epsilon(i)/x0(i))*100; 1 M! t+ D, H& ?; @+ y8 X# u2 n3 Pend   O3 x7 Q% y, P% x0;Hatx0;epsilon;omega;  %必要时去掉％得到各种数据 c=std(epsilon)/std(x0);p=0; for i=1:length(x0) ; V' l1 b: r5 Z    if abs(epsilon(i)-mean(epsilon))<0.6745*std(x0)         p=p+1; $ G- V) }& M! L: Z+ v    end end 5 b0 s. |  b' v" V7 v4 _0 I! zp=p/length(x0) ; L0 X) y, n4 H4 N6 u5 jif p>0.95 & c<0.35 : d& Y7 S9 {. [9 E1 M1 G    disp('The model is good,and the forecast is:'),     disp(Hatx0(length(x0)+T)) ( i' L2 e: Q( ]4 m. yelseif 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:'), . R% N) F( s  j    disp(Hatx0(length(x0)+T)) else p<=0.7 & c>0.65     disp('The model is bad and try again') * l8 J& q. ?0 x9 ]end , P1 X- Y' n7 |: Z$ ifor i=1:length(x0) / f5 |- x6 X2 j# c    Hatx00(i)=Hatx0(i); end z=1:length(x0); / x" }0 |9 ~6 ~! h) q! m; xplot(z,x0,'-',z,Hatx00,':') %将原始数据和模拟值画在一个图上帮助观察 % w* y$ m$ z- N1 ytext(2,x0(2),'History data: real line') text(length(x0)/2,Hatx00(length(x0))/2,'Simulation data:broken line') endT=input('T=');%从键盘输入从最后一个历史数据算起的第T时点????是指什么啊，请大哥们，大姐们教一下，我急用，请快，谢谢我的初始值x0=[1.620938526 0.07925621 0.052318818 2 F2 j' \# l# E5 Z; o) d& U0.041252502 0.021800479 0.053132975 0.089908836 0.109153219 0.079331832 0.342192598 0.099718142 0.135194823 0.109274037 8 H5 ]5 d5 |* H. y# J' E2 }0.08152013 0.067876355 / X  e; F4 u; E+ A, P9 x9 R2 L0.064706843 + H0 q0 C7 g- I2 y7 W0.055562197 6 r/ Q6 f5 Y0 j  M+ m3 V3 T1 m0.050848544 ' f$ T9 ~/ }( h3 ]. U]';

daiqiang5566

楼主莫急..

yysclshi

a= -0.0080
6 z- s+ C0 v9 @9 h' V+ }) p, x" Pu= 0.0713
% k1 e$ @$ H8 ?  P预测值
" X- S& a! C  B1 @, v  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 s$ `: K- f: k0 k! [: p  g6 p0 b; E    0.0924    0.0931    0.0938    0.0946    0.0954    0.0961
初始值
    1.6209    0.0793    0.0523    0.0413    0.0218    0.0531
. P& f, Y* `; [! \! |! N    0.0899    0.1092    0.0793    0.3422    0.0997    0.1352
    0.1093    0.0815    0.0679    0.0647    0.0556    0.0508
残差
# y; y" _. z* ]' M# D' O         0   -0.0053   -0.0329   -0.0447   -0.0648   -0.0342
5 F8 \* p" N1 `8 N& O0 q8 B# e    0.0019    0.0204   -0.0101    0.2520    0.0088    0.0436
# `% V& R7 V0 N. \+ k; a    0.0169   -0.0116   -0.0260   -0.0299   -0.0398   -0.0453
7 o: n% g6 P" M9 _; {6 D相对误差
         0    0.0672    0.6297    1.0835    2.9741    0.6437
" e' D9 \* {0 ^7 ]6 h  ~2 \+ R    0.0209    0.1870    0.1276    0.7365    0.0885    0.3223
! V* Q% v4 i" E    0.1548    0.1420    0.3826    0.4619    0.7162    0.8903
方差比
4 I/ X6 {( }) P, Z    0.1869
3 T7 C0 _- m5 s& ^2 Lp =
9 Q4 v! a$ a6 L" f: V, e     1

yysclshi

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

杨晓敬

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

gxj820

很好啊！！顶顶！

一剑卡卡

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

jshzncd

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