新人求助用龙格库塔法解振动微分方程

jluzpeng · 发表于 2008-7-31 10:12

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

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

x

一个二阶振动微分方程，由于带有大范围运动和振动，因此质量矩阵中的特征值（或者分量）在数值上差别很大，转化为状态方程的过程中需要对质量矩阵求逆，因此就导致一个问题：在转化为状态方程后用ode23或者ode45解状态方程matlab运算速度很慢。请问我该如何解决。举一个简单的列子
M＝【M1，×;×，M2】，其中M1和M2在数值上差别很大
初步计算一下质量矩阵中最大特征值为700多，最小为0.002左右

[ 本帖最后由 jluzpeng 于 2008-7-31 15:11 编辑 ]

alljoyland · 发表于 2008-7-31 14:36

尝试使用ode13s ode113等

jluzpeng · 发表于 2008-7-31 15:07

谢谢，我试试看，如果不行还要来请教

无水1324 · 发表于 2008-8-1 20:13

质量矩阵求逆,本来就是一个非常慢的问题，所以在计算的时候避免出现求逆

[ 本帖最后由无水1324 于 2008-8-2 13:57 编辑 ]

jluzpeng · 发表于 2008-8-3 13:00

为了列状态方程必须求逆啊，我那个要用S－function

jluzpeng · 发表于 2008-8-10 09:04

已解决，原来是我写程序的时候把两项颠倒了

excellentzzg · 发表于 2010-12-20 22:03

请教楼主，怎样用龙格库塔法求解以个四自由度的系统的振动微分方程，能帮忙的话请回复，急求……谢谢

lq271311 · 发表于 2011-9-18 01:21

回复 4 # 无水1324 的帖子

请教一下，解振动微分方程，初值怎么定啊

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

[编程技巧] 新人求助用龙格库塔法解振动微分方程

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

本帖被以下淘专辑推荐:

回复楼主的帖子

评分

回复 2楼的帖子

回复楼主的帖子

评分

回复 4楼的帖子

点评

[编程技巧] 新人求助 用龙格库塔法解振动微分方程

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

本帖被以下淘专辑推荐:

回复 楼主 的帖子

评分

回复 2楼 的帖子

回复 楼主 的帖子

评分

回复 4楼 的帖子

点评

[编程技巧] 新人求助用龙格库塔法解振动微分方程

回复楼主的帖子

回复 2楼的帖子

回复楼主的帖子

回复 4楼的帖子