新人求助 用龙格库塔法解振动微分方程
一个二阶振动微分方程,由于带有大范围运动和振动,因此质量矩阵中的特征值(或者分量)在数值上差别很大,转化为状态方程的过程中需要对质量矩阵求逆,因此就导致一个问题:在转化为状态方程后用ode23或者ode45解状态方程matlab运算速度很慢。请问我该如何解决。举一个简单的列子M=【M1,×;×,M2】,其中M1和M2在数值上差别很大
初步计算一下质量矩阵中最大特征值为700多,最小为0.002左右
[ 本帖最后由 jluzpeng 于 2008-7-31 15:11 编辑 ]
回复 楼主 的帖子
尝试使用ode13s ode113等回复 2楼 的帖子
谢谢,我试试看,如果不行还要来请教回复 楼主 的帖子
质量矩阵求逆,本来就是一个非常慢的问题,所以在计算的时候避免出现求逆[ 本帖最后由 无水1324 于 2008-8-2 13:57 编辑 ]
回复 4楼 的帖子
为了列状态方程必须求逆啊,我那个要用S-function 已解决,原来是我写程序的时候把两项颠倒了 请教楼主,怎样用龙格库塔法求解以个四自由度的系统的振动微分方程,能帮忙的话请回复,急求……谢谢 回复 4 # 无水1324 的帖子请教一下,解振动微分方程,初值怎么定啊
页:
[1]