灰色预测

zhouchang

format long;
: x, O. g5 _$ V' S  S+ h! N: ?X0=[357.35
653.91
1123.29
1601.79
1952.05
1977.59
1964.89
1856.72
, ]9 P9 n3 t2 ]8 ?* q$ t1869.67
8 `7 Q# C0 s/ {9 ?/ F2 L1984.31
2158.41
2828.87
];
n=length(X0);
/ B. y( T6 Y* i! U$ W: y; vX1=[];
. Z  V9 a/ [3 `7 m) [  A' SX1(1)=X0(1);
+ B- V) N' h+ C' l4 {0 ^9 T' [  Ufor i=2:n
    X1(i)=X1(i-1)+X0(i);
: L' k% i9 M+ \2 aend
9 Z8 C3 l$ k3 Z8 L  wfor i=1:n-1
2 Z, B& B( T/ [1 l' b9 d    B(i,1)=-0.5*(X1(i)+X1(i+1));
2 ^2 W* y; R" O; @( G! r    B(i,2)=1;
    Y(i)=X0(i+1);
end
alpha=(B'*B)^(-1)*B'*Y';
a=alpha(1,1);
8 q9 S: B' @$ g1 rb=alpha(2,1);
d=b/a;
/ E8 m  u$ s' {0 Uc=X1(1)-d;
5 D1 i* n. c5 QX2(1)=X0(1);
X(1)=X0(1);
6 F4 c) D6 r# C/ Kfor i=1:n-1
3 T( ]7 G" H3 p, p: Z" k    X2(i+1)=c*exp(-a*i)+d;
    X(i+1)=X2(i+1)-X2(i);
; s% s& [9 ^, x1 l- E; e" b# yend
) @/ H% t4 g6 Q( Z' G: Bfor i=1:n
& m" X' ~6 `8 T' ?' e- Z    error(i)=X(i)-X0(i);
% L) _! v; ?$ [  m    error1(i)=abs(error(i));
& ~1 X  z$ |5 K    error2(i)=error1(i)/X0(i);
end
g=std(error1)/std(X0)
e=c-c*exp(a)
  \# y2 v4 u9 s5 Nfor i=1:18
    X2(i+1)=c*exp(-a*i)+d;
    X(i+1)=X2(i+1)-X2(i);
   
end
X(19)=[];

( H1 s& I! \+ |3 U* T6 k" e8 X
plot(1992:2009,X,'b*-')
hold on
9 m3 X, o7 w- E' s4 `

X0=[
, Q( ]: Q9 d. n: F/ z    357.35
$ q8 @8 z4 I/ s4 i653.91
/ i$ a5 Z0 b9 Y: X" ^2 @5 J! Q1123.29
1601.79
1952.05
3 [/ F) S' ~3 [1 @' d. v1977.59
1964.89
1856.72
1869.67
1984.31
2158.41
0 L' f+ y+ h, U) g/ d, d: ?7 {7 |2828.87
3084.66
3542.55
) V  {7 n+ R6 n; q3925.09
4458.61
7 {! Z5 M( m- D/ A) ^4829.45
8 O7 I* U. c6 b- L5273.3
];

plot(1992:2009,X0,'g+-')
for k=2:18
    l(k)=(X0(k)-X0(k-1))/(X(k)-X(k-1));
5 h% `3 z/ {+ r+ V9 A% Oend
8 f' V2 I( u6 e2 ]2 Y7 l% g% ul
ylabel('投资/亿元');
xlabel('上海       时间/年');
legend('预测值',' 实际值 ' ,'Location','NorthWest')
%legend('预测值')

神秘了一场雨