数学建模社区-数学中国
标题:
求多因子灰色预测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);
0 t# E* z2 \, A* o) U1 c) H" {6 L
x11=cumsum(x10)
4 s; h' J* Q# n& U3 q
x21=cumsum(x20)
( s1 A5 M; R+ z, i/ A
for i=2:n
9 B& g! g+ o! L
z11(i)=0.5*(x11(i)+x11(i-1));
( z) g2 ^2 F- |% ^
end
/ n( J8 V5 j+ Q, ?, C! K4 Z9 [
B=[-z11(2:n)',x21(2:n)'];
% u B" D, @: n' H
Y=x10(2:n)';
$ N7 r( B+ ^: K8 F
u=B\Y
+ I( B3 G5 B' S) f. n
x=dsolve('Dx+a*x=b*x2' ,'x(0)=x0');
' ~2 E) X4 J2 S! ]5 u
x=subs(x,{ 'a' ,'b' ,'x0','x2'},{u(1),u(2),x10(1), 'x21' });
* J! y! \; W3 n/ X: _! Q# w
digits(6),x=vpa(x);x=simple(x)
7 j. M5 _8 l+ n6 D! f# _/ {- e
x=subs(x,{ 't' ,'x21' },{[0:n-1],x21(1:n)})
; W Z t" m6 }& ]
xhat=[x(1),diff(x)]
& K# o( Z* r7 u0 g0 A
epsilon=x10-xhat
- d2 ]; [1 d2 m/ O* `1 D
delta=abs(epsilon./x10)
作者:
lvjianchi
时间:
2013-12-22 22:49
n=length(x10);
$ q( C$ Z+ L9 d( o
x11=cumsum(x10)
7 a& X' ~# W( {
x21=cumsum(x20)
& D2 P/ F4 r7 p% z0 v
for i=2:n
9 {. Y0 h" r6 L4 f. N' P7 x' x
z11(i)=0.5*(x11(i)+x11(i-1));
2 x+ u/ R4 f0 S& H* D# V
end
# f: y" A. l/ `* G5 C
B=[-z11(2:n)',x21(2:n)'];
5 I8 K( i3 G( P4 H0 s
Y=x10(2:n)';
" S" v% K& g) L" `
u=B\Y
+ h6 J# F5 y5 k
x=dsolve('Dx+a*x=b*x2' ,'x(0)=x0');
1 ? @8 @8 a3 P# r) E
x=subs(x,{ 'a' ,'b' ,'x0','x2'},{u(1),u(2),x10(1), 'x21' });
9 Y+ z8 |" n5 W
digits(6),x=vpa(x);x=simple(x)
0 K* R; P& q9 d4 k! d% a( C
x=subs(x,{ 't' ,'x21' },{[0:n-1],x21(1:n)})
( a9 ~, ~5 l+ M$ @3 {
xhat=[x(1),diff(x)]
# y* N7 n( G% [6 o4 N
epsilon=x10-xhat
8 G) B$ D: N3 B" g4 }* @) i
delta=abs(epsilon./x10)
作者:
远行的小船儿666
时间:
2015-12-28 16:00
可以看看的
! A# R. {7 m2 e' E8 {& n
作者:
远行的小船儿666
时间:
2015-12-28 16:14
好的呀，加油
! n0 C! J4 R0 \! l
作者:
千梦君
时间:
2019-1-7 14:14
EHWVFHWDUEWGUEIF
1 h* f# p! X/ d7 W: k- I7 E5 c
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5