clear
close all

load data1_Normalized
y1=data;
time1=linspace(0,length(y1)/fs,length(y1));
Nfft1 =1024*100; %Number of FFT cycles (Points produced)
df1 = fs/(Nfft1-1);
f1 = 0:df1:fs;
fft_Signal1=abs(fft(y1,Nfft1));
fft_Signal1=fft_Signal1/max(fft_Signal1);
load data2_Normalized
y2=data;
df2 = fs/(Nfft1-1);
f2 = 0:df2:fs;
fft_Signal2= abs(fft(y2,Nfft1));
fft_Signal2=fft_Signal2/max(fft_Signal2);
load data
y3=data;
df3 = fs/(Nfft1-1);
f3 = 0:df3:fs;
fft_Signal3= abs(fft(y3,Nfft1));
fft_Signal3=fft_Signal3/max(fft_Signal3);
figure(1)
plot (f1,fft_Signal1,'k'); hold on
plot (f2,fft_Signal2,'r'); 
plot (f3,fft_Signal3,'b'); hold off
ylabel ('fft nomalized amplitude');
xlabel ('Frequency, Hz');
xlim([0 500])