灰色预测

zhouchang

format long;
1 I6 V' o" m6 |6 }; Z; e% i/ NX0=[357.35
653.91
4 ~+ T! Z- A0 |7 e1123.29
1601.79
! I6 H! K: g' Q7 ~2 v7 A- K1952.05
- N/ D+ ?4 b; r: @% ~& E. t- t: J$ Q1977.59
1964.89
1856.72
1869.67
1 V* t2 j, v( G1984.31
' U  B# [: n4 z& ]4 Z2158.41
: o- q7 u5 [: G* }( t2828.87
];
n=length(X0);
) ^. E7 I. ^( @  aX1=[];
' X) l- T' |8 w- oX1(1)=X0(1);
0 f, ]% t3 @  jfor i=2:n
& v. ~) E3 h; P6 y/ c7 z  s2 }    X1(i)=X1(i-1)+X0(i);
) o. _% |6 L% Hend
1 C; ], i$ I2 V  I! M9 Nfor i=1:n-1
    B(i,1)=-0.5*(X1(i)+X1(i+1));
# K% I1 `! `. }8 x* P" p/ J    B(i,2)=1;
    Y(i)=X0(i+1);
end
alpha=(B'*B)^(-1)*B'*Y';
a=alpha(1,1);
( t7 S" X, f8 y2 {6 {2 Ab=alpha(2,1);
d=b/a;
3 J+ B$ y) I" Oc=X1(1)-d;
X2(1)=X0(1);
X(1)=X0(1);
% E7 B7 Z# t8 _; i* c; d& |1 F& afor i=1:n-1
) g' O* T* t! @( ^" d2 J    X2(i+1)=c*exp(-a*i)+d;
0 E/ V' ]- I+ `9 e    X(i+1)=X2(i+1)-X2(i);
end
for i=1:n
    error(i)=X(i)-X0(i);
    error1(i)=abs(error(i));
    error2(i)=error1(i)/X0(i);
+ f  q6 E; K( H4 U% S4 Tend
g=std(error1)/std(X0)
e=c-c*exp(a)
for i=1:18
    X2(i+1)=c*exp(-a*i)+d;
    X(i+1)=X2(i+1)-X2(i);
5 \+ g6 ^. M$ p+ X& l   
end
5 @' H7 V  D$ t4 F  XX(19)=[];
5 A1 c) }1 B+ e
; _: ?8 {4 v. K1 _2 [' B4 m7 k
8 r) ~2 g/ I3 b4 t+ d8 Kplot(1992:2009,X,'b*-')
0 P$ X$ \  I/ L7 K( yhold on

- [7 T% X% S' Z' |  B
2 V+ c# s8 B3 `" d" T1 v6 z9 X. KX0=[
    357.35
* Q/ K# ?) a* t% U653.91
1123.29
3 R' L4 u% V) G: B6 r3 p1601.79
1952.05
; g. o: O: u& N" R# q2 j$ b7 }1977.59
1964.89
1856.72
1869.67
: j& v0 [7 S3 ?- P" k2 G( V1984.31
- b0 Y. @$ c/ L6 i* |/ Y$ I2158.41
2828.87
& \/ ?7 Q& O9 x3 }3084.66
3542.55
3925.09
1 k  ~) X( `4 I" A8 h4458.61
$ f. N( ]- Q6 y6 c( a  e7 u# u& [4829.45
8 u9 _( [% l9 _7 u' [" G8 L5273.3
];

plot(1992:2009,X0,'g+-')
for k=2:18
    l(k)=(X0(k)-X0(k-1))/(X(k)-X(k-1));
end
l
ylabel('投资/亿元');
xlabel('上海       时间/年');
legend('预测值',' 实际值 ' ,'Location','NorthWest')
%legend('预测值')

神秘了一场雨