数学建模社区-数学中国
标题:
求多因子灰色预测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);
% v/ K0 f8 ~7 D& L8 |! l
x11=cumsum(x10)
* f# z2 {4 F" z6 Y! j" ] G
x21=cumsum(x20)
' O. Z& y6 j; V+ R
for i=2:n
2 T& z5 C0 F9 \5 n# w
z11(i)=0.5*(x11(i)+x11(i-1));
. Y5 V, S% T6 A8 f
end
( F6 `; C/ v8 b1 s
B=[-z11(2:n)',x21(2:n)'];
9 O$ t' |) O, j7 m; w2 B
Y=x10(2:n)';
' P, }: g' `5 N C& Z; f
u=B\Y
6 P! R: Z) x7 D& A# S+ p6 H
x=dsolve('Dx+a*x=b*x2' ,'x(0)=x0');
! ~7 I/ x" n( h j M
x=subs(x,{ 'a' ,'b' ,'x0','x2'},{u(1),u(2),x10(1), 'x21' });
4 q) ^9 h. q- T6 A, P) O; s
digits(6),x=vpa(x);x=simple(x)
& b8 E) \/ N8 K. o4 k: _
x=subs(x,{ 't' ,'x21' },{[0:n-1],x21(1:n)})
& J+ k( V4 y7 \- O, |
xhat=[x(1),diff(x)]
0 m1 r$ A v9 g% U5 t* \! F6 f
epsilon=x10-xhat
. j) g0 n0 S4 q% Y
delta=abs(epsilon./x10)
作者:
lvjianchi
时间:
2013-12-22 22:49
n=length(x10);
9 u Z5 z8 R) r4 G! S* m0 L
x11=cumsum(x10)
$ o. B& k% R- X
x21=cumsum(x20)
# O' X: Z2 @$ T# N* a# c( l) `$ c
for i=2:n
" B- ^# Q t8 W5 a% X, t
z11(i)=0.5*(x11(i)+x11(i-1));
3 J* I( X6 P' o# |6 g
end
. j7 O: b; k$ o
B=[-z11(2:n)',x21(2:n)'];
. }) L: D' A* u$ D& n
Y=x10(2:n)';
. o* L" N$ K0 ?$ ?1 G7 p" i
u=B\Y
% ~- m8 d% J" [/ ]$ [0 g' p
x=dsolve('Dx+a*x=b*x2' ,'x(0)=x0');
. B' Q$ z: L* E) x: a! w0 e
x=subs(x,{ 'a' ,'b' ,'x0','x2'},{u(1),u(2),x10(1), 'x21' });
, m) s# M' n" L" s l: \- x5 _ r
digits(6),x=vpa(x);x=simple(x)
8 Z8 x3 J' D2 F) ^
x=subs(x,{ 't' ,'x21' },{[0:n-1],x21(1:n)})
: e% A0 a2 g4 c& _* z
xhat=[x(1),diff(x)]
5 ?; B' r& N* e
epsilon=x10-xhat
! [6 ^( k2 a$ H5 F3 u6 [( C
delta=abs(epsilon./x10)
作者:
远行的小船儿666
时间:
2015-12-28 16:00
可以看看的
# r8 g2 \) c x3 j: T9 u. x# ~6 B
作者:
远行的小船儿666
时间:
2015-12-28 16:14
好的呀，加油
: u" d- p5 |9 y* R% Q8 y3 z7 K+ j
作者:
千梦君
时间:
2019-1-7 14:14
EHWVFHWDUEWGUEIF
4 L+ a5 o9 { E2 m
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5