%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% RESPONSE SPECTRUM 
%
% created by Michalis Fragiadakis, Dec 2013
% mfrag@mail.ntua.gr
%
% please report any bugs
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
includeFile


% assume 5% damping
acc.rec=acc.rec*9.81;
[td,amax,vmax,umax] = spectrum(acc.rec,acc.dtacc,0.05);


% acceleration spectrum
figure(1)
hold on; grid on; box on;
plot(td,amax/9.81,'linewidth',3)
title('acceleration spectrum','FontSize',18)
xlabel('period (sec)','FontSize',18)
ylabel('Spectral acceleration (g)','FontSize',20)

% velocity spectrum
figure(2)
hold on; grid on; box on;
plot(td,vmax,'linewidth',3)
title('velocity spectrum','FontSize',18)
xlabel('period (sec)','FontSize',18)
ylabel('Spectral velocity (m/s)','FontSize',20)

% displacement spectrum
figure(3)
hold on; grid on; box on;
plot(td,umax,'linewidth',3)
title('displacement spectrum','FontSize',18)
xlabel('period (sec)','FontSize',18)
ylabel('Spectral displacement (m)','FontSize',20)


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% GET SPECTRAL VALUES FOR A SINGLE PERIOD VALUE (T=1sec)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

T = 1.5; % sec  <-- period

[umax,vmax,amax]=sdfL(T,acc.rec,0.05,acc.dtacc,acc.nsteps);


% plot
figure(1)
plot(T,amax,'ro','LineWidth',3,'MarkerSize',10)

figure(2)
plot(T,vmax,'ro','LineWidth',3,'MarkerSize',10)

figure(3)
plot(T,umax,'ro','LineWidth',3,'MarkerSize',10)