u_find=find(f1>=3e-3 & f1<=0.3);
f1=f1(u);
aver_FFT=aver_FFT(u);
diff=(psd-aver_FFT)./psd;
red_f=f1*10/u_moy;
