灰色预测

zhouchang

format long;
/ ]) p5 M9 L4 lX0=[357.35
653.91
! T9 F* z6 \% A( X- L- @1123.29
) K0 A6 ]3 J2 \1601.79
. p) ?6 |  p  [  Y4 K1952.05
1977.59
1964.89
1856.72
1869.67
1984.31
' e, ~7 U% r! d; p2158.41
, B& m) d, k$ g2828.87
];
n=length(X0);
X1=[];
X1(1)=X0(1);
$ |( }6 T$ \% Jfor i=2:n
    X1(i)=X1(i-1)+X0(i);
end
/ R' z* O3 S, q: M. mfor i=1:n-1
    B(i,1)=-0.5*(X1(i)+X1(i+1));
0 ~) o9 {, R$ o    B(i,2)=1;
& n7 }" ^/ }% l) R- R% c    Y(i)=X0(i+1);
end
5 c4 S# U9 h) G1 q1 k% }/ Talpha=(B'*B)^(-1)*B'*Y';
a=alpha(1,1);
) a! g* @; y  N9 R: c  R* @b=alpha(2,1);
d=b/a;
c=X1(1)-d;
X2(1)=X0(1);
* j5 p: X$ D4 `X(1)=X0(1);
for i=1:n-1
+ a5 M7 r5 d$ \9 h& L# Y    X2(i+1)=c*exp(-a*i)+d;
    X(i+1)=X2(i+1)-X2(i);
end
  L  f8 m  g& W  R9 Zfor i=1:n
4 [$ e3 M2 R! g9 m1 f    error(i)=X(i)-X0(i);
$ w& e, K/ {' V" m3 l* A- p& c    error1(i)=abs(error(i));
    error2(i)=error1(i)/X0(i);
$ W* o: S! H' _% P9 U$ X6 ^- yend
# X8 W' G% b- G# t0 W) l* e: ~g=std(error1)/std(X0)
e=c-c*exp(a)
  c* R  X* T# C1 J3 Qfor i=1:18
( L  P/ Q- Y9 {% z    X2(i+1)=c*exp(-a*i)+d;
; F& {- O; i4 w* T) w    X(i+1)=X2(i+1)-X2(i);
" i% w9 ?5 Y* V& Y3 m8 a   
9 ~) ?9 H% }3 l/ Uend
; K2 R1 @$ v" p. e8 I( ^X(19)=[];
- E* Q8 t$ a$ O5 c/ l. w$ @
$ d6 p9 ?& g- y: _9 o9 S
0 u# u" T$ A4 b, |6 c) t, splot(1992:2009,X,'b*-')
hold on
9 y  ]8 V9 a1 s0 f/ U0 C

9 e$ {1 c6 c/ ]( H: m2 UX0=[
, M( G2 g1 {  d- s0 ~. m5 A0 f4 ~    357.35
- o5 ?6 Q- k" l0 P+ Y653.91
. i+ r# w+ T2 f8 f  [) L1 X1123.29
1601.79
1952.05
2 }4 I  L4 }# N8 }) S1 q1977.59
1964.89
# M5 H0 X( Z/ _4 S2 ]; i1856.72
5 C7 J; j9 z9 S9 k4 [1869.67
' s5 [* f) r# G  ^6 ~8 Y8 G5 z1984.31
2158.41
2828.87
3084.66
3542.55
$ L8 J+ W9 q, @  z# _  |2 T3925.09
4458.61
* }# C7 v% c2 D: m8 w1 i4829.45
5273.3
% c# h) f" Q8 ^- d0 s];
# D. |2 T9 i, ]  y) U
* r$ M* w9 r0 D  r7 C* F2 O# Xplot(1992:2009,X0,'g+-')
) d! e0 n% d4 u+ T% o1 U9 pfor k=2:18
; V# v7 C4 Z, X5 Y6 ?- \( I4 V    l(k)=(X0(k)-X0(k-1))/(X(k)-X(k-1));
7 I: S9 f" W! _end
l
6 m9 H" E+ c6 k! h- xylabel('投资/亿元');
xlabel('上海       时间/年');
1 {+ M5 Z! e/ b" Alegend('预测值',' 实际值 ' ,'Location','NorthWest')
%legend('预测值')

神秘了一场雨