数学建模社区-数学中国

标题: 求多因子灰色预测GM(1,N)源程序 [打印本页]

作者: hujet    时间: 2012-5-6 15:43
标题: 求多因子灰色预测GM(1,N)源程序
求多因子灰色预测GM(1,N)源程序,不是GM(1,1)的,我要求的是GM(1,4)可是不知道怎么弄程序,谢谢了
作者: chen7705    时间: 2012-8-27 13:13
对,好多都没灰色预测GM(1,N)源程序
作者: 5105    时间: 2013-2-1 19:37
同求啊啊啊……
作者: 栀子流年    时间: 2013-11-18 22:31
+1求程序啊
作者: lvjianchi    时间: 2013-12-22 22:48
n=length(x10); 6 E; E5 \6 d: T5 I3 Z0 @. w
x11=cumsum(x10)  9 ]) r9 X$ ~. ?7 h- _6 i
x21=cumsum(x20)  
8 _/ k0 Z9 o1 C/ J" H4 Ifor  i=2:n
: d/ @' Y/ z* g, t    z11(i)=0.5*(x11(i)+x11(i-1));
) d; X% h) b+ Oend  
1 r& J2 b) x' y* Z7 IB=[-z11(2:n)',x21(2:n)'];  
, Q. R2 W2 W1 \7 G* ?4 ^Y=x10(2:n)'; ) H7 [4 h8 \' ]  O. W; \* z8 r2 q6 q' S
u=B\Y  
# P" q6 @" z. }0 Rx=dsolve('Dx+a*x=b*x2' ,'x(0)=x0');
7 }* r+ f- ~; ?7 V# P- gx=subs(x,{ 'a' ,'b' ,'x0','x2'},{u(1),u(2),x10(1), 'x21' });  8 l8 J! R; u3 K4 U/ n
digits(6),x=vpa(x);x=simple(x)  
- E, |( V6 k8 \% n+ ]x=subs(x,{ 't' ,'x21' },{[0:n-1],x21(1:n)})
" F- f8 E# b5 @; W  F, @+ x7 pxhat=[x(1),diff(x)] 6 U4 g3 {" I" P
epsilon=x10-xhat
& _; Y. B# s! c$ z  U4 ~% zdelta=abs(epsilon./x10)  
作者: lvjianchi    时间: 2013-12-22 22:49
n=length(x10); 6 i8 f& p5 `  `6 a3 ~$ X
x11=cumsum(x10)  $ v9 E6 W" e" Z0 ~! l
x21=cumsum(x20)  
! h! Z7 A* K! j. }% \for  i=2:n
2 e' x) G% m: r! n1 F) S2 \6 V4 y    z11(i)=0.5*(x11(i)+x11(i-1));
! i% [7 R4 c6 K7 ?9 n- U0 send  
" w  o, c3 [4 i; l) RB=[-z11(2:n)',x21(2:n)'];  ' \% t3 q/ C) A
Y=x10(2:n)';
' j0 F  ~0 s* D7 Au=B\Y  3 \+ h8 u9 M" c5 ~, V, g
x=dsolve('Dx+a*x=b*x2' ,'x(0)=x0'); ' e. k3 L. x. m
x=subs(x,{ 'a' ,'b' ,'x0','x2'},{u(1),u(2),x10(1), 'x21' });  ' E# j3 q1 x7 @! h0 _
digits(6),x=vpa(x);x=simple(x)  
% f8 t9 U" T1 w  s! l, t# d9 Nx=subs(x,{ 't' ,'x21' },{[0:n-1],x21(1:n)}) # v& Y& U  z; ^) ]: J( P! m3 T
xhat=[x(1),diff(x)] ; B! n6 S& h; s/ T- s" W
epsilon=x10-xhat
# I1 `# z$ S8 v% d2 ?6 f. S7 ^delta=abs(epsilon./x10)  
作者: 远行的小船儿666    时间: 2015-12-28 16:00
可以看看的                         : D+ h6 U, N3 R2 d, R" @  ^" E

作者: 远行的小船儿666    时间: 2015-12-28 16:14
好的呀,加油                        
3 A6 ?) [. F) f9 J' `" Q# {
作者: 千梦君    时间: 2019-1-7 14:14
EHWVFHWDUEWGUEIF
7 ?) A9 A2 I& y) ]* {# ~( u




欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5