急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
4 h  {- M/ r) z* n- d1 z8 k/ c' gfunction f=gm(x0,m)              %定义为函数gm(x)
1 z1 n" p) B$ j& t2 R" ?n=length(x0);                     
x1=zeros(1,n);                    
x1(1)=x0(1);                  
for i=2:n                         %计算累加序列x1
8 A( c  U9 \/ f- S% c2 r  Y    x1(i)=x1(i-1)+x0(i);         
end
6 D' G% |. {/ T; h7 Mi=2:n;                            %对原始数列平行移位并赋值给y
y(i-1)=x0(i);                     
+ t0 N8 a# B1 p5 Y% ~/ oy=y';                                 
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)];                        
au=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
1 ^) ]% n# k3 F& |- y) u                                            
i=1:m;                             %计算预测累加数列的值
ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
! t6 }! M9 o: V4 H4 C' ?. tyc(1)=ago(1);                           
6 l- y# s/ h8 [  A# Ci=1:m-1;                             %还原数列的值
- v6 }/ L' Q) zyc(i+1)=ago(i+1)-ago(i);              
- u9 u1 t1 a8 `, j" F/ ei=2:n;                              
3 Y' @8 x) b0 Nerror(i)=abs(yc(i)-x0(i));                     %计算残差值
yc(1)=ago(1);                       
3 M  C4 I  A: ^3 M2 c3 ~0 V6 g5 pi=1:m-1;                              %修正还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
& t5 J4 E' h  ^+ H' \c=std(error)/std(x0);                      %计算后验差比 *2
( \, I& `8 Z6 @relerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);                       
p=0;
5 d! @, Q1 f$ c9 w! F% `for i=2:ncol                                    
    if  relerror(1,i)<0.6745*std(x0)               
, b9 T% \& g% j/ D4 F; I        p=p+1;                                 
+ i" ?0 }0 Y3 O" n: q1 d    end
end
p=p/(n-1);                     
w1=min(abs(error));                              
w2=max(abs(error));
" W$ o2 o4 n$ w& oi=1:n;                                     %计算关联度
% |# n# \2 c+ vw(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
3 |( s5 c! H0 [! C, Gw=sum(w)/(n-1);
au                                         %输出参数a,u的值
ago;                                        %输出累加数列ago的值
x0;                                         %输出原始序列值
f=yc;                                        %输出预测的值
error;                                        %输出残差的值
c;                                           %输出后验差比的值
/ ~# k( {  |4 p' A3 |. F  F! M" F. Q6 pp;                                           %输出小误差概率的值
w                                           %输出关联度
end                                   
x0=[31.26
32.09
33.39
+ G  {- M3 }. J# Q- N, {4 ?35.45
, z3 E# _5 K  L5 u# L40.52
43.52
47.86
! i$ j# \2 c% J: _# B) Y$ n/ b51.45
55.6
60.14
5 s5 ?6 G% d- V$ V5 B# Q64.82
68.65
73.22
, B* m) _2 ]- Q- s- o0 f80.22
87.69
, Z* @; j! [; @' Q: P93.97
99.16
103.38
109.46
114.6
119.85
124.92
132.04
139.45
150.93
+ [7 ]. N% G1 Z% q2 W165.13
) Y/ f, B+ Y7 W( I; P181.93
196.83
" p3 l0 ]$ V+ w, N4 a212.38
; `% |, t) v- \' P228.07
" J; O5 w4 w+ [1 O. _6 x& Q: V241.45
* w% A& ~8 v# u- m' U$ }+ L251.03
4 A$ N& ~  F1 A! K: r" W5 h]
我在窗口输入：gm(x0) 为啥会出现
Error in ==> gm at 22
) g8 o* G$ P7 G; I/ t1 t; E; aerror(i)=abs(yc(i)-x0(i));                     %计算残差值
3 J- [+ x, j% V0 P那位高手帮忙 谢谢

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉