% clear;
close all;

% load enregistrement;

% Calcul de la transformée de Fourier à court terme :
n_fenetre = 2^10; % Largeur de la fenêtre (en nombre d'échantillons)
n_decalage = n_fenetre; % Décalage entre positions successives de la fenêtre (en nombre d'échantillons)
fenetre = 'rect'; % Type de la fenêtre : 'rect' ou 'hann'

[Y, valeurs_tau, valeurs_f] = TFCT(y, f_ech, n_fenetre, n_decalage, fenetre);
S = 20 * log10(abs(Y) + eps);

% Affichage du module de la transformée de Fourier à court terme :
figure('Name', 'Transformée de  Fourier');
imagesc(valeurs_tau, valeurs_f, S);
axis xy;
xlabel('Temps (s)');
ylabel('Frequence (Hz)');
drawnow;

[Y, valeurs_tau, valeurs_f] = TFCTST(y, f_ech, n_fenetre, n_decalage, fenetre);
S = 20 * log10(abs(Y) + eps);

% Affichage du module de la transformée de Fourier à court terme :
figure('Name', 'Transformée de  Fourier steg');
imagesc(valeurs_tau, valeurs_f, S);
axis xy;
xlabel('Temps (s)');
ylabel('Frequence (Hz)');
drawnow;

[signal_restitue, ~] = ITFCT(Y, f_ech, n_decalage, fenetre);
audiowrite("test.wav", signal_restitue, f_ech);

% save exercice_1;