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

请问怎么解决？

回复此楼

» 猜你喜欢

回收溶剂求助已经有7人回复
职称评审没过，求安慰已经有40人回复
硝基苯如何除去已经有3人回复
A期刊撤稿已经有4人回复
垃圾破二本职称评审标准已经有17人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
申请26博士已经有5人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
求助文献已经有3人回复

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

一阶微分方程求解，并优化参数值已经有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的回帖

回帖支持 ( 显示支持度最高的前 50 名 )

WF1987

铜虫 (小有名气)

应助: 1 (幼儿园)
金币: 87.3
帖子: 90
在线: 20.6小时
虫号: 1338224
注册: 2011-07-05
性别: GG
专业: 职业卫生

★ ★
小木虫(金币+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];

赞一下(3人)

回复此楼

8楼2011-12-16 16:36:38

已阅回复此楼关注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 三个参数吧呵呵这个看得挺累的.....

赞一下(1人)

回复此楼

工作了,偶尔会上来~可以关注新浪微博 @云是风的梦_Matlab