commit 61a45d296bbd81a9d9c09c3dfca62aefba97e5a4
Author: Laureηt <laurent@fainsin.bzh>
Date:   Sat Jun 10 21:09:28 2023 +0200

    init

diff --git a/TP1/main.m b/TP1/main.m
new file mode 100755
index 0000000..0bc71b7
--- /dev/null
+++ b/TP1/main.m
@@ -0,0 +1,36 @@
+clear;
+
+f0 = 1100;
+N = 90;
+
+Fe1 = 10000;
+Te1 = 1/Fe1;
+T1 = [0:Te1:(N-1)*Te1];
+Y1 = cos(2*pi*f0 * T1);
+
+% beau cosinus
+% plot(T1, Y1);
+
+N_prime = 2^10;
+% fft discrète + zero padding
+% semilogy([0:N_prime-1]/N_prime*Fe1, abs(fft(Y1, N_prime)));
+
+% chelou
+% periodogram(Y1);
+
+% periodogram maison
+% plot( [0:N_prime-1]/N_prime*Fe1, abs(fft(Y1, N_prime)).^2/N_prime );
+
+% broken
+% plot(xcorr(Y1));
+
+
+
+Fe2 = 1000;
+Te2 = 1/Fe2;
+T2 = [0:Te2:(N-1)*Te2];
+Y2 = cos(2*pi*f0 * T2);
+
+% plot(T2, Y2);
+% semilogy([0:N-1]/N*Fe2, abs(fft(Y2)));
+
diff --git a/TP2/main.m b/TP2/main.m
new file mode 100755
index 0000000..5024976
--- /dev/null
+++ b/TP2/main.m
@@ -0,0 +1,61 @@
+clear;
+
+% exercice 1
+
+f1 = 1000;
+f2 = 3000;
+N = 100;
+
+Fe = 10000;
+Te = 1/Fe;
+
+T = [0:Te:(N-1)*Te];
+
+Y1 = cos(2*pi*f1 * T);
+Y2 = cos(2*pi*f2 * T);
+
+Ys = Y1 + Y2;
+% plot(T, Ys);
+% hold;
+
+% fft discrète
+% semilogy([0:N-1]/N*Fe, abs(fft(Ys, N)));
+
+
+
+% exercice 2
+
+fc_tild = 1000/Fe;
+A = 1;
+
+ordre = 10;
+B = 2*fc_tild * sinc(2*fc_tild * [-ordre/2:1:ordre/2]);
+%plot(B);
+%hold;
+
+ordre = 60;
+B = 2*fc_tild * sinc(2*fc_tild * [-ordre/2:1:ordre/2]);
+%plot(B);
+
+
+
+% exercice 3
+
+fc_tild = 1000/Fe;
+A = 1;
+
+ordre = 10;
+B = 2*fc_tild*sinc(2*fc_tild * [-ordre/2:1:ordre/2]);
+Yf = filter( B, A, Ys);
+semilogy([0:N-1]/N*Fe, abs(fft(Yf, N)));
+
+hold;
+
+ordre = 60;
+B = 2*fc_tild*sinc(2*fc_tild * [-ordre/2:1:ordre/2]);
+Yf = filter( B, A, Ys);
+semilogy([0:N-1]/N*Fe, abs(fft(Ys, N)));
+
+
+
+legend;
\ No newline at end of file