function [HF,F] = FourierTransform(y, dt)

HF = fftshift(fft(y));
F = (-0.5:1/(numel(HF)):0.5-1/(numel(HF)))/dt;

end