灰色预测

zhouchang

format long;
% ?% L) Y$ I; F$ EX0=[357.35
( ]8 f5 t, p# Q2 i: B653.91
- R4 c: W* _/ y' [$ b9 l3 Y1123.29
, J3 i, J: N2 f* e* i3 d1601.79
1952.05
1977.59
$ |! }& X  O7 R' ^7 v1964.89
1856.72
: J8 Y  e* O% U9 V: R( Y7 T1869.67
  v. \  |2 g6 K6 i. p, {. J) v; _1984.31
2158.41
2828.87
: W7 L% z6 [: R4 W ];
. B/ x0 r  p+ V+ e. @7 M5 F7 Vn=length(X0);
X1=[];
X1(1)=X0(1);
# N+ U4 r/ n2 J2 ^3 s$ V  w# h6 Kfor i=2:n
+ W$ f& Y3 u. ]; Q* M    X1(i)=X1(i-1)+X0(i);
end
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
2 U+ G- y( ]4 i8 i7 p2 Balpha=(B'*B)^(-1)*B'*Y';
a=alpha(1,1);
b=alpha(2,1);
d=b/a;
+ w% @5 y- E3 ]+ o! Q  [! uc=X1(1)-d;
" U. P1 w9 `9 E* K* V$ zX2(1)=X0(1);
X(1)=X0(1);
1 j4 n9 i; x; n6 a) ^for i=1:n-1
2 ]- L( [3 w8 J# V' V    X2(i+1)=c*exp(-a*i)+d;
+ n' {" b3 I* M& l& G: ?    X(i+1)=X2(i+1)-X2(i);
end
/ U' z( P1 j4 F( ~) E# g: qfor i=1:n
9 T3 q) n3 V2 z) m, f! b/ D    error(i)=X(i)-X0(i);
/ w8 C9 U/ W9 X6 ]( O5 _    error1(i)=abs(error(i));
    error2(i)=error1(i)/X0(i);
) t: ?  @( O5 D% y4 h2 pend
8 D8 x( |$ W. _* P9 @1 d. o6 s5 Rg=std(error1)/std(X0)
e=c-c*exp(a)
for i=1:18
8 ^4 @( [+ ?3 \" ^& @% f6 g) O, r    X2(i+1)=c*exp(-a*i)+d;
    X(i+1)=X2(i+1)-X2(i);
   
end
" Q  c, D! ?; ~: I2 R/ C- OX(19)=[];
* v% U/ S( {: m1 m' [8 F

9 w7 K  N6 a) lplot(1992:2009,X,'b*-')
hold on
1 [9 ^* w' }& @( I- G

X0=[
1 `( e( M( M# J# i9 w/ Z/ {    357.35
653.91
1123.29
1601.79
7 a$ s) c: o5 C1952.05
' {: I" R* I; Z/ [6 s. r1977.59
1964.89
1856.72
1869.67
1984.31
7 n- v+ u7 J7 C/ X* I! p: N& _2158.41
; l5 W$ q8 t% c! e2828.87
- f7 H7 ]6 O, ~% z) {3 D. d3084.66
3542.55
9 b6 N7 K% y) e6 w- [7 T3925.09
6 @8 i1 H7 o+ i, ]7 X  O9 Q4 e2 a4458.61
/ ^3 v- m7 S# P2 d5 j, Q# R4829.45
5273.3
3 m* m! i: J: E( h7 G];
& W$ S  [1 |* Z. M& g' r
plot(1992:2009,X0,'g+-')
for k=2:18
3 j" b. y* l) d9 `8 k1 U    l(k)=(X0(k)-X0(k-1))/(X(k)-X(k-1));
end
l
& M, s: @" t1 r1 yylabel('投资/亿元');
& O5 }$ a/ P3 u8 {: t7 uxlabel('上海       时间/年');
& e! ]# s3 B  J* A/ d3 l6 N  Y6 llegend('预测值',' 实际值 ' ,'Location','NorthWest')
% x' _, w. ~* l1 U%legend('预测值')

; `; ?: ?6 p, n, b" y5 {# E

神秘了一场雨