沖擊響應(yīng)譜計(jì)算的matlab程序

1、disp(' ') disp(' srs.m ver 2.0 July 3, 2006') disp(' by Tom Irvine Email: tomirvine@aol.com') disp(' ') disp(' This program calculates the shock response spectrum')disp('

2、 of an acceleration time history, which is pre-loaded into Matlab.')disp(' The time history must have two columns: time(sec) clear y;clear yy;clear n; clear fn;clear a1;clear a2clear b1;clear b2

3、; 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 ASC

4、II file ');disp(' 2=file preloaded into Matlab '); file_choice = input(''); %if(file_choice==1)[filename, pathname] = uigetfile('*.*');filename = fullfile(pathname, filename); % fid = f

5、open(filename,'r');THM = fscanf(fid,'%g %g',[2 inf]);THM=THM';elseTHM = input(' Enter the matrix name: '); end% t=double(THM(:,1));disp(' Calculating response..... ') % % SRS engine

6、% for j=1:1000%omega=2.*pi*fn(j);omegad=omega*sqrt(1.-(damp^2));cosd=cos(omegad*dt);sind=sin(omegad*dt);domegadt=damp*omega*dt;%if(ialgorithm==1)a1(j)=2.*exp(-domegadt)*cosd;a2(j)=-exp(-2.*domegadt);b1(j)=2.*domegadt;b2(

7、j)=omega*dt*exp(-domegadt);b2(j)=b2(j)*( (omega/omegad)*(1.-2.*(damp^2))*sind -2.*damp*cosd );b3(j)=0;%elseE=exp(-damp*omega*dt);K=omegad*dt;C=E*cos(K);S=E*sin(K);Sp=S/K;%a1(j)=2*C;a2(j)=-E^2; b1(j)=1.-Sp; b2(j)=2.*(Sp-C