急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
function f=gm(x0,m)              %定义为函数gm(x)
0 ^9 H# r  Z; H6 {# f! wn=length(x0);                     
$ e0 b3 }' _; i% Rx1=zeros(1,n);                    
4 w8 ]' Z! L2 rx1(1)=x0(1);                  
0 \; r$ t6 \0 k7 i" _for i=2:n                         %计算累加序列x1
! c7 o, W) i; N' B    x1(i)=x1(i-1)+x0(i);         
end
  f+ V, C, ?" c1 |1 {i=2:n;                            %对原始数列平行移位并赋值给y
4 B" L% y+ ^3 By(i-1)=x0(i);                     
* W/ I  ~" o$ Qy=y';                                 
i=1:n-1;                             
c(i)=-0.5*(x1(i)+x1(i+1));                 c(i)=-0.5*(x1(i)+x1(i+1))
+ L; c) k9 ]" r: f7 `B=[c' ones(n-1,1)];                        
& L6 G% O7 Q6 I* G/ X, gau=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
4 Z: p7 @8 i' f: O& K                                            
" p; W' c: ~+ C' H9 R& ki=1:m;                             %计算预测累加数列的值
ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
yc(1)=ago(1);                           
i=1:m-1;                             %还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
i=2:n;                              
error(i)=abs(yc(i)-x0(i));                     %计算残差值
& F7 H* p* P  B2 yyc(1)=ago(1);                       
i=1:m-1;                              %修正还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
c=std(error)/std(x0);                      %计算后验差比 *2
relerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);                       
9 p1 Z. T7 \/ W5 Fp=0;
# p8 |1 _* B7 `for i=2:ncol                                    
    if  relerror(1,i)<0.6745*std(x0)               
        p=p+1;                                 
    end
$ W( @, Q) ?: B; b' [4 s6 Bend
p=p/(n-1);                     
w1=min(abs(error));                              
% z6 A: s3 a' ]/ p+ s) q3 Aw2=max(abs(error));
i=1:n;                                     %计算关联度
4 x( c* ^, ]( W8 Y; yw(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
$ p+ d6 r. G% T! j0 fw=sum(w)/(n-1);
au                                         %输出参数a,u的值
ago;                                        %输出累加数列ago的值
0 w) l( I* _1 Z/ H; `4 V3 Wx0;                                         %输出原始序列值
& R& N' J0 m" O- u& M5 |, Hf=yc;                                        %输出预测的值
! L/ Q' q6 j" ~error;                                        %输出残差的值
) {/ k1 q0 t& ]3 Gc;                                           %输出后验差比的值
p;                                           %输出小误差概率的值
w                                           %输出关联度
" l! \$ n- q3 v  mend                                   
x0=[31.26
32.09
) H9 p$ }: r1 w' T6 v9 S. C- S! c7 I33.39
35.45
& |4 D5 f2 ]9 I1 a40.52
43.52
47.86
! E8 }$ [1 f$ B% A* I0 l51.45
55.6
$ \. t7 f: X+ ~7 |  t, d- t60.14
64.82
& q( x; f6 m  ~% r# R& N" s68.65
73.22
. K& U3 X" e+ K5 V80.22
87.69
" t2 I9 |3 F5 g4 V: h0 @, y# Q4 g93.97
99.16
103.38
' r1 u* m3 a4 a2 z. N109.46
- R2 G9 s1 y( N! [114.6
119.85
124.92
% ~; C0 [* I" I/ g" Y132.04
! ?8 h$ `9 ]  {( q139.45
150.93
8 J- |9 b. ^/ v5 j5 G! U6 K165.13
181.93
! l' ~& L: U% _( c3 F196.83
212.38
& |7 _1 x( O& C* q, N$ {9 Z228.07
5 @$ o% I9 m; a) ]241.45
/ _3 G5 F5 k9 l. o: {. l# K251.03
]
0 |% P0 R: l* `7 Z1 z我在窗口输入：gm(x0) 为啥会出现
7 I! g: n: ?9 x; ~8 Q  k  }$ Y7 bError in ==> gm at 22
6 o2 ?* |$ c: b; ferror(i)=abs(yc(i)-x0(i));                     %计算残差值
) e' H& i6 v; y) {7 C- [那位高手帮忙 谢谢
6 B# m( ]. K3 Z1 r: \7 S

wuzhenhua

把abs（）去掉

神经病个人组

我也想知道