灰色预测

zhouchang

format long;
+ ^5 U3 p* `& y5 T4 ^X0=[357.35
! }' e. b+ O" }( @2 O653.91
& j" a8 r) ~& X. X1123.29
+ D& ~6 }/ H) L1601.79
1952.05
1977.59
1964.89
1856.72
1869.67
& U/ n6 c& u+ O2 b% u1 @- B7 `1984.31
2158.41
2828.87
$ {! C$ g' j' R! N+ P: x ];
3 ]( T; W2 i' E- Q2 x; R3 \n=length(X0);
+ @, }# _* V) _3 j- d) }4 iX1=[];
X1(1)=X0(1);
6 I& E* M- u9 zfor i=2:n
  S; a; T- `* @! C% k: }  ?    X1(i)=X1(i-1)+X0(i);
% ^) r& B: R2 p0 I$ b. Eend
for i=1:n-1
    B(i,1)=-0.5*(X1(i)+X1(i+1));
    B(i,2)=1;
    Y(i)=X0(i+1);
end
alpha=(B'*B)^(-1)*B'*Y';
! K! f! D; o/ r8 P6 A+ V8 Ga=alpha(1,1);
b=alpha(2,1);
d=b/a;
c=X1(1)-d;
0 l/ v7 [  a1 q: s1 C% D+ @X2(1)=X0(1);
X(1)=X0(1);
for i=1:n-1
  q3 ^+ {1 q+ G( S/ J    X2(i+1)=c*exp(-a*i)+d;
    X(i+1)=X2(i+1)-X2(i);
end
for i=1:n
    error(i)=X(i)-X0(i);
    error1(i)=abs(error(i));
* E; U5 k  w; {    error2(i)=error1(i)/X0(i);
end
. k# j1 s  Z8 f( k, Z( ng=std(error1)/std(X0)
5 V% b# S- E- R: s8 {e=c-c*exp(a)
for i=1:18
    X2(i+1)=c*exp(-a*i)+d;
    X(i+1)=X2(i+1)-X2(i);
   
end
X(19)=[];
& i1 y1 x. l: T% m5 w1 y* |5 Y8 x

plot(1992:2009,X,'b*-')
9 E4 x- n/ w3 k5 }$ Hhold on


X0=[
    357.35
653.91
1123.29
1601.79
4 P3 K7 R9 v1 f* `- J" @: u8 w1952.05
0 l, g1 z, j, w" m* p3 ^1 T1977.59
6 M2 h* Y, g, K- O' z1 X1964.89
( h2 j! f* Y+ n- Y" k1856.72
/ [: R6 ?% W% k5 a; r1869.67
1984.31
+ g5 L0 D! \! D$ B0 i, E2158.41
8 W. Z" ?# G8 b4 }6 @7 t" ]1 S2828.87
7 H" o' d* j) @  b3084.66
+ p! h! s* Q/ q: a' B3542.55
, E6 G( T1 j" Z6 p  V8 C; `& X3925.09
* O8 i# U1 L4 g4458.61
4829.45
5273.3
];
2 Z. x5 y) d" r5 O' ~1 S/ W% u# b
( F! b7 G0 }0 q1 c% Qplot(1992:2009,X0,'g+-')
for k=2:18
& m: {) W) Q' z- I" }! M    l(k)=(X0(k)-X0(k-1))/(X(k)-X(k-1));
- D: ^# |% N* t9 B& R  qend
9 }& Z9 g! }2 e& d' d$ T% E3 ]  Gl
ylabel('投资/亿元');
/ v  X: ?! g& G- x& Uxlabel('上海       时间/年');
legend('预测值',' 实际值 ' ,'Location','NorthWest')
" T- a0 H/ g6 V: X%legend('预测值')
  A% l2 l- n  A
, I$ y  d/ ~* b4 f/ r* U' |1 B

神秘了一场雨