马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。
您需要 登录 才可以下载或查看，没有账号？我要加入
x
看到很多同学问有阻尼的系统方程如何解，这是我在书上看到的，希望对大家有帮助！！
下面是matlab程序，输入M，D，K，即可求得：
function [v,w,zeta]=vbr_sf(m,d,k)
%vbr_sf vbr_sf(m,d,k)
% [v,w,zeta]=vbr4(m,d,k)
% function vbr_sf finds the mode shapes and natural frequencies of
% a linear second order matrix equation.
% [v,w]=vbr_sf(m,k) finds the mode shapes and natural frequencies
% for the undamped case.
if nargin==2
k=d;
[v,w]=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.')
[v,w]=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=[0*k eye(length(k));-m\k -m\d];
[v,w1]=eig(a);
w=abs(w1);
zeta=-real(w1)/w;
end
end
w=diag(w);zeta=diag(zeta);
[ 本帖最后由 gxc8208 于 2007-5-24 16:41 编辑 ] |
发表于 2007-5-24 16:34
提升卡
变色卡
显身卡
发表于 2013-11-25 23:46
