v2_009a.m

% v2_009a.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 rot=input('Rotation View Angle: (60)'); if isempty(rot),rot=60;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)))); view=[rot,30]; mesh(real(H),view) % title('Transfer Function (SDOF): Real Part') if plt==1,print -f1 -deps v2_009a,end; pause mesh(imag(H),view) % title('Transfer Function (SDOF): Imaginary Part') if plt==1,print -f1 -deps v2_009b,end; pause mesh(abs(H),view) % title('Transfer Function (SDOF): Magnitude') if plt==1,print -f1 -deps v2_009c,end; pause mesh(angle(H),view) % title('Transfer Function (SDOF): Phase') if plt==1,print -f1 -deps v2_009d,end; pause mesh(log(abs(H)),view) % title('Transfer Function (SDOF): Log Magnitude') if plt==1,print -f1 -deps v2_009e,end;