% 
% GETFFTFREQ: get frequencies for fft
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 
% Luigi Ferraioli 04-02-2011
% 
% % $Id: getfftfreq.m,v 1.1 2011/02/04 11:11:29 luigi Exp $
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function f = getfftfreq(nfft,fs,type)
  switch lower(type)
    case 'plain'
      f = (0:nfft-1).*fs./nfft;
    case 'one'
      f = (0:floor(nfft/2)).*fs./nfft;
    case 'two'
      if rem(nfft,2) % odd number of data
        f = fftshift([0:floor(nfft/2) -(floor(nfft/2)):-1].*fs./nfft);
      else % even number of data
        f = fftshift([0:floor(nfft/2)-1 -floor(nfft/2):-1].*fs./nfft);
      end
  end
end