amloved 发表于 2008-10-14 09:23

求助多自由度线性阻尼系统的频响函数

比如一个四自由度线性阻尼系统，m1=0.4,m2=0.8,m3=1.2,m4=1.6,c1=c2=c3=c4=12,k1=k2=k3=k4=36000,如何用matlab求出它的频响函数理论值，画出它的频响函数幅值曲线。有类似的程序吗，谢谢了！！！

yejet 发表于 2008-10-14 20:14

要用matlab实现的话
先用ode程序求解得到响应
然后fft

amloved 发表于 2008-10-15 08:53

好的，我试试。

咕噜噜 发表于 2008-10-15 09:47

线性系统啊，线性系统的话是比较简单的，可以正则化为单自由度系统之后按照单自由度的方法来画频域曲线，这个matlab里面有现成命令

ChaChing 发表于 2008-10-18 01:20

回复 楼主 amloved 的帖子

总感觉楼主给的条件不足
给了m1..m4/c1../c4/k1..k4, 但未给连接状况是列不出唯一运动方程的

amloved 发表于 2008-10-19 11:17

MDOF约束系统力学模型

vib 发表于 2008-10-19 13:21

你的结构图为什莫都没有给出呢？你的c，k都在哪些质量块之间？没有这些信息怎末确定你这个系统呢？怎末确定频响函数呢？

amloved 发表于 2008-10-19 15:20

=

=

=

这是它的刚度矩阵，质量矩阵，阻尼矩阵，求它的原点导纳H11的理论值，并画出它的伯德图

[ 本帖最后由 amloved 于 2008-10-19 15:28 编辑 ]

amloved 发表于 2008-10-20 12:18

大家帮帮忙吧，谢谢了。

ChaChing 发表于 2008-10-20 14:28

已经有m, c, k 矩阵, 就可以得出状态方程(state space), 再以bode求解
Help ss; Help bode ...很多, 查查找你需要的

咕噜噜 发表于 2008-10-20 14:33

是的，有了质量、刚度、阻尼矩阵就可以计算频响函数了，用bode命令可以
对于一些知识点还是自己掌握的好，找一下书看看，编出程序来这样以后就不会出现问题了
假如编程序过程中出现错误可以上来讨论

[ 本帖最后由 咕噜噜 于 2008-10-20 14:35 编辑 ]

amloved 发表于 2008-10-21 22:38

m1=;
c=;
k=;
num=1;
den=;
sys=tf(num,den);
bode(sys(1,1));
grid on
不对，分子分母应该指定为非空的列向量，什么意思？

ChaChing 发表于 2008-10-22 08:56

运动方程Mx"+Cx'+Kx=f, 设q1=x, q2=x' => '=A*+B*f
其中A即为特徵矩阵 A=; B=inv(M)
y=C*x+D*f;sys = ss(a,b,c,d); Bode(sys)
细节我没试, 自己试试看吧

amloved 发表于 2008-10-23 13:23

现在很晕
M=;
C=;
K=;
A=;
B=;
C=;
D=0;
sys=ss(A,B,C,D);
bode(sys)
应该怎么改？:@L

[ 本帖最后由 amloved 于 2008-10-23 14:16 编辑 ]

咕噜噜 发表于 2008-10-23 14:43

bode命令不是这样简单用的，给你推荐一本书看看吧
自己灵活运用一下
MATLAB原理与工程应用一个美国人写的
高会生等人翻译   电子工业出版社
页: [1] 2
查看完整版本: 求助多自由度线性阻尼系统的频响函数