查看: 11413|回复: 8

[问题求助] 求多因子灰色预测GM(1,N)源程序

[复制链接]

字体大小: 正常放大

hujet

8 主题	4 听众	99 积分

升级 98.95%

TA的每日心情

	擦汗 2013-4-5 11:27

签到天数: 38 天

[LV.5]常住居民I

群组: 2012第三期美赛培训

群组: 数学建摸协会

群组: 学术交流A

群组: 学术交流B

群组: 学术交流C

电梯直达

^#

发表于 2012-5-6 15:43 |只看该作者 |正序浏览

|招呼Ta 关注Ta |邮箱已经成功绑定

30体力

求多因子灰色预测GM(1,N)源程序，不是GM(1,1）的，我要求的是GM(1,4)可是不知道怎么弄程序，谢谢了

zan

源程序

转播0 淘帖0 分享0 收藏0 支持0 反对0 微信

相关帖子

使用道具举报

千梦君

0 主题	2 听众	40 积分

升级 36.84%

TA的每日心情

	开心 2018-3-7 21:12

签到天数: 1 天

[LV.1]初来乍到

自我介绍: 暂无

群组: 2018MCM科技论文写作

8^#

发表于 2019-1-7 14:14 |只看该作者 |招呼Ta 关注Ta

EHWVFHWDUEWGUEIF

我也说一句

发表

使用道具举报

远行的小船儿666

14 主题	10 听众	1262 积分

升级 26.2%

TA的每日心情

	开心 2017-11-9 21:17

签到天数: 198 天

[LV.7]常住居民III

网络挑战赛参赛者

群组: 大学生软件交流群

群组: 2016国赛备战群组

7^#

发表于 2015-12-28 16:14 |只看该作者 |招呼Ta 关注Ta |邮箱已经成功绑定

好的呀，加油

我也说一句

发表

使用道具举报

远行的小船儿666

14 主题	10 听众	1262 积分

升级 26.2%

TA的每日心情

	开心 2017-11-9 21:17

签到天数: 198 天

[LV.7]常住居民III

网络挑战赛参赛者

群组: 大学生软件交流群

群组: 2016国赛备战群组

6^#

发表于 2015-12-28 16:00 |只看该作者 |招呼Ta 关注Ta |邮箱已经成功绑定

可以看看的

我也说一句

发表

使用道具举报

lvjianchi

2 主题	4 听众	235 积分

升级 67.5%

TA的每日心情

	怒 2014-9-12 23:04

签到天数: 83 天

[LV.6]常住居民II

群组: 西安交大数学建模

群组: 中国矿业大学数模培训

群组: 2012第三期美赛培训

群组: 学术交流A

群组: 2014年网络挑战赛交流

5^#

发表于 2013-12-22 22:49 |只看该作者 |招呼Ta 关注Ta

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' },{[0:n-1],x21(1:n)})
xhat=[x(1),diff(x)]
epsilon=x10-xhat
delta=abs(epsilon./x10)

我也说一句

发表

使用道具举报

lvjianchi

2 主题	4 听众	235 积分

升级 67.5%

TA的每日心情

	怒 2014-9-12 23:04

签到天数: 83 天

[LV.6]常住居民II

群组: 西安交大数学建模

群组: 中国矿业大学数模培训

群组: 2012第三期美赛培训

群组: 学术交流A

群组: 2014年网络挑战赛交流

4^#

发表于 2013-12-22 22:48 |只看该作者 |招呼Ta 关注Ta

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' },{[0:n-1],x21(1:n)})
xhat=[x(1),diff(x)]
epsilon=x10-xhat
delta=abs(epsilon./x10)