v2_020.m

美国辛辛那提大学的振动教材

% v2_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:  [phone]
% FAX:  [phone]
% E-MAIL: [e-mail]
%*********************************************************************
%
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_020a,end;
pause
mesh(imag(H),view)
%     title('Transfer Function (SDOF):  Imaginary Part')
if plt==1,print -f1 -deps v2_020b,end;
pause
mesh(abs(H),view)
%     title('Transfer Function (SDOF):  Magnitude')
if plt==1,print -f1 -deps v2_020c,end;
pause
mesh(angle(H),view)
%     title('Transfer Function (SDOF):  Phase')
if plt==1,print -f1 -deps v2_020d,end;
pause
mesh(log(abs(H)),view)
%     title('Transfer Function (SDOF):  Log Magnitude')
if plt==1,print -f1 -deps v2_020e,end;