灰色预测

zhouchang

format long;
X0=[357.35
653.91
1123.29
1601.79
1952.05
1977.59
% \! m, Q5 ?; w; T# l/ c- j1964.89
7 e0 `" W" M% [$ m1856.72
1869.67
9 R9 l0 r# I7 w1984.31
! u7 B* |" K) w6 U2158.41
3 z/ h% x( D9 X# F2828.87
];
n=length(X0);
X1=[];
X1(1)=X0(1);
for i=2:n
' I+ \8 f( W' C. {" d$ k4 |7 ~    X1(i)=X1(i-1)+X0(i);
end
for i=1:n-1
/ T( J- P6 ?' |0 Z    B(i,1)=-0.5*(X1(i)+X1(i+1));
! L: R' A) B- J7 Y( F    B(i,2)=1;
    Y(i)=X0(i+1);
end
alpha=(B'*B)^(-1)*B'*Y';
a=alpha(1,1);
( Y5 V% w7 @2 b% o, V; Vb=alpha(2,1);
" d* T6 t" B3 {: f  [7 k9 A% Hd=b/a;
c=X1(1)-d;
9 `# C# J& X7 D) r( cX2(1)=X0(1);
) x- c/ _2 v" @( j; [! v3 M' e, L# pX(1)=X0(1);
for i=1:n-1
% e- Z9 V% P9 t- z( O2 L; W# M; g    X2(i+1)=c*exp(-a*i)+d;
    X(i+1)=X2(i+1)-X2(i);
: N0 B& m# k4 g7 M4 }end
- L* o1 o% q- X, L2 K& bfor i=1:n
    error(i)=X(i)-X0(i);
, ?6 A- z4 T5 G7 t    error1(i)=abs(error(i));
    error2(i)=error1(i)/X0(i);
$ |7 D) F" a; lend
" ]1 d: f1 h( `6 g3 Y4 }g=std(error1)/std(X0)
7 w; C7 ]% S- Y% pe=c-c*exp(a)
for i=1:18
1 s/ w* p3 C4 S7 C    X2(i+1)=c*exp(-a*i)+d;
2 q, A0 D$ _, R+ q' t    X(i+1)=X2(i+1)-X2(i);
9 K9 d3 i" t4 B   
end
X(19)=[];

% S9 E' K; T6 i9 v" I
: s3 K& G; P& u6 w9 D# N9 Yplot(1992:2009,X,'b*-')
hold on
1 p& V$ Q7 T% w$ O$ M
6 B* M3 l$ L8 F& W$ C" i
1 i) H- c* n. A$ ^8 f0 NX0=[
( a3 _! R( O0 `& n7 Y    357.35
653.91
1123.29
' ]1 s& y% `2 q6 o. C% m; e! X1601.79
1952.05
1977.59
1964.89
0 h7 \# m' b8 n& y7 h1856.72
1869.67
1984.31
% u+ G6 b  ^$ ?: f2158.41
* i+ A& o# C. W2828.87
; e( F" {  \( H  a# X2 r3084.66
3542.55
3925.09
4458.61
4829.45
5273.3
];
' N+ ?2 s( }- y2 T- U
3 D) ~5 A1 z4 U% z/ \* ^plot(1992:2009,X0,'g+-')
) l8 n0 b3 C& u3 ofor k=2:18
" [1 k. S6 G9 _    l(k)=(X0(k)-X0(k-1))/(X(k)-X(k-1));
& c2 Q; y: e& e/ uend
5 G) }# T+ ~  m4 [0 s$ Al
ylabel('投资/亿元');
( s' H& B8 D5 J& v! `5 Hxlabel('上海       时间/年');
legend('预测值',' 实际值 ' ,'Location','NorthWest')
%legend('预测值')

神秘了一场雨