急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
function f=gm(x0,m)              %定义为函数gm(x)
n=length(x0);                     
, S3 N6 X* ]. G8 F3 sx1=zeros(1,n);                    
x1(1)=x0(1);                  
for i=2:n                         %计算累加序列x1
: I4 n* C2 E$ U9 u$ p. _$ h, l    x1(i)=x1(i-1)+x0(i);         
. \( D& @5 R% v; e9 a2 D8 \end
. U. }6 E9 ]; T" u1 Q5 T; mi=2:n;                            %对原始数列平行移位并赋值给y
y(i-1)=x0(i);                     
) d% V1 Q2 ?/ j$ Q: vy=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)];                        
% M5 U0 d5 q! m$ c: p0 Xau=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
; v( P/ d6 q4 P1 B. f/ V                                            
% D( {2 E" @9 d+ m* n: h  v( P1 qi=1:m;                             %计算预测累加数列的值
ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
yc(1)=ago(1);                           
$ o. L$ |/ T$ l9 i' [$ Li=1:m-1;                             %还原数列的值
, ?6 ~- w$ e1 N* ^$ b1 {- |+ byc(i+1)=ago(i+1)-ago(i);              
i=2:n;                              
error(i)=abs(yc(i)-x0(i));                     %计算残差值
" I  O* I% J* ?2 _  q$ k  H0 Kyc(1)=ago(1);                       
& b8 @. o. b( B- Pi=1:m-1;                              %修正还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
c=std(error)/std(x0);                      %计算后验差比 *2
# l4 \" l5 j5 [" {8 z) f/ Z& Hrelerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);                       
2 i7 j  ?/ v* U) u1 p. H% hp=0;
for i=2:ncol                                    
    if  relerror(1,i)<0.6745*std(x0)               
        p=p+1;                                 
    end
" O# s2 ]" G2 k4 T  R" {end
! J0 t8 b- f! n  d; Fp=p/(n-1);                     
w1=min(abs(error));                              
- H' V& |* U+ P' P! M& c) y/ {# ~w2=max(abs(error));
& C3 S( U/ h: S+ C! u0 Y1 Di=1:n;                                     %计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
w=sum(w)/(n-1);
- g: n* d# X( |, [) y( [' W7 j& ]au                                         %输出参数a,u的值
* D- a+ k- _+ E7 `  G$ I6 P9 p/ kago;                                        %输出累加数列ago的值
x0;                                         %输出原始序列值
7 d. `9 y& m" if=yc;                                        %输出预测的值
error;                                        %输出残差的值
c;                                           %输出后验差比的值
7 W+ g9 J" C+ j8 h: v2 Y( D8 P( tp;                                           %输出小误差概率的值
w                                           %输出关联度
; |. u* R, }. w9 nend                                   
! i1 W* E% w/ `x0=[31.26
32.09
33.39
8 C: K8 C6 u* h; @# r35.45
40.52
43.52
. e+ i- O7 e0 f47.86
$ v) z% C" ?6 g6 r5 ~9 Z51.45
( Z2 a3 l* A9 b- ]55.6
6 g) H7 C, s4 a& C* Z! B1 v3 G' H60.14
- ~( G3 p; x) [& X! g7 ?9 ?64.82
68.65
' t" x+ x8 F! A3 p& i6 R4 v! t73.22
- H8 B" _3 L  i+ F9 i0 d80.22
87.69
8 `# ^* U5 H% C0 @8 e! T93.97
9 @0 K& _+ O; o. e- @99.16
103.38
109.46
1 I" \4 B( Z# T; t114.6
; o8 ]( O; o  T' z( l- u119.85
1 W/ c  v7 ^  ?2 x6 _( m3 Z  r124.92
132.04
139.45
150.93
- `% H* @5 R! ?' U: n9 S* W, C! ^165.13
181.93
  K2 E/ K) l3 r! ~196.83
212.38
228.07
4 E. l2 {0 \! m! P2 }& }% ~- A( ^241.45
251.03
]
. x' k- ]1 E- R( O  D我在窗口输入：gm(x0) 为啥会出现
6 P3 _& a3 P+ n; Y6 yError in ==> gm at 22
( \9 Y) e' o/ e1 Verror(i)=abs(yc(i)-x0(i));                     %计算残差值
那位高手帮忙 谢谢
6 C! j" e) U; t7 W( M. ^3 M

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉