急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
- q' \8 X( P& j/ i2 V8 z1 Hfunction f=gm(x0,m)              %定义为函数gm(x)
n=length(x0);                     
x1=zeros(1,n);                    
x1(1)=x0(1);                  
8 c, W& h) D, ~0 J! y  rfor i=2:n                         %计算累加序列x1
4 b4 W+ V8 T& ]    x1(i)=x1(i-1)+x0(i);         
# C5 T# W0 e; v, I. r# nend
i=2:n;                            %对原始数列平行移位并赋值给y
* j+ G" e+ u* i9 R5 @- T! hy(i-1)=x0(i);                     
y=y';                                 
# g1 H& ^& _/ I0 P7 ei=1:n-1;                             
( j3 O- p5 i# J1 G# M: P* Pc(i)=-0.5*(x1(i)+x1(i+1));                 c(i)=-0.5*(x1(i)+x1(i+1))
B=[c' ones(n-1,1)];                        
+ p; f9 T* _' A( j$ j8 H+ V, Aau=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
1 X4 ^1 P- [7 x7 G: K6 M6 H4 T                                            
+ T3 b, ?2 g. J  c& o5 {i=1:m;                             %计算预测累加数列的值
2 N( Y* |7 H" rago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
yc(1)=ago(1);                           
9 u8 @$ P- u1 c. |! Ri=1:m-1;                             %还原数列的值
5 `) h: a/ J' a& h- ^) L( ^7 pyc(i+1)=ago(i+1)-ago(i);              
i=2:n;                              
, Q/ S  [6 y. Derror(i)=abs(yc(i)-x0(i));                     %计算残差值
yc(1)=ago(1);                       
0 f" \1 p3 A0 ^) b/ w7 r- t! Ui=1:m-1;                              %修正还原数列的值
4 t' |9 f. c6 }% E, g7 Uyc(i+1)=ago(i+1)-ago(i);              
c=std(error)/std(x0);                      %计算后验差比 *2
. m4 V3 ?5 Q% B% n& urelerror=abs((error)-mean(error)*ones(size(error)));
. E. E+ x- @/ q- ~[nrow,ncol]=size(relerror);                       
p=0;
9 R4 F" ?0 y. o( g6 hfor i=2:ncol                                    
    if  relerror(1,i)<0.6745*std(x0)               
$ J8 {  C' s( n! H( C# {+ N$ I( H        p=p+1;                                 
" t+ E9 O! b% b% \& a2 k    end
# I8 r# m( p& B1 a, Y4 h4 t3 qend
  W3 C3 }' F4 _3 ]! s% b4 ^  Yp=p/(n-1);                     
w1=min(abs(error));                              
w2=max(abs(error));
# q4 C2 P* e" d7 p( P' vi=1:n;                                     %计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
w=sum(w)/(n-1);
au                                         %输出参数a,u的值
8 A0 d: E4 b  y* aago;                                        %输出累加数列ago的值
* T5 u. x- R. ?* i) A7 \  Fx0;                                         %输出原始序列值
f=yc;                                        %输出预测的值
error;                                        %输出残差的值
c;                                           %输出后验差比的值
0 C9 ~  v" W& {2 _7 Mp;                                           %输出小误差概率的值
w                                           %输出关联度
, _4 i8 Q5 c( h! V$ Nend                                   
' }2 v% Z+ E( T$ n! x7 hx0=[31.26
3 E. r2 [0 e4 x5 v$ Y$ G! N: \( p32.09
6 l* M: H" n- S( O33.39
35.45
$ Z# z# Y. T% y: _0 x0 b" ?0 T* [# h/ q40.52
) Y% o, h$ d' x6 L43.52
8 i5 o1 I" d8 n6 J: y0 |" e47.86
* T! y* c2 j8 e% z51.45
55.6
5 Q& H5 ]4 I; c60.14
64.82
68.65
6 t6 `( j1 t; b' J73.22
80.22
87.69
93.97
99.16
& m+ {3 }$ }( y$ O103.38
109.46
# L$ b, ~9 @, T114.6
' \6 A/ L3 f9 |8 p/ u4 o& U  |119.85
2 G7 Q7 K; b0 q4 q5 b9 Z124.92
) w5 t8 W- e% f. Y132.04
' d) m; f6 L5 S7 P$ I4 f- ]% T139.45
150.93
165.13
181.93
9 p8 ]3 U  g% |0 a196.83
; s; p3 w) ~1 z% E; }' V212.38
$ R8 G; g" r% P3 V- s3 q( w" d228.07
' H% _2 r$ ?) ?241.45
251.03
]
1 D# T. \0 h% W% U& v/ [我在窗口输入：gm(x0) 为啥会出现
Error in ==> gm at 22
+ J3 V! s; K8 Y" |  Rerror(i)=abs(yc(i)-x0(i));                     %计算残差值
那位高手帮忙 谢谢
, F! t% A& ?8 C+ T$ O: i4 E( V: l

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉