急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
function f=gm(x0,m)              %定义为函数gm(x)
n=length(x0);                     
x1=zeros(1,n);                    
: z6 t( {8 l  E/ _x1(1)=x0(1);                  
% V1 A& k9 U  r' C) u$ |, l. {for i=2:n                         %计算累加序列x1
    x1(i)=x1(i-1)+x0(i);         
, F* t0 ]0 H7 N" p5 m& |, mend
i=2:n;                            %对原始数列平行移位并赋值给y
y(i-1)=x0(i);                     
* H# g) k* T1 I  Fy=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)];                        
. o* N! p  B5 H% r. R2 tau=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
1 h3 ^' @0 }: L) b6 v3 n                                            
3 g* S5 C! @; M/ }; t8 n/ A% ci=1:m;                             %计算预测累加数列的值
ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
: q' \4 d7 S& P4 ?3 c/ Y$ z) myc(1)=ago(1);                           
: \3 w  m/ H' e" y$ [8 Ni=1:m-1;                             %还原数列的值
: O4 i' z6 T) J* y: D4 C  Z4 Xyc(i+1)=ago(i+1)-ago(i);              
i=2:n;                              
- O1 J0 |: ]3 H% h1 jerror(i)=abs(yc(i)-x0(i));                     %计算残差值
7 S! G. e( s$ y6 Wyc(1)=ago(1);                       
i=1:m-1;                              %修正还原数列的值
* n0 Y" Y; C7 g. R4 D6 Zyc(i+1)=ago(i+1)-ago(i);              
/ z, }9 s0 g5 l3 t6 `c=std(error)/std(x0);                      %计算后验差比 *2
) `+ |8 y' m! v" I! jrelerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);                       
' }  T! s2 [+ \p=0;
5 j& w' L0 j7 n% L- Tfor i=2:ncol                                    
2 I4 [/ S! k  ~; J+ @7 g% M    if  relerror(1,i)<0.6745*std(x0)               
        p=p+1;                                 
    end
$ k' U9 Z) K) e. qend
p=p/(n-1);                     
w1=min(abs(error));                              
8 E7 |! Y& x8 B9 A1 Fw2=max(abs(error));
i=1:n;                                     %计算关联度
, K7 a! I% I- ew(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
& Z2 n5 M- {7 X. g4 b; j* Dw=sum(w)/(n-1);
au                                         %输出参数a,u的值
1 R  \8 l! H3 R# Uago;                                        %输出累加数列ago的值
x0;                                         %输出原始序列值
; N+ @( \0 M- j" ~f=yc;                                        %输出预测的值
' ?: ]  l) U+ }: Ierror;                                        %输出残差的值
1 {4 A" c5 O! O# {. \2 t0 L  hc;                                           %输出后验差比的值
# _7 |) S0 ]' t5 jp;                                           %输出小误差概率的值
" e! _: |# e3 T  y: G4 M9 dw                                           %输出关联度
& \$ G2 z- {9 A% o' Rend                                   
% z1 B. w8 S8 |* W: V- f, ^x0=[31.26
32.09
33.39
35.45
8 [2 X% R8 o- j40.52
43.52
47.86
51.45
7 {1 [- ~- F1 D+ Y+ a- I" H; W4 W9 j; g55.6
2 ~( Z+ j) e: l7 F" h0 W60.14
/ v2 D1 }- h  |64.82
. T* ]% M& O, y68.65
73.22
3 s! G3 ~" K; K6 p3 v2 R+ G80.22
! h4 m: m- w$ G! ]( n" {( b3 ?87.69
93.97
99.16
& `9 m& m( u7 g0 C& a$ U" O103.38
. c+ a" O9 _, U! Q8 {& D6 W* h109.46
  e% w* u" ^' y% Q114.6
119.85
9 v* a% o$ M6 Y" P: f2 r1 ?124.92
; y1 J0 l, E1 M132.04
139.45
150.93
, u4 c8 x- Q- Z  o9 }! A1 a1 M165.13
) ^) d9 P1 Q2 }5 \3 G181.93
& C$ B! K1 n' L! T- m) C% `8 U. ~' I196.83
9 @! [% B( a+ P2 i) Z212.38
228.07
; {7 ^: {. P3 O1 i1 \241.45
251.03
" @  e. Q" m( ]]
我在窗口输入：gm(x0) 为啥会出现
Error in ==> gm at 22
% k! i3 ]! Q: \$ ]: _error(i)=abs(yc(i)-x0(i));                     %计算残差值
$ N/ C2 u; T4 Z" A  B& z" \& N那位高手帮忙 谢谢
( W, Z  C$ t% c/ M7 R1 ~0 K

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉