急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
8 A' s) A  f1 A; Yfunction f=gm(x0,m)              %定义为函数gm(x)
n=length(x0);                     
x1=zeros(1,n);                    
x1(1)=x0(1);                  
2 u6 P$ k( p5 I4 ^for i=2:n                         %计算累加序列x1
    x1(i)=x1(i-1)+x0(i);         
% V% b  t) a+ |5 V! ]0 _end
i=2:n;                            %对原始数列平行移位并赋值给y
1 h( r$ r& W+ h  Iy(i-1)=x0(i);                     
y=y';                                 
9 E; p  S7 a" w; z- C( w. e9 xi=1:n-1;                             
+ e, B2 r. B# L; g) X, dc(i)=-0.5*(x1(i)+x1(i+1));                 c(i)=-0.5*(x1(i)+x1(i+1))
B=[c' ones(n-1,1)];                        
* V7 _6 ?7 h: Uau=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)   
; Z/ J6 k& P0 H) Q# [yc(1)=ago(1);                           
i=1:m-1;                             %还原数列的值
# i, H# ]: x& p2 {yc(i+1)=ago(i+1)-ago(i);              
% P! ?7 m! S7 |2 A: |i=2:n;                              
4 w4 M  _. w4 n5 {# Aerror(i)=abs(yc(i)-x0(i));                     %计算残差值
yc(1)=ago(1);                       
i=1:m-1;                              %修正还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
c=std(error)/std(x0);                      %计算后验差比 *2
, `& b  C6 f! O5 Hrelerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);                       
1 i7 _) ?! n$ x7 X: Dp=0;
' O- a$ t6 p% S2 Nfor i=2:ncol                                    
( A+ D3 R- L" l; p% v7 ]    if  relerror(1,i)<0.6745*std(x0)               
        p=p+1;                                 
: k$ y4 D5 e. X, b. M3 M7 [7 z. p    end
end
p=p/(n-1);                     
w1=min(abs(error));                              
w2=max(abs(error));
i=1:n;                                     %计算关联度
' B) J7 L: Q. R2 C. Pw(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
& u. _3 ]. S( ^# W  ~0 m7 }w=sum(w)/(n-1);
' p6 u! u5 t1 r. T0 |au                                         %输出参数a,u的值
ago;                                        %输出累加数列ago的值
x0;                                         %输出原始序列值
' {% x& d: }% a- G1 c9 z* }f=yc;                                        %输出预测的值
4 |3 b4 O6 S' a* M2 l# K# ^error;                                        %输出残差的值
1 d4 r4 a0 d9 o3 P4 z3 dc;                                           %输出后验差比的值
# Q+ B+ e9 X% g- i+ m2 yp;                                           %输出小误差概率的值
* G) ]$ V5 V6 l, e, q7 _+ ]' x+ }w                                           %输出关联度
end                                   
/ R; o7 O# `4 L+ L& ^* yx0=[31.26
32.09
! L7 d4 |$ C5 z1 J# Y, J6 D33.39
% M; o9 G; d7 L# Y35.45
0 Q; W* @! c' X40.52
43.52
- {9 M& s& }8 n47.86
51.45
4 e: X3 b, @+ j6 m4 n55.6
* y( [/ F- D) q60.14
64.82
. D* n% E' F% F: I5 \68.65
7 p9 P4 b% D! z& W0 l, m: D73.22
80.22
! I6 \7 u& ~! J0 w% N, l3 a( Q87.69
! K  d' h% m) O0 N3 M93.97
. G! t  V0 i, D9 V0 c99.16
  O9 B( _2 |0 ^9 P1 v5 u+ H) D& h% i$ o103.38
1 g9 w* Y4 a/ P7 R" Z+ d109.46
& I+ P0 f3 Z' y* \3 s! u, b' z114.6
119.85
124.92
+ q, i( t$ }5 A; Y, |132.04
139.45
# V6 b' Z) ]/ B! y, O150.93
9 }# `: J1 g! a5 \: ~+ z7 l' X165.13
181.93
6 J# R9 i0 P- T1 f1 |) U5 s196.83
212.38
6 `! F1 c2 W# D8 x228.07
241.45
251.03
]
我在窗口输入：gm(x0) 为啥会出现
9 e" ?4 L* O5 m) @/ J/ _3 Z! a/ k" YError in ==> gm at 22
error(i)=abs(yc(i)-x0(i));                     %计算残差值
那位高手帮忙 谢谢

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉