| 查看: 625 | 回复: 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 ] |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有14人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有16人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有15人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有17人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有21人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有21人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
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
50