%f=k*fe/N
%fmin = f(k=0) = 0
%fmax = f(k=N-1) = (N-1)/N * fe ~= fe
%k=0 correspond à la fréquence f=0
%k=N-1 """""""""" max f= (N-1)/N * fe

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

%la fréquence est 20Hz
f=20
%Le théorème de shanon est respecté, car 128 est plus grand que 2x20

k=0:N+nz-1
%z = zeros(1,nz);
%x=[x z]
X=fft(x,448+64);
%X=fftshift(X);
magX=abs(X);
f=k/(N+nz)/Te %pour remplacer k par des Hz
plot(f,magX),title('module de la TF'), xlabel('k'), grid

