24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

beysem

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 11455.2
散金: 3
红花: 1
帖子: 842
在线: 93.5小时
虫号: 831634
注册: 2009-08-21
性别: GG
专业: 安全科学与工程

[交流] 【求助】常微分方程组求解中系数与某变量值关联的问题

以常见的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的值第一次就判断好，后来都不会跳转，很头疼。

请问怎么解决？

回复此楼

» 猜你喜欢

北京211副教授，35岁，想重新出发，去国外做博后，怎么样？已经有3人回复
自荐读博已经有3人回复
求助:我三月中下旬出站，青基依托单位怎么办？已经有5人回复
论文终于录用啦！满足毕业条件了已经有22人回复
不自信的我已经有5人回复
磺酰氟产物，毕不了业了！已经有4人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复

» 本主题相关价值贴推荐，对您同样有帮助:

一阶微分方程求解，并优化参数值已经有14人回复
高手进来用matlab求解常微分方程，十分感谢！已经有17人回复
matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较已经有4人回复
matlab拟合方程参数时初值的选择已经有15人回复
隐式欧拉法求解一阶常微分方程已经有7人回复
求回归方程参数估计的最优化算法matlab代码已经有9人回复
【求助】matlab 求解微分方程中的未知参数已经有20人回复
【求助】请教非齐次常微分方程组的解析解法已经有4人回复
【求助】求助一个一阶常微分方程的初值问题已经有16人回复
【求助】偏微分方程和常微分方程组成的方程组怎么求，大家帮帮忙啊已经有3人回复
λh方程的参数如何求解（一天之内回答加40个金币）已经有2人回复
【求助】如何从Murnaghan方程的结果中得到优化的晶格参数？已经有11人回复
【求助】不仅控制方程里含有未知参数，边界条件里也有已经有5人回复
【求助】如何用Runge-Kutta迭代求解二阶常微分方程组【已解决】已经有9人回复

1楼 2010-08-28 08:59:20

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lijinfeng042

木虫 (小有名气)

Matlab

仿真EPI: 2
应助: 1 (幼儿园)
金币: 2156.1
散金: 115
帖子: 291
在线: 31.5小时
虫号: 1019062
注册: 2010-05-15
性别: GG
专业: 通信理论与系统

★
ctgu_zheng(金币+1):鼓励应助 2010-08-28 16:48:10
beysem(金币+1):感谢关注，但我的问题还是没有解决哦 2010-08-29 07:48:36

引用回帖:

Originally posted by beysem at 2010-08-28 08:59:20:
以常见的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. ...

这是你转换好的方程式吧那原来的方程组呢 u r1 r2 三个参数吧呵呵这个看得挺累的.....