急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
function f=gm(x0,m)              %定义为函数gm(x)
n=length(x0);                     
# q; }7 Y2 B' |* G) I" lx1=zeros(1,n);                    
" L1 S; f4 o$ @5 w" P& {+ Ex1(1)=x0(1);                  
for i=2:n                         %计算累加序列x1
! c5 j) G& [) K4 p2 S# `    x1(i)=x1(i-1)+x0(i);         
end
% Q3 ^% T, y- S8 ]i=2:n;                            %对原始数列平行移位并赋值给y
  C. z+ W9 B; k  ^9 ly(i-1)=x0(i);                     
1 o8 l! ?) L. z( [8 U7 Ty=y';                                 
8 ~3 P' a  i; v0 J6 bi=1:n-1;                             
' c3 i/ z( L5 ]# g% I. b) xc(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为弧度
                                            
i=1:m;                             %计算预测累加数列的值
0 y3 Z5 Z  [9 x% D, r5 Z+ p4 ]ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
yc(1)=ago(1);                           
, H8 ^' r! Z) j& J- N4 H' \i=1:m-1;                             %还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
0 i( R% f7 ]# J6 N- O/ ^i=2:n;                              
2 \8 I0 f9 K  g8 M- Oerror(i)=abs(yc(i)-x0(i));                     %计算残差值
yc(1)=ago(1);                       
4 h2 b% @; ]# c' wi=1:m-1;                              %修正还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
' r  E' g$ A: U8 Y3 [c=std(error)/std(x0);                      %计算后验差比 *2
* ?3 Z$ e/ v, v, w* I6 srelerror=abs((error)-mean(error)*ones(size(error)));
/ m/ V2 S& C; j. M3 F& ~[nrow,ncol]=size(relerror);                       
2 [. b0 R8 H0 w' x  R8 Up=0;
+ t" v/ i# {1 U. }5 \  {* Lfor i=2:ncol                                    
2 A9 d! }, `9 ^2 [: W1 ]; s    if  relerror(1,i)<0.6745*std(x0)               
        p=p+1;                                 
    end
end
p=p/(n-1);                     
0 ]. c- f1 v/ O0 y6 d( I- cw1=min(abs(error));                              
8 H+ |9 E% G* X3 @1 L: H3 U! D. Xw2=max(abs(error));
5 _& X) D5 o9 D3 M! {# Hi=1:n;                                     %计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
% s8 d5 [( p/ {9 ~/ cw=sum(w)/(n-1);
. _$ `/ S/ ?7 F# dau                                         %输出参数a,u的值
ago;                                        %输出累加数列ago的值
x0;                                         %输出原始序列值
f=yc;                                        %输出预测的值
error;                                        %输出残差的值
" O# G8 D5 I  q( G, V, {2 Bc;                                           %输出后验差比的值
+ l! b; l& _5 I* _" lp;                                           %输出小误差概率的值
w                                           %输出关联度
5 f# G' Q. I7 V& Z' }0 s5 Eend                                   
x0=[31.26
, h; R+ ^* w: y* m# Y32.09
33.39
35.45
" |6 I) e: B$ l, h. g40.52
- X9 c8 p  j* J7 l# x43.52
47.86
# ]5 }0 k, b9 b2 ^1 p7 F( b3 v51.45
55.6
) ^) U+ C# ]. Q! s7 K- t60.14
/ {1 k, D, T/ s9 x64.82
9 S! [# x1 g# |# Z/ p  y1 p/ H68.65
73.22
80.22
87.69
93.97
99.16
103.38
/ W9 r$ q- `& q# e109.46
114.6
119.85
2 D, J9 n  W1 V+ Y) {% _& i124.92
/ D  P. P/ o, a! c* ^; W" R0 k132.04
139.45
150.93
165.13
181.93
# y; r# h  m8 ^196.83
212.38
+ m- a& O8 \; ~) k3 K228.07
* P) S0 {" o8 z* z$ e241.45
251.03
]
4 y% ]9 K3 Z" D1 u4 q4 @; w1 a# ^我在窗口输入：gm(x0) 为啥会出现
2 d+ C+ s- \4 {5 Y( I; \2 zError in ==> gm at 22
error(i)=abs(yc(i)-x0(i));                     %计算残差值
5 s) X/ b$ c+ w+ s  F, ?& V那位高手帮忙 谢谢
" R3 k( u  m! j6 f; x/ Z( ~

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉