急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
function f=gm(x0,m)              %定义为函数gm(x)
: b7 u- A; {9 Nn=length(x0);                     
x1=zeros(1,n);                    
! p5 w; \* |$ F1 g5 E2 C: J- Qx1(1)=x0(1);                  
for i=2:n                         %计算累加序列x1
    x1(i)=x1(i-1)+x0(i);         
3 x0 X7 y; T) vend
8 o# _/ h# S8 hi=2:n;                            %对原始数列平行移位并赋值给y
y(i-1)=x0(i);                     
0 Y& V( `7 M, C( Iy=y';                                 
i=1:n-1;                             
; `; L& ]5 M9 |6 ec(i)=-0.5*(x1(i)+x1(i+1));                 c(i)=-0.5*(x1(i)+x1(i+1))
! D" ^" v2 N& Y' r/ r0 xB=[c' ones(n-1,1)];                        
) `% ~0 r+ E% t8 r. ~au=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
) o9 Z9 Y8 S/ O" d2 c                                            
; u/ J$ a8 k5 W/ g& |2 Z; a- Ji=1:m;                             %计算预测累加数列的值
: L$ t1 O2 b4 e) _0 T" j0 ~6 }ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
9 R" G* C" e/ {# X! lyc(1)=ago(1);                           
, f- j8 p/ A0 V( u9 Hi=1:m-1;                             %还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
i=2:n;                              
; D: ?: g! a+ f/ i* {# @$ Cerror(i)=abs(yc(i)-x0(i));                     %计算残差值
yc(1)=ago(1);                       
i=1:m-1;                              %修正还原数列的值
' W5 o, }: n& N/ ?& n  I/ n9 ^: ayc(i+1)=ago(i+1)-ago(i);              
c=std(error)/std(x0);                      %计算后验差比 *2
relerror=abs((error)-mean(error)*ones(size(error)));
: W4 u- l4 K' o3 Y- j[nrow,ncol]=size(relerror);                       
p=0;
for i=2:ncol                                    
) O$ \5 L3 i# G/ e( m6 x    if  relerror(1,i)<0.6745*std(x0)               
        p=p+1;                                 
$ o) @$ D6 z; Q9 P    end
+ t8 D* n! v! M5 M5 \end
5 P2 Z) u4 _1 K: t6 P6 Up=p/(n-1);                     
' t% N! p) W% `+ n9 n7 \, i; v, \w1=min(abs(error));                              
; ?; a* ]0 \, x$ cw2=max(abs(error));
, D, t1 h. P! g  t0 j6 U1 u# M' i, Xi=1:n;                                     %计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
5 t7 w1 J2 Y' ~6 j( Mw=sum(w)/(n-1);
au                                         %输出参数a,u的值
ago;                                        %输出累加数列ago的值
; {& Q8 E' E  a" {x0;                                         %输出原始序列值
$ p8 t7 F) b* o  x+ hf=yc;                                        %输出预测的值
error;                                        %输出残差的值
c;                                           %输出后验差比的值
p;                                           %输出小误差概率的值
w                                           %输出关联度
- a9 y8 i1 P  }( gend                                   
x0=[31.26
* n/ u" K( Q+ N% g/ A* M32.09
33.39
35.45
40.52
43.52
4 T( a3 W' A9 h% A; O- H8 ?47.86
, y2 h$ Z) ^! q& h) S- j- \8 u51.45
55.6
60.14
64.82
68.65
73.22
80.22
$ i7 @+ r. y( |" z* j87.69
0 Q* g+ i! C: J- T93.97
# R- q9 t+ Y; V99.16
103.38
109.46
$ G* [  V+ \* ?* W) e9 T114.6
119.85
+ Q8 C& W% [5 a1 X0 `! h2 g124.92
2 [9 c4 a! ~. A3 A" N132.04
139.45
% q; i6 g9 _6 O0 I$ P; z150.93
$ i+ a9 y$ L5 y7 p; H% B8 q$ j165.13
1 F3 G" W: W5 r9 M% Z: \# K4 f) r181.93
196.83
3 Y% i1 T' |) _212.38
228.07
; Q/ v! Q0 T. Z2 {: W% p241.45
) e4 t9 S( a5 L7 s0 i251.03
- z' G: K& Z/ Z' O! s' k$ v6 c/ K/ |]
我在窗口输入：gm(x0) 为啥会出现
6 p+ w0 K7 f* _; u9 p5 s4 p& fError in ==> gm at 22
1 h7 _# V4 j, P# I, |error(i)=abs(yc(i)-x0(i));                     %计算残差值
! d( o& i( t$ Z* U5 |- W( v那位高手帮忙 谢谢
' b6 ]! m. E5 U9 i8 s2 @

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉