| 查看: 497 | 回复: 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 ] |
回复此楼» 猜你喜欢
复试调剂
已经有7人回复
-
一志愿北京化工085600 310分求调剂
已经有13人回复
-
085600，专业课化工原理，320分求调剂
已经有10人回复
-
288求调剂
已经有5人回复
-
0703化学321分求调剂
已经有16人回复
-
318求调剂
已经有13人回复
-
一志愿南航，数一英一学硕317求调剂！！
已经有5人回复
-
一志愿9材料学硕297已过六级求调剂推荐
已经有11人回复
-
085602调剂 初试总分335
已经有5人回复
-
302分 085601求调剂推荐
已经有8人回复
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
3楼2009-12-29 19:57:11
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
2楼2009-12-29 17:18: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
