急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
function f=gm(x0,m)              %定义为函数gm(x)
n=length(x0);                     
x1=zeros(1,n);                    
x1(1)=x0(1);                  
for i=2:n                         %计算累加序列x1
    x1(i)=x1(i-1)+x0(i);         
end
. [' l% z$ z3 ^1 Ti=2:n;                            %对原始数列平行移位并赋值给y
8 C# C8 I" R0 `7 Sy(i-1)=x0(i);                     
+ S) x) a7 i. B% q' _! ^) Jy=y';                                 
! A1 s& Q, t+ D6 q" {" l* Ai=1:n-1;                             
c(i)=-0.5*(x1(i)+x1(i+1));                 c(i)=-0.5*(x1(i)+x1(i+1))
4 N1 P. k! X: oB=[c' ones(n-1,1)];                        
. b1 ]' |3 H, ~9 C5 H* eau=inv(B'*B)*B'*y;              %计算参数a,u矩阵inv 渐开线函数inv（a）=tan(a)-a a为弧度
2 e, ^2 q5 F" g, Y                                            
7 Y  Z( e4 S2 N: O+ _6 R1 zi=1:m;                             %计算预测累加数列的值
ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
. ^! c1 U0 ^- e- @yc(1)=ago(1);                           
i=1:m-1;                             %还原数列的值
: Q( o! ~4 ~+ j" W- J; c" Byc(i+1)=ago(i+1)-ago(i);              
i=2:n;                              
error(i)=abs(yc(i)-x0(i));                     %计算残差值
0 M3 d, t+ E8 v9 `- i) x- L" [yc(1)=ago(1);                       
( }3 ]0 v7 i- M+ W: ~: `( Pi=1:m-1;                              %修正还原数列的值
2 J6 L7 P) W- [) S3 L8 f: V- Vyc(i+1)=ago(i+1)-ago(i);              
# e* b0 y9 Z3 \c=std(error)/std(x0);                      %计算后验差比 *2
relerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);                       
; t; @6 T& t. e8 Ep=0;
for i=2:ncol                                    
    if  relerror(1,i)<0.6745*std(x0)               
        p=p+1;                                 
+ x7 U/ K* u) J# _( @" q& Y    end
! [* c! M' c1 t0 wend
5 a% o( e% a! d" V4 p& cp=p/(n-1);                     
5 c3 Z; T7 w3 F4 ~: `8 [w1=min(abs(error));                              
w2=max(abs(error));
i=1:n;                                     %计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
w=sum(w)/(n-1);
) k* i, v; _# X7 Oau                                         %输出参数a,u的值
0 t" z, w/ p4 Sago;                                        %输出累加数列ago的值
- f3 }1 d& z* q7 @: f  s" Dx0;                                         %输出原始序列值
& @' [. @; ~0 s  [' }( a& j# mf=yc;                                        %输出预测的值
error;                                        %输出残差的值
c;                                           %输出后验差比的值
p;                                           %输出小误差概率的值
% q" N( @- N$ y- q3 `+ _6 Bw                                           %输出关联度
end                                   
2 P: n7 t; t; Fx0=[31.26
32.09
33.39
* e5 o9 Y  C. G" [, l9 u% e35.45
9 d" n" L# q3 j- T40.52
43.52
47.86
+ X+ B( V1 u: a  N3 o51.45
6 W- f5 B: Y0 @% y* O- Q55.6
60.14
+ o. T! |+ G0 E" w64.82
68.65
73.22
80.22
/ w9 T, ^7 F/ g; O9 e. Q87.69
93.97
99.16
# f5 {: W, l5 T103.38
' V' Z4 o& ^8 B. S0 u5 o) D7 G109.46
114.6
119.85
124.92
. G1 b8 V" A" d' i0 `% E132.04
: r2 t" D! v8 ]6 X! w139.45
150.93
165.13
181.93
1 l5 @# r2 U) W, F3 S196.83
6 V' X" Z* `& J0 y3 n8 H5 _3 h212.38
228.07
. K2 ]1 Y0 r( F241.45
251.03
6 h; f# |/ G7 m6 Q5 x0 N1 P]
2 o2 ]6 r8 o% p. d7 z# c我在窗口输入：gm(x0) 为啥会出现
/ v3 P4 w' ^& X) b" VError in ==> gm at 22
error(i)=abs(yc(i)-x0(i));                     %计算残差值
那位高手帮忙 谢谢

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉