请教高手此微分方程组如何求解

sswwss · 发表于 2008-12-7 16:51

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

先谢谢大家！！

方程中θ和φ是变量, 如果它们的系数是常数我会解，但是方程中有θ 和φ各阶导数的乘积项
向各位请教如何求解？
谢谢！！！

ch_j1985 · 发表于 2008-12-7 20:50

应该要进行解耦吧

whatisi123 · 发表于 2008-12-8 14:27

你解决了吗？我也是遇到了类似的问题？不知可否留下联系方式？
我的QQ：277930401:loveliness:

sswwss · 发表于 2008-12-8 15:05

还没解决
好像也不好解耦

sogooda · 发表于 2008-12-8 16:48

先变量代换，x1=θ,x2=θ',x3=φ,x4=φ',之后就可以写出x'=f(x)的形式，然后就可以用ode45求解了。

whatisi123 · 发表于 2008-12-8 16:57

我是按照下面的高手做的但是在解决问题的时候会出现这样的警告：
Warning: Failure at t=5.780352e+000. Unable to meet integration tolerances without
reducing the step size below the smallest value allowed (1.421085e-014) at time t.
不是很明白？:@Q

whatisi123 · 发表于 2008-12-8 17:01

如果上面方程的系数是关于θ,θ',φ,φ',的函数，那么“先变量代换，x1=θ,x2=θ',x3=φ,x4=φ',之后就可以写出x'=f(x)的形式，然后就可以用ode45求解了。”可以这样解决吗？

sogooda · 发表于 2008-12-8 21:13

如果是显函数的话应该就可以。如果不行，可以再试试其他求解器，比如ode15s什么的。

[ 本帖最后由 sogooda 于 2008-12-8 21:14 编辑 ]

yufeng · 发表于 2008-12-9 08:55

我只能得到数值解

sogooda · 发表于 2008-12-9 09:35

楼主应该是要求数值解吧？虽然也可以先用dsolve试一试，不过感觉上能够求出解析解的可能性不大啊。

sswwss · 发表于 2008-12-9 15:04

谢谢各位！！
我用ode45得到了数值解

账号		自动登录	找回密码
密码			我要加入

[综合讨论] 请教高手此微分方程组如何求解

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复楼主 sswwss 的帖子

回复楼主 sswwss 的帖子

回复楼主 sswwss 的帖子

回复地板 sswwss 的帖子

回复 7楼 whatisi123 的帖子

回复 9楼 yufeng 的帖子

浏览过的版块

[综合讨论] 请教高手此微分方程组如何求解

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复 楼主 sswwss 的帖子

回复 楼主 sswwss 的帖子

回复 楼主 sswwss 的帖子

回复 地板 sswwss 的帖子

回复 7楼 whatisi123 的帖子

回复 9楼 yufeng 的帖子

浏览过的版块

回复楼主 sswwss 的帖子

回复楼主 sswwss 的帖子

回复楼主 sswwss 的帖子

回复地板 sswwss 的帖子