急需求助 灰色模型

1530568658

我要预测10年之后的人口，用灰色模型
function f=gm(x0,m)              %定义为函数gm(x)
n=length(x0);                     
# ]- X/ h- a2 p$ V, vx1=zeros(1,n);                    
. z4 S& Z) q& D! ^& F: mx1(1)=x0(1);                  
; w! r( g2 q1 c% X5 C+ lfor i=2:n                         %计算累加序列x1
    x1(i)=x1(i-1)+x0(i);         
end
/ G+ {; G$ b; M5 ^$ si=2:n;                            %对原始数列平行移位并赋值给y
) s- S% ~* y5 m% ay(i-1)=x0(i);                     
- N% h7 q0 J( Z5 Ay=y';                                 
8 ^) c, B8 {3 s, C2 B$ p2 h, hi=1:n-1;                             
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为弧度
9 {" ]7 P, c+ y$ t8 G" Y0 T4 u8 G2 L                                            
i=1:m;                             %计算预测累加数列的值
ago(i)=(x0(1)-au(2)/au(1))*exp(-au(1)*(i-1))+au(2)/au(1)   
yc(1)=ago(1);                           
6 \  _+ p! s+ i1 C; H0 r6 e( xi=1:m-1;                             %还原数列的值
/ l* O3 |" \6 \7 ^/ @$ Yyc(i+1)=ago(i+1)-ago(i);              
" `7 t' Z) K7 [" b' r# Oi=2:n;                              
- o0 f! b& b- T, N% S& a' P3 o9 Rerror(i)=abs(yc(i)-x0(i));                     %计算残差值
1 F# K& f8 [: \3 t! tyc(1)=ago(1);                       
i=1:m-1;                              %修正还原数列的值
4 F* W3 k( T! F8 f* s' cyc(i+1)=ago(i+1)-ago(i);              
6 I& L0 \! r5 Ac=std(error)/std(x0);                      %计算后验差比 *2
relerror=abs((error)-mean(error)*ones(size(error)));
[nrow,ncol]=size(relerror);                       
7 F; e. M/ Y) ?( K; E0 l! `p=0;
, A) h5 B- v: I0 Hfor i=2:ncol                                    
    if  relerror(1,i)<0.6745*std(x0)               
1 [, d6 F8 B" s/ e. y$ h! V        p=p+1;                                 
" V$ J: n# t# w) ]  `3 c    end
end
p=p/(n-1);                     
: i8 W! \9 K  O+ c# {  M" X# Bw1=min(abs(error));                              
w2=max(abs(error));
i=1:n;                                     %计算关联度
w(i)=(w1+0.5*w2)./(abs(error(i))+0.5*w2);
) f, H0 {" F" S( h  Uw=sum(w)/(n-1);
: j6 I/ i, Q% l1 S# P) L; Kau                                         %输出参数a,u的值
& L, Z$ Z4 {" ]' |: f, a4 Jago;                                        %输出累加数列ago的值
$ [  L/ @% V# C8 J$ W" Y+ Ax0;                                         %输出原始序列值
7 x/ t& T$ @& a- [9 Q$ D% pf=yc;                                        %输出预测的值
* c/ R" [# `* K, |; @5 T6 Cerror;                                        %输出残差的值
/ i6 ^" P' C# L- J  L- H7 H* Cc;                                           %输出后验差比的值
  E* z( n8 N5 k+ p6 I: Zp;                                           %输出小误差概率的值
w                                           %输出关联度
6 d4 E) F: [# V5 Cend                                   
x0=[31.26
" c' K% }( |, o' e, }32.09
% j3 @# |3 H- W3 A6 w  J& F" |33.39
35.45
1 H9 A+ R& j5 u! t40.52
43.52
7 \4 [* B- J) k6 D# y7 h5 l47.86
9 @; @/ p4 H) g! N5 l6 ^0 _1 l51.45
55.6
60.14
64.82
68.65
73.22
80.22
* ~: d9 A5 o4 x8 ~& }* s87.69
- c  ^# c0 I) h2 x/ b: {4 R1 q93.97
; D  w5 ~/ A0 [9 E' R' E99.16
" n7 ?6 I7 r3 S; n% o( z103.38
5 W7 `/ X# j' }- W109.46
114.6
119.85
1 e& d1 [0 k3 k# L( {124.92
132.04
139.45
  B* j+ t& V5 }% I4 y; [: x150.93
165.13
# |9 E2 d5 C: l. M. [181.93
196.83
5 H, K" O; \3 L) B$ v; P  B. H212.38
228.07
! x: n; l: O5 S241.45
251.03
1 P4 D3 j7 e9 E- s; }) S]
  [" n: {" G' L- h' W, F0 `我在窗口输入：gm(x0) 为啥会出现
# Q3 C& @" ~* U5 B6 [. i+ n% rError in ==> gm at 22
9 p" r3 n9 G7 s) ]7 H# Serror(i)=abs(yc(i)-x0(i));                     %计算残差值
9 Y; X5 R" r+ j! x那位高手帮忙 谢谢

神经病个人组

我也想知道

wuzhenhua

把abs（）去掉