数学建模社区-数学中国
标题:
求多因子灰色预测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 I
for i=2:n
: d/ @' Y/ z* g, t
z11(i)=0.5*(x11(i)+x11(i-1));
) d; X% h) b+ O
end
1 r& J2 b) x' y* Z7 I
B=[-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 R
x=dsolve('Dx+a*x=b*x2' ,'x(0)=x0');
7 }* r+ f- ~; ?7 V# P- g
x=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 p
xhat=[x(1),diff(x)]
6 U4 g3 {" I" P
epsilon=x10-xhat
& _; Y. B# s! c$ z U4 ~% z
delta=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 s
end
" w o, c3 [4 i; l) R
B=[-z11(2:n)',x21(2:n)'];
' \% t3 q/ C) A
Y=x10(2:n)';
' j0 F ~0 s* D7 A
u=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 N
x=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