数学建模社区-数学中国»论坛 › 【数模论坛】资源分类 › 数模问题互助 › 求解一个五阶微分方程组，用ode命令失效怎么办？

查看: 1361|回复: 2

[问题求助] 求解一个五阶微分方程组，用ode命令失效怎么办？

[复制链接]

字体大小: 正常放大

Reader_Founder

52 主题	6 听众	495 积分

升级 65%

TA的每日心情

	开心 2014-7-12 17:09

签到天数: 115 天

[LV.6]常住居民II

自我介绍: 数学本科

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

电梯直达

1^#

发表于 2013-8-27 13:11 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

画的图显示出各因变量随时间不变，怎么处理?
在线等，谢谢~

zan

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

使用道具举报

甘乐桑

3 主题	6 听众	163 积分

升级 31.5%

TA的每日心情

	开心 2013-11-29 16:30

签到天数: 51 天

[LV.5]常住居民I

2^#

发表于 2013-8-27 15:44 |只看该作者 |招呼Ta 关注Ta

本帖最后由甘乐桑于 2013-8-27 15:45 编辑

方程贴出来。
可能是初始条件给的不好，会求出常数也很正常。。

点评

Reader_Founder function dy=cb(t,y); dy=[-0.01*y(1)*y(4); 0.01*y(1)*y(4)-0.021*y(2); 0.021*y(2)-0.55*y(3); 0.5*y(3)-0.1*y(4); 0.55*y(3)+0.1*y(4)]; end clear; ts=0:3:72; x0=[100000000, 0, 0, 100, 0]; [t, 详情回复发表于 2013-8-27 18:47

我也说一句

发表

使用道具举报

Reader_Founder

52 主题	6 听众	495 积分

升级 65%

TA的每日心情

	开心 2014-7-12 17:09

签到天数: 115 天

[LV.6]常住居民II

自我介绍: 数学本科

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

3^#

发表于 2013-8-27 18:47 |只看该作者 |招呼Ta 关注Ta

甘乐桑发表于 2013-8-27 15:44
方程贴出来。
可能是初始条件给的不好，会求出常数也很正常。。

function dy=cb(t,y);

dy=[-0.01*y(1)*y(4); 0.01*y(1)*y(4)-0.021*y(2); 0.021*y(2)-0.55*y(3); 0.5*y(3)-0.1*y(4); 0.55*y(3)+0.1*y(4)];

end

clear;
ts=0:3:72;
x0=[100000000, 0, 0, 100, 0];
[t,y]=ode15s('cb',ts,x0);

我也说一句

发表

使用道具举报

返回列表

帐号		密码		只需要一步，快速开始		注册地址	找回密码

[问题求助] 求解一个五阶微分方程组，用ode命令失效怎么办？

点评

浏览过的版块

QQ

电话咨询

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