灰色预测

zhouchang

format long;
X0=[357.35
653.91
1123.29
+ A0 [  ?- S: u9 m+ T1601.79
0 t+ N. p8 J  `7 U* I7 G: q1952.05
1977.59
1964.89
  @% V- O' V/ ~" ?; Y1856.72
8 m* [+ |7 V) k$ {5 a8 V1869.67
4 r% Z" p+ F# U& i1984.31
' |3 h  r2 C/ k2158.41
2828.87
' ^6 b) M+ [" |: O' }. m1 B ];
  {7 G6 D$ h9 @0 A( p, qn=length(X0);
X1=[];
( C; D. y2 M. M- c. @X1(1)=X0(1);
for i=2:n
0 @: `# g, [6 N6 V0 `: G$ i    X1(i)=X1(i-1)+X0(i);
8 u( A# O# \# ^2 q: R0 a6 @3 send
* q& X& N! ^0 n6 Vfor i=1:n-1
5 J, `4 B2 Y, P: w/ H) V    B(i,1)=-0.5*(X1(i)+X1(i+1));
! U/ |5 Y1 n* p; K7 R  }! f    B(i,2)=1;
    Y(i)=X0(i+1);
end
alpha=(B'*B)^(-1)*B'*Y';
- t3 s1 ^3 v+ J0 R5 _* Qa=alpha(1,1);
8 d1 t. J' L8 z. nb=alpha(2,1);
4 J" k) K: e: h/ {1 id=b/a;
. Q& q6 A" ~  [) N, [# u0 E! ?c=X1(1)-d;
X2(1)=X0(1);
X(1)=X0(1);
for 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
4 k+ A. b, D+ n, X    error(i)=X(i)-X0(i);
; z8 J- \% R/ Y; L- G; I    error1(i)=abs(error(i));
    error2(i)=error1(i)/X0(i);
end
" \- Z4 L; k. |2 y: h9 wg=std(error1)/std(X0)
2 s  `' }8 M7 J( W; r5 n5 M" }e=c-c*exp(a)
3 ]0 `+ S: o. b0 E, zfor i=1:18
7 Y7 Y/ j/ C4 o2 `$ x/ ^    X2(i+1)=c*exp(-a*i)+d;
    X(i+1)=X2(i+1)-X2(i);
# K# [) s% B4 A2 x   
- A% u+ S% i) H, mend
* |% l7 T# m2 a2 {X(19)=[];
3 M4 x; J. U8 w& W' b/ c

plot(1992:2009,X,'b*-')
1 K$ M2 @6 Q. A3 v/ Phold on
& j3 g4 O' a# }: ~* W0 a) g. z) h

X0=[
/ F" X; T6 O: \7 @9 B4 k    357.35
653.91
1123.29
% Y# A( S& p' y1601.79
1952.05
1977.59
1964.89
& n. I) j3 O2 y! f1856.72
# ?$ \! m' V4 l1869.67
; N" J" v  p* x) b1984.31
2158.41
# x+ o+ {& @  z; R1 F0 j% d( D2828.87
3084.66
3 g$ \1 C2 ~. M3542.55
: E' n" x5 D) J+ M5 Z3925.09
4458.61
4829.45
5273.3
];
) s; q& |3 P2 o  ^
plot(1992:2009,X0,'g+-')
/ v; X. B9 @; }% e( G- t4 Afor k=2:18
6 {* t3 R1 c7 A3 l4 O3 N- e    l(k)=(X0(k)-X0(k-1))/(X(k)-X(k-1));
2 z% V  H" F) c7 V# r* lend
1 Y' w: @# d& W- u! r, Tl
ylabel('投资/亿元');
xlabel('上海       时间/年');
legend('预测值',' 实际值 ' ,'Location','NorthWest')
%legend('预测值')
% \, ]6 O1 B+ X9 c+ g  s" T

神秘了一场雨