r=randn(16,1000)+1;
rf=fft(r);
m=mean(rf')
s=std(rf')
abs(m)