急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
) ?. \" u, C9 m% U3 e' Yfunction f=gm(x0,m)              %定义为函数gm(x)
n=length(x0);                     
' Z% q0 t) v! v& Rx1=zeros(1,n);                    
. X1 x) o5 T' Y2 E+ p- z5 V; Z. kx1(1)=x0(1);                  
for i=2:n                         %计算累加序列x1
    x1(i)=x1(i-1)+x0(i);         
end
i=2:n;                            %对原始数列平行移位并赋值给y
5 b5 I7 T) X1 k9 A1 j% X9 h# ~1 Fy(i-1)=x0(i);                     
y=y';                                 
i=1:n-1;                             
c(i)=-0.5*(x1(i)+x1(i+1));                 c(i)=-0.5*(x1(i)+x1(i+1))
" {( {( {" C! ~, @& A: BB=[c' ones(n-1,1)];                        
au=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
                                            
/ ~, \: ^2 G7 |) S4 X2 ]i=1:m;                             %计算预测累加数列的值
ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
4 A. D) A6 i9 I* I7 Qyc(1)=ago(1);                           
1 @3 D; H) z2 g! {( K/ O2 C# `i=1:m-1;                             %还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
0 @; F& X1 d' h. b9 d: |i=2:n;                              
5 C6 }, h2 J8 ]7 ^8 Xerror(i)=abs(yc(i)-x0(i));                     %计算残差值
, _" ?, t8 N: ^, W* n2 d) x6 O% Dyc(1)=ago(1);                       
2 W- {! H1 u. A3 ii=1:m-1;                              %修正还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
, M/ ]% F8 i/ l5 G4 O/ n2 W; lc=std(error)/std(x0);                      %计算后验差比 *2
4 d2 }' z0 N+ w  E" ], A4 H- R$ grelerror=abs((error)-mean(error)*ones(size(error)));
6 V+ t, }, s1 u- |[nrow,ncol]=size(relerror);                       
+ i: D: R6 \. |: F/ \. Tp=0;
1 X$ N2 _0 B9 K7 }$ @( Jfor i=2:ncol                                    
( B/ G! j# X2 P5 d    if  relerror(1,i)<0.6745*std(x0)               
        p=p+1;                                 
6 w& k, h. F1 Q0 R9 X* o0 K( j    end
, d) n6 R: L9 i+ _0 G* E. W7 L; D. Vend
p=p/(n-1);                     
w1=min(abs(error));                              
2 j$ |: d2 R2 @3 C  A. |0 E2 Bw2=max(abs(error));
i=1:n;                                     %计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
! x6 @$ b& t0 {" Dw=sum(w)/(n-1);
+ @. _  o, ^% o/ t! |au                                         %输出参数a,u的值
% S( s& h1 ]1 p3 z7 s" K/ qago;                                        %输出累加数列ago的值
$ \# K* z# a: R5 K' nx0;                                         %输出原始序列值
f=yc;                                        %输出预测的值
0 r( q# J9 G* J) p5 q7 {error;                                        %输出残差的值
" ~4 m/ f' ?- M+ a* E+ bc;                                           %输出后验差比的值
p;                                           %输出小误差概率的值
1 [3 z1 D4 u  z: a- aw                                           %输出关联度
end                                   
x0=[31.26
32.09
33.39
, a9 _- Y/ |% y: T' ]35.45
40.52
43.52
! `9 @7 y" Y; K7 q5 c9 F% w47.86
51.45
55.6
& d. S. l; M9 b4 P7 ~60.14
64.82
( B! B6 A) M/ D. H! V! G! R68.65
73.22
2 f0 ]8 p' i' A0 g# d) u" `$ p80.22
. t' O( {! p3 ?! a* d* c: c  L" T7 V87.69
0 {: [0 C* b  @2 K93.97
99.16
103.38
109.46
' g4 c6 @+ {3 X. [+ ~114.6
119.85
124.92
$ P* }1 |3 z! l2 O  T9 h% r: Z9 B6 @( _132.04
1 x; Q" Q$ R; S139.45
150.93
0 u/ N4 D" ~5 _+ A$ O0 E7 e165.13
181.93
196.83
% |; l- p/ z* F1 a% k212.38
228.07
241.45
( _+ _7 V, e  J+ w& H$ O251.03
]
我在窗口输入：gm(x0) 为啥会出现
Error in ==> gm at 22
error(i)=abs(yc(i)-x0(i));                     %计算残差值
; N/ y' b0 r' E6 l  v那位高手帮忙 谢谢

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉