close all;
T = 0.0002;
Fe = 320000;
Te = 1/Fe;
Tsim = T-Te;
M = 64;
deltaF = 5000;

Res=sim('simul_carre');

figure(1);
plot(Res.Sinus_Echanti);
grid;
title('tracé temporel');

freq =deltaF*[0:1:M-1];
yt = fft(Res.Sinus_Echanti.Data);
figure(2);
stem(freq,abs(yt)/M,'r*');
grid;

figure(3);
num = [1];
den = [1.783*10.^(-23) 7.6663*10.^(-18) 1.162*10.^(-11) 3.0332*10.^(-6) 1];
h = tf(num,den);
bode(h);