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); 4 s" `% `9 f0 d7 d2 e* \9 Ryn=zeros(length(x0)-1,1);Hatx0=zeros(1,length(x0)+T); : A( g+ a* {4 {0 \7 H3 V6 mHatx00=zeros(1,length(x0));Hatx1=zeros(1,length(x0)+T); epsilon=zeros(length(x0),1);omega=zeros(length(x0),1); 8 S4 O* P$ f- U! I& @for i=1:length(x0) - N* c5 P0 @- i0 T* Q* u    for j=1:i 9 V) K1 l. {) A, X4 `# d+ A! P6 D        x1(i)=x1(i)+x0(j);     end end for i=1:length(x0)-1     B(i,1)=(-1/2)*(x1(i)+x1(i+1)); . r# A5 L4 r0 f9 u* K* c    B(i,2)=1;     yn(i)=x0(i+1); end / P2 l; i2 j! x$ X6 w' MHatA=(inv(B'*B))*B'*yn % GM(1,1)模型参数估计 / q: i1 E' i! yfor k=1:length(x0)+T $ P3 `9 b5 k6 k6 D! M4 @    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);%累计还原得到历史数据的模拟值 . F, O5 c6 O4 V3 rend for i=1:length(x0) %开始模型检验     epsilon(i)=x0(i)-Hatx0(i);     omega(i)=(epsilon(i)/x0(i))*100; end % x0;Hatx0;epsilon;omega;  %必要时去掉％得到各种数据 , ]! j. S( G% b* @6 Nc=std(epsilon)/std(x0);p=0; + e0 }6 V. z) N" G! K7 X' b" w, L0 afor i=1:length(x0)     if abs(epsilon(i)-mean(epsilon))<0.6745*std(x0) 2 |: n; z, ]/ z6 [        p=p+1;     end end 4 `$ c( o! V& k3 o3 m4 yp=p/length(x0) if p>0.95 & c<0.35 + ^2 \" B1 S, Y( ^2 p    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:'),   N3 r# K6 j" K- b& T    disp(Hatx0(length(x0)+T)) elseif p>0.7 & c>0.65 + a" |1 k' D1 w" {    disp('The model is not good,and the forecast is:'),     disp(Hatx0(length(x0)+T)) else p<=0.7 & c>0.65 8 |8 P0 F) ~/ @9 \* A7 p8 F    disp('The model is bad and try again') end $ b, D7 \0 z3 A! R9 B& bfor i=1:length(x0)     Hatx00(i)=Hatx0(i); end z=1:length(x0); ' R6 x5 E3 n; F5 fplot(z,x0,'-',z,Hatx00,':') %将原始数据和模拟值画在一个图上帮助观察 text(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 ( t  A7 b$ ^0 W  {0.041252502 0.021800479 / A- _$ Z) k9 S2 W( {3 C0.053132975 0.089908836 8 }# `2 E* X$ u/ S1 B9 R7 F/ O0.109153219 0.079331832 0.342192598 . U/ `+ {% k4 M; ^; L+ ~; y! M: y; {* ]$ ~0.099718142 4 E' u" \5 {& G) R- k/ o( t0.135194823 0.109274037 0.08152013 0.067876355 * m9 G* A: m$ C' s, r% J6 W0 W0.064706843 $ Q% i+ \% H$ L5 A5 {0.055562197 9 x. A; [+ L8 {5 Z7 x' c0.050848544 ]';

daiqiang5566

楼主莫急..

yysclshi

a= -0.0080
u= 0.0713
预测值
    1.6209    0.0846    0.0853    0.0859    0.0866    0.0873
6 C6 N4 `* O  P) ~5 @( G" R    0.0880    0.0887    0.0895    0.0902    0.0909    0.0916
( d1 m5 Y9 A4 F6 l; m    0.0924    0.0931    0.0938    0.0946    0.0954    0.0961
初始值
' n% \. d& q% C2 W; L2 }8 ^2 y8 o! S    1.6209    0.0793    0.0523    0.0413    0.0218    0.0531
    0.0899    0.1092    0.0793    0.3422    0.0997    0.1352
9 r# k& a$ p4 i% s    0.1093    0.0815    0.0679    0.0647    0.0556    0.0508
残差
  d7 c, t2 P* q# {5 y         0   -0.0053   -0.0329   -0.0447   -0.0648   -0.0342
. O, z4 O/ E3 k4 k0 d' v; n7 L+ [    0.0019    0.0204   -0.0101    0.2520    0.0088    0.0436
    0.0169   -0.0116   -0.0260   -0.0299   -0.0398   -0.0453
相对误差
) w3 J7 R3 p6 ~  R4 i         0    0.0672    0.6297    1.0835    2.9741    0.6437
    0.0209    0.1870    0.1276    0.7365    0.0885    0.3223
2 |) B" a8 _$ O" K* i1 h2 M$ o) P9 }5 E0 Q    0.1548    0.1420    0.3826    0.4619    0.7162    0.8903
0 c0 X4 w8 ]/ G$ t1 W+ {- s" G方差比
) U- f# K' A" S: {    0.1869
, }1 F8 N9 _7 X, Tp =
     1

yysclshi

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

杨晓敬

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

gxj820

很好啊！！顶顶！

一剑卡卡

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

jshzncd

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