声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1483|回复: 5

[绘图技巧] 时变微分方程组的幅频曲线如何画?

[复制链接]
发表于 2007-9-6 23:55 | 显示全部楼层 |阅读模式

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

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

x
我浏览了许多关于时频曲线画法的帖子,但还是不得要领。特别是我还不会使用maple。
哪位高手能给我一个matlab算法的具体思路(最好有个例子什么的)。
问题如下:AX‘’+B(t)X‘+C(t)=0,其中X为4*1矩阵,A为4*4矩阵,B(t)和C(t)都是含有sin(t)的矩阵,
我现在想画出X中两个变量的时频曲线,该怎么办?
初次拜访,请多指教!

[ 本帖最后由 花如月 于 2007-9-7 12:56 编辑 ]
回复
分享到:

使用道具 举报

发表于 2007-9-7 08:26 | 显示全部楼层

回复 #1 IceHill 的帖子

X中的两个变量是什么意思?
是什么变量?
 楼主| 发表于 2007-9-7 10:30 | 显示全部楼层
不好意思,没说清楚。
就是说X是4*1列向量,2个代表位移,两个代表受力
我想知道头两个变量的幅频曲线。
我查了一些资料,用谐波平衡法求解,可是我所知道的谐波法
都是一些1元或者2元的,都是先手算,然后在代入matlab
现在是四元的,而且是变系数的,手工明显不可以了,所以
我想知道matlab里面有没有一些库函数可以直接算。
谢谢大家。
发表于 2007-9-7 12:32 | 显示全部楼层

同样的疑问!

我也在这个问题上打转呢,我是12自由度的二阶微分方程组,刚度矩阵是时变的,方程组矩阵形式为:
[M]X''+CX'+K(t)X=F(t),现在只是刚刚解出时域数值解,我还想得出频域解,做出各变量随频率的变化情况曲线图。请做过这方面研究的朋友,指点一下,非常着急,非常感谢!!!
发表于 2007-9-7 12:55 | 显示全部楼层
请楼上的和楼主把问题描述清楚,新来的先阅读一下置顶帖子了解一下版规。这样可以节省你们的时间,也可以节省我们的时间。不要让我们带着疑问去探讨你们的问题。

最后说明一下,本版区关于函数绘图的例子和资料非常丰富。应该可以解决你们的问题,从你们的描述来看解都有了。就是一个数据绘图的问题,能用简单的例子描述清楚就好了
发表于 2007-9-7 13:08 | 显示全部楼层

数值解出来了,不过是时域解!

方程形如:[m]{X''}+[C]{X'}+[K(t)]{X}={P}
其中[m]=diag{m1,m2,m3,m4,m5,m6,m7,m8,m9,m10,m11,m12}
    {X}=(x1,x2,x3,x4,x5,x6,x7,x8,x9,x0,x11,x12)
    [C]为12×12阶的常数矩阵
    [K]为12×12阶的时变矩阵,其中每个元素均是时间t的函数,Kij=Kij(t)
    {P}={p1(t),p2(t),p3(t),p4(t),p5(t),p6(t),p7(t),p8(t),p9(t),p10(t),p11(t),p12(t)}’
用数值解法(四阶龙格库塔法),借助matlab进行求解,解出了每个自变量不同时刻的解,但是怎么得到频域解,输出频率特性曲线呢?非常感谢,请明示!

[ 本帖最后由 eight 于 2007-9-7 13:51 编辑 ]
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-14 04:18 , Processed in 0.063491 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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