灰色预测

zhouchang

format long;
/ j- C7 ?7 O1 p8 n  I2 J1 }X0=[357.35
1 b& F" {9 S5 f( R* @' ]653.91
* M6 @: B2 t& C4 o( v1123.29
1601.79
1952.05
& O: K* j8 V+ w1977.59
1964.89
' u! I! P$ S, p3 c; p( o2 G1856.72
* r( n" h# \& Z% W1869.67
1984.31
2158.41
2828.87
6 h& k  J) }: d ];
* c* T0 F3 K7 X) d- C6 ^+ ~. in=length(X0);
0 o7 _& M2 z1 _7 X' mX1=[];
X1(1)=X0(1);
for i=2:n
) Q! T3 n) w: U# y' G) D1 O! ]    X1(i)=X1(i-1)+X0(i);
end
4 p9 S7 [6 a& U0 tfor i=1:n-1
    B(i,1)=-0.5*(X1(i)+X1(i+1));
    B(i,2)=1;
    Y(i)=X0(i+1);
end
3 d- P9 N3 O# r- p/ Ialpha=(B'*B)^(-1)*B'*Y';
5 F* I- ~' p& n. W9 K& c, _4 Za=alpha(1,1);
3 f5 q- J) R/ H4 ~b=alpha(2,1);
) y5 r4 H6 g: ?) _4 |% |) U& Dd=b/a;
7 V7 u4 [/ j# q, C% ec=X1(1)-d;
X2(1)=X0(1);
& |  H/ g8 \/ h* y% k; g( q2 g2 NX(1)=X0(1);
for i=1:n-1
- G9 Y8 O- {# L- e# a7 Z    X2(i+1)=c*exp(-a*i)+d;
5 T# \/ ?4 G5 n/ q# R0 j  |9 ~( s# G    X(i+1)=X2(i+1)-X2(i);
' V2 i0 n( U3 T( D+ Yend
3 y. d7 K% L7 Q, B+ v/ Z5 kfor i=1:n
, s( z1 A% p( `# M' ^    error(i)=X(i)-X0(i);
    error1(i)=abs(error(i));
    error2(i)=error1(i)/X0(i);
end
g=std(error1)/std(X0)
8 i/ P! i/ E. {1 w9 M! z8 T; pe=c-c*exp(a)
for i=1:18
9 o+ W) t9 x9 h, e6 ?    X2(i+1)=c*exp(-a*i)+d;
8 Q# Y* b7 S0 j7 R1 b3 ^0 V' V    X(i+1)=X2(i+1)-X2(i);
   
8 S7 B# f! N* k" c& _) oend
/ X# Z  C, v. H, {% |" nX(19)=[];

3 ^3 E! q0 A9 d, f8 p  n
plot(1992:2009,X,'b*-')
hold on

& l1 W5 C$ p2 k' Z. Z6 |
9 F1 b2 O# `( Y; {8 Y& sX0=[
; w! \8 o1 n( P3 v$ U9 \    357.35
# t; O0 n- p, i: [653.91
7 w3 i9 w6 Z9 ?4 ^1123.29
, ~: G' k+ v/ d9 N- V  Z1601.79
1952.05
; A* N5 X4 Y# K9 b* r( w0 l+ J1977.59
1964.89
6 ]# e/ `' G" n6 |% M2 t1856.72
& q) g! X4 y) k+ K9 W2 i1869.67
1984.31
2158.41
7 T! y9 V0 C1 T/ V* o8 @6 I  L5 I3 o) o2828.87
3084.66
# G5 b: P8 q; `! s. y  H+ k& D3542.55
9 e2 k& W  P5 @3 |9 [  G3925.09
& z# ]% T& T/ N  q4458.61
4829.45
5273.3
: o( _9 R8 `0 O( @];
. x8 D* N" k% |. h: i
2 T3 i, r9 P$ A) ~/ Bplot(1992:2009,X0,'g+-')
8 \! {$ J  z5 K; Ffor k=2:18
    l(k)=(X0(k)-X0(k-1))/(X(k)-X(k-1));
; I1 c. o- j7 G+ Gend
l
ylabel('投资/亿元');
xlabel('上海       时间/年');
legend('预测值',' 实际值 ' ,'Location','NorthWest')
- s9 J1 x1 Q5 i6 @! u4 g0 |%legend('预测值')
2 `0 p+ M) F& w2 ~

神秘了一场雨