niluohe 发表于 2008-12-29 00:02

画加速度相频曲线时遇到的问题

Kr、Mr、Cr分别为已知的模态刚度、质量、阻尼矩阵
v为已知的固有振型
H为加速度频响函数
w为频率
要画H11(即频响函数矩阵的第一个元素)的相频曲线,下面的语句有什么问题么?

syms w
H=-w^2*v*inv(Kr-w^2*Mr+j*w*Cr)*v';
H11=H(1,1);
lim=;
ezplot(180*angle(H11)/pi,lim)

报错——??? Attempt to execute SCRIPT angle as a function.

ChaChing 发表于 2008-12-29 09:51

symbolic我不熟
angle函数支援symbolic吗?

niluohe 发表于 2008-12-29 11:06

原帖由 ChaChing 于 2008-12-29 09:51 发表 http://www.chinavib.com/forum/images/common/back.gif
symbolic我不熟
angle函数支援symbolic吗?
我是菜鸟,symbolic也没怎么用过
如果用atan(imag(H11)/real(H11))来替换angle(H11)
没有报错,但得出的图像明显不对,会是什么问题呢,相频曲线到底怎么画呢?
页: [1]
查看完整版本: 画加速度相频曲线时遇到的问题