clc;
clear all;
close all;

T = 1/5000;
M = 64;
Te = T/M;
Tsim = T-Te;
Fe = 1/Te;

%%Fsin = 4;
%%Fsin = 60;
%%Fsin = 64; shannon < Fe/2
%%Fsin = 100;
Fsin=85000;
Res = sim('SimulDFT');

plot(Res.Sinus_Continu);
hold on; % permet de superposer la courbe à suivre
plot(Res.Sinus_Echanti,'o');

%figure
f = fft(Res.Sinus_Echanti.data, M)/M;
%stem(linspace(0, M-1, M), abs(f),'o');
figure
stem(linspace(0, M-1, M) / (M*Te), abs(f),'o');

%%figure
%%plot(linspace(0, M-1, M) / (M*Te/100), abs(fft(Res.Sinus_Continu.data, M)/M));
%%hold on;
%%plot(linspace(0, M-1, M) / (M*Te), abs(fft(Res.Sinus_Echanti.data, M)/M),'o');