matlab 冲击响应,冲击响应谱计算的matlab程序.doc

disp(' ')

disp(' srs.m ver 2.0 July 3, 2006')

disp(' by Tom Irvine Email: tomirvine@')

disp(' ')

disp(' This program calculates the shock response spectrum')

disp(' of an acceleration time history, which is pre-loaded into Matlab.')

disp(' The time history must have two columns: time(sec) & acceleration')

disp(' ')

%

clear t;

clear y;

clear yy;

clear n;

clear fn;

clear a1;

clear a2

clear b1;

clear b2;

clear jnum;

clear THM;

clear resp;

clear x_pos;

clear x_neg;

%

iunit=input(' Enter acceleration unit: 1= G 2= m/sec^2 ');

%

disp(' ')

disp(' Select file input method ');

disp(' 1=external ASCII file ');

disp(' 2=file preloaded into Matlab ');

file_choice = input('');

%

if(file_choice==1)

[fil