急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
# \. j# l9 G% ^) Tfunction f=gm(x0,m)              %定义为函数gm(x)
, T+ ^& Q+ u" g; V, |+ [& @n=length(x0);                     
x1=zeros(1,n);                    
x1(1)=x0(1);                  
+ }+ A8 V7 Y5 I  u* wfor i=2:n                         %计算累加序列x1
& H* u! N% _7 }3 X' B; r    x1(i)=x1(i-1)+x0(i);         
end
i=2:n;                            %对原始数列平行移位并赋值给y
y(i-1)=x0(i);                     
y=y';                                 
i=1:n-1;                             
2 P' {6 h8 E, y' ~c(i)=-0.5*(x1(i)+x1(i+1));                 c(i)=-0.5*(x1(i)+x1(i+1))
B=[c' ones(n-1,1)];                        
au=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
8 @2 \' w) B. S                                            
* S* B, Q& k9 f9 Ai=1:m;                             %计算预测累加数列的值
ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
- z/ u6 K; m* P" q3 A2 eyc(1)=ago(1);                           
i=1:m-1;                             %还原数列的值
, [+ g, r, P  e; j% K% U$ ^1 e' nyc(i+1)=ago(i+1)-ago(i);              
  [5 u8 q& U, H' X! u  C) M$ _i=2:n;                              
error(i)=abs(yc(i)-x0(i));                     %计算残差值
* u7 s& N' A5 syc(1)=ago(1);                       
i=1:m-1;                              %修正还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
c=std(error)/std(x0);                      %计算后验差比 *2
: }/ \: c% W& s3 H; b) Xrelerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);                       
p=0;
for i=2:ncol                                    
    if  relerror(1,i)<0.6745*std(x0)               
        p=p+1;                                 
    end
end
p=p/(n-1);                     
w1=min(abs(error));                              
( ^$ L0 Z, i! Y1 B1 j1 |: xw2=max(abs(error));
! C* [" h7 {' j3 n/ Z: Ni=1:n;                                     %计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
$ E6 U2 V; {/ G1 \% Nw=sum(w)/(n-1);
% X/ u) y0 Z4 f0 k/ r0 oau                                         %输出参数a,u的值
7 k* v7 s) ]/ \8 Z/ s6 m4 qago;                                        %输出累加数列ago的值
4 ~) x$ f/ y. d" E) ~x0;                                         %输出原始序列值
& r* a9 e3 M8 F" z3 j( _0 w+ {f=yc;                                        %输出预测的值
error;                                        %输出残差的值
c;                                           %输出后验差比的值
$ i' F) A! B! D# [p;                                           %输出小误差概率的值
+ I1 Z+ `! M) F6 X/ T8 Uw                                           %输出关联度
end                                   
x0=[31.26
, A9 ?) R. @- K1 N7 W32.09
33.39
35.45
40.52
43.52
47.86
( L1 b$ Q, k5 V% s- K" |3 }51.45
9 N/ d' d$ }3 P0 M# @55.6
60.14
64.82
68.65
' u* l. G3 |3 T8 ?4 \6 e: l73.22
80.22
87.69
93.97
' l! t+ Y% b. ]$ H6 [99.16
103.38
, I$ t" J( h% N: D! f. I+ O2 x109.46
114.6
119.85
124.92
132.04
/ X( D" W" @5 ]; B1 ]139.45
150.93
! f, a/ y% [; L8 j/ h9 }165.13
, ?/ [; s' |  H, W5 A6 ?181.93
196.83
! L; f2 r$ T, x212.38
228.07
9 L  h  Y! J0 i5 V  L241.45
251.03
* G' z, S1 {! A+ M* L8 \9 g. h3 A]
我在窗口输入：gm(x0) 为啥会出现
Error in ==> gm at 22
error(i)=abs(yc(i)-x0(i));                     %计算残差值
4 E' E7 l! m& o那位高手帮忙 谢谢

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉