灰色预测

zhouchang

format long;
5 O2 o) r  D' b3 S* N! DX0=[357.35
653.91
1123.29
- f+ z# t7 Z9 T. t# P- }- o7 M1601.79
1952.05
1977.59
$ {+ r* T- u6 d1 L- U1964.89
6 C. Y1 [$ G( F$ U2 z( E- d# ~2 h1856.72
1869.67
1984.31
2 r5 r$ l* U# e, i2158.41
2828.87
- j% S+ H6 B6 w8 g% Z' P; E. [ ];
4 @" x4 ?$ O& Dn=length(X0);
7 K9 e8 u, H: I& mX1=[];
6 Q0 C, `3 g5 n6 I5 f8 m# tX1(1)=X0(1);
for i=2:n
    X1(i)=X1(i-1)+X0(i);
end
' l2 _6 }' g% E; s2 Tfor i=1:n-1
    B(i,1)=-0.5*(X1(i)+X1(i+1));
# Z# B0 W, R/ a( H    B(i,2)=1;
    Y(i)=X0(i+1);
end
. ~- v/ N. N  Y; zalpha=(B'*B)^(-1)*B'*Y';
! G5 J2 `. S7 ya=alpha(1,1);
b=alpha(2,1);
d=b/a;
c=X1(1)-d;
X2(1)=X0(1);
% [5 L+ k/ B8 }- L" k. p) M6 mX(1)=X0(1);
for i=1:n-1
    X2(i+1)=c*exp(-a*i)+d;
! q) v! M& S; k( {    X(i+1)=X2(i+1)-X2(i);
end
1 z( Q7 y- R! {- q: B+ ^for i=1:n
    error(i)=X(i)-X0(i);
    error1(i)=abs(error(i));
    error2(i)=error1(i)/X0(i);
- F7 W6 B$ X* H& k+ I; Vend
g=std(error1)/std(X0)
3 s8 W4 m9 R' m0 o, _: j) De=c-c*exp(a)
for i=1:18
* h% X/ B" g* d# k! D2 Y+ n    X2(i+1)=c*exp(-a*i)+d;
    X(i+1)=X2(i+1)-X2(i);
$ Z2 F3 Y* a1 ~7 E$ H& L   
end
X(19)=[];

1 B6 n9 y6 @9 j3 Q8 Z1 C+ d- Y" k1 I
- x: W$ w( [4 C3 D, fplot(1992:2009,X,'b*-')
3 K& P9 h  w% M8 e7 t7 {hold on

+ X# |  ~3 y3 {& ?: ^2 e: ~
- U) E3 E" N# Y. `) F& ?$ EX0=[
9 P4 O& q" [2 K7 g1 M    357.35
# M! G! Z6 J; K' W& \# I653.91
3 ?0 @$ |% `+ C8 U: I1123.29
1601.79
) S* `6 O3 e( W+ O+ F" i& E0 u* t5 ]1952.05
1977.59
1964.89
6 h0 Z) \9 Y2 G. f: Y3 O. ?1856.72
1869.67
1984.31
, I3 R( q% Q( F1 F  [3 _2158.41
2828.87
3084.66
. }) p8 H) q" t- l) c9 K% C3542.55
: s8 E; `% ~9 n) O# K1 K, u3925.09
4458.61
4829.45
5273.3
* H2 T* K4 t3 T];
/ ?0 _. I/ e1 X0 ^. C
3 U1 e+ @1 E( p  o% N8 Dplot(1992:2009,X0,'g+-')
- r& @- h# z* S) s. ]7 ~1 Ofor k=2:18
    l(k)=(X0(k)-X0(k-1))/(X(k)-X(k-1));
end
l
, N; ?7 C( f: w  |: Z; Rylabel('投资/亿元');
xlabel('上海       时间/年');
7 Q8 x" y+ P# g$ [- `8 Flegend('预测值',' 实际值 ' ,'Location','NorthWest')
%legend('预测值')

6 ^4 N2 @& |3 m$ I

神秘了一场雨