clear all;
close all;

%
%% Paramêtres utilisatrice
%

T = 0.25;               % Fenêtre d'observation [s]
M=32;                   % Nombre d'échantillons dans la fenêtre
Fsin = 4 ;              % Fréquence d'un sinus [Hz]

%
%% Paramêtres déduits

Te = T/M;               % Période d'échantillonnage [s]
Tsim = T-Te;            % Instant final de simulation avec M points [s] 
fprintf("Fréquence = %f Hz", 1/Te)

%
%% Calculs préliminaires
%% Lancement de simulink

tic
out=sim("simDFT"); 
toc

%
%% Visualisations temporelle

plot(out.sin_continu,'.-');
hold on;

t =   out.sin_echanti.Time;
sig = out.sin_echanti.Data;

stem(t, sig , 'o');
grid on;

%
%% Calcul DFT

Tfd = fft(sig);
abscisses=1:length(Tfd); 

%% visualisation DFT

figure()
plot(abscisses,abs(Tfd),'o-');
ylabel("Module Transformée de Fourier")
xlabel("???  indices de tableau Matlab de 1 à 32 ???")