数学建模社区-数学中国
标题:
求多因子灰色预测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);
$ [ B5 Q. o- B V/ t# o
x11=cumsum(x10)
( I h& [0 V3 K
x21=cumsum(x20)
; c+ A: j) ?9 I g! k. d' ?
for i=2:n
- S2 a9 f- O0 N
z11(i)=0.5*(x11(i)+x11(i-1));
$ `$ {9 r& K5 h
end
) i2 }- U+ C: `8 k% g. G
B=[-z11(2:n)',x21(2:n)'];
& M2 B' g- ^) ~; h; j+ f
Y=x10(2:n)';
0 Z5 [0 b. o: M7 H3 ]' `* c: b- H
u=B\Y
! ` u! ~. E1 W
x=dsolve('Dx+a*x=b*x2' ,'x(0)=x0');
* {9 v0 R$ D0 V: q9 H
x=subs(x,{ 'a' ,'b' ,'x0','x2'},{u(1),u(2),x10(1), 'x21' });
$ \7 x# j) a6 ^
digits(6),x=vpa(x);x=simple(x)
, W: z4 m1 h- d9 G
x=subs(x,{ 't' ,'x21' },{[0:n-1],x21(1:n)})
$ q" [) }) L8 x2 L ^( D; r
xhat=[x(1),diff(x)]
6 o2 @' z1 V) o& h1 ?3 G8 P# m
epsilon=x10-xhat
; X( f9 Y" y3 i; V+ D* w2 R7 n
delta=abs(epsilon./x10)
作者:
lvjianchi
时间:
2013-12-22 22:49
n=length(x10);
! {9 n1 Y7 H, z1 U) K7 J
x11=cumsum(x10)
3 `8 W4 i- K! s- E/ `* D
x21=cumsum(x20)
. Q6 G+ @7 ?# x) o1 I
for i=2:n
1 R' G, m. Y$ m- e
z11(i)=0.5*(x11(i)+x11(i-1));
; [; ]1 [" o5 l8 O2 W- \8 B
end
! x! @# u' `& N( [2 o7 T
B=[-z11(2:n)',x21(2:n)'];
- w5 W: x" I M9 v8 W [* F
Y=x10(2:n)';
' H- K0 `9 P2 f0 Y( B' `
u=B\Y
1 j0 Z1 h# m9 b
x=dsolve('Dx+a*x=b*x2' ,'x(0)=x0');
. Z4 P- h; L0 h
x=subs(x,{ 'a' ,'b' ,'x0','x2'},{u(1),u(2),x10(1), 'x21' });
; y0 q. M3 e; m+ ]/ ?& P: V5 V; ^
digits(6),x=vpa(x);x=simple(x)
* m P8 S/ I" T9 q# a
x=subs(x,{ 't' ,'x21' },{[0:n-1],x21(1:n)})
0 j0 I* C) P. b) ~
xhat=[x(1),diff(x)]
6 y5 f) f8 T) U. S5 p7 n" C
epsilon=x10-xhat
0 L) ^. X7 j9 g, d5 t
delta=abs(epsilon./x10)
作者:
远行的小船儿666
时间:
2015-12-28 16:00
可以看看的
( f8 U5 T P) t: y. B+ T5 g
作者:
远行的小船儿666
时间:
2015-12-28 16:14
好的呀，加油
/ X+ D- P) }! B' w; A
作者:
千梦君
时间:
2019-1-7 14:14
EHWVFHWDUEWGUEIF
! m2 M$ s* _6 _( u) j* `+ s+ u6 l, e) w
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5