灰色预测

zhouchang

format long;
! [) C! M2 O0 dX0=[357.35
0 h4 I; ?9 V$ r# b  `653.91
* c3 ?2 T1 G/ f* ]1123.29
9 s; G0 s1 F1 S* U( W% Z1601.79
/ l/ x- v- G6 c1952.05
9 k# n4 l3 L4 [8 [4 o' p1977.59
1964.89
1856.72
1869.67
7 Y: N1 z. v4 O8 B, I1984.31
$ w8 Y7 l# n! ?* F) M  C: R2158.41
7 v9 V6 ?" l! j5 A* J7 Z+ ?2828.87
];
/ T, G: H+ a3 L, q8 {$ ^: gn=length(X0);
X1=[];
: c2 y, ^* S, Y  M+ F+ lX1(1)=X0(1);
# o! c6 K' ]$ C6 }for i=2:n
; ^* D) u1 {7 Y; ]% b. `0 O# v    X1(i)=X1(i-1)+X0(i);
end
for i=1:n-1
    B(i,1)=-0.5*(X1(i)+X1(i+1));
* W" u: Q) n& Y) j; B    B(i,2)=1;
: b. l' g! G1 e& H9 M    Y(i)=X0(i+1);
end
1 W: f3 T8 K$ M$ m5 v8 u1 ualpha=(B'*B)^(-1)*B'*Y';
+ M" h" A' x: m/ ^a=alpha(1,1);
6 b( y; l, w/ U5 A0 Y& Zb=alpha(2,1);
d=b/a;
( B, R7 J) i& Q% ]- V' Zc=X1(1)-d;
X2(1)=X0(1);
X(1)=X0(1);
for i=1:n-1
3 s3 O! G( ]( i- V9 P- F    X2(i+1)=c*exp(-a*i)+d;
7 E+ R" X% {9 J/ I1 r* E9 @$ s, j1 g    X(i+1)=X2(i+1)-X2(i);
end
& m$ ~# {6 b+ D4 e% d4 `for i=1:n
    error(i)=X(i)-X0(i);
: K1 x0 G5 `  b+ I! `& u    error1(i)=abs(error(i));
/ j: `, L% q9 [    error2(i)=error1(i)/X0(i);
end
g=std(error1)/std(X0)
e=c-c*exp(a)
0 U6 n4 |* m. _for i=1:18
    X2(i+1)=c*exp(-a*i)+d;
    X(i+1)=X2(i+1)-X2(i);
' w5 W3 C; X6 \1 ~. W   
end
X(19)=[];

9 `4 ?* l; M- U4 H  {! s
plot(1992:2009,X,'b*-')
# _# t9 C; \5 y! j( i& J7 K3 Ehold on


) @8 I! H5 l- y9 n8 A- q) ?X0=[
    357.35
653.91
1123.29
1601.79
) i5 N% ^2 }# `; L1952.05
1977.59
2 X7 [* x$ f3 h8 h9 s1964.89
4 L; L7 w0 i$ j# {% J) i* c) k1856.72
  H" ]" y" N- L+ I( [1869.67
1984.31
$ e6 _/ K: n( r8 t" U5 P; p2158.41
2828.87
3084.66
3542.55
& l6 [: e- ^, \+ ?; b- }3925.09
5 }3 \$ Q' P5 W' U' y" }3 @4458.61
4829.45
5273.3
( {- W2 R( K+ T% e. J$ P+ J+ T];

plot(1992:2009,X0,'g+-')
for k=2:18
" h# l3 _7 b$ k6 B1 \( K' G- J    l(k)=(X0(k)-X0(k-1))/(X(k)-X(k-1));
end
l
ylabel('投资/亿元');
3 O( l/ T2 k- Z; a3 L# `) nxlabel('上海       时间/年');
legend('预测值',' 实际值 ' ,'Location','NorthWest')
%legend('预测值')

神秘了一场雨