灰色预测

zhouchang

format long;
X0=[357.35
- n. y# O! L! P0 y653.91
1123.29
1 _; m& C# X3 L9 U* W1601.79
1952.05
5 \% m- J2 P& b% v) J" U+ {* j1977.59
: C, u) l) l$ J/ z) _6 `: U8 K1964.89
1856.72
1869.67
1984.31
" F* i  a0 n$ K& i2158.41
7 G  ^/ u! i3 ]2 P6 {! d# J8 s2828.87
];
, y+ G& i! c/ q& H& u7 c7 ]8 j7 An=length(X0);
X1=[];
X1(1)=X0(1);
1 a1 m' [  T" v, }; Z% wfor i=2:n
    X1(i)=X1(i-1)+X0(i);
" s1 P0 X) T4 j" |& k' w0 i7 tend
for i=1:n-1
3 P; W7 a7 T8 h: w  n/ O5 a' r    B(i,1)=-0.5*(X1(i)+X1(i+1));
    B(i,2)=1;
* }  Q, q( I3 y: k3 [    Y(i)=X0(i+1);
* k* Z5 t! S3 R8 D: }end
1 l4 D4 \' |* X9 D/ R; F8 valpha=(B'*B)^(-1)*B'*Y';
a=alpha(1,1);
b=alpha(2,1);
d=b/a;
c=X1(1)-d;
. y* r: W0 W0 N5 j: O" NX2(1)=X0(1);
* X9 k* c: [  d2 U4 D' h3 eX(1)=X0(1);
for i=1:n-1
. f- L0 x/ D9 a( ]6 y; p  q    X2(i+1)=c*exp(-a*i)+d;
    X(i+1)=X2(i+1)-X2(i);
end
7 k- V: [- n/ C3 v+ Pfor i=1:n
. w: x1 M5 b' G% h/ w    error(i)=X(i)-X0(i);
  n/ ^/ W0 _7 S- j1 O7 U/ w    error1(i)=abs(error(i));
1 k. p3 m6 K, H, Z    error2(i)=error1(i)/X0(i);
end
g=std(error1)/std(X0)
9 p. _7 A1 R1 Je=c-c*exp(a)
  C0 u: V+ B. E* b! S( W& l, ~  Pfor i=1:18
    X2(i+1)=c*exp(-a*i)+d;
    X(i+1)=X2(i+1)-X2(i);
   
end
4 P2 x: G7 V1 T! w  KX(19)=[];

. ]+ N$ u. p9 Z- i/ h
plot(1992:2009,X,'b*-')
hold on
8 d) l' D+ a' ?& x2 q$ }  C  }
! ?. f- T3 V6 l' J$ T; G8 J
X0=[
1 ^: ?1 o, T, {" Z! A: w    357.35
( u8 V& \- }2 C9 F% g2 G653.91
! x5 U8 V8 V( e+ q& i5 _1123.29
1601.79
. N  p' J+ E3 J( |1 }3 J7 F1952.05
5 A. M% E) {5 W& N1 R1977.59
4 Z5 V2 q/ ?+ {! W$ e# a  d3 l4 S& Y1964.89
1856.72
1869.67
( ~; `, j/ z. U7 ~, u( I4 x! G0 ?  |: Q1984.31
' ^7 ]! d1 s1 e/ a2158.41
2828.87
3084.66
. {) p; B2 b9 f3542.55
3925.09
4458.61
8 e: R2 U6 C. j. Y4 s. k4829.45
  _! k2 y# U0 {% [# Y" j5273.3
: k8 W- t2 I% t$ N$ G# e; w; ?];

plot(1992:2009,X0,'g+-')
for k=2:18
4 S% J: D9 u/ K! H. X1 `7 l" y3 o' \    l(k)=(X0(k)-X0(k-1))/(X(k)-X(k-1));
end
: l/ u+ q: O2 `l
ylabel('投资/亿元');
xlabel('上海       时间/年');
legend('预测值',' 实际值 ' ,'Location','NorthWest')
%legend('预测值')
3 m( f% g% |# X, @2 u+ h& t2 p
7 |% T6 |# y1 [% H: N: n) q' |" [0 w

神秘了一场雨