ENIG TP Traitement de signal AU 2023/2024
GCR2
TP3
Traitement Statistique des signaux discrets
I. But
Le but de ce TP est d'analyser des estimateurs de la fonction de corrélation et de la densité
spectrale de puissance (DSP).
1. Estimation des fonctions de corrélation
la fonction de corrélation se définit de différentes façons suivant la classe de signaux à
laquelle on s'adresse:
signaux déterministes
+∞
Energie finie: 𝑅𝑥,𝑦 (𝜏) = ∫−∞ 𝑥(𝑡) 𝑦 ∗ (𝑡 − 𝜏)𝑑𝑡
1 𝑇
Puissance finie: 𝑅𝑥,𝑦 (𝜏) = lim𝑇→∞ 𝑇 ∫0 𝑥(𝑡) 𝑦 ∗ (𝑡 − 𝜏)𝑑𝑡
signaux aléatoires
𝑅𝑥,𝑦 (𝜏) = 𝐸[𝑥(𝑡)𝑦 ∗ (𝑡 − 𝜏)]
Lorsque x=y, on parle de la fonction d’autocorrélation ( 𝑅𝑥𝑥 (𝜏) = 𝑅𝑥 (𝜏) ).
2. Propriétés de la fonction d'autocorrélation
Parité: 𝑅𝑥 (−𝜏) = 𝑅𝑥 (𝜏)
Maximum en zéro: |𝑅𝑥 (𝜏)| ≤ 𝑅𝑥 (0)
Puissance moyenne du signal = 𝑅𝑥 (0)
3. Propriétés de la fonction d'inter corrélation
∗ (−𝜏)
Symétrie hermitienne : 𝑅𝑥,𝑦 (𝜏) = 𝑅𝑥,𝑦
1
Majoration: |𝑅𝑥,𝑦 (𝜏)| ≤ 2 (𝑅𝑥 (0) + 𝑅𝑦 (0))
4. Algorithme de calcul
L'estimation de la fonction de corrélation fait appel aux algorithmes temporels et
fréquentiels.
Dans le calcul temporel, 𝑅𝑥,𝑦 (𝑘) est estimé à partir de 𝑁 − 𝑘 valeurs:
1
ENIG TP Traitement de signal AU 2023/2024
GCR2
1
𝑅̂𝑥𝑦 (𝑘) = 𝑁 ∑𝑁−𝑘−1
𝑛=0 𝑥(𝑛)𝑦(𝑛 + 𝑘) 0≤ 𝑘 ≤ 𝑁−1
1
Lorsque 𝑁 → 𝑁 − 1, le terme de normalisation est égal à 𝑁 . Par conséquent, un biais est
introduit dans l'estimation: la corrélation est pondérée par une fenêtre triangulaire.
Pour éliminer ce biais, un second estimateur peut être défini par:
1
𝑅̂𝑥𝑦 (𝑘) = 𝑁−𝑘 ∑𝑁−𝑘−1
𝑛=0 𝑥(𝑛)𝑦(𝑛 + 𝑘) 0≤ 𝑘 ≤ 𝑁−1
Par calcul fréquentiel, le calcul est basé sur une fonction de corrélation. Les estimateurs
deviennent:
1
Estimateur biaisé: 𝑅̂𝑥𝑦 (𝑘) = 𝑁 𝐹𝐹𝑇 −1 (𝑋(𝑓)𝑌 ∗ (𝑓))
1
Estimateur non biaisé: 𝑅̂𝑥𝑦 (𝑘) = 𝑁−𝑘 𝐹𝐹𝑇 −1 (𝑋(𝑓)𝑌 ∗ (𝑓))
où 𝑋(𝑓) = 𝐹𝐹𝑇(𝑥(𝑘)) et 𝑌(𝑓) = 𝐹𝐹𝑇(𝑦(𝑘)), avec 𝑛 = 0, . . . , 𝑁 où 𝑁 représente le
nombre de points calculés en fréquence.
II. Estimation de la densité spectrale de puissance
Pour l'estimation du spectre, il existe deux techniques classiques :
Le PERIODOGRAMME
Le CORRELOGRAMME
FFT
𝒙(𝒌) 𝑿(𝒇) ̂𝒙 (𝒌) = 𝟏 |𝑿(𝒇)|𝟐
𝑺
𝑵
signal temporel PERIODOGRAMME
FFT
𝑪𝒙 (𝒌) ̂𝒙 (𝒏) = 𝑭𝑭𝑻(𝑪𝒙 (𝒌))
𝑺
Autocorrélation CORRELOGRAMME
2
ENIG TP Traitement de signal AU 2023/2024
GCR2
III. Travail demandé :
Autocorrélation
1) Autocorrélation d'un sinus
a) Générer 128 échantillons d'une sinusoïde de fréquence 15𝐻𝑧 sur un intervalle
[0𝑠, 1𝑠]. On la note 𝑥1. Observer son spectre X1.
b) Observer les estimations de son autocorrélation.
c) Quelle est l'allure du biais ?
2) Autocorrélation d'un sinus bruité
a) Générer un bruit blanc avec la fonction 𝑟𝑎𝑛𝑑𝑛.
b) Observer les estimations de son autocorrélation. Laquelle parait la plus satisfaisante ?
Retrouver les caractéristiques du bruit.
c) Ajouter au sinusoïde le bruit blanc généré.
d) Observer le signal bruité et l'estimation biaisée de son autocorrélation.
Estimations spectrales
1) Périodogramme Simple
a) Calculer et comparer les périodogrammes des fonctions sinus et sinus bruité.
b) Que représente la valeur du périodogrammes pour la fréquence 0 ?
c) Trouver la variance du signal.
2) Corrélogramme
a) A l'aide de l'estimateur biaisé de l'autocorrélation, comparer le périodogramme et le
corrélogramme de la fonction sinus bruité générée.
Expliquer.
b) A l'aide de l'estimateur non biaisé de l'autocorrélation, représenter le corrélogramme
de la fonction sinus bruité générée. Quel est le principal inconvénient ?
Limitations de Fourier: Etude de la résolution de fréquence
a) Générer un sinusoïde 𝑥2 de fréquence 𝑓2 = 16𝐻𝑧. Observer 𝑥2 ainsi que son spectre
𝑋2.
b) Ajouter le signal 𝑥2 au signal 𝑥1. On le note 𝑦1. Observer son spectre 𝑌1.
3
ENIG TP Traitement de signal AU 2023/2024
GCR2
c) Observer les estimations de corrélation possibles.
d) Observer les différentes fenêtres temporelles d'apodisation (traingulaire, Hamming,
Hanning).
e) Pondérer le signal sinusoïde généré par ces fenêtres. Observer leur spectres centrés à
zéro.
f) Comparer les effets des différentes fenêtres en fonction de leur pouvoir à réduire
l'amplitude des lobes secondaires.
g) Qu’elle est la résolution de fréquence ?
h) Ajouter un signal 𝑥3 de fréquence 𝑓3 = 55𝐻𝑧 au signal 𝑥1. On le note 𝑦2
i) Recommencer les questions b) et d) avec le signal 𝑦2.
j) Que devient la résolution de fréquence ?