急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
& A3 A' |- ~/ ~) B, v2 Mfunction f=gm(x0,m)              %定义为函数gm(x)
n=length(x0);                     
  g1 W) K/ _( Q! zx1=zeros(1,n);                    
' j" E0 {) L0 H7 Nx1(1)=x0(1);                  
for i=2:n                         %计算累加序列x1
    x1(i)=x1(i-1)+x0(i);         
end
+ v& F& k) O2 ^6 ?- n& U2 ?i=2:n;                            %对原始数列平行移位并赋值给y
. z( \7 r. D; x: m5 ]+ Z: P+ Yy(i-1)=x0(i);                     
y=y';                                 
$ h9 e4 K( Z: b! Si=1:n-1;                             
5 h/ P; \. Y1 v, N0 Hc(i)=-0.5*(x1(i)+x1(i+1));                 c(i)=-0.5*(x1(i)+x1(i+1))
B=[c' ones(n-1,1)];                        
2 n( `, Y$ C0 z8 f4 ^- e  c' xau=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
                                            
i=1:m;                             %计算预测累加数列的值
ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
0 g7 n  I8 v+ v" t9 g0 Lyc(1)=ago(1);                           
i=1:m-1;                             %还原数列的值
9 x$ p( C2 r- ?yc(i+1)=ago(i+1)-ago(i);              
( Y) ~( v1 G5 O1 e& D* Pi=2:n;                              
error(i)=abs(yc(i)-x0(i));                     %计算残差值
1 L; d, G, o* b4 o  j  Nyc(1)=ago(1);                       
i=1:m-1;                              %修正还原数列的值
, b/ d) t4 t, K! ryc(i+1)=ago(i+1)-ago(i);              
8 }- Y" S. e' t% f1 Q% yc=std(error)/std(x0);                      %计算后验差比 *2
. S) T& V, g. N& j2 H( Erelerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);                       
4 ]2 a) }. M4 Np=0;
, o0 `1 x, P3 X& Afor i=2:ncol                                    
8 @5 A7 X% y4 c    if  relerror(1,i)<0.6745*std(x0)               
- n8 ~. \' M- R; U        p=p+1;                                 
) I) _5 O3 P# `+ @5 c8 O/ j    end
end
) z5 |; \6 p& b4 s6 S7 vp=p/(n-1);                     
w1=min(abs(error));                              
( \- |; P: M5 y1 \% J" p+ Aw2=max(abs(error));
i=1:n;                                     %计算关联度
9 Z& r  f3 @8 ^% ~: `! F4 `: Tw(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
) b. `& Y) e3 `8 ]! ^w=sum(w)/(n-1);
) K+ F3 [3 x8 ?' hau                                         %输出参数a,u的值
+ }- W% H0 y2 ]/ m2 K7 W* X' _; Sago;                                        %输出累加数列ago的值
x0;                                         %输出原始序列值
7 g! M" d' `% _! _4 f+ [$ ?f=yc;                                        %输出预测的值
: |! K  C- ^" e1 @  ~error;                                        %输出残差的值
c;                                           %输出后验差比的值
8 b2 M- D8 {+ T9 Zp;                                           %输出小误差概率的值
& E0 E9 i) Q6 N4 A' kw                                           %输出关联度
end                                   
x0=[31.26
) L& y  y; R" M# t32.09
33.39
, r8 s6 G% Z8 H+ x3 F+ z35.45
! v( \/ o5 g9 i6 O  x3 i1 \4 S% ~. a40.52
43.52
47.86
51.45
55.6
60.14
" Q0 \0 I& k1 ^64.82
68.65
/ y% v$ y; d+ g$ t73.22
# O% i5 K' Y& Z/ {6 R4 ^80.22
) x, q- B, V) A6 j# \, K1 a! l6 o/ i7 [87.69
93.97
99.16
103.38
. t2 [- r: G9 F: W109.46
114.6
119.85
; f$ k9 a' n: Q" ~+ [9 }124.92
1 {: J. m7 e9 a7 L" G132.04
* Q; j: s* G* @139.45
150.93
* O* |6 }. _- m; e2 L165.13
# k9 G2 x; |" T; z181.93
196.83
212.38
228.07
% a* m' o" \7 b( k7 {7 i241.45
251.03
( U5 m  G  Q* M, d]
我在窗口输入：gm(x0) 为啥会出现
7 y" v9 L5 x! EError in ==> gm at 22
error(i)=abs(yc(i)-x0(i));                     %计算残差值
! u$ o) N# F; q" [& H6 C  f那位高手帮忙 谢谢

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉