samer@32: function f=freqscale(N,m,fs)
samer@32: % freqscale - return frequency scale for FFT analysis
samer@32: %
samer@32: %    FREQSCALE(N,m,fs): frequency scale for the first m bins of an N-pt FT, 
samer@32: %			with a sampling rate=fs, defaults to 1.
samer@32: 
samer@32: if nargin<2, fs=1; end
samer@32: f=fs*(0:(m-1))/N;