急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
+ r/ P1 e8 y( ]  n  Q4 `function f=gm(x0,m)              %定义为函数gm(x)
n=length(x0);                     
x1=zeros(1,n);                    
. R9 M: x; s: H1 H" Mx1(1)=x0(1);                  
- Z; z8 @& k+ W! r: W2 s0 Dfor i=2:n                         %计算累加序列x1
) I$ G- I/ S0 @0 i4 G    x1(i)=x1(i-1)+x0(i);         
end
i=2:n;                            %对原始数列平行移位并赋值给y
y(i-1)=x0(i);                     
y=y';                                 
i=1:n-1;                             
4 I- L* v0 ?1 P8 qc(i)=-0.5*(x1(i)+x1(i+1));                 c(i)=-0.5*(x1(i)+x1(i+1))
. X" _1 S, F! jB=[c' ones(n-1,1)];                        
au=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
0 u' o% U" f& `+ r% [- n+ t                                            
- H* B! g) X  Di=1:m;                             %计算预测累加数列的值
ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
" D/ |  H3 P1 T, K( ~) X5 Eyc(1)=ago(1);                           
6 l' O6 [* N  I4 a! D2 Xi=1:m-1;                             %还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
i=2:n;                              
error(i)=abs(yc(i)-x0(i));                     %计算残差值
yc(1)=ago(1);                       
+ ?1 [7 T$ S5 o& x" J/ E* Ii=1:m-1;                              %修正还原数列的值
3 I5 A5 q! n  d* A) `) Xyc(i+1)=ago(i+1)-ago(i);              
' i! ~& i1 `9 _: }+ S$ S1 M  hc=std(error)/std(x0);                      %计算后验差比 *2
' A- M; e  M* x5 F; ~relerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);                       
p=0;
for i=2:ncol                                    
/ X1 W+ r' I6 u4 I! z8 p    if  relerror(1,i)<0.6745*std(x0)               
        p=p+1;                                 
! d5 }9 T. i5 Q) T+ D    end
( D3 x- F' c6 A# Zend
# z9 d6 a6 R0 C0 H) pp=p/(n-1);                     
$ u% \# N' R% L) I* d* P; Mw1=min(abs(error));                              
7 J. T% t$ z, T+ H4 u* ^( ]* a$ mw2=max(abs(error));
i=1:n;                                     %计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
w=sum(w)/(n-1);
au                                         %输出参数a,u的值
ago;                                        %输出累加数列ago的值
x0;                                         %输出原始序列值
, v% W- `4 ~( k2 A! Ef=yc;                                        %输出预测的值
error;                                        %输出残差的值
c;                                           %输出后验差比的值
: P# }# B8 y+ [* D! `p;                                           %输出小误差概率的值
w                                           %输出关联度
8 ]) @( m& w0 f' `  k5 E6 dend                                   
x0=[31.26
32.09
- |) L: C# l# I) t# C33.39
35.45
7 t9 i  x+ ]; y* ]$ O; v4 O/ v40.52
43.52
47.86
8 a) ]2 ?% g8 F5 ?51.45
55.6
60.14
7 H9 q4 d4 f, P, {3 w64.82
68.65
0 l+ }/ |0 N3 m  ?% ~2 M4 |73.22
80.22
, w" l, c- L6 o; E5 `87.69
2 x$ T8 \0 r$ r; |0 C" R0 X93.97
# b, `; U# L3 M1 h- {4 f99.16
103.38
0 J+ C- u1 q% h+ B109.46
( k2 H' ]2 r2 p114.6
/ w& b; {& M' v119.85
% ?% ~4 K7 G, o/ L9 s+ P124.92
132.04
* C) t. D, d. }+ V* I7 f* _. e7 c" l139.45
150.93
165.13
181.93
3 e9 m( i1 @# e3 }0 y196.83
' Y$ H+ F* |" R* S9 Z  H( S, o212.38
! y7 @' S2 ?4 }& F  X  Y$ D- W228.07
& H  }: v2 S1 k( @241.45
6 M/ J/ N- O6 |3 y- j251.03
]
我在窗口输入：gm(x0) 为啥会出现
Error in ==> gm at 22
error(i)=abs(yc(i)-x0(i));                     %计算残差值
, @" w6 M! ]* G# ?那位高手帮忙 谢谢

wuzhenhua

把abs（）去掉

神经病个人组

我也想知道