声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1881|回复: 7

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

[复制链接]
发表于 2008-7-31 10:12 | 显示全部楼层 |阅读模式

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

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

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

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

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2008-7-31 14:36 | 显示全部楼层

回复 楼主 的帖子

尝试使用ode13s ode113等

评分

1

查看全部评分

 楼主| 发表于 2008-7-31 15:07 | 显示全部楼层

回复 2楼 的帖子

谢谢,我试试看,如果不行还要来请教
发表于 2008-8-1 20:13 | 显示全部楼层

回复 楼主 的帖子

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

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

评分

1

查看全部评分

 楼主| 发表于 2008-8-3 13:00 | 显示全部楼层

回复 4楼 的帖子

为了列状态方程必须求逆啊,我那个要用S-function
 楼主| 发表于 2008-8-10 09:04 | 显示全部楼层
已解决,原来是我写程序的时候把两项颠倒了
发表于 2010-12-20 22:03 | 显示全部楼层
请教楼主,怎样用龙格库塔法求解以个四自由度的系统的振动微分方程,能帮忙的话请回复,急求……谢谢

点评

动手搜一下  发表于 2010-12-21 00:08
发表于 2011-9-18 01:21 | 显示全部楼层
回复 4 # 无水1324 的帖子

请教一下,解振动微分方程,初值怎么定啊
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-13 15:08 , Processed in 0.078571 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表