close all
clear all
clc
 
load frf_3497_1;
dims = size(frf_3497);
clear frf_3497;

total_resp = zeros(dims(1),dims(2));

tic
for k=1:18

    eval(sprintf('load frf_3497_%d;',k)); %

    total_resp = total_resp + frf_3497;
k 
end
toc


%IFFT section 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%
  
freq = 0:0.5:2000;

response = zeros(dims(1), 2*dims(2));

c = length(freq);

for k=1:dims(1)
[response(k,:), time] = IFFT_time_resp(total_resp(k,1:c),freq(1:c));
k
end

save hoop_strain_response response;
save time time;