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]
查看完整版本: 新人求助 用龙格库塔法解振动微分方程