close all;
Flasers = [85.005E3,90E3,95E3,100E3,115E3,120E3];
Fsin = Flasers;
deltaF = 5E3;
T = 1/deltaF;
Fe = 320E3;%kHz
Te = 1/Fe;
M= Fe*T;
Tsim = T-Te;
frequence = deltaF*(0:1:M-1);


Res=sim('SimuliDFT');



figure
plot(Res.temps,Res.Sinus_Continu);
hold on;
plot(Res.temps_Echanti,Res.Sinus_Echanti, 'o');
sig= Res.Sinus_Echanti;
%sig(1:20,4)=0;
echantFft = fft(sig,M);


figure
stem(frequence,abs(echantFft)/M);


%carreeee-----------------------

% figure(1)
% plot(Res.temps,Res.Carre_Continu1);
% hold on;
% plot(Res.temps,Res.Sinus_Carre_Continu);
% hold on;
% plot(Res.temps_Echanti,Res.Sinus_Carre_Echanti, 'o');
% echantFft_Carre = fft(Res.Sinus_Carre_Echanti,M);
% 
% 
% 
% figure(2)
% stem(frequence,abs(echantFft_Carre)/M);
% 
% Fp = tf(1,[1.7483*10^(-23), 7.6663*10^(-18), 1.162*10^(-11), 3.0332*10^(-6),1]);
% figure(3)
% bode(Fp);