急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
" `+ V9 r' L* V! _function f=gm(x0,m)              %定义为函数gm(x)
& F7 P8 G- c: w- Bn=length(x0);                     
1 b' o8 h; p; e0 }; y2 ^x1=zeros(1,n);                    
x1(1)=x0(1);                  
: H3 \6 e9 _3 |7 }* U& n% R/ lfor i=2:n                         %计算累加序列x1
9 @! `6 M6 G6 x    x1(i)=x1(i-1)+x0(i);         
end
* j. q# y9 f+ X% |4 a9 Vi=2:n;                            %对原始数列平行移位并赋值给y
- y( e% R. \, E# {y(i-1)=x0(i);                     
' k5 [2 n( \9 Z# C9 D* e2 g& Ty=y';                                 
i=1:n-1;                             
, C+ U: f1 C1 e# q+ s" ^0 L; Pc(i)=-0.5*(x1(i)+x1(i+1));                 c(i)=-0.5*(x1(i)+x1(i+1))
7 L. e# U. C( M- yB=[c' ones(n-1,1)];                        
/ n& f) |9 f4 l' {7 f( ?au=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
$ C1 k& D$ I; g9 l! h3 {+ O' X                                            
i=1:m;                             %计算预测累加数列的值
ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
yc(1)=ago(1);                           
i=1:m-1;                             %还原数列的值
7 X0 y4 p8 G& r9 u1 I' X7 zyc(i+1)=ago(i+1)-ago(i);              
i=2:n;                              
/ c+ k! m: D6 x! C  _  Serror(i)=abs(yc(i)-x0(i));                     %计算残差值
: ]+ ~) @* O7 ]9 }1 Syc(1)=ago(1);                       
i=1:m-1;                              %修正还原数列的值
" C! V  m/ ]: _& ]yc(i+1)=ago(i+1)-ago(i);              
c=std(error)/std(x0);                      %计算后验差比 *2
relerror=abs((error)-mean(error)*ones(size(error)));
$ w- r9 R" H5 \* c7 [; t; c[nrow,ncol]=size(relerror);                       
p=0;
for i=2:ncol                                    
    if  relerror(1,i)<0.6745*std(x0)               
) Q# E: {( N( g7 H9 `. H. G        p=p+1;                                 
0 Q& B# J" R0 P8 V/ P# l    end
, z. @5 a  ~: F; O* cend
- y6 T, E5 K3 ~0 N2 Kp=p/(n-1);                     
$ M& c3 y# Y3 E7 z1 b! e% g* E2 g" @w1=min(abs(error));                              
6 T2 H! [' |3 |/ b8 L: S4 R+ Aw2=max(abs(error));
i=1:n;                                     %计算关联度
4 _5 g3 V' J7 n4 p$ h2 bw(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
' @$ ?3 G! U5 D) g. bw=sum(w)/(n-1);
au                                         %输出参数a,u的值
ago;                                        %输出累加数列ago的值
x0;                                         %输出原始序列值
$ f: t3 ~& w; K. p! H6 ^f=yc;                                        %输出预测的值
, S* i+ U9 F% I' V. l8 A( n! q3 verror;                                        %输出残差的值
c;                                           %输出后验差比的值
4 F9 t0 |, j% ip;                                           %输出小误差概率的值
$ O7 A1 Q8 g" a6 P/ `w                                           %输出关联度
end                                   
x0=[31.26
32.09
% r  {8 y4 Q0 z( V33.39
35.45
40.52
9 v2 P* Y/ w# Z3 e& ^+ V43.52
3 f$ }# U6 Q) ~/ z$ I47.86
! m2 h+ ~$ x5 I51.45
  q, A" C3 b& E$ K/ L; w+ _$ a1 a55.6
5 c' }, ~2 L' C. L60.14
8 ]" U+ A/ R, K" j! B5 k( B64.82
68.65
; [8 X* C% e- ?( M2 X73.22
7 B% u. B* K) Z* O3 }80.22
87.69
93.97
" Z2 n2 U! m* \8 q99.16
; T2 J4 j  C/ x8 D1 W2 y* n! I( e103.38
% J/ u/ s( ?' s% G9 {9 {109.46
114.6
$ A3 l; r9 u3 z2 K7 ]119.85
124.92
132.04
$ s( x1 v+ g+ _$ I139.45
150.93
& `( E: r! b* G0 l. C( h4 O/ p165.13
7 k$ ?& z4 W4 y; Y181.93
196.83
! d' @( u3 g! m% f* T' ]' t212.38
4 f5 Y' D2 w# R; O' h' q( E228.07
, I& B5 L* B/ ~2 K  C$ n241.45
8 {& H! K% B: ~: ^% Q3 J251.03
" Q8 ^% @/ }( t- m$ s1 \0 z, }9 }]
我在窗口输入：gm(x0) 为啥会出现
& O- D- H4 W  ~Error in ==> gm at 22
4 `* u6 F+ v% q. X" L) X% c  merror(i)=abs(yc(i)-x0(i));                     %计算残差值
2 r+ r/ v6 h% N3 C# t: i7 P那位高手帮忙 谢谢
) O9 g  R( y3 R& J" X! \

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉