急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
function f=gm(x0,m)              %定义为函数gm(x)
. {% X2 _" K3 c  L/ T. g2 Dn=length(x0);                     
x1=zeros(1,n);                    
x1(1)=x0(1);                  
$ ]" s* z1 t; e! C; d) \! [for i=2:n                         %计算累加序列x1
    x1(i)=x1(i-1)+x0(i);         
: s# Z+ M& V0 C) Qend
i=2:n;                            %对原始数列平行移位并赋值给y
y(i-1)=x0(i);                     
y=y';                                 
i=1:n-1;                             
$ q* O4 l" U! o3 mc(i)=-0.5*(x1(i)+x1(i+1));                 c(i)=-0.5*(x1(i)+x1(i+1))
2 G  W2 Y& T6 p  Y+ R# o9 p/ |B=[c' ones(n-1,1)];                        
6 g% N% `5 B6 G  L' nau=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
                                            
# I+ h! \$ J+ w) r; k1 ?6 o) k4 G: ci=1:m;                             %计算预测累加数列的值
ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
5 E; V" q+ X% Y4 eyc(1)=ago(1);                           
) A6 w! E  _. i8 r- g  Ci=1:m-1;                             %还原数列的值
' J- L) X1 R: G( k! hyc(i+1)=ago(i+1)-ago(i);              
i=2:n;                              
" ~9 f& Q8 q( \+ y1 _& s- O; E5 y; berror(i)=abs(yc(i)-x0(i));                     %计算残差值
yc(1)=ago(1);                       
i=1:m-1;                              %修正还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
c=std(error)/std(x0);                      %计算后验差比 *2
relerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);                       
, K) X, b  e7 H8 ~; `* j9 Sp=0;
6 i  w1 M- l! p3 f' W& C9 O' ]$ ifor i=2:ncol                                    
) g4 O9 O4 v. E8 Q) u    if  relerror(1,i)<0.6745*std(x0)               
! v3 e2 e6 W% F4 K. k! m        p=p+1;                                 
    end
7 q* J' f1 D8 Mend
& [# T$ V# O9 D8 Wp=p/(n-1);                     
w1=min(abs(error));                              
* `- k7 B" d/ }- S% E$ lw2=max(abs(error));
i=1:n;                                     %计算关联度
( }! I! z- }# |w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
w=sum(w)/(n-1);
  ?& \% Z$ m) M, [2 E- iau                                         %输出参数a,u的值
ago;                                        %输出累加数列ago的值
x0;                                         %输出原始序列值
6 f1 {7 H% N! zf=yc;                                        %输出预测的值
6 J- l+ \/ v$ \5 Gerror;                                        %输出残差的值
: \; L& u- `  M2 N! i5 V+ j3 H, F( cc;                                           %输出后验差比的值
p;                                           %输出小误差概率的值
w                                           %输出关联度
end                                   
2 r: x1 V1 T8 L! \x0=[31.26
32.09
33.39
35.45
40.52
. [9 L; Y$ w+ L4 p" s( F# o43.52
0 e2 y& ~$ I7 n* R$ h47.86
3 u! z/ P' b% O' J4 a51.45
) b7 ?4 f9 Y, F" Y" |  m6 N" i3 s55.6
60.14
64.82
68.65
: J7 ?6 r% u( A) J73.22
80.22
87.69
93.97
99.16
% D6 f6 E2 C; O% n# [0 h/ t' }& @# _103.38
1 P3 W9 W6 W4 F' B+ G109.46
114.6
119.85
124.92
2 t6 f+ l0 v# I' X, ~132.04
) v0 f1 i% Z$ n, d139.45
150.93
: I% v+ P1 h5 E# ]165.13
181.93
$ K9 t/ c. A' j' Z* h  _# h196.83
212.38
228.07
5 x  }, N. q; d) t6 y" e241.45
251.03
]
我在窗口输入：gm(x0) 为啥会出现
" {1 {* j9 H% Z% h, TError in ==> gm at 22
* S/ D0 k* N: p* c  |error(i)=abs(yc(i)-x0(i));                     %计算残差值
那位高手帮忙 谢谢
! u3 [7 x; W9 E2 J$ f3 b

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉