求多因子灰色预测GM(1,N)源程序
求多因子灰色预测GM(1,N)源程序,不是GM(1,1)的,我要求的是GM(1,4)可是不知道怎么弄程序,谢谢了 对,好多都没灰色预测GM(1,N)源程序 同求啊啊啊…… +1求程序啊 n=length(x10);x11=cumsum(x10)
x21=cumsum(x20)
for i=2:n
z11(i)=0.5*(x11(i)+x11(i-1));
end
B=[-z11(2:n)',x21(2:n)'];
Y=x10(2:n)';
u=B\Y
x=dsolve('Dx+a*x=b*x2' ,'x(0)=x0');
x=subs(x,{ 'a' ,'b' ,'x0','x2'},{u(1),u(2),x10(1), 'x21' });
digits(6),x=vpa(x);x=simple(x)
x=subs(x,{ 't' ,'x21' },{,x21(1:n)})
xhat=
epsilon=x10-xhat
delta=abs(epsilon./x10) n=length(x10);
x11=cumsum(x10)
x21=cumsum(x20)
for i=2:n
z11(i)=0.5*(x11(i)+x11(i-1));
end
B=[-z11(2:n)',x21(2:n)'];
Y=x10(2:n)';
u=B\Y
x=dsolve('Dx+a*x=b*x2' ,'x(0)=x0');
x=subs(x,{ 'a' ,'b' ,'x0','x2'},{u(1),u(2),x10(1), 'x21' });
digits(6),x=vpa(x);x=simple(x)
x=subs(x,{ 't' ,'x21' },{,x21(1:n)})
xhat=
epsilon=x10-xhat
delta=abs(epsilon./x10) 可以看看的
好的呀,加油
EHWVFHWDUEWGUEIF
页:
[1]