灰色预测

zhouchang

format long;
X0=[357.35
653.91
7 k/ ~4 f( e8 o1 q1123.29
; c) V( F! n, P2 G$ H: d1601.79
1952.05
1977.59
4 v  D) t3 k' ]1964.89
+ u* a- d# L/ V$ ?2 p1856.72
1869.67
! _' M0 X0 \* i- v1984.31
$ I9 b6 }) u+ _9 ^) c2158.41
2828.87
% @3 b9 m; U* r ];
, ^( _4 v- v: p- l: cn=length(X0);
- n4 T5 A# M6 gX1=[];
; a" r7 y1 [1 D5 e3 Y% u, M3 cX1(1)=X0(1);
! q$ h: b" m6 j, r8 q* [$ N! p* Lfor i=2:n
; L2 G+ j4 p* Z# v    X1(i)=X1(i-1)+X0(i);
end
$ g8 L3 w  K1 Y- T2 Q( s5 n6 J. ifor i=1:n-1
    B(i,1)=-0.5*(X1(i)+X1(i+1));
    B(i,2)=1;
    Y(i)=X0(i+1);
0 M' D; j, L5 O/ G. P+ a- fend
alpha=(B'*B)^(-1)*B'*Y';
  c! x/ N  B& v( J: S: x& Ka=alpha(1,1);
b=alpha(2,1);
& o% ]- g8 O* ]9 ]/ ~$ Ud=b/a;
c=X1(1)-d;
X2(1)=X0(1);
7 S. E- r0 C! T' b6 r3 {9 f0 zX(1)=X0(1);
for i=1:n-1
    X2(i+1)=c*exp(-a*i)+d;
; u, j. T2 A: ~. e    X(i+1)=X2(i+1)-X2(i);
end
+ ~' f& G( L/ `for i=1:n
    error(i)=X(i)-X0(i);
* Z! [& _; d) R2 n3 _: O! N) k    error1(i)=abs(error(i));
- |+ b/ }9 V; G; ]& Y5 Y9 ~/ Z    error2(i)=error1(i)/X0(i);
8 P% p0 l- z/ i% i/ L$ s/ g7 Eend
4 T' Q' ~' E# t) k  g/ Ag=std(error1)/std(X0)
e=c-c*exp(a)
for i=1:18
    X2(i+1)=c*exp(-a*i)+d;
    X(i+1)=X2(i+1)-X2(i);
* l3 O4 D- m9 E   
end
0 Q6 N. @* L  fX(19)=[];
5 Y% o1 o8 U; a2 ^- x8 Y+ m7 H

% \4 m1 q0 H, Y; g/ xplot(1992:2009,X,'b*-')
hold on


* z+ J& _1 P; t# U* P) M  I" eX0=[
    357.35
653.91
/ X6 v( D4 `0 k) i# }& u1123.29
" y' ]9 u1 [' I% O3 y1601.79
( V' L2 J! b7 f: _# v3 v$ c1952.05
1977.59
1964.89
" h- U' E" |- k4 R; R1856.72
1869.67
# g9 {* g: n* l+ ]  i7 m$ E1984.31
2158.41
% |$ L: U5 q! G4 Z' a$ D6 h# A5 S2828.87
# Y/ C9 \: P( g& T: l* F! [3084.66
3542.55
; l% X( t4 A# v+ A; Y) `  [3925.09
4458.61
4829.45
/ \& R$ P( d9 o! I& h5273.3
];

% t1 }! t/ J  j3 p" V* yplot(1992:2009,X0,'g+-')
, \/ F; n2 i4 K% `for k=2:18
    l(k)=(X0(k)-X0(k-1))/(X(k)-X(k-1));
end
% [* F- I6 X& E/ ~- Tl
ylabel('投资/亿元');
  Z/ b- N7 b1 s8 Lxlabel('上海       时间/年');
legend('预测值',' 实际值 ' ,'Location','NorthWest')
%legend('预测值')

: c$ L; |7 B6 J

神秘了一场雨