// Q 2.2 

// pour effacer les fenetres de commande et les espaces mémoires
clear all;
//***********************************************************

// génération d'un signal sinusoidal de fréquence :
fs=6000;  // Hz
fe=16000;
Te=1/fe;
N=32;
n=0:N-1;
t=n*Te;   // en s

// domaine temporel
xs=sin(2*%pi*fs*t);

// domaine fréquentiel
Xs=fft(xs.*window('tr',N),-1);
Xs=fftshift(Xs);


xset("window",0);xbasc(); // selectionne (ou crée) la fenetre 0 et l'efface
titlepage("signal sinusoidal, frequence 6 kHz, Fe = 16 kHz, FFT sur 32 points")

// visualisation des 20 premiers points
subplot(211);plot2d(t(1:20),xs(1:20));
xtitle("Representation temporelle","t (s)")

// Spectre 
f=(-.5:1/N:.5-1/N)*(fe*0.001);   // en kHz (dénormalisation de la fréquence normalisée)
subplot(212);plot2d(f,abs(Xs/length(Xs)),1) // NB : normalisation / NFFT
xtitle("Representation frequentielle","f (kHz)")