关于MATLAB中求固有频率以及振型的问题
本帖最后由 牛小贱 于 2014-12-17 09:56 编辑在论坛里看到前辈给出的求固有频率以及振型的程序,但是运行不成功,不知道为什么。我才开始学习matlab,所以请各位前辈指点一下,多谢。程序如下:
输入M,D,K
function =vbr_sf(m,d,k)
%vbr_sfvbr_sf(m,d,k)
% =vbr4(m,d,k)
% function vbr_sf finds the mode shapes and natural frequencies of
% a linear second order matrix equation.
% =vbr_sf(m,k) finds the mode shapes and natural frequencies
% for the undamped case.
if nargin==2
k=d;
=eig(m\k);
w=sqrt(w);
end
if nargin==3
if norm(d/m*k-k/m*d) < 1e-8*norm(k/m*d)
%disp('Damping is proportional, eigenvectors are real.')
=eig(m\k);
w=sqrt(w);
zeta=(v'*m*v)\(v'*d*v)/2/w;
else
%disp('Damping is non-proportional, eigenvectors are complex.')
a=;
=eig(a);
w=abs(w1);
zeta=-real(w1)/w;
end
end
w=diag(w);zeta=diag(zeta);
求助完整格式:出错代码和出错提示 ChaChing 发表于 2014-11-30 20:17
求助完整格式:出错代码和出错提示
多谢,我已经解决了 ChaChing 发表于 2014-11-30 20:17
求助完整格式:出错代码和出错提示
前辈您好,我还有一个问题想请教您。我在论坛里问了一个问题,您能不能帮我看看如何解决呢http://forum.vibunion.com/thread-134903-1-1.html 十分感谢!
页:
[1]