灰色预测

zhouchang

format long;
X0=[357.35
653.91
9 @* c' ?3 Y- s+ h; h* b# c' R1123.29
1601.79
8 ?: I7 o$ b9 ]+ `# R. f  O! Y1952.05
1977.59
1964.89
. [7 G$ q& \" ~7 X6 @; Q1856.72
# O- X- p9 U# F# O5 h1869.67
3 z8 g) B" P, d5 I" a1984.31
2158.41
2828.87
];
n=length(X0);
X1=[];
X1(1)=X0(1);
for i=2:n
    X1(i)=X1(i-1)+X0(i);
. L8 O# s& K* ^$ w, d- _& [: wend
* W. t/ j. U: G' rfor i=1:n-1
: J& v9 h. H5 w& P; r    B(i,1)=-0.5*(X1(i)+X1(i+1));
  N# z. \5 i  I+ F* ^) x    B(i,2)=1;
, M' y1 w' P$ x    Y(i)=X0(i+1);
end
alpha=(B'*B)^(-1)*B'*Y';
* `1 H4 c& n2 l  q: i. Ea=alpha(1,1);
b=alpha(2,1);
d=b/a;
6 \! s  W7 x1 [$ ]1 v. L. [c=X1(1)-d;
X2(1)=X0(1);
2 h( y+ ]! q9 L/ ]X(1)=X0(1);
for i=1:n-1
) n5 N( Z4 X5 z, \9 `5 b    X2(i+1)=c*exp(-a*i)+d;
    X(i+1)=X2(i+1)-X2(i);
$ M. m3 k7 `$ F& A- uend
( o) L1 F& J( @: J# }: I% J$ sfor i=1:n
    error(i)=X(i)-X0(i);
* A; T' q* Z: ?+ }# U' b/ A9 N    error1(i)=abs(error(i));
3 h8 U3 b$ f. G    error2(i)=error1(i)/X0(i);
( ?2 e, s/ v' A- Lend
g=std(error1)/std(X0)
6 P- s$ `" c7 n+ q4 i0 B& T" L4 xe=c-c*exp(a)
1 ~3 _8 G, A: }for i=1:18
    X2(i+1)=c*exp(-a*i)+d;
% X+ a0 K: X7 r    X(i+1)=X2(i+1)-X2(i);
2 L1 I9 B& t: x, X6 k: e; A   
end
X(19)=[];
/ }) R2 p. N6 {  r+ k: \
2 A( Z6 e& t+ j8 h
plot(1992:2009,X,'b*-')
hold on
$ D) h+ H' H, }, n" l6 D0 c
7 v, |* s2 j/ P+ y9 l8 {8 y* k8 r! \  L
X0=[
6 f9 A  w5 r9 a) w    357.35
653.91
1123.29
1601.79
# w! _1 y5 D* @- k1 W  S1952.05
% m9 n" R; n! U, y% p1977.59
3 g4 P8 T' }$ \1964.89
, h2 e. ]0 |# S6 _% y3 ~9 f1856.72
1869.67
1984.31
2158.41
5 Q( [3 i  c2 H. O2 }; b& x2828.87
: A6 A0 N/ B' ^+ r9 Z. [3084.66
8 }, _& i! ]8 C7 U4 L! L3542.55
3925.09
4458.61
, ?! F# _3 g" r/ |4829.45
5273.3
* k& ]& Z% F7 j  h];

plot(1992:2009,X0,'g+-')
for k=2:18
  B: Q$ Y5 z# G; M- r4 Y    l(k)=(X0(k)-X0(k-1))/(X(k)-X(k-1));
end
l
ylabel('投资/亿元');
. j8 m: u, w0 A6 n/ z; `, Exlabel('上海       时间/年');
$ K; o; J- Q* e. r7 y4 {0 ^legend('预测值',' 实际值 ' ,'Location','NorthWest')
%legend('预测值')
+ F: _) n; w/ [! J* i+ A" N- |
" a4 U/ q, R, O1 W5 {+ N

神秘了一场雨