灰色预测

zhouchang

format long;
X0=[357.35
4 o/ y9 t( P: x! _5 T/ J  {; c653.91
% B5 K0 k# f$ a  k1123.29
  G" j# Q' j& n6 G( z( g1601.79
1952.05
1977.59
1964.89
3 d1 l' {: e2 }8 Z9 ?1856.72
* I$ M; F1 n$ R) E1869.67
$ R& A: `% [) d0 l% e) F1984.31
8 k' v) U; k2 a! h& [2158.41
2828.87
4 d% Z4 r. E9 X2 H& X3 {5 ^+ \ ];
) S$ Q$ X3 V% o0 ~8 Pn=length(X0);
X1=[];
X1(1)=X0(1);
, [% Y" ?. D) ~( kfor i=2:n
: C. ?9 X8 X, s    X1(i)=X1(i-1)+X0(i);
end
9 z, L3 S! B& ?! Z- D" ^& s$ b2 i9 S  `for i=1:n-1
& L8 B3 m/ A, f5 T- Y    B(i,1)=-0.5*(X1(i)+X1(i+1));
8 _, s. {1 k) {0 F, Q* D' r' v    B(i,2)=1;
    Y(i)=X0(i+1);
end
% |& b" ]$ {& q: }8 calpha=(B'*B)^(-1)*B'*Y';
0 m0 [8 p* E5 v8 F  @7 f* ba=alpha(1,1);
b=alpha(2,1);
7 q& k4 n) f- ]" {6 x5 C6 h- xd=b/a;
c=X1(1)-d;
X2(1)=X0(1);
0 z& J" m' {7 ]: g# u/ Z1 ^X(1)=X0(1);
+ N( P+ r9 d4 a0 Z3 u% Vfor i=1:n-1
9 Y$ z; h. T) S9 }  l    X2(i+1)=c*exp(-a*i)+d;
- Y6 t& k& @0 d( w4 \    X(i+1)=X2(i+1)-X2(i);
end
for i=1:n
1 |5 D9 u" R5 F1 g# g+ _: E6 @+ [    error(i)=X(i)-X0(i);
& S2 v1 `4 E" ?! E: p' I    error1(i)=abs(error(i));
    error2(i)=error1(i)/X0(i);
' \( K* u$ m4 k+ M) eend
2 i' h' M$ g) D2 b' p* J8 Z: Y1 ]# }g=std(error1)/std(X0)
5 E3 q! x, a9 a  h! T5 a' Be=c-c*exp(a)
- I1 a7 I4 I& _$ L0 \for i=1:18
    X2(i+1)=c*exp(-a*i)+d;
  U- Z4 O( r4 A- [  ]/ t% o4 e    X(i+1)=X2(i+1)-X2(i);
8 b: X# Q+ T3 T2 N3 n9 }9 T   
0 l4 ^" [7 k( J5 d; dend
  o& w' [& [' Q+ AX(19)=[];

9 o# E  W3 ~$ `' B- M
0 c/ y5 f! s* m% Y6 `/ Dplot(1992:2009,X,'b*-')
hold on

1 ]) ^5 m4 X% a% ~$ ?9 A5 U" [
X0=[
$ V0 M8 ]8 r1 a8 f8 h# ^4 D    357.35
* ^  ?) Y: K0 `, G, B653.91
- Z: v$ ], N) V  B% B/ Y1123.29
1 o. T6 ^( L: ?/ W* o1601.79
1952.05
1977.59
4 [( N$ W; f; V# Y; D/ p3 a/ h1964.89
$ d% M$ x1 r* I6 R! y9 M- r1856.72
7 i" N6 t9 Q. {; C& j1869.67
1984.31
2158.41
5 Q. t& S! _  K, ]2828.87
3084.66
) p! z3 d3 A7 B6 z2 m2 i3542.55
3925.09
! ^& [) |8 I8 B, q+ H( {4458.61
4829.45
5273.3
' |/ s: U! Q2 g9 V: m];
2 u" m* N- D5 i/ V6 U$ Y  |7 M
plot(1992:2009,X0,'g+-')
, i1 i: y( A2 n) rfor k=2:18
    l(k)=(X0(k)-X0(k-1))/(X(k)-X(k-1));
* s% y# v5 [/ u+ s, i; nend
' P6 B7 ^% a& P+ }( `l
ylabel('投资/亿元');
xlabel('上海       时间/年');
( Q  v/ V$ \% Y& @. ulegend('预测值',' 实际值 ' ,'Location','NorthWest')
; S7 {! Z* o  t; N: v/ D  d' ^%legend('预测值')

" e5 J" E+ q1 O: g

神秘了一场雨