灰色预测

zhouchang

format long;
/ E9 F, P( J" W# _+ K1 d7 ^X0=[357.35
653.91
1123.29
1601.79
1952.05
1977.59
1964.89
/ c. b7 W9 W6 e, D1856.72
0 O% Y/ r8 K$ i* h1869.67
1984.31
2158.41
2828.87
2 G$ S$ ~& t, J ];
n=length(X0);
/ n3 n7 K6 _# P/ ]! L) u" RX1=[];
, Q* p) F: v/ |4 Z1 |X1(1)=X0(1);
0 Q: L* ^8 X0 E4 ~; V' U5 cfor i=2:n
4 b9 J, [6 g' q    X1(i)=X1(i-1)+X0(i);
end
% f- w6 m- u9 C0 b+ b/ N6 x3 s; Gfor i=1:n-1
4 G* h4 l1 O& C: h7 s    B(i,1)=-0.5*(X1(i)+X1(i+1));
8 O, u. o! Z9 V5 B    B(i,2)=1;
* O- D8 @4 P/ `# M  a5 q    Y(i)=X0(i+1);
4 B5 A8 v, ]1 I6 eend
" ?  X- Z/ {$ Z7 ralpha=(B'*B)^(-1)*B'*Y';
0 u5 c8 ?. ?- C) U6 G& Aa=alpha(1,1);
b=alpha(2,1);
( L; s; r! m7 ~1 ud=b/a;
8 t- S: B! J) ]2 f4 M. l) J/ D) }c=X1(1)-d;
X2(1)=X0(1);
" q& m# S- C( k% YX(1)=X0(1);
  _' G( A/ P) y6 y* o, o1 Qfor i=1:n-1
    X2(i+1)=c*exp(-a*i)+d;
    X(i+1)=X2(i+1)-X2(i);
) O( {- U! l8 Xend
8 |& I3 m# Q- T, K5 m- ]9 a7 r" bfor i=1:n
5 y! E' q4 N# _8 `0 o0 S. `3 O3 ~8 O    error(i)=X(i)-X0(i);
    error1(i)=abs(error(i));
    error2(i)=error1(i)/X0(i);
$ v% O) X8 \" |& w2 J, bend
3 N/ W. k. {) c6 {, Z* z0 jg=std(error1)/std(X0)
e=c-c*exp(a)
for i=1:18
    X2(i+1)=c*exp(-a*i)+d;
# h# k# o% ^/ Y' w6 U4 k    X(i+1)=X2(i+1)-X2(i);
   
end
+ f& \+ Q- O6 S  IX(19)=[];
8 _) ]4 V' J# k- V$ `  F
( j8 s$ ~# ]+ X( o) m. s
- W  S/ L3 V$ ^- mplot(1992:2009,X,'b*-')
hold on
3 e0 G- c! f7 V$ Z# C- F4 `/ j1 E
% ?4 d* c! E, z; c! I
/ F6 ~! V1 Z# b9 h3 bX0=[
    357.35
653.91
$ L) e& J7 Q  M4 l( s. g1123.29
6 C, O& ^( n; ~7 a: {+ |1601.79
9 v9 ^% B, r. S) V8 g6 y: Q5 N) r7 w, j1952.05
& P, Y) Y% f. ^! ?1977.59
1964.89
% S( {0 }- k; @0 J1856.72
1869.67
' v# h3 h' z" |& c) S  N  Q# p* y1984.31
6 T4 ~1 C1 G" b2158.41
2828.87
3084.66
4 E6 e7 u- O1 K8 ^  Q: ?3542.55
& `. F7 C3 l: X  D" N! p3925.09
4458.61
4829.45
2 Y2 v+ n- T+ x$ }7 t9 Y5273.3
5 j. p# Z- V8 N, V, C];

2 k2 m# i. B& o; g2 jplot(1992:2009,X0,'g+-')
: e. V& g2 H3 Q  d* [+ Yfor k=2:18
' g' C5 ~3 c3 d1 P$ R) d: ?  X2 T6 K4 E    l(k)=(X0(k)-X0(k-1))/(X(k)-X(k-1));
0 g, d1 N2 m  \' F+ |2 u) g6 qend
) g, G  {& V9 }" e( V' il
ylabel('投资/亿元');
' r1 h( V3 z0 K$ |xlabel('上海       时间/年');
" m, d& C' Z8 @( zlegend('预测值',' 实际值 ' ,'Location','NorthWest')
%legend('预测值')
2 k# O8 R; P) x) ^! |

神秘了一场雨