| 查看: 577 | 回复: 4 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 15322002 的 5 个金币 | |||
15322002银虫 (小有名气)
|
[交流]
【求助】matlab中的ode45
|
||
|
ode45的求解过程究竟是怎样的呢？我的问题如下，比如matlab的帮助中ode45的例子： function dy = rigid(t,y) dy = zeros(3,1); % a column vector dy(1) = y(2) * y(3); dy(2) = -y(1) * y(3); dy(3) = -0.51 * y(1) * y(2); vector [0 1 1] at time 0.options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]); [T,Y] = ode45(@rigid,[0 12],[0 1 1],options); 如果我想控制Y（1）不能小于0，判断如果小于0就赋值为0，应该怎么实现呢？按照下面的方式改过后没有效果，请大家指教 function dy = rigid(t,y) dy = zeros(3,1); % a column vector dy(1) = y(2) * y(3); if y(1)<0 y(1)=0 end dy(2) = -y(1) * y(3); dy(3) = -0.51 * y(1) * y(2 [ Last edited by nono2009 on 2009-12-29 at 16:28 ] |
回复此楼» 猜你喜欢
面上项目filecode邪修
已经有6人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有14人回复
-
2026国自然放榜时间
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有22人回复
-
这样的filecode谁见过
已经有21人回复
-
filecode与中标关系的预测
已经有5人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
2楼2009-12-29 17:18:11
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
3楼2009-12-29 19:57:11
15322002
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 472.4
- 散金: 429
- 帖子: 139
- 在线: 281.4小时
- 虫号: 709669
- 注册: 2009-02-26
- 专业: 工程热物理相关交叉领域
4楼2009-12-29 20:26:52
15322002
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 472.4
- 散金: 429
- 帖子: 139
- 在线: 281.4小时
- 虫号: 709669
- 注册: 2009-02-26
- 专业: 工程热物理相关交叉领域
5楼2009-12-29 20:28:03
100