FORMULAIRE POUR LA THÉORIE DU SIGNAL

Mathematics, Signal Processing · course

Voir tous les documents en électronique et automatique

FORMULAIRE POUR LA THÉORIE DU SIGNAL

Christian JUTTEN

Laboratoire GIPSA, Département Signal et Images (CNRS, INPG, UJF), Grenoble, France

[email protected]

1. FORMULES TRIGONOMÉTRIQUES

2. MODÈLES MATHÉMATIQUES DE SIGNAUX

1.1. Somme et différence d’angles

2.1. Signaux simples

• sin(α + β) = sin α cos β + cos α sin β,

• sin(α − β) = sin α cos β − cos α sin β,

• rect[(t−τ )/T ] : fonction rectangle centrée sur τ , d’amplitude

1 et de largeur T (l’aire vaut T ),

• cos(α + β) = cos α cos β − sin α sin β,

• tri[(t−τ )/T ] : fonction triangle centrée sur τ , d’amplitude

• cos(α − β) = cos α cos β + sin α sin β.

1.2. Expression en fonction des angles doubles

• 2 sin2 α = 1 − cos 2α,

• 2 cos2 α = 1 + cos 2α,

• 2 sin α cos α = sin 2α.

1.3. Produits de fonctions trigonométriques

• 2 sin α sin β = cos(α − β) − cos(α + β),

• 2 cos α cos β = cos(α − β) + cos(α + β),

• 2 sin α cos β = sin(α − β) + sin(α + β).

1.4. Somme et différence de fonctions trigonométriques

1 et de largeur 2T (l’aire vaut T ),

• δ(t) : distribution de Dirac, satisfaisant δ(t) = 0 pour

t 6= 0, et

+

−∞

R

δ(t)dt = 1,

• sinc(u) = sin(πu)

πu

1).

: sinus cardinal de u (son aire vaut

2.2. Opérateur de répétition

• repT (.) : répétition de . avec une période T ,

• repT (x(t)) =

+

∞k=

−∞

P

x(t − kT ),

• repT (δ(t)) = δT (t) =

Dirac.

+

∞k=

−∞

P

δ(t − kT ) : peigne de

• sin α + sin β = 2 sin[(α + β)/2] cos[(α − β)/2],

2.3. Opérateur de convolution

• sin α − sin β = 2 sin[(α − β)/2] cos[(α + β)/2],

Définition

• cos α + cos β = 2 cos[(α + β)/2] cos[(α − β)/2],

• cos α − cos β = −2 sin[(α + β)/2] sin[(α − β)/2].

• x(t) ∗ y(t) = (x ∗ y)(t) =

x(t − v)y(v)dv,

+

−∞

R

x(u)y(t − u)du =

+

−∞

R

1.5. Formes exponentielles

• exp(jα) = cos α + j sin α, avec j2 = −1,

• sin α = exp(jα)

exp(

2j

jα)

• cos α = exp(jα)+exp(

2

jα)

,

.

• x(t) ∗ y(t − t0) =

y)(t − t0),

• x(t) ∗ y(at + b) =

+

−∞

R

+

−∞

R

Propriétés

x(u)y(t − u − t0)du = (x ∗

x(u)y(a(t − u) + b)du

1.6. Développements limités

• Commutativité : x(t) ∗ y(t) = y(t) ∗ x(t),

• sin α = α − α3/3! + α5/5! . . . + (−1)pα2p+1/(2p +

1)! + . . .,

• cos α = 1− α2/2! + α4/4! . . . + (−1)pα2p/(2p)! + . . .,

• Associativité : [x(t)∗y(t)]∗z(t) = x(t)∗[y(t)∗z(t)] =

x(t) ∗ y(t) ∗ z(t),

• Distributivité par rapport à l’addition : x(t) ∗ [y(t) +

• exp α = 1 + α + α2/2! + . . . + αk/k! + . . ..

z(t)] = x(t) ∗ y(t) + x(t) ∗ z(t)

2.4. Propriétés de la distribution de Dirac δ(t)

+

• x(t0) =

−∞

R

• x(t) ∗ δ(t) = x(t),

x(t)δ(t − t0)dt

• x(t) ∗ δ(t − t0) = x(t − t0),

• x(t − t1) ∗ δ(t − t2) = x(t − t1 − t2),

• δ(at) = |a|−

1δ(t).

3. VALEUR MOYENNE TEMPORELLE, ÉNERGIE

ET PUISSANCE

• ¯x = limT

+

1

T

+T /2

T /2 x(t)dt : moyenne.

R

Signaux réels

• Wx =

+

−∞

R

• Px = limT

Signaux complexes

+

1

T

x2(t)dt : énergie,

+T /2

T /2 x2(t)dt : puissance.

R

|x(t)|2dt : énergie,

• Wx =

+

−∞

R

• Px = limT

+

1

T

+T /2

T /2 |x(t)|2dt : puissance.

R

4. REPRÉSENTATION VECTORIELLE

On considère des signaux x(t) et y(t) dans l’espace L2(t1, t2).

4.1. Distance euclidienne

• d(x, y) = kx − yk =

t2

t1 |x(t) − y(t)|2dt

1/2

i

h R

• d(x, y)2 = kxk2 + kyk2 − 2ℜhx, yi

• kxk2 =

• hx, yi =

t2

,

t1 |x(t)|2dt

h R

i

t2

t1 x(t)y∗(t)dt.

R

4.2. Inégalité de Schwartz

• |hx, yi|2 ≤ hx, xi.hy, yi,

t2

t1 x(t)y∗(t)dt

R

2

R

t2

t1 |x(t)|2dt.

t2

t1 |y(t)|2dt.

R

(cid:12)

(cid:12)

(cid:12)

(cid:12)

(cid:12)

(cid:12)

4.3. Orthogonalité de x(t) et y(t)

• hx, yi =

R

t2

t1 x(t)y∗(t)dt = 0.

4.4. Développement sur une base orthogonale {ψk(t)}

• x(t) =

P

Publicité

• αk = 1

λkk

+

∞k=1 αkψk(t),

hx, ψki = 1

λkk

• λkk = hψk, ψki = kψkk2 =

t2

t1 |x(t)|2dt =

R

P

R

+

∞k=1 |αk|2λkk.

t2

t1 x(t)ψ∗k(t)dt,

R

t2

t1 |ψk(t)|2dt,

4.5. Développement en séries de Fourier (DSF) dans L2(t1, t1+

T )

On utilise ψk(t) = exp

j2π k

T t

.

(cid:0)

+

∞k=1 αk exp

(cid:1)

j2π k

T t

,

| exp

(cid:0)

j2π k

T t

(cid:1)

|2dt = T ,

hx, exp

(cid:0)

j2π k

T t

(cid:1)

i = 1

T

(cid:0)

(cid:1)

• x(t) =

• λkk =

P

t1+T

t1

R

• αk = 1

λkk

j2π k

dt.

T t

(cid:1)

t1+T

t1

R

x(t) exp

(cid:0)

5. TRANSFORMÉES DE FOURIER

5.1. Transformée de fonctions sommables

Pour des fonctions x(t) bornées et absolument sommables

(

|x(t)|dt < +∞), on définit la transformée de Fourier (TF),

R

notée X(f ) :

x(t) ⇋ X(f ) =

+

Z

−∞

x(t) exp(−j2πf t)dt.

De même, étant donnée X(f ), on peut calculer x(t) par trans-

formée de Fourier inverse :

X(f ) ⇌ x(t) =

+

Z

−∞

X(f ) exp(+j2πf t)df.

5.2. Notations

On notera

• x(t) ⇋ X(f ),

• X(f ) = F{x(t)} : transformée de Fourier,

• x(t) = F −

1{X(f )} : transformée de Fourier inverse.

5.3. Propriétés

• La TF est linéaire : ax(t) + by(t) ⇋ aX(f ) + bY (f )

• La TF conserve la parité.

• Si la fonction x(t) est réelle (imaginaire, resp.)

paire, X(f ) est imaginaire (réelle, resp.) impaire.

im-

5.4. Règles de calculs

Soit une fonction x(t), telle que x(t) ⇋ X(f ).

• x(−t) ⇋ X(−f ) : retournement temporel,

• x∗(t) ⇋ X ∗(−f ) : complexe conjuguée,

• x(at) ⇋ |a|−

1X(f /a) : effet Doppler,

• x(t − t0) ⇋ X(f ) exp(−j2πf t0) : translation tem-

porelle,

• x(t) exp(j2πf0t) ⇋ X(f −f0) : translation en fréquence

ou modulation,

• dn

dtn x(t) ⇋ (j2πf )nX(f ) : dérivation sur t.

Soient deux fonctions x(t) et y(t), telles que x(t) ⇋ X(f ) et

y(t) ⇋ Y (f ), on montre le Théorème de Plancherel :

• x(t) ∗ y(t) ⇋ X(f )Y (f ),

• x(t)y(t) ⇋ X(f ) ∗ Y (f ).

5.5. Transformées de Fourier de quelques signaux usuels

• rect(t/T ) ⇋ T sinc(f T ) : fonction rectangle,

• tri(t/T ) ⇋ T sinc2(f T ) : fonction triangle,

5.8.2. Transformées de Fourier

Les transformées de Fourier de signaux périodiques sont des

spectres de raies, aux fréquences n/∆, dont l’enveloppe spec-

trale est la transformée de Fourier d’une période du signal di-

visée par la période ∆.

• X(f ) =

dessus,

P

+

∞n=

−∞

• Y (f ) =

dessus.

P

+

∞n=

−∞

Xnδ(f − n/∆), avec Xn défini ci-

Ynδ(f − n/∆), avec Yn défini ci-

• exp(−at)ǫ(t) ⇋ 1

a+j2πf :

impulsion exponentielle

5.8.3. Transformées de signaux périodiques usuels

unilatérale,

• exp(−a|t|) ⇋

double,

2a

a2+(2πf )2 :

impulsion exponentielle

• cos(2πf0t) ⇋ 1

2 [δ(f + f0) + δ(f − f0)] : signal cosi-

nusoïdal,

• exp(−πt2) ⇋ exp(−πf 2) : impulsion gaussienne,

• sin(2πf0t) ⇋ j

2 [δ(f + f0) − δ(f − f0)] : signal sinu-

soïdal,

• exp(−at) sin(2πf0t)ǫ(t) ⇋

soïde amortie,

2πf0

(a+j2πf )2+(2πf0)2 : sinu-

• exp(j2πf0t) ⇋ δ(f − f0) : phaseur à la fréquence f0,

• exp(−at) cos(2πf0t)ǫ(t) ⇋

nusoïde amortie.

a+j2πf0

(a+j2πf )2+(2πf0)2 : cosi-

5.6. Transformées de Fourier de distributions

• ǫ(t) ⇋ 1

j2πf + δ(f )

2

: échelon unité,

• sgn(t) ⇋ 1

jπf : fonction signe,

• δ(t) ⇋ 1 : impulsion de Dirac,

• k ⇋ kδ(f ) : constante.

5.7. Transformée de Fourier de signaux à moyenne non

nulle

• Soit x(t) = ¯x + x0(t) où ¯x est la moyenne de x(t), on

a F{x(t)} = F {

0(t)

x

j2πf + ¯xδ(f ),

}

• ǫ(t) ⇋ 1

2 δ(f ) + 1

j2πf ,

• sgn(t) ⇋ 1

jπf .

k δ(t−k∆) ⇋ 1

n δ(f −n/∆) : peigne de Dirac,

P

• δ∆(t) ⇋ 1

tions abrégées,

P

∆ δ1/∆(f ) : peigne de Dirac avec les nota-

+

∞n=

Xn exp(j2πnt/∆) ⇋

signal x(t) périodique,

P

−∞

+

∞n=

−∞

P

Xnδ(f −n/∆):

• Arep∆{rect(t/θ} ⇋

∆ sinc(nθ/∆).

P

+

Publicité

∞n=

−∞

Xnδ(f −n/∆) avec Xn =

6. CORRÉLATIONS ET DENSITÉS SPECTRALES

6.1. Signaux à énergie finie

6.1.1. Auto- et inter-corrélation

• Auto-corrélation : Γxx(τ ) =

R

• Inter-corrélation : Γxy(τ ) =

R

+

−∞

+

−∞

x(t)x∗(t − τ )dt,

x(t)y∗(t − τ )dt.

6.1.2. Relation entre corrélation et convolution

• Γxy(τ ) = x(τ ) ∗ y∗(−τ ).

5.8. Transformées de Fourier de signaux périodiques

6.1.3. Propriétés des fonctions de corrélation

Soient deux signaux x(t) et y(t), périodiques de même péri-

ode ∆.

5.8.1. Développement en séries de Fourier (DSF)

Les signaux étant périodiques, si les conditions de Diriclet

sont satisfaites, on peut calculer les DSF de x(t) et y(t) :

• x(t) =

+

∞n=

avec Xn = [

P

−∞

Xn exp(j2πnt/∆),

+∆/2

∆/2 x(t) exp(−j2πnt/∆)dt]/∆,

• y(t) =

Yn exp(j2πnt/∆),

avec Yn = [

P

−∞

+∆/2

∆/2 y(t) exp(−j2πnt/∆)dt]/∆.

R

+

∞n=

R

• Symétrie hermitienne :

◦ Γxy(τ ) = Γ∗yx(−τ ),

◦ Γxx(τ ) = Γ∗xx(−τ ),

• Bornes

◦ Cas général : |Γxy(τ )|2 ≤ Γxx(0)Γyy(0),

◦ Pour x(t) = y(t) : |Γxx(τ )| ≤ Γxx(0),

◦ Γxx(0) ∈ R est l’énergie du signal.

• Dans le cas de signaux réels, la fonction d’autocorrélation

est paire et maximale en 0, où Γxx(0) est l’énergie du

signal.

6.1.4. Densité spectrale d’énergie

• Γxx(τ ) ⇋ Sxx(f ) =

Γxx(τ ) exp(−j2πf τ )dτ ,

R

• Sxx(f ) = |X(f )|2 : densité spectrale d’énergie,

• Γxy(τ ) ⇋ Sxy(f ) =

Γxy(τ ) exp(−j2πf τ )dτ :

densité inter-spectrale d’énergie,

+

−∞

+

−∞

R

• Sxy(f ) = X(f )Y (f )∗,

6.1.5. Relations de Parseval

x(t)y∗(t)dt =

X(f )Y ∗(f )df ,

+

−∞

+

R

−∞

R

0,

Γxy(0) =

+

−∞

R

+

−∞

R

Sxx(f )df : énergie du signal.

+

−∞

R

|x(t)|2dt =

|X(f )|2df , ou, pour le cas τ =

6.1.6. Dérivation de la fonction de corrélation

• dΓxy(τ )

dτ = Γx′y(τ ) = −Γxy′ (τ )

6.2. Signaux à puissance moyenne finie

6.2.1. Fonctions de corrélation

• Γxx(τ ) = limT

corrélation,

→∞

• Γxy(τ ) = limT

corrélation.

→∞

1

T

1

T

+T /2

T /2 x(t)x∗(t − τ )dt : auto-

R

+T /2

T /2 x(t)y∗(t − τ )dt : inter-

R

6.2.2. Densités spectrale et interspectrale de puissance

• Sxx(f ) = F{Γxx(τ )} : densité spectrale de puissance,

• Sxy(f ) = F{Γxy(τ )} : densité inter-spectrale de puis-

sance,

• Sxx(f ) = limT

X(f, T ) = F{x(t, T )} = F{rect(t/T )x(t)}.

1

T |X(f, T )|2, avec

→∞

6.2.3. Relation de Parseval

• Px = Γxx(0) =

Sxx(f ).

+

−∞

R

6.3. Signaux périodiques

• Inter-corrélation

◦ Γxy(τ ) = 1

+∆/2

∆/2 x(t)y∗(t − τ )dt,

t0+∆

t0

x(t)y∗(t − τ )dt,

R

◦ Γxy(τ ) = 1

◦ Γxy(τ ) =

R

+

∞n=

isant le DSF de x(t) et y(t).

−∞

P

XnY ∗n exp(j2πnτ /∆), en util-

6.3.2. Densités spectrale et inter-spectrale de puissance

Les densités spectrale et inter-spectrale de puissance sont les

transformées de Fourier de fonctions périodiques de périodes

∆. On obtient donc des spectres de raies, aux fréquences

k/∆, dont l’enveloppe spectrale est la transformée de Fourier

d’une période de la fonction de corrélation divisée par la péri-

ode ∆.

• Densité spectrale

◦ Sxx(f ) =

+

∞n=

|Xn|2δ(f − n/∆), en util-

−∞

isant le DSF de x(t),

P

∆2 |X(f, ∆)|2δ1/∆(f ), avec

X(f, ∆) = F{x(t, ∆)} = F{rect((t)/∆)x(t)}.

◦ Sxx(f ) = 1

• Densité inter-spectrale

◦ Sxy(f ) =

+

∞n=

XnY ∗n δ(f − n/∆), en util-

−∞

isant le DSF de x(t),

P

◦ Sxy(f ) = 1

∆2 X(f, ∆)Y ∗(f, ∆)δ1/∆(f ), avec

X(f, ∆) = F{x(t, ∆)} = F{rect((t)/∆)x(t)}

et Y (f, ∆) = F{y(t, ∆)} = F{rect((t)/∆)y(t)}.

6.3.3. Relations de Parseval

• Px = Γxx(0) =

+

−∞

R

Sxx(f )df =

+

∞n=

−∞

|Xn|2.

P

7. SIGNAUX ALÉATOIRES

7.1. Variables aléatoires

7.1.1. Densités de probabilité

On note pX (u) la densité d’une variable aléatoire (VA) X.

• Si X et Y sont deux VA indépendantes, la densité de la

Publicité

VA Z = X + Y est : pZ(u) = (pX ∗ py)(u),

Soient deux signaux x(t) et y(t), périodiques de même péri-

ode ∆.

• Si Y = f (X), avec f inversible,

pY (u) = pX (f −

1(u))/|f ′(f −

1(u))|

6.3.1. Fonctions de corrélation

7.1.2. Moments

Les auto et inter-corrélations sont périodiques de période ∆,

et on peut les calculer par intégration sur une période ∆.

• Auto-corrélation

◦ Γxx(τ ) = 1

◦ Γxx(τ ) = 1

◦ Γxx(τ ) =

R

+

∞n=

isant le DSF de x(t).

P

−∞

+∆/2

∆/2 x(t)x∗(t − τ )dt,

t0+∆

t0

x(t)x∗(t − τ )dt,

R

|Xn|2 exp(j2πnτ /∆), en util-

Moyenne d’une VA continue ou discrète

+

• µ1 = E[X] =

• µ1 = E[X] =

R

upX (u)du,

−∞

i uiPr(X = ui).

P

Moment d’ordre k d’une VA continue ou discrète

• µk = E[X k] =

• µk = E[X k] =

R

+

−∞

i uk

ukpX (u)du,

i Pr(X = ui).

P

Moment centré d’ordre k d’une VA continue ou discrète

7.3. Auto-corrélation

• µ′k = E[(X − µ1)k] =

• µ′k = E[(X − µ1)k] =

R

(u − µ1)kpX (u)du,

−∞

i(ui − µ1)kPr(X = ui).

+

P

Variance : c’est le moment centré d’ordre 2. On la note

fréquemment σ2 plutôt que µ′2.

• σ2 = µ′2 = E[(X − µ1)2] =

• σ2 = µ′2 = E[(X −µ1)2] =

R

(u − µ1)2pX (u)du,

−∞

i(ui−µ1)2Pr(X = ui).

+

Relation entre variance et moments d’ordre 1 et 2

P

• σ2 = E[(X − µ1)2] = E[X 2] − E[X]2 = µ2 − µ2

1.

7.1.3. Distributions de quelques lois usuelles

• Pr(k) = C k

npk(1 − p)n

k : loi Binomiale, de la VA

somme de n VA binaires, avec p = Pr(X = x1) et

1 − p = Pr(X = x0),

• Pr(n) = an

n! exp(−a) :

a > 0,

loi de Poisson de paramètre

• p(u) = 1

− (u

2σ2

moyenne m et de variance σ2, notée N (m, σ2).

√2πσ

exp

m)

(cid:16)

(cid:17)

:

loi gaussienne de

2

7.2. Vecteurs aléatoires

7.2.1. Densités de probabilité

On note pX(u1, . . . , un) la densité d’un vecteur aléatoire (VcA)

X de dimension n.

Soit Y = f (X), avec f inversible. On note f −

mation inverse, |Jf −1| son jacobien et (f −

i de la transformation f −

1 la transfor-

1)i la composante

1 :

• pY(u1, . . . , un) =

pX((f −

1)1(u1), . . . , (f −

1)n(un))|Jf −1|.

Dans le cas d’une transformation linéaire Y = AX où A est

une matrice régulière et u = (u1, . . . , un)T :

• pY(u1, . . . , un) =

1u)| det A−

pX(A−

1|.

Dans le cas de la somme Z = X + Y de deux variables

aléatoires, la densité pZ(u) vaut :

• pZ(u) =

pXY (x, u − x)dx,

R

• pZ(u) =

pX (x)pY (u − x)dx = (pX ∗ pY )(u), si X

et Y sont des VA indépendantes.

R

7.2.2. Théorème de la limite centrale

Théorème 7.2.1 La distribution statistique d’une somme de

n variables aléatoires indépendantes, possédant la même loi

tend asymptotiquement vers une distribution gaussienne quelle

que soit la distribution des termes individuels.

On note de façon générale Xi = x(ti) la VA associée au

processus aléatoire x(t) pour t = ti.

7.3.1. Auto-corrélation statistique

Cas général

• RXX (t1, t2) = E[X1X2] =

x1x2p(x1, x2)dx1dx2.

Cas d’un signal stationnaire, en notant τ = t1 − t2

R R

• RXX (τ ) = E[X(t)X(t−τ )] =

x1x2p(x1, x2)dx1dx2.

Les VA sont orthogonales ssi RXX (τ ) = 0.

R R

7.3.2. Auto-covariance statistique

C’est l’auto-corrélation des processus aléatoires centrés. On

la note CXX (τ ).

Cas général

• CXX (t1, t2) = E[(X1 − E(X1))(X2 − E(X2))]

(x1 − E(X1))(x2 − E(X2))p(x1, x2)dx1dx2

=

= RXX (t1, t2) − E(X1)E(X2).

R R

Cas d’un signal stationnaire, en notant τ = t1 − t2

• CXX (τ ) = E[(X(t) − E(X))(X(t − τ ) − E(X))]

x1x2p(x1, x2)dx1dx2

=

= RXX (τ ) − E(X)2.

R R

Les VA sont non corrélées ssi CXX (τ ) = 0.

7.3.3. Cas de processus ergodiques et stationnaires

Dans ce cas, on peut remplacer les moyennes statistiques par

des moyennes temporelles, car on a asymptotiquement :

• RXX (τ ) = limT

→∞

1

T

R

T /2

T /2 x(t)x(t−τ )dt = ΓXX (τ ).

7.3.4. Coefficient de corrélation

Par définition, c’est l’auto-covariance normalisée, notée ρX (τ ) :

• ρX (τ ) = CXX (τ )

CXX (0) = CXX (τ )

σ2

X

.

7.3.5. Propriétés pour les signaux réels

Les fonction d’auto-corrélation et d’auto-covariance ont les

propriétés intéressantes pour des processus aléatoires réels.

De plus, en τ = 0, on a les relations énergétiques.

• elles sont paires,

• le maximum est atteint en τ = 0 :

◦ |CXX (τ )| ≤ CXX (0)

◦ |RXX (τ )| ≤ RXX (0)

◦ CXX (0) = σ2

X ,

◦ RXX (0) = σ2

X + µ2

X .

7.4. Densité spectrale de puissance (DSP)

• RXY (τ ) = E[X(t)Y (t − τ )] (cas stationnaire).

Pour un processus aléatoire x(t), on ne peut pas calculer di-

rectement la transformée de Fourier, puisque x(t) ne peut pas

être décrit. Par conséquent, on ne peut pas calculer directe-

ment sa DSP. On peut en revanche calculer sa fonction d’auto-

corrélation par une moyenne statistique ou temporelle, et en

déduire sa DSP en appliquant le théorème ci-dessous.

et d’inter-covariance :

• CXY (t1, t2) = E[(X(t1) − µX )(Y (t2) − µY )]

= RXY (t1, t2) − E[X(t1))E(Y (t2)] (cas général),

• CXY (τ ) = E[(X(t) − µX )(Y (t − τ ) − µY )]

= RXY (τ ) − µX µY (cas stationnaire).

7.4.1. Théorème de Wiener-Kintchine

Le théorème de Wiener-Kintchine établit le lien entre la fonc-

tion d’auto-corrélation et la DSP pour des signaux aléatoires.

7.5.2. Coefficient d’inter-corrélation

Le coefficient d’inter-corrélation est défini par :

Publicité

Théorème 7.4.1 La DSP d’un processus aléatoire station-

naire au sens large est le transformée de Fourier de sa fonc-

tion d’auto-corrélation :

• ρXY (τ ) = CXY (τ )

σX σY

.

7.5.3. Densité inter-spectrale

RXX (τ ) exp(−j2πf τ )dτ.

La densité inter-spectrale de puissance est la transformée de

Fourier de la fonction d’inter-corrélation :

SXX (f ) =

+

Z

−∞

7.4.2. Conséquence

Si on connaît la DSP SXX (f ) d’un processus aléatoire x(t),

on peut déduire la fonction d’auto-corrélation par transformée

de Fourier inverse :

RXX (τ ) =

+

Z

−∞

SXX (f ) exp(j2πf τ )df.

On peut aussi l’exprimer en fonction de l’auto-covariance, car

RXX (τ ) = CXX (τ ) + µ2

X :

• SXY (f ) = F{RXY (τ )}.

7.5.4. Cohérence

• γXY (f ) = |

2

SXY (f )

SXX (f )SY Y (f ) .

|

La cohérence est proche de 1 si les signaux x(t) et y(t) sont

liés par une relation linéaire (filtre). Dans le cas d’une relation

non linéaire ou en présence d’un bruit additif, la cohérence

devient très inférieure à 1.

SXX (f ) = F{CXX (τ )} + µ2

X δ(f ).

8. FILTRES LINÉAIRES INVARIANTS

En τ = 0, on a :

• RXX (0) = PX =

SXX (f )df .

On considère un filtre de réponse impulsionnelle g(t) et en

fréquence G(f ), dont les signaux d’entrée et de sortie sont

x(t) et y(t), respectivement.

+

−∞

R

Si le processus est aussi ergodique, on a

• RXX (0) = limT

→∞

SXX (f )df .

=

+

−∞

R

1

T

T /2

T /2 |x(t)|2dt

R

7.4.3. Bruit blanc

Définition 7.4.2 On appelle bruit blanc un processus aléa-

toire b(t) dont la densité spectrale de puissance est constante,

∀f :

SBB(f ) = η/2.

Par transformée de Fourier inverse, on déduit la fonction d’auto-

corrélation d’un bruit blanc b(t) :

RBB(τ ) = F −

1{SXX (f )} = δ(τ )η/2.

7.5. Inter-corrélation et densité inter-spectrale

7.5.1. Inter-corrélation et inter-covariance

On définit les fonctions d’inter-corrélation :

• RXY (t1, t2) = E[X(t1)Y (t2)] (cas général)

• Y (f ) = G(f )X(f ) = X(f )G(f ),

• y(t) = (g ∗ x)(t) = (x ∗ g)(t).

Pour une cascade de filtres, on a g(t) = (g1 ∗ g2 ∗ . . . ∗ gn)(t)

ou, en fréquence, G(f ) = G1(f )G2(f ) . . . Gn(f ).

8.1. Formule des interférences

On considère deux filtres de réponses impulsionnelles g1(t) et

g2(t) et de réponses en fréquence (TF) G1(f ) et G2(f ), dont

les signaux d’entrée sont x1(t) et x2(t), et de sortie y1(t) et

y2(t).

• SY1Y2 (f ) = SX1X2(f )G1(f )G∗2(f ).

8.2. Relation entrée-sortie des DSP

En conséquence de la formule des interférences, on a pour les

DSP :

• SY Y (f ) = |G(f )|2SXX (f ),

• SY X (f ) = G(f )SXX (f ).

8.3. Relation entrée-sortie entre corrélations

9.2. Opérateur de moyenne temporelle

En conséquence des formules entre DSP, par transformée de

Fourier inverse, on a pour les auto-corrélations :

• RY Y (τ ) = Γgg(τ ) ∗ RXX (τ ) pour des signaux aléa-

toires, et

• ΓY Y (τ ) = Γgg(τ ) ∗ ΓXX (τ ) si les signaux sont er-

godiques, ou certains à puissance moyenne finie ou à

énergie finie (attention aux définitions de Γ).

et pour les inter-corrélations :

• g(t) = rect((t − T /2)/T ),

• G(f ) = sinc(T f ) exp(−jπf T ),

• µy = G(0)µX = µX (moyenne),

• Γgg(τ ) = tri(τ /T )/T (auto-corrélation).

9.3. Opérateur de filtrage passe-bas idéal

Le filtre causal de largeur de bande B n’étant pas réalisable,

on décale g(t) de t0.

• RY X (τ ) = g(τ ) ∗ RXX (τ ), pour des signaux aléa-

toires, et

• G(f ) = rect[f /(2B)] exp(−j2πf t0),

avec |G(f )| = rect[f /(2B)] et φ(G(f )) = −2πf t0,

• ΓY X (τ ) = g(τ )∗ΓXX (τ ) si les signaux sont ergodiques,

ou certains à puissance moyenne finie ou à énergie finie

(attention aux définitions de Γ).

• g(t) = 2Bsinc[2B(t − t0)],

• Γgg(τ ) = 2Bsinc[2Bτ ] (auto-corrélation).

8.4. Statistiques en sortie d’un opérateur de filtrage

On considère y(t) = (g ∗ x)(t). Les statistiques de y(t) sup-

posé stationnaire et ergodique sont :

• µY = µX G(0) (moyenne),

• PY = RY Y (0)

=

=

RXX (τ )Γgg(τ )dτ

SXX (f )|G(f )|2df (moment d’ordre 2),

R

R

Y = CY Y (0) = PY − µ2

• σ2

X

=

CXX (τ )Γgg(τ )dτ (variance)

R

• SY X (f ) = G(f )SXX (f ) (auto-corrélation).

8.4.1. Filtre passe-bas sans perte

G(0) = 1 et la moyenne en sortie : µY = µX .

8.4.2. Filtre passe-haut

G(0) = 0 et la moyenne en sortie µY = 0.

8.4.3. Signal d’entrée blanc, centré et de DSP η/2

moyenne en sortie : µY = 0,

auto-corrélation de l’entrée : ΓXX (τ ) = δ(τ )η/2,

auto-corrélation en sortie : RY Y (τ ) = Γgg(τ )η/2,

variance : σ2

Y = Γgg(0)η/2.

9.4. Opérateur de filtrage passe-bande idéal

Le filtre causal de largeur de bande B, centré sur les fréquences

±f0 n’étant pas réalisable, on décale g(t) de t0.

• G(f ) =

rect[f /B] exp(−j2πf t0) ∗ [δ(f − f0) + δ(f + f0)],

avec |G(f )| = rect[f /(2B)] et φ(G(f )) = −2πf t0,

• g(t) = 2Bsinc[2B(t − t0)] cos(2πf0t),

• Γgg(τ ) = 2Bsinc[Bτ ] cos(2πf0τ ) (auto-corrélation).

10. FILTRE ADAPTÉ

Problème : concevoir un filtre h(t) qui détecte un signal x(t)

connu dans une observation bruitée, x(t) + n(t), où n(t) est

un bruit de DSP SN N (f ), avec le meilleur rapport signal/bruit

au temps t0. On note y(t) = (x + n) ∗ h(t) la sortie du filtre

h(t).

Le filtre optimal est le filtre adapté (au signal x(t)) :

H(f ) = k

X ∗(f )

SN N (f )

exp(−j2πf t0),

(1)

avec le rapport S/B :

S

B (cid:17)

(cid:16)

+

2

=

Z

|X(f )|2

SN N (f )

df.

(2)

−∞

Dans le cas où le bruit n(t) est un bruit blanc de DSP égale à

N0 :

9. AUTRES OPÉRATEURS

ou, dans le domaine temporel :

H(f ) = k′X ∗(f ) exp(−j2πf t0),

9.1. Opérateur de retard t0

• g(t) = δ(t − t0),

• G(f ) = exp(−j2πf t0), soit |G(f )| = 1

et φ(G(f )) = −2πf t0 (phase linéaire),

h(t) = k′x∗(t0 − t),

avec le rapport S/B :

2

S

B (cid:17)

(cid:16)

=

1

N0 Z

+

−∞

|X(f )|2df.

(3)

(4)

(5)

• Γgg(τ ) = δ(τ ) (auto-corrélation).

Pour un signal x(t) réel, on a simplement h(t) = k′x(t0 − t).