+
−% 
+
−% 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