急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
# ?% m% o7 U; t  |function f=gm(x0,m)              %定义为函数gm(x)
2 h' Z* V6 l1 I, Cn=length(x0);                     
3 t" z' ]% u1 m7 L# E- v+ }x1=zeros(1,n);                    
x1(1)=x0(1);                  
for i=2:n                         %计算累加序列x1
    x1(i)=x1(i-1)+x0(i);         
end
, t4 P% s9 [1 x& `- Di=2:n;                            %对原始数列平行移位并赋值给y
8 [9 o7 [  l4 L& Y# U$ My(i-1)=x0(i);                     
+ _9 F  o( y& E) b$ s$ Ly=y';                                 
. o9 |9 A* V5 a" ?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)];                        
1 S3 b- l$ _+ pau=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
6 l" o0 v1 C1 r1 [                                            
2 ?" Z/ y! q7 i- j' ^i=1:m;                             %计算预测累加数列的值
6 o) q; c: \9 G2 B- s1 P- Rago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
yc(1)=ago(1);                           
i=1:m-1;                             %还原数列的值
7 I4 w2 b+ e; s5 I) S+ u! ?yc(i+1)=ago(i+1)-ago(i);              
i=2:n;                              
+ d& w+ e8 A7 g# `& ?0 @7 Uerror(i)=abs(yc(i)-x0(i));                     %计算残差值
  }! x1 ~0 a: M6 _; x; ^yc(1)=ago(1);                       
% L! O. k4 G5 ]) w* t9 G# d3 ii=1:m-1;                              %修正还原数列的值
: @, r$ }: b( i! M! G, t0 pyc(i+1)=ago(i+1)-ago(i);              
0 a3 z- h: q3 @  T4 o- D1 d) Hc=std(error)/std(x0);                      %计算后验差比 *2
, B4 M! e4 O' i( e4 a4 lrelerror=abs((error)-mean(error)*ones(size(error)));
' V* c/ P! s7 L/ v[nrow,ncol]=size(relerror);                       
& d# c! t" A/ w" O/ k4 N5 Y. e/ k& r* Z5 Pp=0;
1 Y% k- `# F7 p! j8 }3 B: Jfor i=2:ncol                                    
3 s% I2 F0 z, g% E* w7 W    if  relerror(1,i)<0.6745*std(x0)               
2 j1 @( J. B" S) z/ w, U        p=p+1;                                 
9 y! e7 a* F1 ~% l% X7 I  E2 K! i    end
+ M4 {4 Y" H, Q' X" F1 X: Tend
p=p/(n-1);                     
/ f/ r% E( X# N% Kw1=min(abs(error));                              
w2=max(abs(error));
% f: s) H9 [' M- C% ~i=1:n;                                     %计算关联度
" G# P" M' Y1 xw(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
# G7 F, w7 x8 P/ H+ _w=sum(w)/(n-1);
) B+ e1 R- Y: @8 s  C" Vau                                         %输出参数a,u的值
; `4 m7 T" i* C% R2 tago;                                        %输出累加数列ago的值
x0;                                         %输出原始序列值
f=yc;                                        %输出预测的值
error;                                        %输出残差的值
! b) G7 p6 G( yc;                                           %输出后验差比的值
p;                                           %输出小误差概率的值
w                                           %输出关联度
end                                   
* a/ [8 c( K2 J; cx0=[31.26
32.09
! a* S* \1 s+ p- p! Z4 y7 t33.39
35.45
40.52
43.52
47.86
2 t' w* C6 K4 j8 {  w- B0 f51.45
( U' s% H/ R! U6 K3 L55.6
60.14
+ w- l* l' @2 b" W+ A64.82
" E0 L8 i# Q3 }( @* ^  D. F6 j68.65
  [- [) m8 l: l+ k73.22
80.22
% {* s/ F. d& r# [: f8 ^87.69
93.97
. T7 [9 o- R( L6 Z, {+ z99.16
# _4 D5 p4 u# S0 _. x103.38
109.46
114.6
119.85
% ]% w$ [$ N4 ~3 y8 G124.92
132.04
139.45
+ E% c4 u' r2 l; ]1 l. e150.93
  n: m3 W2 C4 B' r165.13
9 ?# \' U! h# Y- z! E181.93
196.83
9 ]) q0 |& Y2 e9 W' @6 z! b: m' Q212.38
228.07
241.45
# `2 w1 m/ o* y: a& m! A3 u' Q251.03
2 O4 L; C2 h. ~) [0 J6 V8 A5 j]
我在窗口输入：gm(x0) 为啥会出现
Error in ==> gm at 22
error(i)=abs(yc(i)-x0(i));                     %计算残差值
) [4 t- A! J9 t* V8 u& i那位高手帮忙 谢谢

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉