clc;
clear all;
close all;

F1 = 85005.9;
F2= 90000;
F3 = 94986.8;
F4 = 100000;
F5 = 115015.9;
F6 = 120000;
Fe = 320000;
F = 5000;

Te = 1/Fe;
T = 1/F ;
M = T/Te;
Tsim = T-Te;
Res = sim('LG');

plot(Res.Sinus_Continu);
hold on;
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
bode(tf(1,[1.7483E-23 7.6663E-18 1.162E-11 3.0332E-6 1]))