function PlotFFT(x, Fs);
%% PlotFFT
%% ==============
%% Function to plot (scaled) fft magnitude of a signal.
%%
%% Code fragment adapted from:-
%% http://www.mathworks.com/support/tech-notes/1700/1702.html
%%
%% x is input signal
%% Fs is sampling frequency
%%
%% Created  : 23 Aug 2005
%% Modified : 23 Aug 2005
%%
%% 2005 Salman Durrani .
%% ------------------------------------------------------------

Fn=Fs/2;                 
NFFT=2.^(ceil(log(length(x))/log(2)));
FFTX=fft(x,NFFT);
NumUniquePts = ceil((NFFT+1)/2);
FFTX=FFTX(1:NumUniquePts);
MX=abs(FFTX);            
MX=MX*2;
MX(1)=MX(1)/2;   
MX(length(MX))=MX(length(MX))/2;  
MX=MX/length(x);
f=(0:NumUniquePts-1)*2*Fn/NFFT;
plot(f,MX);
xlabel('Frequency (Hz)')
ylabel('Magnitude')