灰色预测

zhouchang

format long;
. t* B! l& q/ p' _8 O$ Q0 T& vX0=[357.35
. G% I. D8 ~: @. B; p653.91
1123.29
  ^) I, N- _9 S: E) |1601.79
1952.05
* w0 J" r8 X5 P& k1977.59
1964.89
1856.72
1869.67
: M7 w: p  i: l3 Z% U9 G) Z5 |- w1984.31
& i- X; s4 f# I1 k: s0 U" T# \* o, w2158.41
2828.87
];
5 ^5 h3 P4 o+ C# _! n3 xn=length(X0);
' u; L$ g/ Y; I9 s/ y) ~X1=[];
X1(1)=X0(1);
/ J* ?  t- N5 p! Tfor i=2:n
9 ^3 ]7 J) O/ k' O% ^# T, ?; x% n    X1(i)=X1(i-1)+X0(i);
0 ?* a2 N- e8 Dend
2 G* l) `( I* l2 O$ {for i=1:n-1
    B(i,1)=-0.5*(X1(i)+X1(i+1));
6 s( m6 t% b, g# a0 v+ F6 }5 G    B(i,2)=1;
. b4 y) ^5 i$ b; W- E' c- H  Q    Y(i)=X0(i+1);
end
alpha=(B'*B)^(-1)*B'*Y';
: R9 ^3 J! x6 O; u) z  G% }a=alpha(1,1);
, \: B' ]0 N7 l; C) [+ g: Kb=alpha(2,1);
; P. {% N* i! I) d" g3 td=b/a;
c=X1(1)-d;
; y6 a( i8 V# e* ^" g6 UX2(1)=X0(1);
7 k5 N( i% H1 f. @X(1)=X0(1);
for i=1:n-1
    X2(i+1)=c*exp(-a*i)+d;
8 p6 K, H& H  E* a, |2 A1 {    X(i+1)=X2(i+1)-X2(i);
end
for i=1:n
) {& R1 i0 O% n    error(i)=X(i)-X0(i);
$ v5 z& r% s! X1 Z. b    error1(i)=abs(error(i));
/ j# E4 O' F8 S" N    error2(i)=error1(i)/X0(i);
end
4 ^3 R% g+ m0 z( |! r2 @* H) \g=std(error1)/std(X0)
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
/ ]  p: J) A4 q' m0 JX(19)=[];
  f) {2 k- T* Z4 P# o- u
1 {4 G  D# ]  M: z
plot(1992:2009,X,'b*-')
hold on
- k! v* C: U8 _* @1 X8 W

X0=[
    357.35
653.91
1123.29
1601.79
1952.05
9 J# |1 I- |) a- C+ x1977.59
1964.89
1 B& ^. o) d( ]5 T* C9 }8 i1856.72
. n4 R! N5 K4 D6 G% h9 B$ D9 C! H1869.67
3 e: S: E& P* S8 Z7 b2 E1984.31
; ~7 X$ z- r' M' \: ^7 q/ Y$ u6 R. f$ s2158.41
2828.87
, w3 r8 `% i. ~. S+ ^. a& L3084.66
3542.55
3925.09
8 v' |6 u" ^, e  a4458.61
4 d( D7 s& @) _: F" v/ w! ^9 v" O4829.45
5273.3
];
: W( b  x. x- `
plot(1992:2009,X0,'g+-')
for k=2:18
    l(k)=(X0(k)-X0(k-1))/(X(k)-X(k-1));
* r3 P5 U- r  J: dend
l
( [8 I/ e" T5 T5 q* o! Oylabel('投资/亿元');
; A; m3 S# a+ X6 [xlabel('上海       时间/年');
: _% V7 \3 d  M1 Elegend('预测值',' 实际值 ' ,'Location','NorthWest')
5 B8 z! g; J+ U: l2 ]4 b%legend('预测值')

" p6 d( i/ t' h+ x1 ~* ]

神秘了一场雨