灰色预测

zhouchang

format long;
X0=[357.35
653.91
# ?( P4 m3 {; r4 u. _: h1123.29
1601.79
1952.05
1977.59
% x5 }  J2 o1 p' h6 f- U1964.89
  V, z+ a, N  m8 M" h$ w$ o1856.72
7 k2 j8 c( O& j. v9 }1869.67
+ f. u) t  y$ |# w( t1984.31
! {" x; D4 Z  |% y" H5 P( K2158.41
* f+ X  |$ @9 \% [+ R+ z  b0 W; q2828.87
4 \# J: y/ {9 m ];
n=length(X0);
4 P+ f# M! s1 p- cX1=[];
- A3 T# I% ]$ p2 U7 cX1(1)=X0(1);
, Q' V$ @. }/ a$ O2 R! F" h: Lfor i=2:n
6 I/ V" H& c, e" f    X1(i)=X1(i-1)+X0(i);
end
% @. x1 S0 `/ F/ H* j" ^6 Jfor i=1:n-1
    B(i,1)=-0.5*(X1(i)+X1(i+1));
$ |1 t  |# }  o5 n" F    B(i,2)=1;
    Y(i)=X0(i+1);
end
alpha=(B'*B)^(-1)*B'*Y';
a=alpha(1,1);
b=alpha(2,1);
d=b/a;
c=X1(1)-d;
X2(1)=X0(1);
X(1)=X0(1);
for i=1:n-1
( |7 h- V8 f2 L    X2(i+1)=c*exp(-a*i)+d;
& j! l+ I4 S9 M    X(i+1)=X2(i+1)-X2(i);
end
for i=1:n
3 K" ?- W4 w$ b  T* {    error(i)=X(i)-X0(i);
    error1(i)=abs(error(i));
    error2(i)=error1(i)/X0(i);
end
/ s+ c9 L9 ]6 T9 ~: {g=std(error1)/std(X0)
2 B$ S3 s$ {- ~# ?. g7 be=c-c*exp(a)
for i=1:18
    X2(i+1)=c*exp(-a*i)+d;
    X(i+1)=X2(i+1)-X2(i);
6 z0 _" _; |; d* ^3 n2 {   
end
X(19)=[];


plot(1992:2009,X,'b*-')
+ m$ {% T, d- a( Q# Bhold on


2 `. U6 w7 g3 ?+ r* L1 T5 CX0=[
- m0 V; D) b( A, D1 i    357.35
& J2 R  g& ^2 D- N8 q653.91
1123.29
1601.79
1952.05
1977.59
1964.89
  l4 Y; \" S3 q8 Q: i3 I1856.72
1869.67
) i6 y+ L% p8 V9 G5 o' o1984.31
2158.41
2828.87
4 u* U0 n4 I  x: ~3084.66
3542.55
3925.09
4458.61
& y* r0 f) o$ x0 [+ z4829.45
0 h9 h8 f, E  o$ l: E+ {5273.3
];

' m$ z, r8 N' ~plot(1992:2009,X0,'g+-')
( W8 g2 o) t1 B) Z) lfor k=2:18
5 X& j4 t4 g4 f' C4 U# U5 i/ a    l(k)=(X0(k)-X0(k-1))/(X(k)-X(k-1));
end
l
ylabel('投资/亿元');
) A) b+ k( v7 _xlabel('上海       时间/年');
% Z  G/ B1 z0 ^9 i$ Q* Slegend('预测值',' 实际值 ' ,'Location','NorthWest')
%legend('预测值')

神秘了一场雨