Document Actions
sdof3.m
% SDOF3.M This is a script file to solve a sdof system
% given the mass, damping and stiffness terms
% in dimensionless units when the mass is varied.
% The output includes
% poles, residues (modal coefficients) and
% frequency domain plots of the
% frequency response functions.
damp=10
stiff=1600
f=linspace(0,50,500);
hold off
clg
for i=1:10
mass=2 + (i-1).*1
a=[mass,damp,stiff];
b=[0,0,1];
[r,p,k]=residue(b,a);
residue = r(1);
lambda=p(1);
xf=residue./(j.*f-lambda) + residue'./(j.*f-lambda');
semilogy(f,abs(xf))
hold on
end
title('Effect of Changing Mass')
xlabel('Frequency (Hz)')
ylabel('Magnitude')
grid
% meta plot1
hold off