急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
% W; p4 Q0 N7 Z+ p/ T, v# X$ jfunction f=gm(x0,m)              %定义为函数gm(x)
/ P+ e3 F2 a- J! X4 C/ qn=length(x0);                     
x1=zeros(1,n);                    
x1(1)=x0(1);                  
for i=2:n                         %计算累加序列x1
; n/ A) q2 j+ \! w2 E; m    x1(i)=x1(i-1)+x0(i);         
& R- U% J+ i9 o+ s5 ]. U3 U9 oend
7 Q2 J( {% N! n+ h1 V2 `# \2 Ji=2:n;                            %对原始数列平行移位并赋值给y
y(i-1)=x0(i);                     
* d! g) L" o. O4 y- t% C1 A6 n+ Vy=y';                                 
i=1:n-1;                             
: d( k- F5 u) K2 U' uc(i)=-0.5*(x1(i)+x1(i+1));                 c(i)=-0.5*(x1(i)+x1(i+1))
+ ^/ b7 m# y3 CB=[c' ones(n-1,1)];                        
" W- t* M& W# v. {au=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
1 K5 S2 i7 n8 \' ^5 i. \( W2 x* }                                            
i=1:m;                             %计算预测累加数列的值
ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
5 [  e) e, e' L) M$ F9 {yc(1)=ago(1);                           
9 q5 W7 `2 e/ J) w  A. @2 W, t# E. K" Vi=1:m-1;                             %还原数列的值
( r+ L+ y" K6 a9 I" Iyc(i+1)=ago(i+1)-ago(i);              
i=2:n;                              
error(i)=abs(yc(i)-x0(i));                     %计算残差值
yc(1)=ago(1);                       
7 _  f- V) j5 ^9 G# k; li=1:m-1;                              %修正还原数列的值
; Y& T7 S( \) G, pyc(i+1)=ago(i+1)-ago(i);              
& e. h& Q; n7 V* Q4 a$ Y8 y5 qc=std(error)/std(x0);                      %计算后验差比 *2
% t( u0 z% M& B. f) Xrelerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);                       
1 ~3 Y, I  f3 Z4 _  sp=0;
for i=2:ncol                                    
  u9 s7 @' F& Q4 A0 H7 G    if  relerror(1,i)<0.6745*std(x0)               
        p=p+1;                                 
+ f5 Q% L2 g3 a) C    end
* n4 u  x7 H0 `, Cend
* P% N9 ^% h$ A" Wp=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);
7 f5 D8 y0 _4 Gw=sum(w)/(n-1);
au                                         %输出参数a,u的值
ago;                                        %输出累加数列ago的值
. T) h6 q3 }/ }* |& _x0;                                         %输出原始序列值
f=yc;                                        %输出预测的值
error;                                        %输出残差的值
# \8 Q3 S. j$ k  A0 ]! Jc;                                           %输出后验差比的值
; @: N) N0 Y1 W1 H& |7 Mp;                                           %输出小误差概率的值
- _$ g/ l8 u. Lw                                           %输出关联度
0 `. ~1 B* L; pend                                   
& r" A0 w6 T. Ax0=[31.26
32.09
/ D2 W$ k" U/ G. J/ H, l5 O33.39
35.45
40.52
2 m; w9 ]5 k$ {$ @8 }8 P# Y43.52
6 J9 u- L4 G' b8 B- l) _47.86
51.45
55.6
60.14
) a5 H+ X2 z: j- j8 s  H' W64.82
& P$ p) f0 e/ {68.65
73.22
80.22
87.69
- c. l+ M6 |6 _9 [: R. C0 P93.97
99.16
103.38
109.46
114.6
119.85
124.92
132.04
139.45
150.93
  c! H2 ?2 g8 E. u6 n165.13
. J, \* m7 Y3 ]1 ?  [0 X4 m3 e; d181.93
196.83
212.38
6 y! F  h) T' ]% e* t( D0 p0 \228.07
0 [2 A8 C1 z6 d  {5 i' T241.45
+ c. {2 P  p2 Q  S8 k8 y7 o251.03
]
我在窗口输入：gm(x0) 为啥会出现
' h+ `6 n3 ?9 c1 a- }  F0 VError in ==> gm at 22
- o. i$ F- k& K/ K* Verror(i)=abs(yc(i)-x0(i));                     %计算残差值
那位高手帮忙 谢谢
: T; J# [8 r; }8 W

wuzhenhua

把abs（）去掉

神经病个人组

我也想知道