v1_020.m

% v1_020.m % % This is a script file to solve a sdof system % given the mass, damping and stiffness terms % in dimensionless units. The output is a three % dimensional plot in the s domain (complex independent % variable. % % SDOF System, Laplace Domain (3D) plot % Figures for UC-SDRL-CN-20-263-662, Chapter 2 %********************************************************************** % Author: Randall J. Allemang % Date: 18-Apr-94 % Structural Dynamics Research Lab % University of Cincinnati % Cincinnati, Ohio 45221-0072 % TEL: 513-556-2725 % FAX: 513-556-3390 % E-MAIL: randy.allemang@uc.edu %********************************************************************* % clear,clg plt=input('Store plots to file (Yes=1): (0)');if isempty(plt),plt=0;end mass=10; a=input('Real Part of Pole: (-0.0625)'); if isempty(a),a=-0.0625;end b=input('Imaginary Part of Pole: (0.51)'); if isempty(b),b=0.51;end lambda(1)=a+j*b; lambda(2)=a-j*b; % % Set up mesh for only quadrants two and three % [sigma,omega]=meshdom(-0.2:0.005:0,-1:0.02:1); s=sigma+j*omega; H=(1.0./mass).*(1.0./((s-lambda(1)).*(s-lambda(2)))); fig1=figure(1); mesh(sigma,omega,real(H)) view([1 1.5 1]); title('Transfer Function (SDOF): Real Part') pause fig2=figure(2); mesh(sigma,omega,imag(H)) view([1 1.5 1]); title('Transfer Function (SDOF): Imaginary Part') pause fig3=figure(3); mesh(sigma,omega,abs(H)) view([1 1.5 1]); title('Transfer Function (SDOF): Magnitude') pause fig4=figure(4); mesh(sigma,omega,angle(H)) view([1 1.5 1]); title('Transfer Function (SDOF): Phase') pause fig5=figure(5); mesh(sigma,omega,log(abs(H))) view([1 1.5 1]); title('Transfer Function (SDOF): Log Magnitude') if plt==1 print -deps -f1 v1_020a.eps print -deps -f2 v1_020b.eps print -deps -f3 v1_020c.eps print -deps -f4 v1_020d.eps print -deps -f5 v1_020e.eps end pause