灰色预测

zhouchang

format long;
X0=[357.35
653.91
* s# e2 p3 ?2 d8 T. V  w1123.29
1601.79
1952.05
2 X$ j: v/ a! T4 v1977.59
1964.89
/ b' ^8 M2 C8 O' m  \, O4 T1856.72
1869.67
1984.31
2158.41
4 [/ O3 B3 n) E+ C3 p! @2828.87
- b: I/ n& F! s$ u8 i ];
n=length(X0);
3 [- {" N. Z) T' F6 x: R; B. FX1=[];
5 U: f7 s; ^' A7 `2 \4 Q* e% nX1(1)=X0(1);
, a4 F- V2 e. X0 P5 O& Nfor i=2:n
    X1(i)=X1(i-1)+X0(i);
9 l1 H& B. ?7 K7 b- Fend
0 N# D4 ]8 ?: _: {! Ifor i=1:n-1
    B(i,1)=-0.5*(X1(i)+X1(i+1));
    B(i,2)=1;
    Y(i)=X0(i+1);
  e) n6 g$ j! E$ ^+ a+ X  nend
3 H% {. u! f2 A. M  N, v1 `alpha=(B'*B)^(-1)*B'*Y';
a=alpha(1,1);
b=alpha(2,1);
d=b/a;
c=X1(1)-d;
X2(1)=X0(1);
X(1)=X0(1);
/ i. T! R) k3 K  jfor i=1:n-1
! X0 U  t. d& {5 v$ d6 Z    X2(i+1)=c*exp(-a*i)+d;
    X(i+1)=X2(i+1)-X2(i);
end
for i=1:n
5 p7 ?- m2 J6 b. G" g6 h& y" G    error(i)=X(i)-X0(i);
    error1(i)=abs(error(i));
    error2(i)=error1(i)/X0(i);
$ W# e! a1 x# ]& Iend
g=std(error1)/std(X0)
7 ^% y  W3 {. U  g5 Q2 E! De=c-c*exp(a)
for i=1:18
6 n1 m2 R7 t  P2 u; G0 H/ [    X2(i+1)=c*exp(-a*i)+d;
6 t9 ]* S" `; k/ N$ M    X(i+1)=X2(i+1)-X2(i);
. v* ^7 K) @7 I3 ?# ^' {! W3 N3 }   
$ h- C9 P+ H+ J7 E, Yend
X(19)=[];
& M( S3 U8 t& O( ?  F' N& P6 a& \

( @7 m+ w! f$ j- R0 E; dplot(1992:2009,X,'b*-')
hold on
. I1 T- m; b1 l% \2 g- P
& p5 X7 {6 u. j
' O, @' a* h6 D# a, m- wX0=[
    357.35
653.91
2 ^) ?+ [  P9 q2 O# _6 u+ l! }1123.29
" i& ?) f) N7 ]* u1601.79
1952.05
- B( p1 d3 j3 M- \1977.59
! C1 h; o  E+ m% N& [8 {1964.89
1856.72
7 q) T* Y# h5 \8 u8 d1869.67
1984.31
2158.41
7 k! G7 ^! X% y1 q. r3 m2 B' \3 f2828.87
1 H$ u8 g# |* t3084.66
- `$ @; r3 f8 D2 f% r3542.55
3925.09
5 w* w* B( |' Q$ e, d/ v4458.61
4829.45
5273.3
];

plot(1992:2009,X0,'g+-')
for k=2:18
    l(k)=(X0(k)-X0(k-1))/(X(k)-X(k-1));
end
6 ~$ U' a- \4 a5 m: _l
ylabel('投资/亿元');
( e! ^4 d" O2 F' Hxlabel('上海       时间/年');
( _5 G( k, X' _) V3 M/ N9 ~5 Elegend('预测值',' 实际值 ' ,'Location','NorthWest')
& m. o6 H: E# `( m! c% Q%legend('预测值')

神秘了一场雨