function[S,f]=fouriert_b(s,Dt)
Df=(1/((length(s)-1)*Dt));
f=(-0.5/Dt):Df:(0.5/Dt);
S=fftshift(fft(fftshift(s)));




