vtb5_1.m

来自「振动工具箱」· M 代码 · 共 32 行

function vtb5_1(m,c,k)

%VTB5_1 Transmissibility ratio of a SDOF system.
%  VTB5_1(m,c,k) plots the tranmissibility ratio for the system
%  described by the mass m, damping coefficient c, and stiffness k.
%  VTB5_1(zeta,w) plots the transmissibility ratio for the system
%  described by the damping ratio zeta and undamped natural
%  frequency w.
% VTB5_1(23,0.9,34);

if nargin==2
   z=m;
   w=c;
else
   w=sqrt(k/m);
   z=c/(2*m*w);
end
r=[0:.01:2]';
%Calculate various parts of equation 5.7 for TR
TR1=ones(201,1)+(2*z.*r).^2;
TR2=(ones(201,1)-r.^2).^2+(2*z.*r).^2;
TR=sqrt(TR1./TR2);

aa=version;ll=length(aa);

plot(r,TR)
title(['Transmissibility plot for \zeta = ',num2str(z),' \omega = ',num2str(w),' rad/s'])
grid on
ylabel('Transmissibility Ratio')
xlabel('Dimensionless Frequency')