灰色预测

zhouchang

format long;
( o  q' y0 R1 z+ |X0=[357.35
) R; W/ D' E, E% G4 s6 q653.91
1123.29
2 t: F3 v! H0 i/ w' i1601.79
5 s: n% K8 S/ D& V0 `) g" Y1952.05
  `4 A0 h, T  e  T- M! W0 h3 v' a) ?1977.59
1964.89
1856.72
" d9 e, y# T; r$ }) m/ G9 E1869.67
2 Q# C2 W2 K3 Q$ S1 ?( c1984.31
2158.41
: x3 X- U; Q$ o, X3 I7 |  ^2828.87
];
  ]# \* r% e$ ]3 D$ A; t. g  qn=length(X0);
X1=[];
$ ~8 N) ]/ N( N; K/ m2 oX1(1)=X0(1);
: c+ m3 n$ |, D4 [for i=2:n
    X1(i)=X1(i-1)+X0(i);
$ V7 u1 N1 W# b; p' r, S7 Bend
0 M" u- [2 M4 {* h5 l# Mfor i=1:n-1
: q+ \' S- S" U+ [    B(i,1)=-0.5*(X1(i)+X1(i+1));
    B(i,2)=1;
    Y(i)=X0(i+1);
; d: q1 }2 ^$ tend
alpha=(B'*B)^(-1)*B'*Y';
a=alpha(1,1);
  U7 ~5 |- L( R" O/ r1 K* u4 O: Ab=alpha(2,1);
d=b/a;
c=X1(1)-d;
+ }$ ], n6 i- H1 ]& `, K/ DX2(1)=X0(1);
X(1)=X0(1);
for i=1:n-1
9 n: `2 I" h& m    X2(i+1)=c*exp(-a*i)+d;
    X(i+1)=X2(i+1)-X2(i);
end
# Z2 {( x# o5 X  {/ Ofor i=1:n
    error(i)=X(i)-X0(i);
    error1(i)=abs(error(i));
    error2(i)=error1(i)/X0(i);
end
/ Q8 X) C0 J, i2 Jg=std(error1)/std(X0)
e=c-c*exp(a)
for i=1:18
9 c% f5 \, q; N8 i( S' Y    X2(i+1)=c*exp(-a*i)+d;
" c4 N; J# v4 z# F/ n    X(i+1)=X2(i+1)-X2(i);
/ @4 A2 i- j! r* t; P; h   
end
1 @% H. b! m, r' Y" xX(19)=[];
0 k) ]5 I/ o5 F

plot(1992:2009,X,'b*-')
hold on
- J) X. |: u" C4 h* r8 h# i
: W) `$ H3 f1 ~& ^
X0=[
    357.35
0 @% m; F* R' O/ n; N653.91
: S0 l6 s7 {4 V" O" w1123.29
1601.79
. v  e( [, E% }' O! B1952.05
1977.59
1964.89
1856.72
1869.67
, G4 N2 h: ?8 y% I8 B& u1984.31
2158.41
2828.87
3084.66
3542.55
3925.09
4458.61
4829.45
5273.3
! U' k4 K) M* Y1 J];

! }, G- t# I) I3 D7 I  R5 C9 V0 oplot(1992:2009,X0,'g+-')
; s, _; D) l6 rfor k=2:18
    l(k)=(X0(k)-X0(k-1))/(X(k)-X(k-1));
9 F9 I7 \2 P6 C# I; ^& y* J; gend
3 k" X* v4 Q3 K" nl
ylabel('投资/亿元');
xlabel('上海       时间/年');
0 m- R2 f& D7 }( z8 f. flegend('预测值',' 实际值 ' ,'Location','NorthWest')
%legend('预测值')

神秘了一场雨