灰色预测

zhouchang

format long;
X0=[357.35
653.91
8 ?3 j! [, y) _4 a& ?6 t1123.29
; m" O: g: \7 J+ c1601.79
) i" e7 z+ V; z7 O1 _$ p1952.05
% S7 L. ~" X- Z# g1977.59
1964.89
1856.72
2 y( \) Z1 j3 @1869.67
& Q& l0 h/ l+ d' U" ^3 C- Q1984.31
7 A; |; ?4 I, Z" z2158.41
2828.87
];
n=length(X0);
X1=[];
  Q1 p* v9 V- W# i' uX1(1)=X0(1);
0 t0 y* U; ^1 Zfor i=2:n
. ?- e6 q+ n4 F, |) s    X1(i)=X1(i-1)+X0(i);
! K, ?0 h, x4 }! [end
) Q4 N7 ?" s; m1 mfor i=1:n-1
# v2 T+ ^6 P$ A# s4 B+ }    B(i,1)=-0.5*(X1(i)+X1(i+1));
( V7 i) Z) `, E# [    B(i,2)=1;
+ u' y$ W" Q* O7 R# Z    Y(i)=X0(i+1);
5 N' ?& Q) g) F7 K0 U" @end
alpha=(B'*B)^(-1)*B'*Y';
a=alpha(1,1);
! `+ X$ ~# A, I5 m; t: d; I3 Qb=alpha(2,1);
6 Y" ~+ s/ w1 j: d. y  @8 K+ m+ qd=b/a;
* w8 q2 ?- G5 V5 O. I+ yc=X1(1)-d;
. B4 n" \! r7 P3 H1 ]1 s$ JX2(1)=X0(1);
X(1)=X0(1);
' @) C! d! i+ y6 M8 s1 v: Xfor i=1:n-1
: ?9 O2 F' o) B* ]# x# f    X2(i+1)=c*exp(-a*i)+d;
$ Q/ m, B: ~! @    X(i+1)=X2(i+1)-X2(i);
end
for i=1:n
    error(i)=X(i)-X0(i);
    error1(i)=abs(error(i));
    error2(i)=error1(i)/X0(i);
end
1 k4 e  S  }; n- U1 Vg=std(error1)/std(X0)
* }" p. p/ P# ]& A$ `4 X0 V0 ne=c-c*exp(a)
for i=1:18
9 M/ Z8 F! K. l* A& [% ^    X2(i+1)=c*exp(-a*i)+d;
    X(i+1)=X2(i+1)-X2(i);
   
end
- F( v8 e; S' ]' H) _: u  q" M% d/ {X(19)=[];

0 u2 y/ [0 S6 L8 E" s8 p
plot(1992:2009,X,'b*-')
hold on

1 O6 [0 R* I* e8 V' K3 S
X0=[
    357.35
5 h! i- q, u7 y& L2 c) r* H( _653.91
1123.29
1601.79
" ?% w0 Z6 A8 z% ~5 a" @1952.05
- F, k( L! p+ r: i) E; [& T0 M1977.59
' g" u3 i8 [4 q; [3 G8 n. X/ {. Z1964.89
& b/ R" j, V* i1856.72
/ x% C3 s4 [2 z' n1869.67
1984.31
2158.41
2828.87
, R4 @0 \, |" e1 i+ d; T0 ?) b3084.66
6 I1 j# X* e: K7 r' K8 `( p3542.55
3925.09
' X# d  L; e8 v$ N( \4458.61
4829.45
% T3 W6 k, q+ v( T. b5273.3
];
5 r: a7 @6 s! h  y' O
plot(1992:2009,X0,'g+-')
8 f" Z9 c8 Z# \' A% v" ^1 bfor k=2:18
% D9 a) Q/ w9 v+ D. Y    l(k)=(X0(k)-X0(k-1))/(X(k)-X(k-1));
3 t: O9 f, ?2 U2 oend
l
ylabel('投资/亿元');
xlabel('上海       时间/年');
4 g+ Z5 q1 F2 F, T$ blegend('预测值',' 实际值 ' ,'Location','NorthWest')
, _7 n( ^# p% P. |%legend('预测值')

神秘了一场雨