急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
9 x, L; k+ C& N8 Vfunction f=gm(x0,m)              %定义为函数gm(x)
% K3 ]9 b4 b* w2 ~7 c% m6 j+ h, mn=length(x0);                     
9 E* V- b2 {+ a  R1 ~7 Wx1=zeros(1,n);                    
* |/ W9 h! w! Y7 @. j- rx1(1)=x0(1);                  
for i=2:n                         %计算累加序列x1
2 Q% j$ m+ c& M0 D2 @8 D) C. t    x1(i)=x1(i-1)+x0(i);         
end
- e' h5 L, s4 l- j( b4 D( y" @, J/ L) Ji=2:n;                            %对原始数列平行移位并赋值给y
y(i-1)=x0(i);                     
9 K1 i5 g  V& H# e. G1 ], Y8 Oy=y';                                 
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)];                        
, H  j! N) D3 {0 P2 zau=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
                                            
0 \" I# m9 ^: X+ C2 x, |. Y* U; G9 @i=1:m;                             %计算预测累加数列的值
4 K; M! d# F; C, [ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
" l' r  ~- V8 g* ?4 o' [7 Wyc(1)=ago(1);                           
i=1:m-1;                             %还原数列的值
  K; _9 V8 G  E: ?2 }4 L0 v0 vyc(i+1)=ago(i+1)-ago(i);              
6 D  Q, H) E$ h! Ai=2:n;                              
error(i)=abs(yc(i)-x0(i));                     %计算残差值
yc(1)=ago(1);                       
i=1:m-1;                              %修正还原数列的值
; n2 Z+ D$ y2 G2 u, Cyc(i+1)=ago(i+1)-ago(i);              
$ W' f9 S: G: H; Dc=std(error)/std(x0);                      %计算后验差比 *2
% C: w2 P. n% zrelerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);                       
" ^$ Y  U, @9 {( fp=0;
/ _/ b0 o7 ~" h) ?  X: pfor i=2:ncol                                    
    if  relerror(1,i)<0.6745*std(x0)               
        p=p+1;                                 
7 y3 n" n1 s9 R) s  X0 m- K    end
end
2 C7 v& ~* }) v" s" ]7 G0 op=p/(n-1);                     
' Y8 E0 Z7 }0 O( r  r* L/ @w1=min(abs(error));                              
w2=max(abs(error));
3 `( I% w2 y* h( i2 t( Ei=1:n;                                     %计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
w=sum(w)/(n-1);
au                                         %输出参数a,u的值
$ X  h1 q9 G5 s1 i: \) Eago;                                        %输出累加数列ago的值
% \1 r( L/ A: I$ n; \" h! j7 L2 ]x0;                                         %输出原始序列值
' Y' E# r7 _5 y, p- p! Rf=yc;                                        %输出预测的值
4 Z) ]3 F6 K( v6 P8 serror;                                        %输出残差的值
c;                                           %输出后验差比的值
* r! {' ~' C3 U+ ?p;                                           %输出小误差概率的值
w                                           %输出关联度
# K4 k* y: n; G6 Uend                                   
+ f- @7 J) E0 C0 L# w  g8 ox0=[31.26
9 }2 C/ V2 ~- L6 M6 d) ~# n32.09
33.39
35.45
# J1 B; s; }4 i$ I% m1 K40.52
7 P2 x- K7 q2 M& @6 V) g43.52
47.86
; U$ U3 |+ ~+ q9 D51.45
55.6
9 t2 ^# y4 `/ t/ o: m+ N60.14
$ b4 O" R- Q; z' a( m- g; R64.82
68.65
73.22
80.22
0 T$ R, |3 m. U4 ~/ k# ?87.69
93.97
: B, t8 z3 K. ^' Y$ E99.16
( E; m1 b3 r/ ?! Z; F0 t" ~103.38
109.46
114.6
# V! k& I: m& ~$ ]+ d2 n119.85
; C+ y3 ~$ k% v124.92
# @0 \) y7 i' S2 O132.04
9 `4 B3 ^0 E0 e# p/ L139.45
6 d4 S# j( p& J5 C3 h' f150.93
165.13
181.93
" h7 W/ i' A3 {2 ]9 T* |196.83
212.38
8 ^  e2 u* i5 C228.07
241.45
251.03
1 H2 k8 @) Q* t2 p8 R; O+ N]
8 c9 m. j8 H) V1 U我在窗口输入：gm(x0) 为啥会出现
" Z* i+ k6 Z% x% F: \4 h; jError in ==> gm at 22
7 D0 f- ~2 v; D1 n5 n' y$ M* T& derror(i)=abs(yc(i)-x0(i));                     %计算残差值
( f. S' v0 y! O) Y那位高手帮忙 谢谢

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉