急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
function f=gm(x0,m)              %定义为函数gm(x)
0 A3 ]! e2 @( `n=length(x0);                     
# C- D/ S$ `8 ~0 Hx1=zeros(1,n);                    
( h- K/ ^9 `5 l. f* M* U( fx1(1)=x0(1);                  
  C. [. i* G3 a5 ]2 x2 i4 f% X& afor i=2:n                         %计算累加序列x1
: K7 N4 X: p$ Y' @! Z" I# h    x1(i)=x1(i-1)+x0(i);         
! ~/ \) P+ m9 I8 ]end
9 q4 s1 N5 Y# x7 I- wi=2:n;                            %对原始数列平行移位并赋值给y
y(i-1)=x0(i);                     
/ [2 W/ {4 p; n1 q, y  b4 d5 K8 Xy=y';                                 
2 F) V+ z+ v$ i( P: _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为弧度
                                            
# a/ {7 T8 i1 k" Vi=1:m;                             %计算预测累加数列的值
ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
4 h* T7 B" a7 [yc(1)=ago(1);                           
& k0 s* ]3 u) s+ d# ii=1:m-1;                             %还原数列的值
$ u: ^/ S) r5 `/ K/ Kyc(i+1)=ago(i+1)-ago(i);              
( N& S; V6 r" G  O; U) Fi=2:n;                              
6 t0 m' M) j+ P  n8 Zerror(i)=abs(yc(i)-x0(i));                     %计算残差值
, M, \' G1 X9 n+ t( lyc(1)=ago(1);                       
i=1:m-1;                              %修正还原数列的值
' J  \0 b8 T; {' F# o0 Q2 Cyc(i+1)=ago(i+1)-ago(i);              
( s  y7 i) E" Z$ Oc=std(error)/std(x0);                      %计算后验差比 *2
* l2 d/ P- L2 k' I2 c1 ]relerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);                       
* N: s9 N' v. K* m+ w+ p, \6 @p=0;
for i=2:ncol                                    
* |) R& J2 Z! S& @    if  relerror(1,i)<0.6745*std(x0)               
3 ^$ w  J' m% U8 M% R        p=p+1;                                 
    end
end
0 d/ v) P7 F6 H0 o6 n' K% vp=p/(n-1);                     
w1=min(abs(error));                              
w2=max(abs(error));
i=1:n;                                     %计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
w=sum(w)/(n-1);
/ _. `' K! Y) M7 g/ I* Eau                                         %输出参数a,u的值
ago;                                        %输出累加数列ago的值
x0;                                         %输出原始序列值
  P5 `( ~2 [6 j6 g8 m, Gf=yc;                                        %输出预测的值
error;                                        %输出残差的值
c;                                           %输出后验差比的值
  @  y* e% I. g! i' i$ Zp;                                           %输出小误差概率的值
w                                           %输出关联度
3 ~2 O! e# V7 ^& w) x9 \end                                   
x0=[31.26
2 }* _' n3 E' Z1 a( T32.09
2 S6 ?7 k4 `& F0 |* z2 X6 }# g33.39
$ L/ c2 L8 b$ L35.45
& L) ~; r! Z* R# t5 `6 x, a& Z5 n40.52
43.52
$ ~, }' a# r- [- P47.86
* c6 u2 S+ }; }# p- R51.45
5 O7 S  t- J  m( M: c! A  w55.6
60.14
$ M1 e/ D  H, ?4 C; J. n64.82
8 I9 g( G( H1 N: S68.65
73.22
) `2 R6 B+ M, @0 I0 r/ r80.22
87.69
" x8 b4 j! w6 P" @7 j( G93.97
99.16
103.38
109.46
% ^( `* w2 A4 }, J, M114.6
0 G" ~% m5 e, _* ]- e5 Z119.85
' v/ V2 @. S2 p8 t+ G3 e9 [124.92
5 N' k& \# H6 a9 ^# U( P6 }132.04
139.45
4 W+ ]: i; N3 [/ i$ q" `150.93
9 o) c: d( C4 Q# C165.13
; ^" K( Y& q( T, l" r( \4 L181.93
196.83
; B, @/ ?/ J2 z3 o  z212.38
3 v5 g% M! i2 S0 y1 j8 m% k228.07
241.45
( S* Y1 q4 q" ]251.03
]
我在窗口输入：gm(x0) 为啥会出现
Error in ==> gm at 22
9 \+ s9 P! ^) Y: f- Serror(i)=abs(yc(i)-x0(i));                     %计算残差值
5 v: a+ g, C; b& @; z那位高手帮忙 谢谢

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉