急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
function f=gm(x0,m)              %定义为函数gm(x)
n=length(x0);                     
x1=zeros(1,n);                    
x1(1)=x0(1);                  
for i=2:n                         %计算累加序列x1
8 i) D/ l2 |" i' R3 v! y4 L    x1(i)=x1(i-1)+x0(i);         
end
i=2:n;                            %对原始数列平行移位并赋值给y
% Q$ T, l( H! v7 H4 qy(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))
4 V. W0 A/ {8 j7 }" }, YB=[c' ones(n-1,1)];                        
5 `% B" d) g& s" _) v& jau=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
  f0 Y/ z- H! l7 [7 L* q( Q6 |6 U9 p                                            
" v& W( `; z7 {i=1:m;                             %计算预测累加数列的值
ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
4 ?6 ?& ~/ I: v# z+ m& f1 M: g3 gyc(1)=ago(1);                           
7 }% m. j& j$ ^1 h, L+ Z1 yi=1:m-1;                             %还原数列的值
5 |0 j1 c3 T/ v! oyc(i+1)=ago(i+1)-ago(i);              
i=2:n;                              
error(i)=abs(yc(i)-x0(i));                     %计算残差值
yc(1)=ago(1);                       
3 Q# {& q2 F# ii=1:m-1;                              %修正还原数列的值
% S9 P& M" y" O% s2 d$ Vyc(i+1)=ago(i+1)-ago(i);              
c=std(error)/std(x0);                      %计算后验差比 *2
, D* B, O- a3 y$ s- Drelerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);                       
1 Y4 _: E; |# ?p=0;
for i=2:ncol                                    
    if  relerror(1,i)<0.6745*std(x0)               
        p=p+1;                                 
    end
end
5 S( m% ?4 Y+ a$ W- x+ kp=p/(n-1);                     
2 E- [4 s8 H( [" p. `! l  \w1=min(abs(error));                              
* r" }# Q: y* L# k( Dw2=max(abs(error));
i=1:n;                                     %计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
% }. e2 B4 g$ k& Yw=sum(w)/(n-1);
  y! @7 P3 V# K0 v7 fau                                         %输出参数a,u的值
ago;                                        %输出累加数列ago的值
6 i% R9 v, ^0 Kx0;                                         %输出原始序列值
f=yc;                                        %输出预测的值
% J# X. p4 ~- z3 O* kerror;                                        %输出残差的值
c;                                           %输出后验差比的值
p;                                           %输出小误差概率的值
1 v6 ^3 [: P" f, D5 ~) Qw                                           %输出关联度
end                                   
0 G7 |- }1 _- t, ?, dx0=[31.26
# r3 x; S: L! H+ ^32.09
3 S$ f( b) _% U3 U33.39
0 e! d% L& L% K35.45
% o; w3 Z  C- a6 c9 I% O40.52
$ ^3 w) g0 ?1 [- L43.52
7 f9 q: w4 I" t  O2 J& Z$ w0 k3 I47.86
51.45
8 I3 |6 g, ?; a  Y55.6
- W# F) ^8 Z5 l4 w8 G60.14
64.82
/ j2 q' C/ k4 K" Y. V' e68.65
73.22
3 s6 l+ v7 U7 @( P' K0 ?* Z80.22
87.69
! o$ X9 r; L4 @& _1 m% D. H93.97
3 r8 J% p7 Z- _9 d  [$ V/ w: d# P99.16
' _  L( w/ i* T# H+ @  Q103.38
109.46
7 u: Z. x. M/ a2 G- c& k* w114.6
119.85
8 x9 A% P' l0 ]7 t5 f2 y1 i124.92
) C  Y/ I5 q% H) K132.04
' o( j9 b: Y% B139.45
150.93
165.13
181.93
196.83
212.38
6 Z1 N+ m: w! I6 w228.07
241.45
1 Y$ V" b, Y& v6 x; x( Z8 u8 h251.03
$ D# {, `' G9 K" P]
/ ~4 V5 ?/ S2 {2 t, N: b  J+ |我在窗口输入：gm(x0) 为啥会出现
Error in ==> gm at 22
error(i)=abs(yc(i)-x0(i));                     %计算残差值
那位高手帮忙 谢谢
! z0 q! v) a5 u) h: y$ J- h

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉