N=64;
nz=448;
Te=1/128;
k=0:N-1;
x=sin(2*pi*20*k*Te);
plot(k,x);

X=fft(x);
X=fftshift(X)
magX=abs(X);
plot(k, magX), grid;