function FFTejemplo

fs = 100;                          % Frecuencia de muestreo(sample)(Hz)
t = 0:1/fs:10-1/fs;                % Muestra de 10 segundos
x = (1.3)*sin(2*pi*15*t) ...       % 15 Hz 
  + (1.7)*sin(2*pi*40*(t-2)) ...   % 40 Hz 
  + (2.5)*randn(size(t));          % Ruido blanco;

%Calculo de la DFT y espectro de potencia:

m = length(x);          % largo ventana
n = pow2(nextpow2(m));  % largo de la transformada
y = fft(x,n);           % DFT
f = (0:n-1)*(fs/n);     % Rango de Frecuencia
power = y.*conj(y)/n;   % Espectro de frecuencia

%Graficar el espectro
plot(f,power)
xlabel('Frecuencia (Hz)')
ylabel('Potencia')
title('Periodograma')

%Graficar el espectro centrado en 0

y0 = fftshift(y);          % Rearrange y values
f0 = (-n/2:n/2-1)*(fs/n);  % 0-centered frequency range
power0 = y0.*conj(y0)/n;   % 0-centered power

plot(f0,power0)
xlabel('Frecuencia (Hz)')
ylabel('Potencia')
title('Periodograma Centrado en 0')

end
