急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
2 K. ?3 W  n* y6 y3 _function f=gm(x0,m)              %定义为函数gm(x)
n=length(x0);                     
6 p+ X, n; o, Ax1=zeros(1,n);                    
- w+ P6 w, j( G: b4 ex1(1)=x0(1);                  
$ U  G2 Y  N- nfor i=2:n                         %计算累加序列x1
1 U, [+ b9 G6 q7 N! ]    x1(i)=x1(i-1)+x0(i);         
end
9 [  ?6 H) N  o$ ni=2:n;                            %对原始数列平行移位并赋值给y
9 z2 U9 S8 u* t$ `9 yy(i-1)=x0(i);                     
y=y';                                 
) }, Q, T6 J8 \; H9 D! U: Hi=1:n-1;                             
8 ?2 W' [  \# B! ~3 ~8 v: o# Hc(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为弧度
/ m6 i; W6 {6 h  A) o; U2 n                                            
i=1:m;                             %计算预测累加数列的值
3 @- ?! W. p- z2 o9 Dago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
yc(1)=ago(1);                           
4 V: v  T, j6 D# {i=1:m-1;                             %还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
+ w+ ?; q: Y  |i=2:n;                              
error(i)=abs(yc(i)-x0(i));                     %计算残差值
yc(1)=ago(1);                       
1 p, V* e# d- [0 f8 C' r8 Oi=1:m-1;                              %修正还原数列的值
5 }7 l6 ?8 b1 }& ~7 \yc(i+1)=ago(i+1)-ago(i);              
( E) E9 `3 z* i6 s. W2 }% cc=std(error)/std(x0);                      %计算后验差比 *2
. u1 t; d* p' x, {, y" ?relerror=abs((error)-mean(error)*ones(size(error)));
- [( ~% N+ o5 f[nrow,ncol]=size(relerror);                       
p=0;
+ m# U! w, t4 [; Cfor i=2:ncol                                    
    if  relerror(1,i)<0.6745*std(x0)               
        p=p+1;                                 
    end
end
p=p/(n-1);                     
w1=min(abs(error));                              
* c% G5 u4 Y% _w2=max(abs(error));
i=1:n;                                     %计算关联度
/ B+ V# k7 |2 o5 n: N% W( Nw(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
: {* Q; a5 o% gw=sum(w)/(n-1);
au                                         %输出参数a,u的值
3 C9 @& p. P! G3 aago;                                        %输出累加数列ago的值
0 d7 |/ f, `3 f+ K( B3 Xx0;                                         %输出原始序列值
3 R3 I/ a/ n' A7 z$ j4 m2 q1 zf=yc;                                        %输出预测的值
2 H) Q  J* c( Jerror;                                        %输出残差的值
c;                                           %输出后验差比的值
p;                                           %输出小误差概率的值
w                                           %输出关联度
end                                   
x0=[31.26
- r* J4 V5 ]1 Z$ D# S32.09
# B' Q, |- [) u3 F; Q; z33.39
35.45
40.52
* E$ m" d& a" ~7 @& s% I$ g( [0 g43.52
47.86
+ Y' o: F2 a5 v; |* B51.45
; H; C/ |" I9 G0 _55.6
60.14
% B" J  s2 J: [0 y( N$ g# Y64.82
) G; E4 D0 O  x" j4 U/ E! k68.65
73.22
80.22
. v- ^6 e5 E+ S5 ^87.69
6 a- O8 {, n' H2 e1 k93.97
99.16
103.38
109.46
/ J- m& B+ i7 w+ O' V/ j114.6
6 }5 `4 {  A7 g( W! f. U4 g3 K119.85
124.92
132.04
139.45
150.93
165.13
+ r, ?: e' o" j9 r$ E* w3 |181.93
- Q8 _8 x3 N3 G4 g0 K6 h196.83
212.38
228.07
241.45
( f9 {& ]: K5 X9 X) x251.03
]
7 F! Q% z* ?( |+ C3 ^9 V- s我在窗口输入：gm(x0) 为啥会出现
& [( z/ j% H. p! e" G2 M# MError in ==> gm at 22
+ \) |3 `3 j8 i  x) a0 f2 r8 _error(i)=abs(yc(i)-x0(i));                     %计算残差值
那位高手帮忙 谢谢

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉