clc;
clear;
close all;

data1 = importdata ('amort14cort.txt');
tempo1 = data1(:,1);
x1 = data1(:,2);
y1 = data1(:,3);      %aceleração

figure 1
plot (tempo1,y1, 'b');
%xlim ([6.91 7])
xlabel ('Tempo [s]')
ylabel ('Amplitude [g]')
title ('Resposta de vibração ao Impulso')


data2= importdata ('amort222cort.txt');
tempo2= data2 (:,1)
x2= data2(:,2);
y2= data2(:,3);  %aceleração

figure 2
plot(tempo2,y2, 'r');
%xlim([18.4 18.43])

data3= importdata ('amort32cort.txt');
tempo3= data3 (:,1);
x3= data3(:,2); %amplitude
y3= data3(:,3);

figure 3
plot(tempo3,y3, 'g');
%xlim([7.41 7.5])

data4= importdata ('amort431cort.txt');
tempo4= data4 (:,1);
x4= data4(:,2); %amplitude
y4= data4(:,3);

figure 4
plot (tempo4,y4, 'm');
%xlim([11.24 11.36])

figure 5
plot(tempo1-6.9146,y1, 'b',tempo2-18.405,y2, 'r',tempo3-7.4134,y3, 'g',tempo4-11.248,y4, 'm')
xlim([0 0.1])
xlabel('Tempo[s]')
ylabel('Amplitude [g]')
legend ('Resposta ao impulso 1° Amostra', 'Resposta ao Impulso 2° Amostra','Resposta ao Impulso 3° Amostra','Resposta ao Impulso 4° Amostra')
title ('Resposta ao Impulso metodo experimental das amostras')