function [X,freq]=get_fft(x,Fs)
N=(length(x)); k=0:N-1; T=N/Fs; freq=k/T; cutOff = ceil(N/2);
freq = freq(1:cutOff); X=fft(x)/N*2; X = X(1:cutOff); 
end