| 查看: 2953 | 回复: 12 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
beysem金虫 (正式写手)
|
[交流]
【求助】常微分方程组求解中系数与某变量值关联的问题
|
||
|
以常见的Appollo的 微分方程组为例，其function文件如下 11. function dx=appollo(t,x) % 确定几个待定系数 mu/mustar/r1/r2 12. mu=1/82.45; 13. mustar=1-mu; 14. r1=sqrt((x(1)+mu)^2+x(3)^2); 15. r2=sqrt((x(1)-mustar)^2+x(3)^2); % 4个微分方程的表达式 16. dx=[x(2) 17. 2*x(4)+x(1)-mustar*(x(1)+mu)/r1^3-mu*(x(1)-mustar)/r2^3 18. x(4) 19. -2*x(2)+x(3)-mustar*x(3)/r1^3-mu*x(3)/r2^3]; 现在我的问题是 当x(2) 小于 0 时，待定系数为 mu=1/82.45； 现假设： 调用 ode45时，如果 x(2) 的值出现大于0时， 实际需要mu=-1， 我在尝试类似问题时，发现 mu的一旦定义了，哪怕是 if... 语句中， mu的值第一次就判断好，后来都不会跳转，很头疼。 请问怎么解决？ |
回复此楼» 猜你喜欢
现代”学阀”该如何界定
已经有12人回复
-
国社科系统bug了，是不是要放榜了？
已经有10人回复
-
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有4人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有3人回复
-
上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告
已经有11人回复
-
我的奶奶
已经有3人回复
信彼南山
木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 4142.9
- 散金: 1221
- 红花: 16
- 帖子: 1178
- 在线: 233.5小时
- 虫号: 1133529
- 注册: 2010-10-27
- 专业: 导航、制导与传感技术
11楼2013-11-25 20:05:31
lijinfeng042
木虫 (小有名气)
Matlab
- 仿真EPI: 2
- 应助: 1 (幼儿园)
- 金币: 2156.1
- 散金: 115
- 帖子: 291
- 在线: 31.5小时
- 虫号: 1019062
- 注册: 2010-05-15
- 性别: GG
- 专业: 通信理论与系统
2楼2010-08-28 13:26:52
beysem
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 11455.2
- 散金: 3
- 红花: 1
- 帖子: 842
- 在线: 93.5小时
- 虫号: 831634
- 注册: 2009-08-21
- 性别: GG
- 专业: 安全科学与工程
3楼2010-08-28 14:52:46
li2008jie
铁杆木虫 (著名写手)
逍遥孤鹤
- 应助: 10 (幼儿园)
- 金币: 8644.4
- 红花: 7
- 帖子: 1091
- 在线: 532.1小时
- 虫号: 590822
- 注册: 2008-09-01
- 性别: GG
- 专业: 金属功能材料
4楼2010-08-29 07:58:26