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

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

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

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

x

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=[0,50];
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 发表

 登录/注册后可看大图

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

我是菜鸟，symbolic也没怎么用过
如果用atan(imag(H11)/real(H11))来替换angle（H11）
没有报错，但得出的图像明显不对，会是什么问题呢，相频曲线到底怎么画呢？

账号		自动登录	找回密码
密码			我要加入

[编程技巧] 画加速度相频曲线时遇到的问题

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