灰色预测

zhouchang

format long;
X0=[357.35
4 g+ E" n4 Z3 @7 q  q+ [+ n653.91
1123.29
1601.79
4 @0 ]8 \* O) Y) ^7 H% G1952.05
; k1 U# O" f& g1977.59
1964.89
1856.72
# s: S: S9 S* G0 t. F: _. }) A# r1869.67
1984.31
* _* W4 W; U: R4 C2158.41
' y. |8 O/ T5 o2 `. w3 q, |2828.87
];
2 J& @: J! L# p& a- qn=length(X0);
5 Z7 A  f$ t8 a' \0 q# n+ CX1=[];
# S5 X9 V* z. l. q4 A* `8 d- ^X1(1)=X0(1);
for i=2:n
    X1(i)=X1(i-1)+X0(i);
end
for i=1:n-1
    B(i,1)=-0.5*(X1(i)+X1(i+1));
6 y5 e2 ]$ O: a, M    B(i,2)=1;
) _" H1 Q0 E) e* j    Y(i)=X0(i+1);
end
# h  n2 ?, I9 q6 Q$ B* Nalpha=(B'*B)^(-1)*B'*Y';
5 g! Z; Q3 x6 g" F/ h+ [a=alpha(1,1);
* O% H) K/ U1 _b=alpha(2,1);
3 }. b" |& u- G6 Wd=b/a;
9 z/ }. ]; ?+ s4 t7 z& Ic=X1(1)-d;
3 U7 Q5 Z$ u" Z# [/ e3 UX2(1)=X0(1);
* d& P# d, \7 ~1 YX(1)=X0(1);
for i=1:n-1
; Q# T" I& }6 E- v; y4 z1 E( u; y    X2(i+1)=c*exp(-a*i)+d;
2 ^, R. Y+ I( j4 A/ }9 C    X(i+1)=X2(i+1)-X2(i);
. U5 v9 g/ p% y" ?+ Vend
. e, q) F; D, o& y3 p3 dfor i=1:n
    error(i)=X(i)-X0(i);
    error1(i)=abs(error(i));
    error2(i)=error1(i)/X0(i);
; ?/ [- q1 m3 [) L: Zend
g=std(error1)/std(X0)
% p& F8 G2 B$ ^: ~& ue=c-c*exp(a)
for i=1:18
    X2(i+1)=c*exp(-a*i)+d;
$ K+ I  M; F) a9 d4 _    X(i+1)=X2(i+1)-X2(i);
   
end
X(19)=[];


plot(1992:2009,X,'b*-')
& l" C7 j" Q7 ~; |: vhold on


/ S6 f% O0 s: L$ ]1 WX0=[
    357.35
653.91
' e3 ~" Y% }$ k# U8 Y1123.29
  v8 p9 M  X0 G3 k6 j1601.79
, B  u; ?  X* [1952.05
1977.59
: O6 l1 Z8 o+ ~+ S' @, V& K1964.89
1856.72
1869.67
1984.31
2158.41
2828.87
+ u' X7 @* T, ^( D3084.66
3542.55
3925.09
4458.61
, I% h; w, _% w7 d) [  ^2 k, K+ k! ^4829.45
. `" [. A5 ]7 ~3 \5273.3
" t3 j% w  H1 V3 |];

9 M, N3 ]: z- ~! Uplot(1992:2009,X0,'g+-')
/ c$ m) c  x% Gfor k=2:18
    l(k)=(X0(k)-X0(k-1))/(X(k)-X(k-1));
end
, A/ w2 R( K3 U" V- C* Vl
7 |: C& E/ C- _) Yylabel('投资/亿元');
6 n2 {% m3 U" qxlabel('上海       时间/年');
5 U, M, d% K7 jlegend('预测值',' 实际值 ' ,'Location','NorthWest')
%legend('预测值')

& @1 I: Y' u5 r8 G% w8 r: {

神秘了一场雨