急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
; R- e! v* U$ v0 wfunction f=gm(x0,m)              %定义为函数gm(x)
8 I' t  v* Y6 F/ n" R. Rn=length(x0);                     
x1=zeros(1,n);                    
x1(1)=x0(1);                  
" r1 q& D' p2 D  lfor i=2:n                         %计算累加序列x1
7 O6 g/ n) S' o. |3 k    x1(i)=x1(i-1)+x0(i);         
& U% |* M4 O. E' ~- e. S1 dend
0 ~8 G2 c, V7 u& ki=2:n;                            %对原始数列平行移位并赋值给y
3 H5 Y* L9 N: }, q! [; ny(i-1)=x0(i);                     
y=y';                                 
/ V( J6 E! h2 O* `6 G2 B7 }i=1:n-1;                             
" \3 O6 U1 ]  ]$ bc(i)=-0.5*(x1(i)+x1(i+1));                 c(i)=-0.5*(x1(i)+x1(i+1))
B=[c' ones(n-1,1)];                        
3 X6 ?- a4 ]% H2 \; [au=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
                                            
i=1:m;                             %计算预测累加数列的值
ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
  |9 |. T  Z5 I1 `! [" ayc(1)=ago(1);                           
i=1:m-1;                             %还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
; C9 c, d4 M7 Q& Ni=2:n;                              
error(i)=abs(yc(i)-x0(i));                     %计算残差值
4 O2 d8 ]: H" w3 w7 C0 d: Dyc(1)=ago(1);                       
! l' x( D0 R' ?, c* O0 {i=1:m-1;                              %修正还原数列的值
0 U# I$ S3 U/ L: v5 J2 cyc(i+1)=ago(i+1)-ago(i);              
c=std(error)/std(x0);                      %计算后验差比 *2
relerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);                       
p=0;
for i=2:ncol                                    
  I+ P$ g& }4 B) R: a0 H. o    if  relerror(1,i)<0.6745*std(x0)               
  A8 F6 l8 Q2 n! Q/ s        p=p+1;                                 
& t- f, B& I% h# A* V5 w) l    end
end
p=p/(n-1);                     
9 R% }0 p' N$ D: W+ X1 }w1=min(abs(error));                              
- L3 @0 h5 S5 U* t* Q; Y! }w2=max(abs(error));
3 v  Z" T/ m3 n8 Di=1:n;                                     %计算关联度
$ a0 f  h& ^; j$ H* rw(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
w=sum(w)/(n-1);
9 T& e& E! {! W3 V6 ~# ?au                                         %输出参数a,u的值
- |; o' D2 }; y$ t; {: x5 J- K' C+ Dago;                                        %输出累加数列ago的值
( A- L: e6 u9 e* h! Ex0;                                         %输出原始序列值
3 S4 p/ q3 M* Of=yc;                                        %输出预测的值
error;                                        %输出残差的值
c;                                           %输出后验差比的值
p;                                           %输出小误差概率的值
/ U* S+ B- N" e- m3 g% R- w5 Y+ D) [7 Ew                                           %输出关联度
5 p3 M! k6 y+ I& W9 o* S5 Qend                                   
x0=[31.26
32.09
33.39
2 @) Z  G3 P: w8 V+ g- v35.45
* h4 W. G0 i0 ?% x40.52
43.52
7 g( U3 S5 }) \/ {+ a47.86
51.45
% J4 t' W3 G1 Z, e55.6
2 K$ o" z% T+ ?60.14
64.82
' T1 s& _% ~! E2 e0 d! ^5 ^2 l9 A68.65
5 h" I" n& B% _; |73.22
. B4 N6 B: O2 b& v, _80.22
3 \, @; s9 k% c( H/ l( O, C87.69
( {$ W/ y0 Z) H1 @93.97
99.16
$ b' G* X9 `8 k8 Y2 Q9 p103.38
" `" U6 i1 f8 ^1 T8 A" W109.46
114.6
1 {8 L# C/ B* k6 l+ g- H7 a# u5 n119.85
! k" t5 Z0 }3 V5 a124.92
: N5 D! `, s( m132.04
; i: o3 W0 Y4 {8 k) m2 Y6 W139.45
150.93
165.13
7 m6 S) o3 x* i181.93
( L+ i. h! S9 P# ]. G! Q196.83
212.38
228.07
, ^* \3 b( ~7 \9 N; e; N241.45
251.03
6 p0 C' z/ }3 s5 L& b]
1 e. u/ z% y8 w我在窗口输入：gm(x0) 为啥会出现
0 i, r  N5 d! V  J$ L* {; \Error in ==> gm at 22
error(i)=abs(yc(i)-x0(i));                     %计算残差值
! d% g% |7 [  g) X. \9 x那位高手帮忙 谢谢

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉