急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
function f=gm(x0,m)              %定义为函数gm(x)
6 j* |+ J, b8 _n=length(x0);                     
2 A# ^7 ]% N! J" B; t6 Ox1=zeros(1,n);                    
x1(1)=x0(1);                  
  l* r* j4 p, s/ Q7 n. l5 Kfor i=2:n                         %计算累加序列x1
    x1(i)=x1(i-1)+x0(i);         
end
i=2:n;                            %对原始数列平行移位并赋值给y
& L7 i* V; y2 T- by(i-1)=x0(i);                     
3 [- v3 t( _' w! q$ T* Wy=y';                                 
i=1:n-1;                             
, \( ~& e, V' @7 e1 Yc(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为弧度
                                            
  a7 G9 c5 x+ w7 N: `7 pi=1:m;                             %计算预测累加数列的值
ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
3 b% G- Z9 j. I# H& vyc(1)=ago(1);                           
2 r- Y$ \) _, T. G: H( @( W' a1 yi=1:m-1;                             %还原数列的值
  |  N# X/ b5 D2 L  t0 h4 ~yc(i+1)=ago(i+1)-ago(i);              
i=2:n;                              
error(i)=abs(yc(i)-x0(i));                     %计算残差值
, ]3 T- _8 u  M  ?5 R: W* jyc(1)=ago(1);                       
i=1:m-1;                              %修正还原数列的值
$ x" H7 n9 m+ z1 T7 R) qyc(i+1)=ago(i+1)-ago(i);              
c=std(error)/std(x0);                      %计算后验差比 *2
2 [! |, G/ T8 {( N+ m1 u2 jrelerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);                       
8 M! r. ], x8 C! [6 cp=0;
- N1 N! g+ s, X* j' P9 B# Afor i=2:ncol                                    
    if  relerror(1,i)<0.6745*std(x0)               
( Y( x9 M9 a+ O  F: K9 _! A        p=p+1;                                 
& B" e+ v' |+ x* W$ S! J- O+ ^5 a    end
8 T# O0 B. ]2 d6 H/ M) Iend
4 O9 J& l' O5 u9 F( D" D8 h/ zp=p/(n-1);                     
w1=min(abs(error));                              
w2=max(abs(error));
0 m! H& c6 R  w4 g' J# ~i=1:n;                                     %计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
w=sum(w)/(n-1);
" Q% }7 y. E$ d% ^au                                         %输出参数a,u的值
4 c1 x: r8 n) E' K; a( {% ?+ r2 Qago;                                        %输出累加数列ago的值
3 q% m( U0 e# nx0;                                         %输出原始序列值
7 b% d- W3 q$ e* T/ Wf=yc;                                        %输出预测的值
3 h+ z9 K7 l, ferror;                                        %输出残差的值
c;                                           %输出后验差比的值
p;                                           %输出小误差概率的值
w                                           %输出关联度
5 ]1 r+ h# O. J% H/ a$ P! nend                                   
x0=[31.26
32.09
% \- s' @. j" B6 k; `$ U7 Z9 `1 D33.39
35.45
/ q! I" N. `& b2 f40.52
43.52
47.86
51.45
55.6
60.14
% {/ n# z* Z$ W4 y- t. x64.82
68.65
73.22
& i. z* r( L! S6 S6 R2 ^* l80.22
87.69
93.97
) I5 r4 T7 c! p99.16
103.38
: n) t; o. D2 _9 ]" E109.46
2 S( ?$ p! k  b/ I( ^3 f4 H114.6
119.85
' G& i$ x' z: d124.92
132.04
  }) }8 G; r0 e0 @( _139.45
150.93
165.13
" Y% e4 `2 j$ _181.93
' `) m/ C) W5 S* r& T( h0 ^2 Y196.83
212.38
228.07
241.45
251.03
, V; \, z( ~5 J( b9 O]
我在窗口输入：gm(x0) 为啥会出现
Error in ==> gm at 22
- @4 |7 U9 V9 Y8 _. Lerror(i)=abs(yc(i)-x0(i));                     %计算残差值
0 h  }6 E/ ], ~" v那位高手帮忙 谢谢
1 I+ A8 |' a7 D9 i0 H" V0 {- d

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉