求助--怎么用龙格库塔法解多元二阶微分方程组

yjh2911 · 发表于 2006-3-24 08:12

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

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

x

求助--怎么用龙格库塔法解多元二阶微分方程组 有范例更好 谢谢！！

happy · 发表于 2006-3-24 08:19

doc ode45 例子里边也有

siyanger · 发表于 2006-3-30 07:48

我也想知道，例了里的太简单了。

happy · 发表于 2006-3-30 19:16

简单来说就是将方程转化到状态空间 然后根据状态空间中的方程写出function 该function返回值是广义坐标的一阶导数，输入值就是广义坐标，当然还可以有参变量 然后直接调用ode45，ode45的具体参数表示什么含义帮助里说得很清楚了 至于具体例子，论坛中多得很，自己找来看吧，要多复杂都有

siyanger · 发表于 2006-3-30 20:07

<DIV class=quote>以下是引用happy在2006-3-30 19:16:46的发言： 简单来说就是将方程转化到状态空间 然后根据状态空间中的方程写出function 该function返回值是广义坐标的一阶导数，输入值就是广义坐标，当然还可以有参变量 然后直接调用ode45，ode45的具体参数表示什么含义帮助里说得很清楚了 至于具体例子，论坛中多得很，自己找来看吧，要多复杂都有</DIV>
要多复杂都有，不是吧？我另发了一个贴子了。这个微分方程，说来不怕你笑我笨，我解了快一个月了，还没解出来。先是将方程近似线性化，再就无量纲化，推来推去，都不行。结果还是发散哪（Y无限大），但在实际中，这种方程是有解的，而且是有规律的。呈正或余弦变化。

happy · 发表于 2006-3-31 19:45

<DIV class=quote>以下是引用siyanger在2006-3-30 20:07:03的发言： 
要多复杂都有，不是吧？我另发了一个贴子了。这个微分方程，说来不怕你笑我笨，我解了快一个月了，还没解出来。先是将方程近似线性化，再就无量纲化，推来推去，都不行。结果还是发散哪（Y无限大），但在实际中，这种方程是有解的，而且是有规律的。呈正或余弦变化。</DIV>
 呵呵，理论上用rk法求解是没有问题的，但是对于高维系统，rk法效果很差的，发散是常见，换算法吧

siyanger · 发表于 2006-3-31 19:56

能推荐一下什么算法比较好吗？非常感谢！！ 又：我这个方程的自由度说来也不算很多，15个。变成状态方程就是30个了。

happy · 发表于 2006-3-31 19:58

<DIV class=quote>以下是引用siyanger在2006-3-31 19:56:26的发言： 能推荐一下什么算法比较好吗？非常感谢！！ 又：我这个方程的自由度说来也不算很多，15个。变成状态方程就是30个了。</DIV>
 30个求解就非常困难了 曾经用rk法算过16个自由度的，调得只郁闷 虽然最后算出来了，代价太大了 可以试试直接积分法，就是newmark之类的

siyanger · 发表于 2006-3-31 21:11

呵呵，我的方程是15个自由度啊，比起你的16个自由度，还少一个呢。我是说变为状态方程是30个。你的不是8自由度吧？我再试试直接积分法吧，再作不出来，毕业不了啦！谢谢你哟！ 等我解出来，再交流交流！

yjh2911 · 发表于 2006-4-1 16:15

给个实例--谢谢！！

happy · 发表于 2006-4-1 16:31

<DIV class=quote>以下是引用yjh2911在2006-4-1 16:15:09的发言： 给个实例--谢谢！！</DIV>
自己到前面的帖子里找吧

yjh2911 · 发表于 2006-4-12 20:11

我的是6个自由度 状态方程是12个 可是总是提示初始矩阵小，怎么回事？

happy · 发表于 2006-4-14 20:48

<DIV class=quote>以下是引用yjh2911在2006-4-12 20:11:21的发言： 
我的是6个自由度 状态方程是12个 可是总是提示初始矩阵小，怎么回事？</DIV> 把代码给出来，问题应该是你给的初值少于12个量

[此贴子已经被作者于2006-4-14 20:49:30编辑过]

常新宇 · 发表于 2014-5-19 23:32

我也想知道

tnanr · 发表于 2014-7-24 19:53

？？？doc ode45 例子里边也有这个在哪里下载？

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

[编程技巧] 求助--怎么用龙格库塔法解多元二阶微分方程组

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

本帖被以下淘专辑推荐:

回复：（yjh2911）求助--怎么用龙格库塔法解多元二阶微...

回复：（yjh2911）求助--怎么用龙格库塔法解多元二阶微...

回复：（siyanger）以下是引用happy在2006-3-30 19:16...

回复：（siyanger）能推荐一下什么算法比较好吗？非常...

给个实例--谢谢

回复：（yjh2911）给个实例--谢谢

回复：（yjh2911）我的是6个自由度状态方程是12个可是...