数学建模社区-数学中国
标题:
求多因子灰色预测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);
" q7 v9 ]2 g( q% g
x11=cumsum(x10)
; b8 P% G+ b. u: s/ z u) V
x21=cumsum(x20)
# F0 Z+ N1 X: H+ l" U9 R' A$ m
for i=2:n
3 L2 V5 C5 _6 V4 L/ u) n
z11(i)=0.5*(x11(i)+x11(i-1));
$ K/ _ q8 \! r9 P
end
! N6 X2 h0 q4 f( ^& ^) }% P
B=[-z11(2:n)',x21(2:n)'];
3 n% f8 a, |' Y. t
Y=x10(2:n)';
1 ]4 l6 K( e* c$ S) O
u=B\Y
9 x- Y. u' B: q' ^( k
x=dsolve('Dx+a*x=b*x2' ,'x(0)=x0');
) Q: e, @7 f5 Z& ~5 P
x=subs(x,{ 'a' ,'b' ,'x0','x2'},{u(1),u(2),x10(1), 'x21' });
' {5 q ?2 R8 `( z6 G" t! y: ~
digits(6),x=vpa(x);x=simple(x)
+ V; l- G; ]# O
x=subs(x,{ 't' ,'x21' },{[0:n-1],x21(1:n)})
1 @/ \: f: x% o; f; ^+ g% F' o I2 m
xhat=[x(1),diff(x)]
7 |+ W8 z1 _4 M, U# q( D
epsilon=x10-xhat
+ d( @) G- s5 N; ~
delta=abs(epsilon./x10)
作者:
lvjianchi
时间:
2013-12-22 22:49
n=length(x10);
2 p( R3 k% h6 s, A6 I
x11=cumsum(x10)
) G8 l8 g; z) F3 f
x21=cumsum(x20)
- [# Q& L8 B6 K' h' _8 N
for i=2:n
* v+ B- o, L' R ?! ^+ {0 H
z11(i)=0.5*(x11(i)+x11(i-1));
+ K C# {# [$ ]: M2 I. b" m
end
0 u7 s1 \1 j4 @6 h! Q: u4 H% e
B=[-z11(2:n)',x21(2:n)'];
( u4 i" V( s) \; F+ m W; [- Z# A
Y=x10(2:n)';
, H2 S8 ~- n: y9 Z7 O: i
u=B\Y
" ^( G& }- T3 V* {$ j, w
x=dsolve('Dx+a*x=b*x2' ,'x(0)=x0');
4 Z, g; `8 E. M5 L& I4 I$ p
x=subs(x,{ 'a' ,'b' ,'x0','x2'},{u(1),u(2),x10(1), 'x21' });
, }4 Q) W% t& k! k0 `1 E9 D0 t
digits(6),x=vpa(x);x=simple(x)
6 ~' U1 a7 W2 h7 B+ [
x=subs(x,{ 't' ,'x21' },{[0:n-1],x21(1:n)})
: t- `5 e; D- ~! Z+ z
xhat=[x(1),diff(x)]
' z; U+ I4 e1 h" I7 [ D/ T
epsilon=x10-xhat
3 l3 W1 y& A \' \! A4 }
delta=abs(epsilon./x10)
作者:
远行的小船儿666
时间:
2015-12-28 16:00
可以看看的
; L+ \& w L( @) |5 H
作者:
远行的小船儿666
时间:
2015-12-28 16:14
好的呀，加油
# I: J1 e0 Q C- G
作者:
千梦君
时间:
2019-1-7 14:14
EHWVFHWDUEWGUEIF
& j" [* }/ w4 j- g2 W
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5