ENSEEIHT/TP-traitement-signal
1
0
Fork 0
Code Issues Pull requests Packages Projects Releases Wiki Activity

init

Browse source
This commit is contained in:
Laureηt 2023-06-10 21:09:28 +02:00
commit 61a45d296b
Signed by: Laurent
SSH key fingerprint: SHA256:kZEpW8cMJ54PDeCvOhzreNr4FSh6R13CMGH/POoO8DI
2 changed files with 97 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

36
TP1/main.m Executable file
View file

@ -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)));

61
TP2/main.m Executable file
View file

@ -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;