急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
function f=gm(x0,m)              %定义为函数gm(x)
n=length(x0);                     
: V; I, R/ I% x1 ]; M: t, |x1=zeros(1,n);                    
x1(1)=x0(1);                  
for i=2:n                         %计算累加序列x1
    x1(i)=x1(i-1)+x0(i);         
end
i=2:n;                            %对原始数列平行移位并赋值给y
y(i-1)=x0(i);                     
y=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)];                        
4 [1 c1 T9 ]; B0 c$ w% k1 T1 uau=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
$ A: x! \; }/ ?0 F+ y                                            
i=1:m;                             %计算预测累加数列的值
ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
4 v3 m( V4 k! y! Pyc(1)=ago(1);                           
3 Q5 \; ^( `! k4 si=1:m-1;                             %还原数列的值
1 ?1 x( M' ^# |* s$ vyc(i+1)=ago(i+1)-ago(i);              
i=2:n;                              
error(i)=abs(yc(i)-x0(i));                     %计算残差值
5 S+ s: m) P7 p7 Tyc(1)=ago(1);                       
1 q/ ]- x6 Y% mi=1:m-1;                              %修正还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
  Q. X% M) M" m, s) P9 kc=std(error)/std(x0);                      %计算后验差比 *2
relerror=abs((error)-mean(error)*ones(size(error)));
& x( W) R% V# s' P[nrow,ncol]=size(relerror);                       
p=0;
for i=2:ncol                                    
$ Y: P( W$ P! A1 g    if  relerror(1,i)<0.6745*std(x0)               
        p=p+1;                                 
    end
; u3 k- p. l& A, _9 c. R1 o6 Dend
) K, j4 x4 ]" @4 Hp=p/(n-1);                     
w1=min(abs(error));                              
w2=max(abs(error));
i=1:n;                                     %计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
- ~$ M: G" E" bw=sum(w)/(n-1);
3 ]0 j8 ~3 C1 z4 wau                                         %输出参数a,u的值
  l- h3 y8 R" F& I' V, c" l; |ago;                                        %输出累加数列ago的值
x0;                                         %输出原始序列值
f=yc;                                        %输出预测的值
error;                                        %输出残差的值
3 O) I1 T! |- o! h/ ^0 jc;                                           %输出后验差比的值
p;                                           %输出小误差概率的值
. I' f' \1 }6 n# _( }w                                           %输出关联度
- j$ }/ D+ F- b, O# Pend                                   
: C: o! _* p6 bx0=[31.26
. l6 a) z/ Y- `2 X& b  k( c32.09
$ f" o2 ~% A2 W! P* u33.39
35.45
' U4 Y! K3 C* ^* ]: A40.52
43.52
. V! R8 h& w9 {: O47.86
51.45
' o; U' P4 ^1 x/ {5 O55.6
+ ^  ?2 ?- V' t& w8 n60.14
0 K2 w" z2 G5 M8 ?64.82
' F- J& K: W2 W6 O68.65
3 h: W4 K. r% T) x: ^3 A  g73.22
3 D, k3 o  w+ I! d, w) Z80.22
87.69
. v$ r/ p9 V) ?4 E; [; Q93.97
99.16
103.38
" e7 p: Q3 b4 ?- Y. ?1 o- f109.46
114.6
119.85
; \$ P' R& H. X3 \: j3 y/ a  g- q124.92
132.04
139.45
, [* Z3 j) \1 E  `, o150.93
7 L4 ^$ u9 g1 o$ C/ V. {165.13
181.93
; T( e) {' B2 w- a; A8 [0 _6 N! \196.83
212.38
228.07
241.45
+ v& j& v; F3 \( B251.03
]
2 k/ R7 b5 c# z1 \# |5 [: k6 L! F  T我在窗口输入：gm(x0) 为啥会出现
, ?$ {# A7 d  J1 e. i: AError in ==> gm at 22
, k" \6 L5 w: \7 }: }% Z# m% Ierror(i)=abs(yc(i)-x0(i));                     %计算残差值
那位高手帮忙 谢谢
5 j# d' H! j4 K6 I" H9 l0 r4 s

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉