close all;
%% frequences theorique utilisees
F1=85000;
F2=90000;
F3=95000;
F4=100000;
F5=115000;
F6=120000;


%% frequences avec une imprecision
F1 = 85005.9;
F2 = 90000;
F3 = 94986.8;
F4 = 100000;
F5 = 115015.9;
F6 = 120000;

%% calcul et affichage du systeme de son FFT et du spectre de fourrier (Lancer une des sections en haut avant)
Res2 = sim('SimuProject');

figure(5);
plot(Res.Sinus_Continu);
hold on;
plot(Res.Sinus_Echanti,'o');
grid;

figure(6)
Tf = abs(fft(Res2.Sinus_Echanti.Data)/M);
stem(Tf);
title("fft de Sinus");
grid;

figure(7);
df = 1/T;
frequence = df * [0:1:M-1];
semilogy(frequence, Tf,'r*');
grid;
%% cas du pistolet avec des signaux carres

Scarre = sim('Signal_carre');

figure(1);
plot(Scarre.Carre_conti)
hold on;
plot(Scarre.Carre_discret,'o');

Tf = abs(fft(Scarre.Carre_discret.data)/M);

figure(3);
stem(frequence, Tf);

%% affichage bode du filtre de chebychev
F = tf([1],[1.7483E-23,7.6663E-18, 1.162E-11, 3.0332E-6 , 1]);

bode(F);