急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
function f=gm(x0,m)              %定义为函数gm(x)
n=length(x0);                     
x1=zeros(1,n);                    
% h0 d# W. E; Q; W5 r' \x1(1)=x0(1);                  
' [# ]0 x+ r. Z- Q" U6 }for i=2:n                         %计算累加序列x1
    x1(i)=x1(i-1)+x0(i);         
end
i=2:n;                            %对原始数列平行移位并赋值给y
y(i-1)=x0(i);                     
7 P& M$ G# U* o6 D  T6 z' X' Q# Fy=y';                                 
i=1:n-1;                             
! ^/ K) z# w+ wc(i)=-0.5*(x1(i)+x1(i+1));                 c(i)=-0.5*(x1(i)+x1(i+1))
) h5 P( `2 p1 {: h: gB=[c' ones(n-1,1)];                        
au=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
# g+ `) g% R: u0 o0 Y, b                                            
i=1:m;                             %计算预测累加数列的值
ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
- ^1 e$ s( c' n; Y3 l5 Ayc(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));                     %计算残差值
yc(1)=ago(1);                       
- s0 P8 b1 }3 r0 \8 `i=1:m-1;                              %修正还原数列的值
2 v) a! r3 ]+ X+ e+ _! Xyc(i+1)=ago(i+1)-ago(i);              
c=std(error)/std(x0);                      %计算后验差比 *2
5 \: Q% W* c. Y. E8 \& z* }9 L3 n$ yrelerror=abs((error)-mean(error)*ones(size(error)));
% ]# M4 h6 }8 p4 m[nrow,ncol]=size(relerror);                       
p=0;
: k, \1 P$ b, M2 dfor i=2:ncol                                    
    if  relerror(1,i)<0.6745*std(x0)               
1 v% B9 u, S  [% ]; d# V$ K        p=p+1;                                 
    end
end
9 _) D( s+ d" p8 G& o" X- c' [p=p/(n-1);                     
# b% B/ \- l$ I) ~/ uw1=min(abs(error));                              
4 ^) v5 Y- ?! K2 u5 J6 E/ w; A) Fw2=max(abs(error));
i=1:n;                                     %计算关联度
) f. p6 A& e0 [$ ]  t) zw(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
w=sum(w)/(n-1);
au                                         %输出参数a,u的值
ago;                                        %输出累加数列ago的值
  f& q" l8 x$ L% [! c2 Bx0;                                         %输出原始序列值
f=yc;                                        %输出预测的值
6 s! k; k+ q/ ^9 {) V4 yerror;                                        %输出残差的值
& e- B: D7 Q, Yc;                                           %输出后验差比的值
p;                                           %输出小误差概率的值
/ k. P) o1 E2 a- Mw                                           %输出关联度
& p% o2 S0 Q- l# L3 vend                                   
1 B' ]* a1 ]$ xx0=[31.26
' v' l: U5 ~# {2 t6 R& |" {4 s  X32.09
) c8 Z" y) H. ?7 s, b33.39
% Y( M+ x+ w; S8 V7 p: Y9 m4 P, ~& G35.45
! C7 }0 S0 U  X) B$ }" ?40.52
, [: H9 I0 o) C- {( ]43.52
47.86
2 j+ t( \  ]# ]& ?# ^4 u51.45
55.6
$ t9 l& l" A" N7 e+ o+ U7 `  D60.14
64.82
68.65
5 F1 j7 c4 p# r1 }4 H73.22
80.22
; P4 r7 @5 F% g' M% m8 A2 I! Z* U87.69
! o- T* m# B+ w% X' S* a/ K5 n93.97
99.16
103.38
109.46
114.6
119.85
124.92
132.04
5 x8 x8 j8 n' r/ W; P7 P2 g139.45
150.93
4 k9 K+ `( ~6 m165.13
181.93
5 V8 y/ G6 e) u# u196.83
, O/ E: x6 s4 L0 ~$ ^) d8 u4 k212.38
& V2 b, s$ ^$ Z1 ^+ K. J, t228.07
; ?( L0 O2 `" X# [241.45
251.03
( y1 S7 B6 C' |/ |4 ?; w4 [/ {]
我在窗口输入：gm(x0) 为啥会出现
Error in ==> gm at 22
- M1 z7 q2 Q  Q$ N% cerror(i)=abs(yc(i)-x0(i));                     %计算残差值
那位高手帮忙 谢谢

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉