急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
) j( u* N7 N# x7 u# Q9 {+ F5 cfunction f=gm(x0,m)              %定义为函数gm(x)
2 o" r0 B! o+ ?2 u3 p0 sn=length(x0);                     
x1=zeros(1,n);                    
x1(1)=x0(1);                  
for i=2:n                         %计算累加序列x1
7 c2 a) F4 b. J* v1 U. k    x1(i)=x1(i-1)+x0(i);         
end
4 P* l5 C3 m1 `: ?& h$ Ki=2:n;                            %对原始数列平行移位并赋值给y
y(i-1)=x0(i);                     
y=y';                                 
" i& f2 L  H6 h' f3 h, o% ti=1:n-1;                             
c(i)=-0.5*(x1(i)+x1(i+1));                 c(i)=-0.5*(x1(i)+x1(i+1))
7 {  Z2 y* e7 y; q; SB=[c' ones(n-1,1)];                        
- E! Y( n! R9 j) n# u$ A* dau=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
* f4 g2 h0 {3 x. d5 c0 d                                            
3 `! N" ~5 C# c1 M4 pi=1:m;                             %计算预测累加数列的值
ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
yc(1)=ago(1);                           
# ~3 e2 K- _. O$ I* v. a( xi=1:m-1;                             %还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
9 D0 D$ c0 v. D! `( J% li=2:n;                              
' n  O# Q' a$ G- b; @1 ferror(i)=abs(yc(i)-x0(i));                     %计算残差值
6 U$ L5 t8 t% J% u$ {yc(1)=ago(1);                       
i=1:m-1;                              %修正还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
0 |1 k- X$ |' rc=std(error)/std(x0);                      %计算后验差比 *2
relerror=abs((error)-mean(error)*ones(size(error)));
2 w. i% z. E) o* o[nrow,ncol]=size(relerror);                       
p=0;
for i=2:ncol                                    
) N! i& o+ }" J* W: y    if  relerror(1,i)<0.6745*std(x0)               
        p=p+1;                                 
    end
2 G, E1 h- E3 f! ]% \end
' J: I2 U, \3 L1 h& @# m! ^7 Hp=p/(n-1);                     
w1=min(abs(error));                              
w2=max(abs(error));
i=1:n;                                     %计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
7 j9 \6 V% J/ G: l3 |  V* N6 ]w=sum(w)/(n-1);
au                                         %输出参数a,u的值
8 G. t2 [% J6 z7 y/ @' L9 e+ q0 ^& tago;                                        %输出累加数列ago的值
x0;                                         %输出原始序列值
2 S# s- B" \9 m* Jf=yc;                                        %输出预测的值
error;                                        %输出残差的值
8 ~( K' }, M/ ]/ R, x# Dc;                                           %输出后验差比的值
p;                                           %输出小误差概率的值
w                                           %输出关联度
end                                   
& S7 T9 T' H# i7 g: lx0=[31.26
32.09
  Q4 A, V6 ~; `0 z: D! c7 c; V5 d33.39
35.45
40.52
0 J- U; |+ k2 ]4 m* o1 Q43.52
& f! `8 |) z9 e1 l47.86
( c  a( C0 I$ z# y% v  ~51.45
55.6
60.14
64.82
68.65
73.22
9 n" ]/ a' g  U: w% e9 C' I80.22
  y' e- H% J* ^) x. [87.69
93.97
99.16
% W& l( Z8 S' D- g% H9 w103.38
109.46
114.6
/ w( l6 X* V" s5 }1 ?; L& U" b119.85
124.92
0 [8 |7 H+ ]- \6 b$ M132.04
" N9 ?/ S& G2 b4 G# g; J  ?139.45
150.93
: n* n3 g9 }3 e' g4 {: U165.13
181.93
% x% v5 _. t  J, F) U. }196.83
212.38
228.07
241.45
251.03
]
我在窗口输入：gm(x0) 为啥会出现
Error in ==> gm at 22
7 v5 H1 G/ n6 Berror(i)=abs(yc(i)-x0(i));                     %计算残差值
那位高手帮忙 谢谢
  x+ _& i0 g& \! T& @6 k; L6 ?

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉