急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
function f=gm(x0,m)              %定义为函数gm(x)
n=length(x0);                     
x1=zeros(1,n);                    
3 `0 O0 A1 J7 M7 Z2 ]# px1(1)=x0(1);                  
- c3 I2 K$ m9 x% R, t$ ^* Afor i=2:n                         %计算累加序列x1
    x1(i)=x1(i-1)+x0(i);         
) p, B% i% q0 z  j/ V& m! v0 r$ e5 J+ Mend
5 m) Q$ B6 n) c3 Z; Xi=2:n;                            %对原始数列平行移位并赋值给y
y(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))
B=[c' ones(n-1,1)];                        
+ z7 k6 s2 v; b3 u/ ?. Cau=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
                                            
i=1:m;                             %计算预测累加数列的值
" g( T* x  `) _! m) x8 w" h$ Zago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
5 ]2 l" T, Y! g$ uyc(1)=ago(1);                           
  P% C2 I7 g$ M: P% y1 L/ Ti=1:m-1;                             %还原数列的值
4 D: [: `* W4 y. G% n/ Jyc(i+1)=ago(i+1)-ago(i);              
i=2:n;                              
  Z9 f0 A2 \* F+ M; herror(i)=abs(yc(i)-x0(i));                     %计算残差值
5 Z6 Z: t9 S7 Oyc(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)));
- I* U& T# h1 |$ {6 o, X7 K[nrow,ncol]=size(relerror);                       
) G( o9 H5 S3 R/ zp=0;
* d$ q5 F* z# C( z1 gfor i=2:ncol                                    
  a! o" l* C& o* V( p/ c3 I1 X    if  relerror(1,i)<0.6745*std(x0)               
        p=p+1;                                 
    end
4 u: _9 e) {9 w+ E5 J, J8 }end
' _6 n, s, Z3 L3 U2 ~! k. rp=p/(n-1);                     
w1=min(abs(error));                              
w2=max(abs(error));
" g$ j9 _2 N3 `( K& r. oi=1:n;                                     %计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
w=sum(w)/(n-1);
au                                         %输出参数a,u的值
- b; @; }: D# H% Xago;                                        %输出累加数列ago的值
x0;                                         %输出原始序列值
3 T5 ^4 k% J3 o; I, b8 i6 Gf=yc;                                        %输出预测的值
error;                                        %输出残差的值
" W( Q* K( F$ }* A) Ic;                                           %输出后验差比的值
* ?9 r# B: g6 }, up;                                           %输出小误差概率的值
/ t; w, |, z9 h* l0 gw                                           %输出关联度
& {/ [) P( C" d. ~- h1 v) G* Hend                                   
, g) {. M* |; m) I# Ax0=[31.26
32.09
33.39
8 {6 H  C/ k! t2 c+ S& {. w# _35.45
40.52
1 k0 M' K' }* ^* _# L# ]8 N% ^3 ]" {43.52
) U' b% G2 Z8 |- q47.86
5 `+ o$ ]- F$ {. Y7 O$ D) b: g51.45
55.6
4 `  t# Y' |* e/ s6 A9 [' ?2 ~60.14
/ D+ d/ z4 o' A6 X, [- C9 ^8 ]8 ?64.82
68.65
73.22
80.22
87.69
1 v) {( i1 J8 {7 d93.97
4 N, T/ e+ D0 j  O+ E99.16
103.38
109.46
3 ~& M" M) u$ _/ V" d- [114.6
119.85
124.92
5 B9 {1 o& b+ y1 W" E$ N$ h132.04
139.45
150.93
165.13
181.93
. L8 @: D8 _2 M8 z" r196.83
  x3 p) c! `  T212.38
228.07
241.45
251.03
]
. R+ L, ~# V) z2 n6 [9 k我在窗口输入：gm(x0) 为啥会出现
Error in ==> gm at 22
error(i)=abs(yc(i)-x0(i));                     %计算残差值
那位高手帮忙 谢谢
  E/ E8 L& D# H5 I. F$ b# P& d

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉