| 查看: 2819 | 回复: 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的值第一次就判断好，后来都不会跳转，很头疼。 请问怎么解决？ |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有17人回复
-
为什么 国际(地区)合作与交流项目 没有放榜？
已经有11人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有4人回复
-
让我中一个面上吧！
已经有17人回复
-
基金不中，共勉
已经有6人回复
-
国自然面上复盘~欢迎讨论
已经有8人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
我不理解！
已经有23人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 一阶微分方程求解，并优化参数值
已经有14人回复
- 高手进来用matlab求解常微分方程，十分感谢！
已经有17人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
- matlab拟合方程参数时初值的选择
已经有15人回复
- 隐式欧拉法求解一阶常微分方程
已经有7人回复
- 求回归方程参数估计的最优化算法matlab代码
已经有9人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】请教非齐次常微分方程组的解析解法
已经有4人回复
- 【求助】求助一个一阶常微分方程的初值问题
已经有16人回复
- 【求助】偏微分方程和常微分方程组成的方程组怎么求，大家帮帮忙啊
已经有3人回复
- λh方程的参数如何求解（一天之内回答加40个金币）
已经有2人回复
- 【求助】如何从Murnaghan方程的结果中得到优化的晶格参数？
已经有11人回复
- 【求助】不仅控制方程里含有未知参数，边界条件里也有
已经有5人回复
- 【求助】如何用Runge-Kutta迭代求解二阶常微分方程组【已解决】
已经有9人回复
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
xiegangmai(金币+1): 谢谢参与！ 2011-12-16 21:46:57
小木虫(金币+0.5):给个红包，谢谢回帖
xiegangmai(金币+1): 谢谢参与！ 2011-12-16 21:46:57
|
[t,y]=ode45('appollo2',[0,20],[1,-30,1,1]); plot(t,y(:,3));figure;plot(t,y(:,4)) 上为调用程序 %…………………………………………………………………… function dx=appollo2(t,x) mu=1/82.45*(x(2)<0)-1*(x(2)>0); %这是关键 mustar=1-mu; r1=sqrt((x(1)+mu)^2+x(3)^2); r2=sqrt((x(1)-mustar)^2+x(3)^2); dx=[x(2); 2*x(4)+x(1)-mustar*(x(1)+mu)/r1^3-mu*(x(1)-mustar)/r2^3; x(4); -2*x(2)+x(3)-mustar*x(3)/r1^3-mu*x(3)/r2^3]; |
8楼2011-12-16 16:36:38
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 (幼儿园)
- 金币: 8670.4
- 红花: 7
- 帖子: 1090
- 在线: 531.9小时
- 虫号: 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
free--fire
新虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 543.5
- 散金: 570
- 帖子: 199
- 在线: 51.7小时
- 虫号: 1051983
- 注册: 2010-07-05
- 性别: GG
- 专业: 力学中的基本问题和方法
6楼2010-08-30 08:58:31
7楼2011-12-16 15:40:40
9楼2011-12-16 16:37:05
captainmeng
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 291.1
- 散金: 281
- 帖子: 223
- 在线: 66.6小时
- 虫号: 2615493
- 注册: 2013-08-25
- 性别: GG
- 专业: 流体力学
10楼2013-11-22 00:04:25
投票:
