急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
function f=gm(x0,m)              %定义为函数gm(x)
7 {! L, ]3 V: H# b- Q/ y  o+ jn=length(x0);                     
, m- c" K: M/ S2 N& @% V, ox1=zeros(1,n);                    
x1(1)=x0(1);                  
for i=2:n                         %计算累加序列x1
    x1(i)=x1(i-1)+x0(i);         
8 l0 ~  [9 |: F& xend
i=2:n;                            %对原始数列平行移位并赋值给y
2 _- l; g, I7 i% y9 I# qy(i-1)=x0(i);                     
y=y';                                 
i=1:n-1;                             
c(i)=-0.5*(x1(i)+x1(i+1));                 c(i)=-0.5*(x1(i)+x1(i+1))
" \7 e  w. ~! q3 W8 Y, D! M# Z- F& uB=[c' ones(n-1,1)];                        
au=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
                                            
i=1:m;                             %计算预测累加数列的值
- I  F+ j- k" R# r- K4 wago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
yc(1)=ago(1);                           
1 S  L' _: s" R/ ^; Ni=1:m-1;                             %还原数列的值
7 R- V" n0 k& N& r. {yc(i+1)=ago(i+1)-ago(i);              
i=2:n;                              
error(i)=abs(yc(i)-x0(i));                     %计算残差值
: Z: `' l6 d9 u# Eyc(1)=ago(1);                       
i=1:m-1;                              %修正还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
0 ]: ^) h9 _- X6 uc=std(error)/std(x0);                      %计算后验差比 *2
relerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);                       
p=0;
for i=2:ncol                                    
    if  relerror(1,i)<0.6745*std(x0)               
/ d. M: z) i2 y1 r8 Z+ H. m        p=p+1;                                 
    end
end
p=p/(n-1);                     
' D% {) z* u  q+ J! pw1=min(abs(error));                              
% d) v, b3 v# E( jw2=max(abs(error));
i=1:n;                                     %计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
w=sum(w)/(n-1);
au                                         %输出参数a,u的值
7 f- R& f( C) y/ ?+ S$ ~- _, jago;                                        %输出累加数列ago的值
x0;                                         %输出原始序列值
f=yc;                                        %输出预测的值
error;                                        %输出残差的值
7 n  R" F+ g. ^$ ?c;                                           %输出后验差比的值
! z& _- A* m! p! w8 Op;                                           %输出小误差概率的值
w                                           %输出关联度
6 V) f( b5 j2 z/ F: O4 Qend                                   
x0=[31.26
5 R8 F# y8 H3 E! f2 r" G2 T# m32.09
33.39
35.45
40.52
8 @3 \! R: J* K. ^2 L7 _6 T43.52
47.86
51.45
55.6
: J: Q0 |# b* U3 x60.14
64.82
68.65
73.22
80.22
  `" v- g. e. r+ x. N87.69
; Z* c8 L% H9 p6 N- |93.97
8 L1 R& o8 u6 d4 W  I1 Y) f3 M99.16
1 U+ K) e4 y' o$ I$ F+ a7 X5 J103.38
5 j7 @, F- z/ f5 b8 f1 ]109.46
114.6
119.85
124.92
4 U1 H* ?& }) e, x132.04
139.45
5 W% N( X  e) A8 {. [! T7 ~1 j+ F150.93
0 o1 j1 T+ F+ Z4 s165.13
# l2 K* d! R$ L/ S& O( t181.93
196.83
5 N3 Z  B" D% e6 j5 d5 x6 D% c212.38
! e* m* W9 U6 f/ p9 H$ Z4 G228.07
+ W1 x* t  {5 l+ n$ U* k241.45
251.03
$ x, O3 E/ T! ^! ^5 l: S]
我在窗口输入：gm(x0) 为啥会出现
Error in ==> gm at 22
error(i)=abs(yc(i)-x0(i));                     %计算残差值
那位高手帮忙 谢谢

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉