急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
function f=gm(x0,m)              %定义为函数gm(x)
2 i7 f+ K2 q/ o; A: I+ |2 dn=length(x0);                     
x1=zeros(1,n);                    
x1(1)=x0(1);                  
. [6 Y9 F" F$ Gfor i=2:n                         %计算累加序列x1
    x1(i)=x1(i-1)+x0(i);         
end
3 q) b5 ]9 t- l6 _1 Z+ u' m0 v4 pi=2:n;                            %对原始数列平行移位并赋值给y
" y7 ~6 t+ F/ ^4 _! l7 xy(i-1)=x0(i);                     
y=y';                                 
8 \& X) {) P( N( xi=1:n-1;                             
" o- j/ G- s1 p0 v6 A) M' qc(i)=-0.5*(x1(i)+x1(i+1));                 c(i)=-0.5*(x1(i)+x1(i+1))
B=[c' ones(n-1,1)];                        
au=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
                                            
' r$ h7 p$ E& m' z/ X7 q, Gi=1:m;                             %计算预测累加数列的值
  E$ t3 n& e* B# |# p) N0 N  Iago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
yc(1)=ago(1);                           
7 A7 q' I: L3 h6 k8 e2 D+ R6 Mi=1:m-1;                             %还原数列的值
+ w- O6 L! i$ @! ]yc(i+1)=ago(i+1)-ago(i);              
6 {1 y7 ]  b/ a( [* o6 li=2:n;                              
error(i)=abs(yc(i)-x0(i));                     %计算残差值
yc(1)=ago(1);                       
2 C3 d" z# k8 U+ Z+ D& C, I+ C" ui=1:m-1;                              %修正还原数列的值
' d; w- B6 \" Q" h: {9 S; v0 Z: J" vyc(i+1)=ago(i+1)-ago(i);              
+ }, e" P! z; U5 H0 @7 ac=std(error)/std(x0);                      %计算后验差比 *2
relerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);                       
) M; x6 T) Y2 O" dp=0;
( u: j& O' I# k( g/ xfor i=2:ncol                                    
! ^+ H& ^4 \# Y  N: G. V' l    if  relerror(1,i)<0.6745*std(x0)               
5 S4 ~9 `, F' e7 u+ r        p=p+1;                                 
7 T3 C% h+ q9 U6 q3 p  `    end
end
5 G1 I2 x5 p$ ]' q5 l; gp=p/(n-1);                     
) P* ^9 \& z9 }! y. Jw1=min(abs(error));                              
w2=max(abs(error));
i=1:n;                                     %计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
w=sum(w)/(n-1);
$ _' K3 W  s6 ^8 y1 Aau                                         %输出参数a,u的值
7 W1 ?7 y2 |9 c  ~- u  _' X9 s% Gago;                                        %输出累加数列ago的值
+ L4 ]4 n1 N6 [$ vx0;                                         %输出原始序列值
0 g7 _# Q6 _" i% e2 o) ^" Af=yc;                                        %输出预测的值
error;                                        %输出残差的值
c;                                           %输出后验差比的值
p;                                           %输出小误差概率的值
w                                           %输出关联度
end                                   
x0=[31.26
32.09
0 v: o+ B( E" Z1 [0 ^33.39
1 O1 F2 P& M6 h$ g4 D35.45
40.52
43.52
47.86
: }: X3 o: ^( }  v3 F" J; p1 T51.45
55.6
! a& c6 H. Q  R; }4 Y( m0 w8 D6 Y60.14
. g) T* Q0 Y2 N) j6 W( v64.82
# k5 R$ D/ `# h2 \+ S" E68.65
6 A1 [0 E) r" E73.22
5 G7 r  g) R8 ]& z$ t# s' ^- }80.22
87.69
93.97
, [1 L+ d0 x9 m+ e$ I1 [99.16
103.38
- T4 V6 Y0 C, Y4 Y0 Y109.46
114.6
119.85
124.92
132.04
139.45
150.93
# b% M6 l, o2 f' P& y165.13
181.93
196.83
1 z1 p" y/ h1 x1 [! Z' P8 X" k+ I212.38
' P& O: z& w0 C8 o, u' T. D) Q$ ~228.07
; {% E& E! F" N5 S" P' A241.45
251.03
# u- P8 }" M; I' y/ @5 a7 f]
我在窗口输入：gm(x0) 为啥会出现
Error in ==> gm at 22
error(i)=abs(yc(i)-x0(i));                     %计算残差值
那位高手帮忙 谢谢

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉