急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
function f=gm(x0,m)              %定义为函数gm(x)
) @$ I5 t' s9 [7 u8 s% S/ Zn=length(x0);                     
, }  A0 L2 S. z( k" qx1=zeros(1,n);                    
6 H& f) b% a6 }% Nx1(1)=x0(1);                  
; ]7 R. t% N$ C- c6 {for i=2:n                         %计算累加序列x1
    x1(i)=x1(i-1)+x0(i);         
7 V8 P+ c) T) l# E7 R! v0 cend
i=2:n;                            %对原始数列平行移位并赋值给y
1 v4 A1 h9 b0 d8 \y(i-1)=x0(i);                     
: |4 m0 u3 P7 b$ Vy=y';                                 
i=1:n-1;                             
c(i)=-0.5*(x1(i)+x1(i+1));                 c(i)=-0.5*(x1(i)+x1(i+1))
B=[c' ones(n-1,1)];                        
au=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
, ]4 E1 o2 e$ T                                            
i=1:m;                             %计算预测累加数列的值
7 n1 V+ d( z$ Dago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
% q- a, J( M. x1 E. {5 ayc(1)=ago(1);                           
i=1:m-1;                             %还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
/ B* E5 }( U7 t7 mi=2:n;                              
error(i)=abs(yc(i)-x0(i));                     %计算残差值
( B9 m4 {) S$ S% E, S" C- ?yc(1)=ago(1);                       
i=1:m-1;                              %修正还原数列的值
/ t' m4 T' T7 S, y& s1 Y* }yc(i+1)=ago(i+1)-ago(i);              
c=std(error)/std(x0);                      %计算后验差比 *2
8 z; F+ K- R& Yrelerror=abs((error)-mean(error)*ones(size(error)));
' j- M1 v, b" F) H[nrow,ncol]=size(relerror);                       
p=0;
for i=2:ncol                                    
    if  relerror(1,i)<0.6745*std(x0)               
7 T' J. u3 `7 i; N        p=p+1;                                 
    end
end
p=p/(n-1);                     
' s$ R0 i8 b7 K% h+ x" f3 yw1=min(abs(error));                              
w2=max(abs(error));
  P( d9 i+ Q- a; b. Pi=1:n;                                     %计算关联度
8 j; X! @8 u, b  Hw(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
  N" ?# Y+ n0 H, nw=sum(w)/(n-1);
au                                         %输出参数a,u的值
ago;                                        %输出累加数列ago的值
, f; r, m/ O( cx0;                                         %输出原始序列值
- Q5 c  g  C6 l7 I0 Yf=yc;                                        %输出预测的值
error;                                        %输出残差的值
  V' X0 K+ Z2 O3 D, oc;                                           %输出后验差比的值
p;                                           %输出小误差概率的值
w                                           %输出关联度
3 {1 V  l5 X9 cend                                   
, g. E: t# g. I7 C1 ^x0=[31.26
4 u4 ?  _/ X5 q$ n9 ~32.09
33.39
5 n" M- k# |1 k8 d: t35.45
40.52
43.52
47.86
( ^# m. _  H0 i  ^9 P' b51.45
8 v; x- u/ _4 b' ~- o5 A- g' ~55.6
( i% ]4 P- G; @9 U- V! ^  [6 V60.14
  Y* k; @0 o% b4 x3 ~. G7 ?; i64.82
/ ^5 ^. E3 M+ H' r68.65
73.22
7 _- e3 q7 |0 h1 K- B0 O80.22
! k- t5 B- p' e( F0 v87.69
93.97
2 I) ?+ w" f  Z' M( W0 X$ O99.16
  {; X  E6 I$ Z1 Y103.38
6 |" @# [% `# W  N0 D9 B109.46
5 _6 j2 o8 t) V; f8 m114.6
119.85
124.92
132.04
& h8 l8 g; z! P; I5 j139.45
% D2 b, I& m! |' Q& J* J2 y150.93
165.13
181.93
196.83
212.38
228.07
241.45
251.03
' U5 b% t% o' [  k+ `$ ?! s]
我在窗口输入：gm(x0) 为啥会出现
Error in ==> gm at 22
error(i)=abs(yc(i)-x0(i));                     %计算残差值
" M( B6 S. _2 t, i# Z" [3 s$ s那位高手帮忙 谢谢

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉