急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
1 ~+ k3 N5 h& x9 [) `* o0 }function f=gm(x0,m)              %定义为函数gm(x)
! B# V/ M$ K* z0 j% ?4 N# s) gn=length(x0);                     
! m' D7 G$ Q  d# t1 Ux1=zeros(1,n);                    
x1(1)=x0(1);                  
* R- D& V: T( _1 r. gfor i=2:n                         %计算累加序列x1
* D/ p. x0 R6 Q" p: W/ G    x1(i)=x1(i-1)+x0(i);         
end
i=2:n;                            %对原始数列平行移位并赋值给y
2 V- ^/ m, \, x* jy(i-1)=x0(i);                     
9 i" l% e2 R! s2 F, j3 Fy=y';                                 
i=1:n-1;                             
c(i)=-0.5*(x1(i)+x1(i+1));                 c(i)=-0.5*(x1(i)+x1(i+1))
, r+ M, H% ^: p8 j7 A8 S2 nB=[c' ones(n-1,1)];                        
+ K9 X: ~. b+ J; ?: u1 ?3 cau=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
                                            
, ^8 N8 f/ [& M0 K. e% I# Q/ Vi=1:m;                             %计算预测累加数列的值
ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
yc(1)=ago(1);                           
i=1:m-1;                             %还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
i=2:n;                              
4 G8 G& M' U) I5 A# m; d5 K7 L2 derror(i)=abs(yc(i)-x0(i));                     %计算残差值
9 Q" w  I/ ?, Q0 fyc(1)=ago(1);                       
+ B! e4 G# i5 {: n. g; E( l6 di=1:m-1;                              %修正还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
- _9 D- V' N( ^) t$ }1 K6 V8 _c=std(error)/std(x0);                      %计算后验差比 *2
3 S* {% c4 a3 L, @/ P# `1 Hrelerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);                       
; ]0 o$ j$ n9 Z2 V7 k7 Gp=0;
* K* H6 M; l, O# a( wfor i=2:ncol                                    
    if  relerror(1,i)<0.6745*std(x0)               
        p=p+1;                                 
/ ^3 R/ K  q. t- F& k' i  [7 I    end
6 F5 V9 |+ k0 v! t/ X$ O# kend
p=p/(n-1);                     
5 U! S' L/ T9 v0 D. z# Nw1=min(abs(error));                              
w2=max(abs(error));
0 _: J2 z; q0 Y' P7 Ai=1:n;                                     %计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
w=sum(w)/(n-1);
au                                         %输出参数a,u的值
ago;                                        %输出累加数列ago的值
x0;                                         %输出原始序列值
) i0 J4 w5 {/ ?2 Y, i; lf=yc;                                        %输出预测的值
7 W" C$ i: t1 v; Qerror;                                        %输出残差的值
c;                                           %输出后验差比的值
* y' T$ r) O. ^4 H: ~3 m) cp;                                           %输出小误差概率的值
w                                           %输出关联度
end                                   
x0=[31.26
/ j) T0 d0 W2 s( q( ~/ r32.09
33.39
35.45
40.52
: z4 S1 u$ ~# J' ^* S3 b" r, y$ a9 u43.52
47.86
" D' V6 t' K' }. L51.45
+ a4 J0 B$ F- j55.6
60.14
$ c2 g: F- _$ D  n+ K64.82
68.65
73.22
& L: u( b1 X% u80.22
5 r/ Y: Y8 a, d- h: K87.69
93.97
- x' M/ n% |% }- l1 D99.16
& X, F4 o# p& j9 r! I' Z% f103.38
109.46
114.6
# A$ t, H4 M3 J# V' U$ H119.85
124.92
6 N4 F2 T6 E( Y$ R132.04
5 `2 B: Z& k# z$ h139.45
150.93
0 P$ }8 |" g! x0 J4 Q! k. `165.13
# l5 w$ f0 m5 O" M181.93
+ ]' r3 o: a  c# R+ ?6 R196.83
: T3 B, ~5 G8 j  D212.38
228.07
: b# |2 D/ P6 H' F' D+ k0 L241.45
8 X8 I6 Y1 D  X0 W. j/ Z3 O- w, e% e251.03
0 J: t9 m9 n( d# f. s]
0 N* q5 ?' P1 u6 b3 @0 w: [+ r" u我在窗口输入：gm(x0) 为啥会出现
Error in ==> gm at 22
error(i)=abs(yc(i)-x0(i));                     %计算残差值
+ a  f, @& X, E& f3 o. M那位高手帮忙 谢谢
4 c6 R. u7 a. I0 X/ w/ f, ?, C

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉