灰色预测

zhouchang

format long;
# I5 z1 F6 y# g8 K" F" @& |X0=[357.35
653.91
1123.29
1601.79
1952.05
1977.59
1964.89
1856.72
' S) B7 s( m7 S- F' B3 _1869.67
" Z5 i, L2 s% Q, D( ~8 u$ ~1984.31
8 E8 H5 X6 K. X2158.41
$ g  c( X6 g' c( G( n- n& o2828.87
];
n=length(X0);
$ T  o; H  Q) n& E0 c' ?X1=[];
6 l, w8 A6 @7 R' X8 k. iX1(1)=X0(1);
for i=2:n
6 m$ h9 S7 j, d2 p2 p    X1(i)=X1(i-1)+X0(i);
/ c' y3 f4 T9 J4 lend
for i=1:n-1
    B(i,1)=-0.5*(X1(i)+X1(i+1));
    B(i,2)=1;
& W' q% q! W, l! t    Y(i)=X0(i+1);
end
alpha=(B'*B)^(-1)*B'*Y';
: m, @+ C7 c7 O0 @( a( Ta=alpha(1,1);
0 r/ R1 o( [1 u. G' O* y  ~5 ab=alpha(2,1);
d=b/a;
c=X1(1)-d;
X2(1)=X0(1);
1 l$ ~2 M+ U1 d' U3 j! _X(1)=X0(1);
for i=1:n-1
' v# y8 ^% T+ E# o) A5 _    X2(i+1)=c*exp(-a*i)+d;
    X(i+1)=X2(i+1)-X2(i);
end
  _. c) D! ]0 D  o$ Z- t; [for i=1:n
8 l6 O3 d+ k  }! l0 Y- a# G    error(i)=X(i)-X0(i);
    error1(i)=abs(error(i));
( }+ H* o* ~5 q    error2(i)=error1(i)/X0(i);
' R! W  M* X) y" a$ ~; Aend
g=std(error1)/std(X0)
e=c-c*exp(a)
% t$ y- r/ P7 f7 t$ H4 f. p  dfor i=1:18
# j9 e3 p( }$ e5 }1 f' c    X2(i+1)=c*exp(-a*i)+d;
    X(i+1)=X2(i+1)-X2(i);
   
end
X(19)=[];
5 v6 u: ~- N0 u  R
2 x0 [/ `, C" Y. Z* o5 ]
# H' ]5 N* y7 t) L3 E/ Q2 Nplot(1992:2009,X,'b*-')
: D& s3 I6 n* _0 m( [hold on
& ^$ M- |& |7 r4 o( d# V: p
( R# ?8 ~% l- t/ }7 ]( s% W
X0=[
    357.35
653.91
" l' W, S7 ]$ P1123.29
1601.79
1952.05
1977.59
1964.89
; X" N0 I3 s" `! [$ }( R1856.72
: z" w3 A, g2 B) ?1869.67
1984.31
2158.41
2828.87
3084.66
3542.55
9 H# h5 U- n- W* L' h& K& |2 W3925.09
4 h6 y) \" |8 M; S- r' H% w  f4458.61
4829.45
6 V0 C& E# C: P7 w4 {5273.3
];

plot(1992:2009,X0,'g+-')
for k=2:18
    l(k)=(X0(k)-X0(k-1))/(X(k)-X(k-1));
1 x0 \( f2 A* t, T) Rend
l
4 l0 n, G. Y& v/ iylabel('投资/亿元');
2 K" C" ^! R& ixlabel('上海       时间/年');
legend('预测值',' 实际值 ' ,'Location','NorthWest')
%legend('预测值')
, Z7 z/ d* e% {% k: Y/ I3 t

神秘了一场雨