急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
: n0 q  T6 F& K7 e) kfunction f=gm(x0,m)              %定义为函数gm(x)
( j+ H6 O4 _+ `! E" _: C; ^n=length(x0);                     
. E. }9 A: w+ G& ]x1=zeros(1,n);                    
x1(1)=x0(1);                  
for i=2:n                         %计算累加序列x1
5 k4 {3 T2 B# q& [- \) }' O3 b    x1(i)=x1(i-1)+x0(i);         
) w9 d2 c; _7 ^% @end
1 W6 Z3 K: @& }i=2:n;                            %对原始数列平行移位并赋值给y
y(i-1)=x0(i);                     
y=y';                                 
i=1:n-1;                             
9 g' F6 p! N& `. Tc(i)=-0.5*(x1(i)+x1(i+1));                 c(i)=-0.5*(x1(i)+x1(i+1))
" s. |/ h- O8 v- s1 y& bB=[c' ones(n-1,1)];                        
au=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
                                            
% ?1 W9 I0 C0 ~, J! A$ Q: k2 ~& ei=1:m;                             %计算预测累加数列的值
7 e) B$ i4 V9 m8 e% [) \ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
yc(1)=ago(1);                           
8 w) {# b3 b1 \$ k5 O! L. @2 ui=1:m-1;                             %还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
3 Y- d9 _. X( y( o! yi=2:n;                              
: M. _# Z  z; N2 Q' [* Jerror(i)=abs(yc(i)-x0(i));                     %计算残差值
yc(1)=ago(1);                       
2 ~  T5 Z& }. T4 \- ei=1:m-1;                              %修正还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
c=std(error)/std(x0);                      %计算后验差比 *2
; A; @4 k  o" P- C$ g% vrelerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);                       
p=0;
3 D, E4 `+ h! c0 Tfor i=2:ncol                                    
    if  relerror(1,i)<0.6745*std(x0)               
+ ]$ k% Y/ X9 d+ }8 }( |        p=p+1;                                 
    end
end
# L( W! }; t7 \# w5 a# K+ `( Pp=p/(n-1);                     
w1=min(abs(error));                              
7 ~) \8 }! q0 H! Rw2=max(abs(error));
i=1:n;                                     %计算关联度
( ^* h' W- F2 b7 h! {. S* |w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
) x4 l/ d5 j5 ~9 R( H) xw=sum(w)/(n-1);
7 j1 s- ?' y& ]0 T5 E! Tau                                         %输出参数a,u的值
% p3 n8 T* v0 k5 y3 F8 C" h( l, ]ago;                                        %输出累加数列ago的值
# e+ n' {6 d) r; `7 kx0;                                         %输出原始序列值
f=yc;                                        %输出预测的值
error;                                        %输出残差的值
# W& O/ H2 }& A$ j$ ec;                                           %输出后验差比的值
- H& L  p1 `* Q' C5 |p;                                           %输出小误差概率的值
9 w4 s* f+ M/ k2 |6 C+ Gw                                           %输出关联度
end                                   
' ^5 b$ D4 s' p* f1 vx0=[31.26
) w4 B6 j/ g8 K: g1 D32.09
33.39
35.45
5 `! P5 O& b9 z  c40.52
) L: H8 w9 h' h  y# A% a43.52
47.86
51.45
55.6
2 r/ c8 u! B9 D5 v- S0 P60.14
# Z$ S" w( S0 B# ~- P64.82
68.65
* h- J( u% F5 }* h2 g. Z73.22
; b0 [. d9 d  j; V+ ^80.22
9 q% s8 s9 H" h" D* h& [. l87.69
93.97
2 Y# R+ Q: U  E7 M99.16
7 c- \# [$ x) ^8 g103.38
109.46
' U1 u# h! t# G& k3 T' s114.6
# e+ [7 |/ H6 Y) `0 ], j119.85
  |5 v7 f( p+ J# x9 E/ t) I124.92
132.04
139.45
1 [# m$ K$ h' E/ W0 Q7 R' U150.93
) _" C' `2 ]% F# g% S  k165.13
3 S: X% f% b( v9 u6 |181.93
196.83
) n3 Z' Z3 p7 t4 b2 W1 Y. q( q0 Q212.38
228.07
241.45
251.03
]
我在窗口输入：gm(x0) 为啥会出现
Error in ==> gm at 22
0 O0 ^$ F6 k/ l" }( \error(i)=abs(yc(i)-x0(i));                     %计算残差值
$ k/ f5 M1 u  p, n那位高手帮忙 谢谢

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉