急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
function f=gm(x0,m)              %定义为函数gm(x)
( I$ G7 Y, @3 l( _1 g3 }0 u$ Jn=length(x0);                     
x1=zeros(1,n);                    
. g7 {& f  e) Q) K! f- Px1(1)=x0(1);                  
for i=2:n                         %计算累加序列x1
6 W$ _2 K. x; ^; O1 E( d9 z4 G    x1(i)=x1(i-1)+x0(i);         
end
i=2:n;                            %对原始数列平行移位并赋值给y
" c3 f' m9 P$ l; By(i-1)=x0(i);                     
( C/ S7 O* n0 K7 z" y4 W5 Fy=y';                                 
0 B; Z3 |( b2 N* ?% q5 C0 E' |i=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)];                        
au=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
  H" u  x' ~- X* k  R8 j* x& u                                            
9 W: ^8 t, f2 T6 w+ @& W3 oi=1:m;                             %计算预测累加数列的值
# r" q) Q1 l& Z  U+ Sago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
yc(1)=ago(1);                           
i=1:m-1;                             %还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
3 a! Z' f) K& i! N3 p( ?/ hi=2:n;                              
+ _: b6 j/ L$ d' N! B! Ferror(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
" g6 M1 H! |6 C; }relerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);                       
p=0;
for i=2:ncol                                    
9 n! d! H) o: o6 A; i/ T' [    if  relerror(1,i)<0.6745*std(x0)               
8 T& X9 n- j' H        p=p+1;                                 
' r, f3 W1 D. X( `! h0 \    end
end
p=p/(n-1);                     
; ~' d( |" C! @' Vw1=min(abs(error));                              
" b6 Q( @0 l  J, S" R/ C  rw2=max(abs(error));
! C# q) V+ J% k% L" T6 Ji=1:n;                                     %计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
w=sum(w)/(n-1);
1 T7 Q1 S; V+ r" y) \au                                         %输出参数a,u的值
ago;                                        %输出累加数列ago的值
x0;                                         %输出原始序列值
f=yc;                                        %输出预测的值
+ m! V& X% T  J; ~) |) v8 cerror;                                        %输出残差的值
: K/ |% F, N6 f5 l$ s4 oc;                                           %输出后验差比的值
p;                                           %输出小误差概率的值
) k) `1 n/ R0 U7 ]6 H1 |) K; v7 z6 ow                                           %输出关联度
+ ~: P$ u% l/ b7 C. a: \2 s, z) eend                                   
x0=[31.26
( J9 g. V0 q# P32.09
33.39
6 }, T, I* }4 e! H; R5 D( h2 Q* w7 j7 R35.45
40.52
43.52
6 p# v( N$ X4 l2 Z& E47.86
51.45
55.6
( \; B# b9 C# ?4 ]1 F0 L& I+ r60.14
3 W, t$ ?6 M" _64.82
. x8 M% E' D" k  y) x2 |68.65
73.22
6 k+ _8 T+ A+ [80.22
87.69
93.97
99.16
103.38
# f) B" t/ O5 X6 i0 Q109.46
114.6
! l0 u4 e6 A+ Q. K2 X+ C( ?- O119.85
124.92
132.04
, k+ y3 q% e9 z- l- L8 C# Y0 S139.45
8 Q% U: R3 v0 M6 }: H! h; R9 f150.93
165.13
181.93
196.83
/ t, B2 P. n$ b: B6 [  l2 J212.38
228.07
7 W2 n4 k* U5 v& j  a241.45
251.03
]
我在窗口输入：gm(x0) 为啥会出现
Error in ==> gm at 22
& X% A; A9 ]1 B8 a3 B8 e4 n: verror(i)=abs(yc(i)-x0(i));                     %计算残差值
  H6 l' E  U2 B2 `+ \: ?那位高手帮忙 谢谢

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉