急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
function f=gm(x0,m)              %定义为函数gm(x)
8 d* {" c% b! _% s# q9 ], K0 tn=length(x0);                     
0 c1 C: I( W( g6 d: m) `3 Rx1=zeros(1,n);                    
x1(1)=x0(1);                  
for i=2:n                         %计算累加序列x1
    x1(i)=x1(i-1)+x0(i);         
: R% h8 M! X' z% |- }& A! send
, m8 Q6 ?% J+ B% d- c; o1 ii=2:n;                            %对原始数列平行移位并赋值给y
y(i-1)=x0(i);                     
1 F% r  H0 B& O8 M' x: Wy=y';                                 
( S; Q4 |/ c  z% P) mi=1:n-1;                             
c(i)=-0.5*(x1(i)+x1(i+1));                 c(i)=-0.5*(x1(i)+x1(i+1))
, F  Y' I, L" ]$ m& {5 R! nB=[c' ones(n-1,1)];                        
0 W$ ?  i5 u6 N! ^$ {au=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
+ F( f8 f& F: ]2 [                                            
& r4 @/ H4 l8 K/ p( Zi=1:m;                             %计算预测累加数列的值
ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
2 k  \" Y  ?; M2 Y, P' C) @# ryc(1)=ago(1);                           
0 G" n+ Q. W# V  o6 q; a6 [: y6 qi=1:m-1;                             %还原数列的值
( n% k  x) |4 ]yc(i+1)=ago(i+1)-ago(i);              
i=2:n;                              
error(i)=abs(yc(i)-x0(i));                     %计算残差值
5 G: N2 V( I  a2 Jyc(1)=ago(1);                       
i=1:m-1;                              %修正还原数列的值
( h6 ]3 ?  l3 ?, @: J. a' pyc(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/ z) K, N0 R2 f! J/ Pp=0;
for i=2:ncol                                    
( h7 ^0 f# l. \8 ~; s    if  relerror(1,i)<0.6745*std(x0)               
* d) L& X' l# m# b. _$ `+ `4 |% [        p=p+1;                                 
7 K2 f# Y/ V6 P" C/ H: u) f9 _- e9 O    end
end
p=p/(n-1);                     
w1=min(abs(error));                              
* g) h$ a" u! Dw2=max(abs(error));
i=1:n;                                     %计算关联度
- F" y# U) {1 O6 L' E0 K0 w/ Aw(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
w=sum(w)/(n-1);
; b% h& V2 K" Y' o' f3 w! \# x) Sau                                         %输出参数a,u的值
' D$ A- P5 Y) J2 ?- {( Aago;                                        %输出累加数列ago的值
9 W$ v7 l/ R! J+ T, A" T7 Fx0;                                         %输出原始序列值
6 X5 \/ O1 @# H# _* \f=yc;                                        %输出预测的值
. [1 x0 o0 Y' z7 w! V' d; a- F+ ?error;                                        %输出残差的值
9 M0 `' N6 O9 [4 W5 Cc;                                           %输出后验差比的值
p;                                           %输出小误差概率的值
w                                           %输出关联度
end                                   
, [- O5 C9 d* m! J" ^( r# w5 Ox0=[31.26
" J  I# m1 ~  e' k32.09
33.39
& R7 d1 e1 T# J' S8 d7 W35.45
; `. b0 d% n& g' y+ ^# x40.52
43.52
# x) i- K3 g( o8 ]8 K% G47.86
51.45
" R5 o- v4 ^& _1 _' A55.6
60.14
64.82
68.65
73.22
80.22
87.69
93.97
99.16
0 q6 h& Y- r2 O. ]0 e* Q0 }103.38
( c  ^3 d3 S. E109.46
114.6
* J) N3 O2 G/ |! o& n1 W& i1 {119.85
124.92
7 r8 A. r( p- O132.04
) \0 [/ m7 N( k* W3 a) b% v3 F; K139.45
150.93
165.13
$ P$ B* u& R( R, q181.93
( O( X, C  f2 e2 K+ M196.83
/ h9 N2 V' N6 D2 Q! p0 h6 R212.38
228.07
241.45
251.03
]
我在窗口输入：gm(x0) 为啥会出现
- R+ A2 S  i0 l# a  zError in ==> gm at 22
error(i)=abs(yc(i)-x0(i));                     %计算残差值
8 Q: y$ H, ]1 ^9 c那位高手帮忙 谢谢
/ z- t3 x! F- V/ p6 \0 S# u

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉