| 查看: 3012 | 回复: 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的值第一次就判断好，后来都不会跳转，很头疼。 请问怎么解决？ |
回复此楼» 猜你喜欢
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
申请科研项目纯靠实力能成功吗？
已经有7人回复
-
极限拉扯，如何抉择
已经有14人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 一阶微分方程求解，并优化参数值
已经有14人回复
- 高手进来用matlab求解常微分方程，十分感谢！
已经有17人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
- matlab拟合方程参数时初值的选择
已经有15人回复
- 隐式欧拉法求解一阶常微分方程
已经有7人回复
- 求回归方程参数估计的最优化算法matlab代码
已经有9人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】请教非齐次常微分方程组的解析解法
已经有4人回复
- 【求助】求助一个一阶常微分方程的初值问题
已经有16人回复
- 【求助】偏微分方程和常微分方程组成的方程组怎么求，大家帮帮忙啊
已经有3人回复
- λh方程的参数如何求解（一天之内回答加40个金币）
已经有2人回复
- 【求助】如何从Murnaghan方程的结果中得到优化的晶格参数？
已经有11人回复
- 【求助】不仅控制方程里含有未知参数，边界条件里也有
已经有5人回复
- 【求助】如何用Runge-Kutta迭代求解二阶常微分方程组【已解决】
已经有9人回复
beysem
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 11455.2
- 散金: 3
- 红花: 1
- 帖子: 842
- 在线: 93.5小时
- 虫号: 831634
- 注册: 2009-08-21
- 性别: GG
- 专业: 安全科学与工程
3楼2010-08-28 14:52:46
lijinfeng042
木虫 (小有名气)
Matlab
- 仿真EPI: 2
- 应助: 1 (幼儿园)
- 金币: 2156.1
- 散金: 115
- 帖子: 291
- 在线: 31.5小时
- 虫号: 1019062
- 注册: 2010-05-15
- 性别: GG
- 专业: 通信理论与系统
2楼2010-08-28 13:26:52
li2008jie
铁杆木虫 (著名写手)
逍遥孤鹤
- 应助: 10 (幼儿园)
- 金币: 8646.9
- 红花: 7
- 帖子: 1091
- 在线: 532.1小时
- 虫号: 590822
- 注册: 2008-09-01
- 性别: GG
- 专业: 金属功能材料
4楼2010-08-29 07:58:26
lijinfeng042
木虫 (小有名气)
Matlab
- 仿真EPI: 2
- 应助: 1 (幼儿园)
- 金币: 2156.1
- 散金: 115
- 帖子: 291
- 在线: 31.5小时
- 虫号: 1019062
- 注册: 2010-05-15
- 性别: GG
- 专业: 通信理论与系统
5楼2010-08-29 19:22:10
10
