查看: 11395|回复: 8

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

[复制链接]

字体大小: 正常放大

hujet

8 主题	4 听众	99 积分

升级 98.95%

TA的每日心情

	擦汗 2013-4-5 11:27

签到天数: 38 天

[LV.5]常住居民I

群组: 2012第三期美赛培训

群组: 数学建摸协会

群组: 学术交流A

群组: 学术交流B

群组: 学术交流C

电梯直达

1^#

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

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

30体力

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

zan

源程序

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

相关帖子

使用道具举报

chen7705

2 主题	5 听众	80 积分

升级 78.95%

TA的每日心情

	开心 2012-10-28 21:16

签到天数: 18 天

[LV.4]偶尔看看III

群组: 学术交流A

2^#

发表于 2012-8-27 13:13 |只看该作者 |招呼Ta 关注Ta

对,好多都没灰色预测GM(1,N)源程序

我也说一句

发表

使用道具举报

5105

0 主题	7 听众	42 积分

升级 38.95%

TA的每日心情

	慵懒 2014-2-10 00:58

签到天数: 12 天

[LV.3]偶尔看看II

3^#

发表于 2013-2-1 19:37 |只看该作者 |招呼Ta 关注Ta |邮箱已经成功绑定

同求啊啊啊……

我也说一句

发表

使用道具举报

栀子流年

2 主题	8 听众	165 积分

升级 32.5%

TA的每日心情

	奋斗 2014-1-10 08:49

签到天数: 46 天

[LV.5]常住居民I

自我介绍: 开朗，勤奋

4^#

发表于 2013-11-18 22:31 |只看该作者 |招呼Ta 关注Ta

+1求程序啊

我也说一句

发表

使用道具举报

lvjianchi

2 主题	4 听众	235 积分

升级 67.5%

TA的每日心情

	怒 2014-9-12 23:04

签到天数: 83 天

[LV.6]常住居民II

群组: 西安交大数学建模

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

群组: 2012第三期美赛培训

群组: 学术交流A

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

5^#

发表于 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)

我也说一句

发表

使用道具举报

lvjianchi

2 主题	4 听众	235 积分

升级 67.5%

TA的每日心情

	怒 2014-9-12 23:04

签到天数: 83 天

[LV.6]常住居民II

群组: 西安交大数学建模

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

群组: 2012第三期美赛培训

群组: 学术交流A

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

6^#

发表于 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)