急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
function f=gm(x0,m)              %定义为函数gm(x)
n=length(x0);                     
x1=zeros(1,n);                    
- t5 s: \" B- L' j( W! Fx1(1)=x0(1);                  
for i=2:n                         %计算累加序列x1
: s* W4 H! a& r9 {& ]0 a( ?    x1(i)=x1(i-1)+x0(i);         
" r7 O: c# O+ J9 [end
i=2:n;                            %对原始数列平行移位并赋值给y
y(i-1)=x0(i);                     
y=y';                                 
* Q) |  `3 @1 ]& di=1:n-1;                             
c(i)=-0.5*(x1(i)+x1(i+1));                 c(i)=-0.5*(x1(i)+x1(i+1))
$ {- q* R4 [6 g% q8 V& f3 u: M5 iB=[c' ones(n-1,1)];                        
) s- R8 J& f# l6 _au=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
! V* ^# ?- [3 c- P. K                                            
i=1:m;                             %计算预测累加数列的值
ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
' C$ D" k$ g9 u+ hyc(1)=ago(1);                           
& t% `  v/ t2 G2 N$ ui=1:m-1;                             %还原数列的值
% ]$ B: `. b: V2 Ayc(i+1)=ago(i+1)-ago(i);              
i=2:n;                              
8 u3 ^3 g! z+ x" h8 Werror(i)=abs(yc(i)-x0(i));                     %计算残差值
yc(1)=ago(1);                       
; S0 q5 a( e- Y8 p$ M8 j+ {- Di=1:m-1;                              %修正还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
c=std(error)/std(x0);                      %计算后验差比 *2
4 a6 X$ n! Q& ?; g! Grelerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);                       
p=0;
( y) G  e; J1 U6 z- Qfor i=2:ncol                                    
    if  relerror(1,i)<0.6745*std(x0)               
        p=p+1;                                 
    end
end
p=p/(n-1);                     
: B' p- D' w' S: W8 hw1=min(abs(error));                              
6 _! q9 L- O* J& r1 a" ew2=max(abs(error));
i=1:n;                                     %计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
6 O; T' d: l8 [9 t) uw=sum(w)/(n-1);
au                                         %输出参数a,u的值
. W; M1 }1 {! _ago;                                        %输出累加数列ago的值
6 d2 k& G; E' U" U; ^$ i5 |7 Hx0;                                         %输出原始序列值
' \1 Y7 T' \/ a& rf=yc;                                        %输出预测的值
6 Q! w3 l0 f) serror;                                        %输出残差的值
c;                                           %输出后验差比的值
" a- {# q* n8 l6 S& V6 \p;                                           %输出小误差概率的值
% `4 r+ y& J$ K8 b8 L$ W( [$ X: yw                                           %输出关联度
end                                   
x0=[31.26
32.09
" O; N' F4 Y+ h8 I3 u33.39
/ M# d5 d. Z4 ?6 p9 G3 P( u35.45
6 T% F% A0 D6 ~. R; v40.52
, X  q2 l, ^: N43.52
! F- R& Z4 Q! q6 |47.86
51.45
55.6
60.14
* d1 Q5 I  d3 T$ Z8 ^64.82
! o( _  M2 v$ W6 l68.65
73.22
" u' P& D" y8 Q/ I3 s" t80.22
7 z" O: g$ ~2 b9 l4 h87.69
/ h. d2 l0 c" A& D3 M+ ]1 z) n93.97
. W0 z: e- T; w7 {99.16
& e0 r( e  K) k9 w/ i5 B* z5 z103.38
109.46
- t) \+ U; [; X, \114.6
119.85
124.92
: x2 l/ E) ?6 B. F6 k132.04
139.45
) A$ f1 K+ j. I1 s# Z* t150.93
165.13
2 [. ?, L8 z: I" I181.93
: K* u+ f! m" H0 x3 E" m1 v196.83
212.38
9 m9 R% ~5 [8 B& z. C) n* v228.07
241.45
4 Z# u$ j; m3 L7 ]( E0 |% J251.03
]
+ _( l8 F% {* }+ j我在窗口输入：gm(x0) 为啥会出现
Error in ==> gm at 22
7 s7 u# l5 }2 z: Herror(i)=abs(yc(i)-x0(i));                     %计算残差值
那位高手帮忙 谢谢

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉