急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
function f=gm(x0,m)              %定义为函数gm(x)
8 E* C( r/ v8 v8 E# F& vn=length(x0);                     
2 `$ ~7 u% j, _/ \x1=zeros(1,n);                    
" U) Z; N5 ]  d; h: w& Z9 u4 Y- [x1(1)=x0(1);                  
for i=2:n                         %计算累加序列x1
    x1(i)=x1(i-1)+x0(i);         
end
% K: b4 g: d/ M! I8 r- Wi=2:n;                            %对原始数列平行移位并赋值给y
y(i-1)=x0(i);                     
( z4 f8 d* P, V: G8 O* Qy=y';                                 
5 s& r* g, p+ U0 D& i' Vi=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)];                        
au=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
! W8 c6 b# u  {/ w                                            
i=1:m;                             %计算预测累加数列的值
ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
yc(1)=ago(1);                           
9 K' k2 \5 L. O3 |3 R, \i=1:m-1;                             %还原数列的值
* B/ c5 b# C* T8 wyc(i+1)=ago(i+1)-ago(i);              
i=2:n;                              
* n5 g# k: I  n4 x- [# Yerror(i)=abs(yc(i)-x0(i));                     %计算残差值
) z/ W0 V$ y7 j( uyc(1)=ago(1);                       
i=1:m-1;                              %修正还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
: w5 E) r* J* @1 ~( \& m  }c=std(error)/std(x0);                      %计算后验差比 *2
/ B  e, \9 e) q8 m6 ?. ]+ [" e8 yrelerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);                       
p=0;
. A6 K, A2 j4 V% A1 h# wfor i=2:ncol                                    
) ~/ g, S: ~8 Y5 d    if  relerror(1,i)<0.6745*std(x0)               
& ]7 A" C7 ]6 F- s: N7 P        p=p+1;                                 
    end
2 i9 e! _) o4 a2 S$ e2 lend
6 ~) Q" S7 [/ h( }7 i) k2 ~& _" b3 n6 Mp=p/(n-1);                     
2 t& d( ]1 T5 Cw1=min(abs(error));                              
$ n" {% x0 ?0 b4 j4 f6 }3 i0 jw2=max(abs(error));
i=1:n;                                     %计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
w=sum(w)/(n-1);
! n$ z7 h" o; y$ @) o& Q0 pau                                         %输出参数a,u的值
ago;                                        %输出累加数列ago的值
- I) J# ~$ G, V1 X. qx0;                                         %输出原始序列值
7 v( L5 q( y9 ]: j8 Lf=yc;                                        %输出预测的值
error;                                        %输出残差的值
c;                                           %输出后验差比的值
p;                                           %输出小误差概率的值
w                                           %输出关联度
' m2 T8 t) x) E& h# G. Send                                   
x0=[31.26
32.09
33.39
# V1 L. R9 [/ t, b$ d: _4 O35.45
40.52
- O0 s. V  ^; v" @43.52
47.86
7 F- U$ o3 p5 F: e" O- {51.45
55.6
60.14
64.82
' G$ x; e0 Y2 K/ [. M68.65
73.22
+ f$ D7 w% i! b: h1 J# H* B' g9 X80.22
2 o8 R( u8 ~! p# i* C  v- D2 W- s87.69
! ^) L( w/ N* Y# D# N/ M" x8 {93.97
99.16
: b$ v, S* D5 X8 N4 Y103.38
109.46
114.6
119.85
124.92
; w  Y& B9 i4 q132.04
0 x, ~* ~0 N8 L* U" e8 L3 z3 |0 v! e139.45
+ [& e+ K% c2 u  `+ S: L! z150.93
165.13
6 }6 ~5 A) p1 ~" e) n7 @181.93
: U3 y/ W7 P) Y0 }196.83
; h0 C0 F& b5 O9 c212.38
228.07
& H* p7 N/ f0 m7 K! H241.45
* {0 h1 q4 ~+ m% v! @- h251.03
]
我在窗口输入：gm(x0) 为啥会出现
Error in ==> gm at 22
4 r% q* I4 Q& S5 d8 lerror(i)=abs(yc(i)-x0(i));                     %计算残差值
" N6 E2 d* Q, [) C0 y那位高手帮忙 谢谢
" V/ n# F: W* Y: v+ {" p

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉