function [S]=ejemplopotencia(w,N);
% entradas
% w ...... pulsacion del tono
% N ...... numero de muestras del proceso x(n)
% salida
% S ...... densidad espectral de potencia
%
M=1024;
if(N>M)
    M=N;
end
x=zeros(N,1);
x=cos(w*[0:N-1]'+2*pi*(rand(1,1)-0.5))+0.001*randn(N,1);
S=abs(fft(x,M)).^2/N;
plot(linspace(0,0.5,M/2),S(1:M/2));
ylabel('Densidad Spectral Potencia')
xlabel('frecuencia');
title(['potencia tono: ' num2str(2*max(S)/N)]);
grid on