急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
9 J% K0 n, o" n4 `: J- N. P: Bfunction f=gm(x0,m)              %定义为函数gm(x)
n=length(x0);                     
, Y- }5 v- m3 |2 [x1=zeros(1,n);                    
2 v8 |- W# ^) |; |x1(1)=x0(1);                  
& M# ~) N7 ~) M8 w2 r& H" E7 ifor i=2:n                         %计算累加序列x1
    x1(i)=x1(i-1)+x0(i);         
+ b; J. G( W' A- l7 O. c* |end
i=2:n;                            %对原始数列平行移位并赋值给y
y(i-1)=x0(i);                     
y=y';                                 
$ S- u! U/ ~" h# w1 U0 ?i=1:n-1;                             
c(i)=-0.5*(x1(i)+x1(i+1));                 c(i)=-0.5*(x1(i)+x1(i+1))
) g. o6 ~: q. ^6 sB=[c' ones(n-1,1)];                        
au=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
: A6 R" C, u1 |' ?                                            
# z5 H7 Y- I$ a2 ]7 Yi=1:m;                             %计算预测累加数列的值
8 y9 u- ^! H% \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);              
, V' C6 Y9 B& v) I1 T9 K7 K5 e4 ai=2:n;                              
error(i)=abs(yc(i)-x0(i));                     %计算残差值
. a. T( b, N; Z) Tyc(1)=ago(1);                       
7 Q' s1 ~+ o# _- Q4 }" o( ni=1:m-1;                              %修正还原数列的值
$ J5 P5 y6 @( g- R( V  v% D9 oyc(i+1)=ago(i+1)-ago(i);              
" p6 @+ y- H% P& i9 z, `c=std(error)/std(x0);                      %计算后验差比 *2
relerror=abs((error)-mean(error)*ones(size(error)));
+ |$ E) x% Y) \; o* \7 R[nrow,ncol]=size(relerror);                       
: _5 D& C) V' U6 p) Gp=0;
for i=2:ncol                                    
    if  relerror(1,i)<0.6745*std(x0)               
  `4 O( y$ J3 B0 l0 b2 j3 c, r        p=p+1;                                 
0 D4 e7 V( O+ v5 K+ b    end
end
p=p/(n-1);                     
! e7 l+ Y% j. {4 |w1=min(abs(error));                              
w2=max(abs(error));
i=1:n;                                     %计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
w=sum(w)/(n-1);
" G6 ^, V& ^2 N6 \' B8 K$ n: Yau                                         %输出参数a,u的值
ago;                                        %输出累加数列ago的值
x0;                                         %输出原始序列值
f=yc;                                        %输出预测的值
( u8 i# Q1 F  }7 S! j0 x" l* h: _error;                                        %输出残差的值
% D' g# A% L+ f4 cc;                                           %输出后验差比的值
p;                                           %输出小误差概率的值
w                                           %输出关联度
% ^$ z; u4 p1 \5 p% u0 |" jend                                   
* N1 G5 h( E1 `3 i$ b7 U* U- N0 xx0=[31.26
32.09
33.39
35.45
40.52
; a5 [  _5 B& ^2 ^8 j; n+ }43.52
0 K$ l  p: Q) [, t0 u. P47.86
51.45
55.6
% e0 G0 Q: w3 R3 t60.14
. K" `! R; N! {6 i  [+ ^  J64.82
- g" I, C+ U0 N- S68.65
1 t  W+ \( z* m* ]4 f73.22
80.22
. l* Y. I4 G; z. p/ I) M$ v. `! S) U87.69
6 [" U2 r- F. F1 o0 U93.97
% }+ u1 @* e  D6 O7 p) W99.16
103.38
# f) H* z( ]8 O9 z2 p7 A109.46
$ u9 m/ w9 @* N0 `3 R" U114.6
9 c2 k2 D# X$ ~6 H1 q3 x  w119.85
124.92
0 Z, R% T- W  H5 P5 r- j132.04
139.45
150.93
165.13
181.93
9 s+ V2 l3 F3 {( i  ^5 ^196.83
6 G- Q. w. {5 y6 b212.38
228.07
241.45
6 V8 b# @9 G+ Q5 i251.03
]
我在窗口输入：gm(x0) 为啥会出现
8 {( q0 D9 Q1 C" _( KError in ==> gm at 22
; D7 a! W  L1 E, z! U+ verror(i)=abs(yc(i)-x0(i));                     %计算残差值
那位高手帮忙 谢谢

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉