TRAITEMENT DE SIGNAL
FICHE DE RAPPORT
SÉANCE #2 (TP2)
FENÊTRES TEMPORELLES
NOM ET PRÉNOM CLASSE
#1
#2
OBJECTIFS
Étudier et comparer l’effet de différentes fenêtres temporelles utilisées pour l’analyse
spectrale de signaux de longues durées.
1
PROGRAMME A COMPLETER
Fe=40000; % Définition de la fréquence d'échantillonnage
Te=1/Fe; % Période d'échantillonnage Te=1/Fe
% Définition des fenêtres de pondération
%*************************************************************
N= % Nombre d'échantillons de la fenêtre
Nfft = % Nombre de points de calcul des différentes FFT
%=============================================================
% Fenêtres dans le domaine temporel
y1= boxcar(N)'; % Fenêtre rectangulaire
y2= blackman(N)'; % Fenêtre de Blackman
%=============================================================
%*************************************************************
% Comparaison des deux fenêtres de pondération
% Deux sinusoïdes de fréquences f1 et f2 et d’amplitude unité
%*************************************************************
f1 = 8000; % Fréquence du premier sinus
f2 = 7200; % Fréquence du deuxième sinus
temps=(0:Te:(N-1)*Te); % les valeurs du temps pour lesquelles
on calculera d(t)
%vecteurs frequences
freq=linspace(0,Fe,Nfft+1);
freq=freq(1:end-1);
freq=freq-Fe/2;
s= ; % Signal à traiter de longueur N échantillons
sig1 = ; Signal tronqué par rectangle
sig_f1=fftshift(abs(fft(sig1,Nfft)));
sig2 = ; Signal tronqué par Blackman
sig_f2=fftshift(abs(fft(sig2,Nfft)));
figure
plot(freq,sig_f1,'b');
hold on
plot(freq,sig_f2,'g');
hold off
grid
title('Deux sinus à 7200Hz et 8000Hz d''amplitude égale à 1')
legend('Rectangulaire','Blackman'),
xlabel('Fréquence en Hz')
ylabel('Module de la TF');
%*************************************************************
% Comparaison des deux fenêtres de pondération
% Deux sinusoïdes de fréquences éloignées mais d'amplitude différente 1 volt et 0.01
volt
%*************************************************************
f1 = 8000; % Fréquence du premier sinus
f2 = 3200; % Fréquence du deuxième sinus
s= ;% Signal à traiter de longueur N échantillons
sig1= ; Signal tronqué par un rectangle
sig_f1=fftshift(abs(fft(sig1,Nfft)));
2
sig2= ; Signal tronqué par Blackman
sig_f2=fftshift(abs(fft(sig,Nfft)));
figure
plot(freq,sig_f1,'b');
hold on
plot(freq,sig_f2,'g');
hold off
grid
title('Deux sinus à 3200Hz et 8000Hz d''amplitudes égale à 0.01 et 1')
legend('Rectangulaire','Blackman'),
xlabel('Fréquence en Hz')
ylabel('Module de la TF');
Contenu du Rapport
Le Rapport doit contenir :
Le programme Matlab ci-dessus complété ;
des figures illustrant la supériorité d’une fenêtre sur l’autre en ce qui a trait à l’étalement
spectral. Quelle fenêtre est meilleure en ce qui concerne l’étalement spectral spectrale ?
des figures illustrant la supériorité d’une fenêtre sur l’autre en ce qui a trait à la
résolution spectrale. Quelle fenêtre est meilleure en ce qui concerne la résolution
spectrale ?