// Q 3.1

fe=8000;
Te=1/fe;
f1=1000; f2=1210; 
N=128;
n=0:N-1;
t=n*Te;   // en s

// domaine temporel
x=sin(2*%pi*f1*t)+sin(2*%pi*f2*t);
// domaine fréquentiel
X=fft(x,-1);
X=fftshift(X);

// dénormalisation de la fréquence normalisées à + ou - 0,5
f=(-.5:1/N:.5-1/N)*(fe*0.001);   // en kHz
xset("window",0);xbasc(); // selectionne (ou crée) la fenetre 0 et l'efface
subplot(212);plot2d(f,20*log10(abs(X/length(X))),1,"011","",[-4 -40 4 0]) // NB : normalisation / NFFT
xtitle("Representation frequentielle (en dB)","f (kHz)")