急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
% h9 O% J* X0 @6 g! Lfunction f=gm(x0,m)              %定义为函数gm(x)
3 a% x: A7 y' L, F4 j  M4 p; ^8 r0 Qn=length(x0);                     
x1=zeros(1,n);                    
7 q5 D- b6 ^: g7 \x1(1)=x0(1);                  
for i=2:n                         %计算累加序列x1
    x1(i)=x1(i-1)+x0(i);         
end
i=2:n;                            %对原始数列平行移位并赋值给y
y(i-1)=x0(i);                     
y=y';                                 
i=1:n-1;                             
  B$ l$ {, k* V+ f/ [6 E* }c(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为弧度
4 Q2 S9 B& {3 n* T8 s                                            
i=1:m;                             %计算预测累加数列的值
' q. y: t; L- D1 i' n# n. \$ B% C/ rago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
# {9 H* J4 s1 O  H+ gyc(1)=ago(1);                           
i=1:m-1;                             %还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
i=2:n;                              
9 U4 J4 _! s/ ^) K0 M1 z% f: werror(i)=abs(yc(i)-x0(i));                     %计算残差值
; D) c. u$ i& Z  lyc(1)=ago(1);                       
i=1:m-1;                              %修正还原数列的值
) [0 e' ~) e% K4 P& j2 K' Q0 F4 _4 zyc(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);                       
7 ~8 s6 m) E8 N! T5 c  R( `p=0;
. r! f! I2 \5 Z/ W: x- Jfor i=2:ncol                                    
    if  relerror(1,i)<0.6745*std(x0)               
% ]2 y# s' I* z8 S/ l) [% [! y. X        p=p+1;                                 
3 j: h, \( A5 C( ~# `( B* j& ]    end
* l8 A: @6 J7 i: I& r; p1 |end
p=p/(n-1);                     
! C, c# Z+ G; h( pw1=min(abs(error));                              
; e5 B5 [* e9 O: Uw2=max(abs(error));
( b& m' |$ v9 N9 U6 qi=1:n;                                     %计算关联度
' G  B3 T; [, E, Tw(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
w=sum(w)/(n-1);
8 e8 k3 x% n0 v* Y" C" h% Gau                                         %输出参数a,u的值
ago;                                        %输出累加数列ago的值
x0;                                         %输出原始序列值
f=yc;                                        %输出预测的值
7 S8 h  f2 S! g: _$ \4 [$ Zerror;                                        %输出残差的值
& ?9 J1 ~1 E* I' Ec;                                           %输出后验差比的值
5 I* |. [0 _, u) L2 }9 Qp;                                           %输出小误差概率的值
w                                           %输出关联度
# v$ S0 ?7 C+ D# ^  u6 Y0 xend                                   
& N9 ?" A5 B( r8 Z/ }9 p( ?x0=[31.26
3 k+ K4 Z7 n4 R2 r, Y% L32.09
+ ~  V# Y- w* o$ ~, N' y  ^7 v0 b2 m33.39
35.45
40.52
43.52
47.86
51.45
2 c/ |$ R  [+ w55.6
7 A4 i3 |/ Z& C5 M; _60.14
64.82
68.65
/ r3 z$ s4 f# e- N, e- O8 h73.22
80.22
: F3 U; y/ M6 a) J) Z87.69
6 M+ e) P( f% B$ I93.97
1 J0 r8 y" p) u6 \9 g) t, C* a# Y( ]2 W99.16
. @6 q% p) S4 D* K/ ~/ t103.38
! S* B& L! y: q5 s& v( ]7 _) G109.46
9 g$ H- D% N. ^2 {* k# V114.6
119.85
) p0 T7 q6 W7 j2 o, r" Y124.92
& m9 n* r3 z. c, i132.04
139.45
150.93
165.13
' a: q" Y2 O$ H  e& x) |181.93
196.83
9 i  R4 ^+ B) s9 \# v% F* h212.38
, r' X% k; ~; d2 Q* ?2 K228.07
241.45
! z. Z/ M; y" C# E251.03
]
我在窗口输入：gm(x0) 为啥会出现
: m4 w- @  o5 g( `9 |7 PError in ==> gm at 22
  c( [5 \2 e/ F' O2 U1 K1 X- m8 K' Cerror(i)=abs(yc(i)-x0(i));                     %计算残差值
7 _+ d$ k& ^8 o: B6 {1 A. |那位高手帮忙 谢谢
& h$ V" h/ E$ r5 e# T6 v

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉