clear
clc
//-----------------------------Calcul des points------------------------

N=64;

f1=50;
T1=1/f1;
Te1=T1/N;
t1=0:Te1:(N-1)*Te1;
s1=sin(2*%pi*f1*t1)+(1/3)*sin(3*2*%pi*f1*t1)+(1/5)*sin(5*2*%pi*f1*t1);

f2=100;
T2=1/f2;
Te2=T2/N;
t2=0:Te2:(N-1)*Te2;
s2=cos(2*%pi*f2*t2)+(1/2)*cos(2*2*%pi*f2*t2)+(1/4)*cos(4*2*%pi*f2*t2);

//----------------------------Calcul de la fft--------------------------

fft_s1=abs(fft(s1));
fft_s2=abs(fft(s2));
