数学建模社区-数学中国»论坛 › 算法论坛 › 灰色算法 › 灰色代码

发新帖

查看: 4014|回复: 5

上一主题

下一主题

[代码资源] 灰色代码

字体大小: 正常放大

pmdswjtu

3 主题	4 听众	25 积分

升级 21.05%

TA的每日心情

	开心 2014-9-22 18:41

签到天数: 5 天

[LV.2]偶尔看看I

电梯直达

跳转到指定楼层

1^#

发表于 2012-2-7 22:52 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

clc,clear
x0=[87 149 100 26 129];
n=length(x0);
x1=cumsum(x0)
a_x0=diff(x0);
a_x0=[0,a_x0]
for i=2:n
z(i)=0.5*(x1(i)+x1(i-1));
end
B=[-x0(2:end)',-z(2:end)',ones(n-1,1)];
Y=a_x0(2:end)';
u=B\Y
x=dsolve('D2x+a1*Dx+a2*x=b','x(0)=c1,x(5)=c2');
x=subs(x,{'a1','a2','b','c1','c2'},{u(1),u(2),u(3),x1(1),x1(5)});
yuce=subs(x,'t',0:n-1);
digits(6),x=vpa(x)
x0_hat=[yuce(1),diff(yuce)]
epsilon=x0-x0_hat
delta=abs(epsilon./x0)

zan

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

回复

使用道具举报

alair006

头像被屏蔽

0 主题	4 听众	558 积分

升级 86%

TA的每日心情

	擦汗 2012-2-8 08:16

签到天数: 4 天

[LV.2]偶尔看看I

2^#

发表于 2012-2-8 08:19 |只看该作者 |招呼Ta 关注Ta

提示: 作者被禁止或删除内容自动屏蔽

回复

使用道具举报

nuanyang

3 主题	5 听众	127 积分

升级 13.5%

TA的每日心情

	奋斗 2012-2-12 14:00

签到天数: 9 天

[LV.3]偶尔看看II

2012挑战赛参赛者

群组: 2012第三期美赛培训

3^#

发表于 2012-2-8 08:34 |只看该作者 |招呼Ta 关注Ta

这是干什么用的？

我也说一句

发表

回复

使用道具举报

alair006

头像被屏蔽

0 主题	4 听众	558 积分

升级 86%

TA的每日心情

	擦汗 2012-2-8 08:16

签到天数: 4 天

[LV.2]偶尔看看I

4^#

发表于 2012-2-8 10:39 |只看该作者 |招呼Ta 关注Ta

提示: 作者被禁止或删除内容自动屏蔽

回复

使用道具举报

467857726de

3 主题	7 听众	179 积分

升级 39.5%

TA的每日心情

	擦汗 2015-1-22 15:49

签到天数: 21 天

[LV.4]偶尔看看III

国际赛参赛者

自我介绍: 我是本科在校学生，从小就喜欢数学，希望可以在这里学到想要的东西。。。

群组: 学术交流B

5^#

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

不错哦谢谢楼主

我也说一句

发表

回复

使用道具举报

实名认证

1 主题	9 听众	1747 积分

TA的每日心情

	开心 2016-7-26 21:58

签到天数: 182 天

[LV.7]常住居民III

群组: 2014年美赛冲刺培训

群组: 数学建模培训课堂1

群组: 物联网工程师培训

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

6^#

发表于 2014-3-9 15:17 |只看该作者 |招呼Ta 关注Ta

非常感谢楼主的福利！

我也说一句

发表

回复

使用道具举报

发新帖

qq

电话咨询
04714969085

fastpost

关于我们| 联系我们| 诚征英才| 对外合作| 产品服务|

手机版|Archiver| |繁體中文手机客户端

蒙公网安备 15010502000194号

Powered by Discuz! X2.5 © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 ) 论坛法律顾问：王兆丰

GMT+8, 2026-4-25 18:07 , Processed in 0.515313 second(s), 78 queries .

回顶部