灰色预测

zhouchang

format long;
, ~. R, |1 I. ]( h& v, ^% Z1 v7 QX0=[357.35
. ?9 L  R' r  A5 J653.91
# S1 d# v2 u2 S3 E+ Q+ G/ T) a* {1 P1123.29
; K) w+ T; r/ t, i" f+ Y1601.79
0 U! n& C6 i; q1952.05
1977.59
2 l0 F4 w* B) g, L. ~  D. @1964.89
1856.72
9 R# `. }6 P4 A( a1869.67
1984.31
2158.41
7 w/ z4 r9 Z7 Q2 L* C* r2828.87
];
4 U5 ~- g/ ^9 U6 n! d- X4 a+ ^1 in=length(X0);
- r& O1 ~- E$ y9 a" [5 }9 ^X1=[];
* ]- a( E1 H! \X1(1)=X0(1);
# A, f6 k5 ?. V9 t$ z; pfor i=2:n
    X1(i)=X1(i-1)+X0(i);
3 \; N, q, c- {0 d: V! Rend
for i=1:n-1
" \% N6 K$ |( y" K9 }    B(i,1)=-0.5*(X1(i)+X1(i+1));
    B(i,2)=1;
    Y(i)=X0(i+1);
end
  P& d% w  a& V6 B3 B  o6 Talpha=(B'*B)^(-1)*B'*Y';
! o2 [! X8 h$ f; @9 ~' ?! T& i: Aa=alpha(1,1);
b=alpha(2,1);
1 C; n, K) m. {& P, g" v/ \d=b/a;
c=X1(1)-d;
/ V% H- \8 \+ IX2(1)=X0(1);
X(1)=X0(1);
/ q+ b, `% h7 }; T* R, ufor i=1:n-1
    X2(i+1)=c*exp(-a*i)+d;
    X(i+1)=X2(i+1)-X2(i);
end
for i=1:n
    error(i)=X(i)-X0(i);
    error1(i)=abs(error(i));
' ?6 ?2 J9 v  V, w# c5 C    error2(i)=error1(i)/X0(i);
end
g=std(error1)/std(X0)
1 p+ u& b: P3 \7 t, ne=c-c*exp(a)
for i=1:18
    X2(i+1)=c*exp(-a*i)+d;
    X(i+1)=X2(i+1)-X2(i);
1 G+ b6 N+ ~) M) W) A   
end
X(19)=[];

$ i7 _/ U& Z# S& C# J
5 ?( F: Q& x1 C) S8 r9 f; rplot(1992:2009,X,'b*-')
hold on
& r7 X8 w5 I0 o5 u( K

0 f; D- z7 [1 x6 t( MX0=[
    357.35
653.91
1123.29
: g4 c  S, {! H1601.79
3 |5 ]3 ~" d# p. K; g& H) Z- r$ v1952.05
1977.59
7 F3 `# r, h$ A; [+ E3 D1964.89
1856.72
9 T0 c8 F. K. p3 d1869.67
1984.31
1 P9 G5 @9 l/ n) c2158.41
) S" b& N- q& \0 N5 N5 I2828.87
) m0 P+ M. k+ V- G" @3084.66
8 _3 R7 a: [; e/ O) K9 ?: F8 t3542.55
3925.09
4458.61
7 B2 C" n/ p0 Q' j* T4829.45
1 c7 u: h8 K, Y  {5273.3
+ X. T& f5 S+ J7 B: N  B# {];

7 F: x: M2 i" l  y4 q+ q( T9 F7 ?plot(1992:2009,X0,'g+-')
for k=2:18
    l(k)=(X0(k)-X0(k-1))/(X(k)-X(k-1));
! E8 J8 K  J0 Q3 lend
7 A7 h1 B! o( ]# L8 zl
! k3 ]. f9 |4 l0 T5 i" E. eylabel('投资/亿元');
xlabel('上海       时间/年');
legend('预测值',' 实际值 ' ,'Location','NorthWest')
%legend('预测值')
  n6 P- l5 |6 i# y: s: a# J  J5 T0 m  v
6 B" H" L! o, O! |2 e2 f9 N

神秘了一场雨