急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
( |8 G9 [1 ~7 Z1 xfunction f=gm(x0,m)              %定义为函数gm(x)
n=length(x0);                     
x1=zeros(1,n);                    
x1(1)=x0(1);                  
for i=2:n                         %计算累加序列x1
0 @  S7 y0 U3 @, W5 W0 C    x1(i)=x1(i-1)+x0(i);         
6 e& F1 I7 b, j% J% m% ^end
  X& O0 C* ~% Y; W$ A1 d3 p" }/ [i=2:n;                            %对原始数列平行移位并赋值给y
y(i-1)=x0(i);                     
* v! M" Y6 ^6 U  Y, \, x3 ly=y';                                 
1 K1 D, F' Q5 I1 ^i=1:n-1;                             
% L' K, [# o! ~) J8 ]c(i)=-0.5*(x1(i)+x1(i+1));                 c(i)=-0.5*(x1(i)+x1(i+1))
B=[c' ones(n-1,1)];                        
, {& }* Z' R0 X6 Q2 mau=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
9 Q' _$ {) C# _+ {; @                                            
6 ]0 }( @+ g3 S3 Fi=1:m;                             %计算预测累加数列的值
1 U! t7 u0 z( d0 Aago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
yc(1)=ago(1);                           
4 @; l5 P6 O' l' D( x/ r' P8 O/ Oi=1:m-1;                             %还原数列的值
& E3 H6 ~2 x$ m( g$ Oyc(i+1)=ago(i+1)-ago(i);              
5 B4 k0 ?; k( c/ B) \- Qi=2:n;                              
, D/ S; F# @7 ~# M1 |error(i)=abs(yc(i)-x0(i));                     %计算残差值
yc(1)=ago(1);                       
' w' [' r8 g. l3 ri=1:m-1;                              %修正还原数列的值
yc(i+1)=ago(i+1)-ago(i);              
) o" d' b2 a0 o% A! dc=std(error)/std(x0);                      %计算后验差比 *2
relerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);                       
7 N/ l; `+ {2 k8 {+ X3 E' kp=0;
/ V& Y% R# T( ?$ I3 Dfor i=2:ncol                                    
    if  relerror(1,i)<0.6745*std(x0)               
        p=p+1;                                 
! k3 a6 h' |: N2 S1 g* Z; S5 k    end
8 m( L3 N$ N/ G+ e! x* K( \end
p=p/(n-1);                     
: T$ m. w7 _# nw1=min(abs(error));                              
w2=max(abs(error));
9 Q. j: M$ ~' r" V1 si=1:n;                                     %计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
. t" G" e2 s7 Fw=sum(w)/(n-1);
au                                         %输出参数a,u的值
# o5 X. _8 Z( o2 y! R2 Q8 m, x$ v+ ?$ x4 Yago;                                        %输出累加数列ago的值
0 ~; u& q2 u& xx0;                                         %输出原始序列值
f=yc;                                        %输出预测的值
error;                                        %输出残差的值
c;                                           %输出后验差比的值
  V# a8 o" K3 T+ Q2 f5 Q9 q6 qp;                                           %输出小误差概率的值
w                                           %输出关联度
- @. `2 }1 N& A! q$ jend                                   
0 l! p6 |: _% G' I( E+ J, a1 s8 Y9 Kx0=[31.26
32.09
' e( t, y; o+ g33.39
: x' ?7 h3 x: J35.45
6 `3 D* i( v- T! o; F& c/ A+ l40.52
43.52
47.86
  \* J, I/ X! `$ \9 d51.45
, ~' a( b" \. t6 ^2 }+ \, L55.6
% ^2 P6 P) n9 u* E60.14
64.82
# l4 Z$ N9 y1 Y68.65
# A# D. S) b5 W1 S: a* u9 f6 T73.22
80.22
. N; f8 c# w) f! G" r2 n+ F( }/ V87.69
93.97
, R% m+ ~9 O. W, Q6 l99.16
103.38
6 ]8 O' p: h$ }6 k' @" z109.46
114.6
119.85
! r; B- W7 @$ C) q$ X124.92
0 d& ~. X8 w3 G6 W132.04
139.45
150.93
165.13
181.93
  ~, b* x" u' v- C! L5 e  N9 i; t196.83
212.38
0 ?. G1 t$ g: N3 |6 D4 x228.07
241.45
251.03
]
我在窗口输入：gm(x0) 为啥会出现
1 X) U; [+ {0 OError in ==> gm at 22
  ~( A0 I( n- r% \# |- W, H; Yerror(i)=abs(yc(i)-x0(i));                     %计算残差值
9 d$ |$ K% F! \- I7 w, x1 S, l那位高手帮忙 谢谢

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉