Traitement numérique du signal
Formulaire
Laurent Oudre
[email protected]
2017-2018
Université Paris 13, Institut Galilée, Ecole d’ingénieurs Sup Galilée
Parcours Informatique et Réseaux Alternance - 1ère année
Etude des signaux dans le domaine temporel
Energies et puissances
Energie totale
Z +∞ +∞
X
2
Ex = |x(t)| dt Ex = |xn |2
−∞ n=−∞
Puissance moyenne totale
τ m
1 1
Z X
Px = lim |x(t)|2 dt Px = lim |xn |2
τ →+∞ 2τ −τ m→+∞ 2m + 1
n=−m
Puissance moyenne totale (cas périodique)
t0 +T n0 +M−1
1 1
Z X
∀t0 ∈ R, Px = |x(t)|2 dt ∀n0 ∈ Z, Px = |xn |2
T t0 M n=n0
Signaux usuels
Signal porte centré de durée L
(
1 si − L2 ≤ t < L
2
ΠL (t) =
0 sinon
Distribution de Dirac (continu)
( Z +∞
+∞ si t = 0
δ(t) = δ(t)dt = 1 x(t) × δ(t − t0 ) = x(t0 )δ(t − t0 ) x(t) ∗ δ(t − t0 ) = x(t − t0 )
0 sinon −∞
Distribution de Dirac (discret) (
1 si n = 0
δn =
0 sinon
Peigne de Dirac
+∞
X
XT (t) = δ(t − nT )
n=−∞
Sinus cardinal (
1 si t = 0
sinc(t) = sin(πt) Annulations pour t = k ∈ Z
πt sinon
Etude des signaux dans le domaine fréquentiel
Transformées de Fourier : définitions
Transformée de Fourier continue
Z +∞ Z +∞
X(f ) = x(t) e−j2πf t dt x(t) = X(f ) ej2πf t df
−∞ −∞
Transformée de Fourier discrète
Fs Fs
N : nombre d’échantillons Fréquences observables : fk = k avec k = 0, · · · , N −1 Résolution en fréquence : ∆f =
N N
N −1 N −1
X kn 1 X kn
Xk = xn e−j2π N xn = Xk ej2π N
n=0
N
k=0
Transformées de Fourier : propriétés
Linéarité
T F {λx(t) + µy(t)} = λT F {x(t)}+µT F {y(t)} T F −1 {λX(f ) + µY (f )} = λT F −1 {X(f )}+µT F −1 {Y (f )}
Produit/Produit de convolution
T F {x(t) × y(t)} = T F {x(t)} ∗ T F {y(t)} T F {x(t) ∗ y(t)} = T F {x(t)} × T F {y(t)}
Translation/Modulation
T F {x(t − t0 )} = e−j2πf t0 T F {x(t)} T F ej2πf0 t x(t) = X(f − f0 )
Divers
Module de la transformée de Fourier pair et argument de la
Signal réel ⇔ transformée de Fourier impair
Transformée de Fourier réelle (et signal réel) ⇔ Signal pair (et réel)
Signal périodique ⇔ Spectre discret
Quelques transformées de Fourier
T F {δ(t − t0 )} = e−j2πf t0 T F ej2πf0 t = δ(f − f0 )
T F {δ(t)} = 1 T F {1} = δ(f )
δ(f − f0 ) + δ(f + f0 ) δ(f − f0 ) − δ(f + f0 )
T F {cos (2πf0 t)} = T F {sin (2πf0 t)} =
2 2j
1
T F {XT (t)} = X T1 (f ) T F {ΠL (t)} = L sinc (Lf )
T
Filtrage des signaux
Produit de convolution
Définitions
Z +∞ Z +∞ +∞
X +∞
X
(f ∗ g)(t) = f (t − τ ) g(τ )dτ = g(t − τ ) f (τ )dτ (f ∗ g)n = fn−m gm = gn−m fm
−∞ −∞ m=−∞ m=−∞
Propriétés
f ∗δ =f f ∗g =g∗f f ∗ (g + h) = f ∗ g + f ∗ h (f ∗ g) ∗ h = f ∗ (g ∗ h)
Propriétés des filtres
Caractérisations
h(t) (ou hn ) : réponse impulsionnelle → sortie du filtre avec un Dirac en entrée
H(f ) = T F {h(t)} : fonction de transfert du filtre
Bandes passantes à -3 dB
1
|H(fc )|2 = max |H(f )|2 |H(fc )|2 dB = 10 log10 max |H(f )|2 − 3
2
Signaux aléatoires
Densité spectrale de puissance
Définition Z +∞
Px = Γx (f )df
−∞
Filtrage
y =h∗x → Γy (f ) = |H(f )|2 Γx (f )
Bruits
Bruit blanc
N0
Γb (f ) = σ 2 =
2
Rapport signal sur bruit
Px Px
SN R = SN RdB = 10 log10
Pb Pb