Single Side Band FFT for Matlab

fftsinglesideband.m

function [resultado n]=fftsinglesideband(senal);

   %% Funcion que retorna una FFT de Single Side
   % ademas retorna n, la cantidad de puntos 
   % considerados para la FFT
   
   % set de frecuencias se puede generar 
   % usando
   % frecuencias       =SAMPLING*(0:(n/2))/n;
   
   L=length(senal);
   n=2^nextpow2(L);
   fftmoduladaFM1            =abs( fft(senal, n)/n );
   fftsinglesideFM1          =fftmoduladaFM1(1:n/2+1);
   fftsinglesideFM1(2:end-1) =2*fftsinglesideFM1(2:end-1);
   resultado=fftsinglesideFM1;
end