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

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

先谢谢大家！！

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

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

回复 楼主 sswwss 的帖子

应该要进行解耦吧

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

回复 楼主 sswwss 的帖子

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

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

还没解决
好像也不好解耦

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

回复 楼主 sswwss 的帖子

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

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

回复 地板 sswwss 的帖子

我是按照下面的高手做的但是在解决问题的时候会出现这样的警告：
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

回复 7楼 whatisi123 的帖子

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

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

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

我只能得到数值解

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

回复 9楼 yufeng 的帖子

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

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

谢谢各位！！
我用ode45得到了数值解
页: [1]
查看完整版本: 请教高手此微分方程组如何求解