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

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

一个二阶振动微分方程，由于带有大范围运动和振动，因此质量矩阵中的特征值（或者分量）在数值上差别很大，转化为状态方程的过程中需要对质量矩阵求逆，因此就导致一个问题：在转化为状态方程后用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

回复 2楼 的帖子

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

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

回复 楼主 的帖子

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

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

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

回复 4楼 的帖子

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

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

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

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

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

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

回复 4 # 无水1324 的帖子

请教一下，解振动微分方程，初值怎么定啊
页: [1]
查看完整版本: 新人求助 用龙格库塔法解振动微分方程