急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
function f=gm(x0,m)              %定义为函数gm(x)
n=length(x0);                     
x1=zeros(1,n);                    
x1(1)=x0(1);                  
8 I! F+ V7 I& P" u; h, F( p- n) Tfor i=2:n                         %计算累加序列x1
    x1(i)=x1(i-1)+x0(i);         
0 n1 Y: I# n$ J7 [3 u. Fend
9 N2 ^+ g$ j9 e. h1 [% x8 A% U' n" Ci=2:n;                            %对原始数列平行移位并赋值给y
y(i-1)=x0(i);                     
5 ?/ R/ [: J4 R, _9 H, By=y';                                 
( d% l: O: Y( Y- [1 I$ Mi=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为弧度
3 R5 A( z( V, ~7 v2 W% u# s                                            
% T! p- ~) a5 s' T* ~* {' ji=1:m;                             %计算预测累加数列的值
7 {+ {& r/ w$ m" {5 H: |ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
" ^! v& x9 B* syc(1)=ago(1);                           
: L5 b& M) S8 Fi=1:m-1;                             %还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
* ~+ F! |9 m9 G& Y1 U( \; pi=2:n;                              
error(i)=abs(yc(i)-x0(i));                     %计算残差值
yc(1)=ago(1);                       
i=1:m-1;                              %修正还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
; H/ }) W  T: M/ L! pc=std(error)/std(x0);                      %计算后验差比 *2
& {+ M, \* t8 W' e/ Nrelerror=abs((error)-mean(error)*ones(size(error)));
! A. n$ w" ^! _9 g[nrow,ncol]=size(relerror);                       
4 W/ p% w% j2 @p=0;
, {+ U5 R# l& s6 m# F: b6 wfor i=2:ncol                                    
0 ^& V! H, F$ {/ o2 x2 F& Q% E    if  relerror(1,i)<0.6745*std(x0)               
! f7 u8 M6 U& X# u. P9 r. f        p=p+1;                                 
    end
end
p=p/(n-1);                     
" E& I- Q  a2 K/ R6 S# y0 u1 `w1=min(abs(error));                              
6 G+ [0 \6 k& z( B. Q6 G+ ew2=max(abs(error));
i=1:n;                                     %计算关联度
  ^, I/ O4 k, d- yw(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
% k7 s6 R3 ^& u$ ]1 cw=sum(w)/(n-1);
au                                         %输出参数a,u的值
ago;                                        %输出累加数列ago的值
x0;                                         %输出原始序列值
$ o8 Z+ R6 z7 K2 B7 B' _3 q* ff=yc;                                        %输出预测的值
error;                                        %输出残差的值
& w* k' |7 k$ Q4 a$ rc;                                           %输出后验差比的值
p;                                           %输出小误差概率的值
w                                           %输出关联度
- O' {4 U( }" k! Y3 X. L5 @end                                   
x0=[31.26
32.09
33.39
35.45
40.52
. h+ i5 g1 U+ D2 l( c7 O43.52
47.86
51.45
55.6
$ z# q' H( f% D  P& V5 M60.14
! z1 M* z# }7 B$ o/ W) _: @5 s4 s64.82
68.65
8 U8 b) O0 g# u4 |) ?& v# T73.22
80.22
87.69
. U% Q; ]. I  _# ?. k/ Z93.97
. _5 A: r6 x, p% U. X4 T99.16
103.38
3 N! ^7 A# k% H  n5 l6 _" b109.46
114.6
( c- B/ P7 R1 C5 o7 g5 |7 X119.85
+ C$ e' `1 T; R' P/ h5 w124.92
$ d7 D2 b+ ^- w4 K( R+ s132.04
139.45
150.93
165.13
2 V# h  p/ H( P! v6 S: V181.93
! t: H8 }2 R8 ]/ v- Z196.83
212.38
228.07
1 {6 D6 C2 e$ y1 q. m. U! h4 e241.45
2 c  ]& ~, G& u! ^251.03
]
我在窗口输入：gm(x0) 为啥会出现
1 K) M$ e! ^* Q2 gError in ==> gm at 22
/ x. t( G' W" }2 L# e$ y; w: p! `error(i)=abs(yc(i)-x0(i));                     %计算残差值
. o1 [: p! G4 b$ y4 {/ _$ l5 o那位高手帮忙 谢谢

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉