灰色预测

zhouchang

format long;
X0=[357.35
) y$ X% w) X8 D1 e653.91
1123.29
3 S+ w8 E, q$ A7 Y1601.79
1952.05
& t: x3 M) f0 O$ Y( O- i1977.59
1964.89
1856.72
1869.67
1984.31
, X# h) A/ z& H: C6 Q2158.41
2828.87
! b. N( p1 d$ C ];
n=length(X0);
2 z% l  K! e# }X1=[];
X1(1)=X0(1);
for i=2:n
    X1(i)=X1(i-1)+X0(i);
# d# }* q+ p' l" g/ P7 r0 @end
for i=1:n-1
; p' y+ `0 U8 _5 P) q    B(i,1)=-0.5*(X1(i)+X1(i+1));
    B(i,2)=1;
1 E( P" t- W% o* F$ t/ \1 |    Y(i)=X0(i+1);
end
alpha=(B'*B)^(-1)*B'*Y';
a=alpha(1,1);
7 L. n" I9 F7 m% t: lb=alpha(2,1);
% M9 J  m: X3 R: s: N$ K* S! Q$ `" j! md=b/a;
c=X1(1)-d;
X2(1)=X0(1);
+ [" {/ K# b2 oX(1)=X0(1);
for i=1:n-1
    X2(i+1)=c*exp(-a*i)+d;
8 m* H- ^' C, u; j    X(i+1)=X2(i+1)-X2(i);
end
for i=1:n
    error(i)=X(i)-X0(i);
    error1(i)=abs(error(i));
    error2(i)=error1(i)/X0(i);
8 X. o  ^0 U6 {) r$ e0 z# J& Pend
& G7 y! l0 g& F2 s8 O" eg=std(error1)/std(X0)
' v- x) N$ u, p+ N+ ~  K! Je=c-c*exp(a)
for i=1:18
" y; m7 I/ Q5 v: K' `* z1 k6 Z    X2(i+1)=c*exp(-a*i)+d;
    X(i+1)=X2(i+1)-X2(i);
" e* y. l( O/ z7 I   
& ^  A5 [/ C" R! k- M+ B7 |end
5 {8 S4 A. ~6 {4 x5 W4 oX(19)=[];
$ _- ?2 T2 L9 c6 G+ E" s

5 t4 m" K/ f* W0 Qplot(1992:2009,X,'b*-')
hold on
2 p* p) o( P9 W
2 ^& r7 o5 g) t' A
% E* m% I, H8 R0 l2 o5 GX0=[
    357.35
) {0 U- K& N: ]& Y: b3 y# J* \653.91
! F8 b/ t0 R. k$ ~* K; c1 s1123.29
1601.79
1952.05
$ D  c1 P4 n4 l' o1977.59
$ {6 C. w0 O% K) n1964.89
2 D! e9 \; g8 C: p1856.72
1869.67
1984.31
2158.41
: E+ y6 P+ S) O' `' a2828.87
3084.66
3542.55
/ i; i+ C& C0 T5 V7 W6 T2 V% `$ P3925.09
6 B& Z5 r! x8 Y9 P+ i% z4458.61
0 G" h. ?& ?* @4829.45
5273.3
  J- U. V6 i- F" q3 C& Q];
) [* X: L( Z4 ?5 ]+ {$ |
plot(1992:2009,X0,'g+-')
1 Z8 u  F7 H; n+ b3 afor k=2:18
, y  Y" J) i% r3 p    l(k)=(X0(k)-X0(k-1))/(X(k)-X(k-1));
3 j  e3 P! N, Cend
* {+ M- ], X& d' L/ C9 k* K2 Ql
+ J8 ^! z& k; g, bylabel('投资/亿元');
, j- b+ I- }; Rxlabel('上海       时间/年');
legend('预测值',' 实际值 ' ,'Location','NorthWest')
- k* W$ L6 S) z" y%legend('预测值')
" ?4 ?- U; J$ x+ g4 G1 p

神秘了一场雨