function [res]=sg(x,N)
x=x.*hamming(length(x))';
y=abs(fft(x,N)).^2;
res=y(1:1:N/2)';
