function PlotFT(sin, fs)
% PlotFT - plot the Fourier transform (FFT) of a signal

ff = (0:length(sin)-1)*fs/length(sin) - fs/2;

Sn = fft(sin);
semilogy(ff,abs(fftshift(Sn)));
xlabel('Frequency (Hz)')
ylabel('|Fourier Transform|')
