急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
- s& I7 I7 G4 Q; r, ?2 |8 i3 _. ^function f=gm(x0,m)              %定义为函数gm(x)
n=length(x0);                     
x1=zeros(1,n);                    
x1(1)=x0(1);                  
& p5 D( {/ W9 V% vfor i=2:n                         %计算累加序列x1
) R% l4 q3 G8 J" o0 c2 g6 H0 z0 w6 C    x1(i)=x1(i-1)+x0(i);         
end
i=2:n;                            %对原始数列平行移位并赋值给y
y(i-1)=x0(i);                     
- Z2 ?2 w7 `  F( E: G0 z. Z# H4 f2 ^' Qy=y';                                 
, N' t  E! x1 J+ R+ T! yi=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)];                        
- r' J# T  J/ oau=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
7 w2 a! K; h* J0 V                                            
i=1:m;                             %计算预测累加数列的值
; C  d7 l2 A; Jago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
yc(1)=ago(1);                           
i=1:m-1;                             %还原数列的值
, w9 `3 Z/ R/ I5 Y! d3 Myc(i+1)=ago(i+1)-ago(i);              
" L" k3 w% q& W$ \( I( Vi=2:n;                              
error(i)=abs(yc(i)-x0(i));                     %计算残差值
yc(1)=ago(1);                       
% |. D- L. t6 J( P- V8 q$ O! wi=1:m-1;                              %修正还原数列的值
1 \: v5 I3 j$ M) H" k! xyc(i+1)=ago(i+1)-ago(i);              
c=std(error)/std(x0);                      %计算后验差比 *2
3 _  N" u* X! U( R6 i6 n8 Rrelerror=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)               
5 e+ i2 |6 H1 |+ X% S4 I6 Z3 k5 d        p=p+1;                                 
6 I3 L3 A: Y/ M  Z    end
4 _) P1 R+ z+ y% d( N8 _end
p=p/(n-1);                     
w1=min(abs(error));                              
0 G1 ]& ~, S2 S4 h  r& @8 jw2=max(abs(error));
i=1:n;                                     %计算关联度
- i5 `/ s8 }: [( ]& Yw(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
- H. G* L' T. |4 iw=sum(w)/(n-1);
au                                         %输出参数a,u的值
ago;                                        %输出累加数列ago的值
x0;                                         %输出原始序列值
( V0 \3 U0 p8 M5 a( Df=yc;                                        %输出预测的值
7 K6 X0 H1 R3 k9 }8 R5 eerror;                                        %输出残差的值
c;                                           %输出后验差比的值
0 q+ K  j5 L2 \p;                                           %输出小误差概率的值
0 [9 U9 ]% [* x5 ow                                           %输出关联度
2 e" t: z8 l! E( Y5 ?& t* o: vend                                   
x0=[31.26
6 v9 ]& ]2 @+ o$ M+ m0 q  g32.09
33.39
35.45
, s3 x" n( `3 f/ z+ @0 l/ y40.52
43.52
47.86
( W$ k7 r1 ^) e7 H  @  r51.45
! c& }' }& Y; w+ i: q55.6
60.14
64.82
6 F. X2 z: y! X6 @, l* W# L0 m' n68.65
2 _2 P7 m* e/ B73.22
80.22
7 H! M% R4 H  C/ f( ]- m87.69
93.97
$ r# ?" z- f$ t& ]" D9 ^99.16
3 E! r. G- H+ u+ J$ l9 M( P103.38
109.46
; y& P, _1 B& n" z! F' }& y114.6
, `+ u6 v7 K% s6 d119.85
124.92
) K) I/ N6 L& M132.04
. q" G% U3 x$ v) j( `& f139.45
- {. T0 q4 w' t( w6 P" v3 M150.93
2 H& G$ [0 i+ x/ s  u165.13
181.93
) }  T7 D; t4 q0 _/ G6 y196.83
212.38
228.07
, V6 O4 t! l7 U% o3 n& e241.45
251.03
]
我在窗口输入：gm(x0) 为啥会出现
Error in ==> gm at 22
0 V* s. Y9 ]/ a) ~6 S% g' xerror(i)=abs(yc(i)-x0(i));                     %计算残差值
那位高手帮忙 谢谢

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉