Communication Num´erique

Filtrage des signaux, Communication numérique · course

Voir tous les documents en électronique et automatique

Communication Num´erique

Communication Num´erique

Filtrage des signaux

Yoann Morel

http://xymaths.free.fr/Signal/Communication-Numerique-cours-TP.php

Communication Num´erique

1 G´en´eralit´es sur les filtres

2 M´ethode de Fourier

3 M´ethode s´equentielle : Filtres RIF et RII

4 Filtrage adapt´e

5 Effet du bruit et probabilit´e d’erreur

6 Application : Crit`ere et Filtre de Nyquist

Interf´erence entre symboles (IES)

Cas id´eal - Filtre de Nyquist

Filtre en cosinus sur´elev´e

Filtre en racine de cosinus sur´elev´e

Filtre en cosinus sur´elev´e avec coefficient de retomb´ee

Communication Num´erique

G´en´eralit´es sur les filtres

1 G´en´eralit´es sur les filtres

2 M´ethode de Fourier

3 M´ethode s´equentielle : Filtres RIF et RII

4 Filtrage adapt´e

5 Effet du bruit et probabilit´e d’erreur

6 Application : Crit`ere et Filtre de Nyquist

Interf´erence entre symboles (IES)

Cas id´eal - Filtre de Nyquist

Filtre en cosinus sur´elev´e

Filtre en racine de cosinus sur´elev´e

Filtre en cosinus sur´elev´e avec coefficient de retomb´ee

Communication Num´erique

G´en´eralit´es sur les filtres

Un canal de transmission r´eel a une bande passante limit´ee :

physiquement (bande passante r´eelle)

par l’utilisateur (limitation de d´ebit, partage entre utilisateurs,

probl`emes de CEM, . . .)

Or la DSP d’un signal r´eel, donc limit´e dans le temps, n’est jamais limit´e.

=⇒ N´ecessit´e de la limiter : recours au filtrage

y = h ∗ x

ˆy = H × ˆx

o`u, H = ˆh est la fonction de transfert du filtre.

Un filtre est caract´eris´e par sa fonction de transfert

(cid:44)→ Modification du contenu spectral du signal

Communication Num´erique

G´en´eralit´es sur les filtres

Rappel sur les filtres

x

(cid:45)

A

R.I. : h

(cid:45)y

Communication Num´erique

G´en´eralit´es sur les filtres

Rappel sur les filtres

x

(cid:45)

A

R.I. : h

(cid:45)y

y = h ∗ x

ˆy = H × ˆx

o`u, H = ˆh est la fonction de transfert du filtre.

Un filtre est caract´eris´e par sa fonction de transfert

(cid:44)→ Modification du contenu spectral du signal

Communication Num´erique

G´en´eralit´es sur les filtres

Principaux filtres

• Filtre passe-bas

|H(f )|

(cid:54)

G

(cid:64)

(cid:64)

(cid:64)

fc

(cid:64)

(cid:64)

(cid:45)

log f

(cid:64)

(cid:64)

(cid:64)

(cid:64)

Communication Num´erique

G´en´eralit´es sur les filtres

Principaux filtres

• Filtre passe-bande

|H(f )|

(cid:54)

G

(cid:0)

(cid:0)

(cid:0)

(cid:0)

(cid:0)

(cid:0)

(cid:0)

(cid:0)

(cid:0)

(cid:64)

(cid:64)

(cid:64)

(cid:64)

(cid:64)

f1

(cid:27)

f2

(cid:45)

∆f

(cid:45)

log f

(cid:64)

(cid:64)

(cid:64)

(cid:64)

Communication Num´erique

G´en´eralit´es sur les filtres

Principaux filtres

• Filtre passe-haut

|H(f )|

(cid:54)

G

(cid:0)

(cid:0)

(cid:0)

(cid:0)

(cid:0)

(cid:0)

(cid:0)

(cid:0)

(cid:0)

fc

(cid:45)

log f

Communication Num´erique

G´en´eralit´es sur les filtres

Synth`ese de filtres

• M´ethode analogique : ´electronique, composants actifs et

passifs

• M´ethode num´erique :

m´ethode de Fourier : ´etude de H(f )

m´ethode s´equentielle : relation xk ↔ yk

Communication Num´erique

M´ethode de Fourier

1 G´en´eralit´es sur les filtres

2 M´ethode de Fourier

3 M´ethode s´equentielle : Filtres RIF et RII

4 Filtrage adapt´e

5 Effet du bruit et probabilit´e d’erreur

6 Application : Crit`ere et Filtre de Nyquist

Interf´erence entre symboles (IES)

Cas id´eal - Filtre de Nyquist

Filtre en cosinus sur´elev´e

Filtre en racine de cosinus sur´elev´e

Filtre en cosinus sur´elev´e avec coefficient de retomb´ee

h(t) = 2aB sinc (2B(t − t0))

Communication Num´erique

M´ethode de Fourier

Filtre passe-bas id´eal

|H(f )|

(cid:54)

−B

B

H(f ) = a Rect2B(f ) e−2iπf t0

(cid:45)

log f

D’o`u la R.I. :

Communication Num´erique

M´ethode de Fourier

Filtre passe-bas id´eal

|H(f )|

(cid:54)

(cid:45)

log f

D’o`u la R.I. :

−B

B

H(f ) = a Rect2B(f ) e−2iπf t0

h(t) = 2aB sinc (2B(t − t0))

R.I. filtre passe-bas

fr´equence porteuse

Communication Num´erique

M´ethode de Fourier

Filtre passe-bande id´eal

|H(f )|

(cid:54)

(cid:27)

(cid:45)

−f0

2B

(cid:45)

log f

(cid:27)

(cid:45)

f0

2B

H(f ) = a [Rect2B(f − f0) + Rect2B(f + f0)] e−2iπf t0

d’o`u la R.I.

h(t) = 4aB sinc (2B(t − t0))

(cid:125)

(cid:123)(cid:122)

(cid:124)

cos(2πf0(t − t0))

(cid:125)

(cid:123)(cid:122)

(cid:124)

Communication Num´erique

M´ethode de Fourier

Filtre passe-bande id´eal

|H(f )|

(cid:54)

(cid:27)

(cid:45)

−f0

2B

(cid:45)

log f

(cid:27)

(cid:45)

f0

2B

H(f ) = a [Rect2B(f − f0) + Rect2B(f + f0)] e−2iπf t0

d’o`u la R.I.

h(t) = 4aB sinc (2B(t − t0))

(cid:125)

(cid:123)(cid:122)

R.I. filtre passe-bas

(cid:124)

cos(2πf0(t − t0))

(cid:124)

(cid:125)

(cid:123)(cid:122)

fr´equence porteuse

Communication Num´erique

M´ethode de Fourier

Exemple de filtrage

D´ebruitage d’un signal

Un bruit blanc additif gaussien (BBAG) est un signal additif

“haute fr´equence”.

=⇒ Filtrage par un filtre passe-bas pour d´ebruiter un signal

(cid:54) |H(f )|

−B

B

(cid:45)

log f

Communication Num´erique

M´ethode de Fourier

DD

Communication Num´erique

M´ethode de Fourier

Propri´et´e : Un filtre est r´ealisable (causal) si et seulement si

Supp h est born´e, i.e. si et seulement si la r´eponse

impulsionnelle h(t) est nulle en dehors d’un intervalle

de IR.

Probl`eme : Ces filtres id´eaux ne sont pas r´ealisables

(cid:54)|H(f )|

−B

B

(cid:45)

log f

par cette troncature, on introduit des oscillations

(Gibbs)

Communication Num´erique

M´ethode de Fourier

On multiplie h(t) par une fenˆetre W telle que

h(t) W (t) = 0 pour t /∈ [t1; t2]

o`u W est une fenˆetre

. . .

Probl`eme :

par cette troncature, on introduit des oscillations

(Gibbs)

Communication Num´erique

M´ethode de Fourier

On multiplie h(t) par une fenˆetre W telle que

h(t) W (t) = 0 pour t /∈ [t1; t2]

o`u W est une fenˆetre

rectangulaire

. . .

Probl`eme :

par cette troncature, on introduit des oscillations

(Gibbs)

Communication Num´erique

M´ethode de Fourier

On multiplie h(t) par une fenˆetre W telle que

h(t) W (t) = 0 pour t /∈ [t1; t2]

o`u W est une fenˆetre

triangulaire

Publicité

. . .

Probl`eme :

par cette troncature, on introduit des oscillations

(Gibbs)

Communication Num´erique

M´ethode de Fourier

On multiplie h(t) par une fenˆetre W telle que

h(t) W (t) = 0 pour t /∈ [t1; t2]

o`u W est une fenˆetre

Hanning (ou plutˆot Hann . . .)

. . .

Probl`eme :

par cette troncature, on introduit des oscillations

(Gibbs)

Communication Num´erique

M´ethode de Fourier

On multiplie h(t) par une fenˆetre W telle que

h(t) W (t) = 0 pour t /∈ [t1; t2]

o`u W est une fenˆetre

Hamming

. . .

Probl`eme :

par cette troncature, on introduit des oscillations

(Gibbs)

Communication Num´erique

M´ethode de Fourier

On multiplie h(t) par une fenˆetre W telle que

h(t) W (t) = 0 pour t /∈ [t1; t2]

o`u W est une fenˆetre

Blackman

. . .

Probl`eme :

par cette troncature, on introduit des oscillations

(Gibbs)

Communication Num´erique

M´ethode de Fourier

On multiplie h(t) par une fenˆetre W telle que

h(t) W (t) = 0 pour t /∈ [t1; t2]

o`u W est une fenˆetre

cosinus sur´elev´e

. . .

Probl`eme :

Communication Num´erique

M´ethode de Fourier

On multiplie h(t) par une fenˆetre W telle que

h(t) W (t) = 0 pour t /∈ [t1; t2]

o`u W est une fenˆetre

. . .

Probl`eme : par cette troncature, on introduit des oscillations

(Gibbs)

Communication Num´erique

M´ethode de Fourier

=⇒ Introduction ondulations :

Gain du filtre n’est plus constant, mais d´epend (l´eg`erement)

de la fr´equence

diminution de la raideur du filtre

Communication Num´erique

M´ethode s´equentielle : Filtres RIF et RII

1 G´en´eralit´es sur les filtres

2 M´ethode de Fourier

3 M´ethode s´equentielle : Filtres RIF et RII

4 Filtrage adapt´e

5 Effet du bruit et probabilit´e d’erreur

6 Application : Crit`ere et Filtre de Nyquist

Interf´erence entre symboles (IES)

Cas id´eal - Filtre de Nyquist

Filtre en cosinus sur´elev´e

Filtre en racine de cosinus sur´elev´e

Filtre en cosinus sur´elev´e avec coefficient de retomb´ee

Communication Num´erique

M´ethode s´equentielle : Filtres RIF et RII

Filtres num´eriques ←→ Signaux ´echantillon´es

(xk)

(cid:45)

A

R.I. : h

(cid:45)(yk)

La relation g´en´erale s’´ecrit alors (convolution discr`ete) :

yk =

l=+∞

(cid:88)

l=−∞

hl xk−l

ou, (hn) est la r´eponse impulsionnelle (discrete) du filtre.

Communication Num´erique

M´ethode s´equentielle : Filtres RIF et RII

Transformation discr`ete de Fourier :

(cid:98)yk =

(cid:88)

n

yk e−2iπf n

(cid:88)

=

n

(cid:88)

=

hl

l

(cid:88)

=

hl

(cid:32)

(cid:88)

(cid:33)

hl xn−l

e−2iπf n

l

(cid:32)

(cid:88)

n

(cid:32)

(cid:88)

(cid:33)

xn−l e−2iπf n

(cid:33)

e−2iπf (n(cid:48)+l)

l

n(cid:48)

(cid:88)

=

hl e−2iπf l

l

(cid:33)

e−2iπf n(cid:48)

(cid:32)

(cid:88)

n(cid:48)

Communication Num´erique

M´ethode s´equentielle : Filtres RIF et RII

On retrouve ainsi la relation,

(cid:98)yk = (cid:99)hk (cid:99)xk

Filtre `a r´eponse impulsionnelle finie : RIF (FIR)

Ce sont des filtres tels que Supp h born´e,

c’est-`a-dire hn = 0 si n /∈ [n1; n2]

On a alors,

et d’autre part,

(cid:98)h(f ) =

n=n2(cid:88)

n=n1

hn e−2iπf n

yk =

l=n2(cid:88)

l=n1

hl xk−l

Communication Num´erique

M´ethode s´equentielle : Filtres RIF et RII

Sch´ema d’un filtre RIF

{xi, xi+1, . . . , xk}

ak

Tk

ak−1

Tk−1

(cid:80) yk

· · ·

· · ·

...

hn d´efinit ainsi un filtre

passe-bas.

Communication Num´erique

M´ethode s´equentielle : Filtres RIF et RII

Exemple simple de filtre num´erique

h tel que h0 = h1 =

1

2

C’est-`a-dire,

, et hn = 0 si n (cid:54)= 0 et n (cid:54)= 1.

(cid:88)

yk =

hl xk−l =

(xk + xk−1)

1

2

l

yk est une moyenne locale du signal xk d’entr´e.

On a alors,

(cid:98)h(f ) =

(cid:88)

n

hn e−2iπf n =

1

2

+

1

2

e−2iπf = e−iπf cos(πf )

et donc,

|(cid:98)h(f )|2 = cos2(πf ) ∼ 1 , f → 0

∼ 0 , f → ±π

Communication Num´erique

M´ethode s´equentielle : Filtres RIF et RII

Exemple simple de filtre num´erique

h tel que h0 = h1 =

1

2

C’est-`a-dire,

, et hn = 0 si n (cid:54)= 0 et n (cid:54)= 1.

(cid:88)

yk =

hl xk−l =

(xk + xk−1)

1

2

l

yk est une moyenne locale du signal xk d’entr´e.

On a alors,

(cid:98)h(f ) =

(cid:88)

n

hn e−2iπf n =

1

2

+

1

2

e−2iπf = e−iπf cos(πf )

et donc,

|(cid:98)h(f )|2 = cos2(πf ) ∼ 1 , f → 0

∼ 0 , f → ±π

hn d´efinit ainsi un filtre

passe-bas.

Communication Num´erique

M´ethode s´equentielle : Filtres RIF et RII

Les filtres FIR ne sont en g´en´eral pas suffisants (trops restrictifs).

On d´efinit alors les filtres `a R´eponse Impulsionnelle Infinie (RII, ou

IIR), pour lesquels on a :

(cid:98)h(f ) =

n=+∞

(cid:88)

n=−∞

hne−2iπf n , hn (cid:54)= 0

En pratique, il est impossible d’impl´ementer des suites de longueur

infinie.

On introduit alors ces filtres de mani`ere r´ecursive :

yk =

M

(cid:88)

n=0

βn xk−n −

N

(cid:88)

n=1

αn yk−n

yk est calcul´e `a partir de l’entr´ee xk et des valeurs ant´erieures de

yk : c’est un calcul it´eratif

Communication Num´erique

M´ethode s´equentielle : Filtres RIF et RII

Sch´ema d’un filtre RII

{xi, . . . , xk}

ak

Tk

bk

Tk

ak−1

Tk−1

bk−1

· · · · · ·

(cid:80)

yk

· · ·

· · ·

...

{yi, · · · , yk}

Communication Num´erique

M´ethode s´equentielle : Filtres RIF et RII

En r´eordonnant, on obtient, en posant α0 = 1,

α0 yk =

M

(cid:88)

n=0

βn xk−n −

N

(cid:88)

n=1

αn yk−n ⇐⇒

N

(cid:88)

n=0

αn yk−n =

M

(cid:88)

n=0

βn xk−n

soit, en prenant la transform´ee de Fourier,

(cid:32) N

(cid:88)

n=0

(cid:33)

Publicité

αn e−2iπf n

(cid:98)y(f ) =

(cid:32) M

(cid:88)

n=0

(cid:33)

βn e−2iπf n

(cid:98)x(f )

d’o`u, la fonction de transfert,

H(f ) = (cid:98)h(f ) = (cid:98)y(f )

(cid:98)x(f )

=

M

(cid:88)

n=0

N

(cid:88)

n=0

βn e−2iπf n

αn e−2iπf n

=⇒ Les fr´equences de coupure sont les racines de P .

=⇒ Probl`eme de stabilit´e du filtre, si le d´enominateur s’an-

nule (racines de Q)

=⇒ En g´en´eral, existence d’un domaine de stabilit´e qu’il faut

respecter.

Communication Num´erique

M´ethode s´equentielle : Filtres RIF et RII

En utilisant une transform´ee en z (z = e−2iπf ),

H(f ) = (cid:101)H(z) =

M

(cid:88)

n=0

N

(cid:88)

n=0

βn zn

αn zn

=

βM zM + · · · + β1z + β0

αN zN + · · · + α1z + α0

=

P (z)

Q(z)

Communication Num´erique

M´ethode s´equentielle : Filtres RIF et RII

En utilisant une transform´ee en z (z = e−2iπf ),

H(f ) = (cid:101)H(z) =

M

(cid:88)

n=0

N

(cid:88)

n=0

βn zn

αn zn

=

βM zM + · · · + β1z + β0

αN zN + · · · + α1z + α0

=

P (z)

Q(z)

=⇒ Les fr´equences de coupure sont les racines de P .

=⇒ Probl`eme de stabilit´e du filtre, si le d´enominateur s’an-

nule (racines de Q)

=⇒ En g´en´eral, existence d’un domaine de stabilit´e qu’il faut

respecter.

Communication Num´erique

Filtrage adapt´e

1 G´en´eralit´es sur les filtres

2 M´ethode de Fourier

3 M´ethode s´equentielle : Filtres RIF et RII

4 Filtrage adapt´e

5 Effet du bruit et probabilit´e d’erreur

6 Application : Crit`ere et Filtre de Nyquist

Interf´erence entre symboles (IES)

Cas id´eal - Filtre de Nyquist

Filtre en cosinus sur´elev´e

Filtre en racine de cosinus sur´elev´e

Filtre en cosinus sur´elev´e avec coefficient de retomb´ee

Communication Num´erique

Filtrage adapt´e

Soit un signal cod´e avec la forme d’onde g(t) (formant du code en

ligne), et bruit´e :

x(t) = g(t) + w(t) , t ∈ [0; Tb]

o`u w est un BBAG de densit´e spectrale W (f ) =

N0

2

On cherche le filtre en reception h(t) :

y(t) = h(t) ∗ x(t) = g0(t) + n(t)

tel que le RSB en t = T (instant de d´ecision) :

η =

|g0(T )|2

E(n2(t))

soit maximal

Communication Num´erique

Filtrage adapt´e

On a :

g0(T ) = h(t) ∗ g(t) =

(cid:90)

IR

H(f ) G(f ) e2iπf T df

E(n2(t)) =

(cid:90)

IR

|H(f )|2 |W (f )|2 df =

N0

2

(cid:90)

IR

|H(f )|2 df

et,

d’o`u,

(cid:12)

(cid:90)

(cid:12)

(cid:12)

(cid:12)

IR

η =

2

(cid:12)

(cid:12)

(cid:12)

(cid:12)

H(f ) G(f ) e2iπf T df

N0

2

(cid:90)

IR

|H(f )|2 df

Communication Num´erique

Filtrage adapt´e

Par Cauchy-Schwarz, on a :

2

N0

avec ´egalit´e si et seulement si,

η ≤

(cid:90)

IR

|G(f )|2 df

H(f ) = k G(f ) e−2iπf T , k ∈ C

soit, pour des filtres r´eels,

h(t) = g(T − t)

Le filtre h(t) = g(T − t) est le filtre optimal, ou filtre adapt´e,

`a la d´etection du signal g(t) bruit´e.

Communication Num´erique

Effet du bruit et probabilit´e d’erreur

1 G´en´eralit´es sur les filtres

2 M´ethode de Fourier

3 M´ethode s´equentielle : Filtres RIF et RII

4 Filtrage adapt´e

5 Effet du bruit et probabilit´e d’erreur

6 Application : Crit`ere et Filtre de Nyquist

Interf´erence entre symboles (IES)

Cas id´eal - Filtre de Nyquist

Filtre en cosinus sur´elev´e

Filtre en racine de cosinus sur´elev´e

Filtre en cosinus sur´elev´e avec coefficient de retomb´ee

Communication Num´erique

Effet du bruit et probabilit´e d’erreur

Effet du bruit et probabilit´e d’erreur

Le filtre adapt´e est le filtre optimal pour la d´etection du signal

utile.

La mesure de la performance du filtrage (adapt´e ou non) est

donn´ee par la probabilit´e d’erreur :

probabilit´e de d´ecider qu’un bit 1 `a ´etait re¸cu alors

qu’un bit 0 a ´et´e envoy´e, ou inversement.

Communication Num´erique

Effet du bruit et probabilit´e d’erreur

Exemple de calcul de la probabilit´e d’erreur

On consid`ere un signal polaire NRZ :

x(t) =

(cid:40) +A + w(t) , si bit ´emis = +1

−A + w(t) , si bit ´emis = +0

D’apres le th´eoreme de Bayes, on a

pe = p(0|1) p1 + p(1|0) p0

avec,

(cid:40) p0, p1 les probabilit´es d’avoir un 0 ou un 1

p(0|1), p(1|0) les probabilit´es conditionnelles

Communication Num´erique

Effet du bruit et probabilit´e d’erreur

Soit Y la sortie du filtre adapt´e :

Y = g(−t) ∗ x(t) =

1

Tb

(cid:90) Tb

0

x(t) dt

et donc, si un 0 a ´et´e envoy´e,

Y = −A +

1

Tb

(cid:90) Tb

0

w(t) dt

d’o`u la densit´e de probabilit´e conditionnelle,

TY |0 =

1

2πσ2

(cid:18)

exp

(cid:19)

(y + A)2

2σ2

avec la variance σ du bruit w(t) gaussien : σ2 =

N0

2Tb

Communication Num´erique

Effet du bruit et probabilit´e d’erreur

On obtient alors, avec un seuil de d´ecision λ,

p(1|0) = pY |0(Y ≥ λ) =

(cid:90) +∞

λ

T (Y |0) dz =

(cid:32)

1

2

erfc

A + λ

(cid:112)N0/Tb

(cid:33)

avec la fonction d’erreur compl´ementaire :

erfc(u) =

2

π

On trouve de mˆeme,

(cid:90) +∞

e−s2

ds

u

(cid:32)

(cid:33)

A − λ

(cid:112)N0/Tb

p(0|1) =

1

2

erfc

Communication Num´erique

Effet du bruit et probabilit´e d’erreur

Au total, on obtient,

(cid:32)

p1

2

erfc

Dans le cas o`u, p0 = p1 =

alors,

, alors pe est minimum pour λ = 0 et

(cid:33)

+

p0

2

erfc

(cid:33)

(cid:32)

A + λ

(cid:112)N0/Tb

A − λ

(cid:112)N0/Tb

1

2

pe =

1

2

erfc

(cid:33)

(cid:32)

A

(cid:112)N0/Tb

et donc, si on pose εb = A2 Tb, l’´energie par bit (impulsion

d’amplitude A de dur´ee Tb), alors,

pe =

1

2

erfc

(cid:19)

(cid:18)(cid:114) εb

N0

Communication Num´erique

Effet du bruit et probabilit´e d’erreur

pe =

1

2

erfc

Publicité

(cid:19)

(cid:18)(cid:114) εb

N0

(cid:114) εb

N0

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

1 G´en´eralit´es sur les filtres

2 M´ethode de Fourier

3 M´ethode s´equentielle : Filtres RIF et RII

4 Filtrage adapt´e

5 Effet du bruit et probabilit´e d’erreur

6 Application : Crit`ere et Filtre de Nyquist

Interf´erence entre symboles (IES)

Cas id´eal - Filtre de Nyquist

Filtre en cosinus sur´elev´e

Filtre en racine de cosinus sur´elev´e

Filtre en cosinus sur´elev´e avec coefficient de retomb´ee

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Interf´erence entre symboles (IES)

1 G´en´eralit´es sur les filtres

2 M´ethode de Fourier

3 M´ethode s´equentielle : Filtres RIF et RII

4 Filtrage adapt´e

5 Effet du bruit et probabilit´e d’erreur

6 Application : Crit`ere et Filtre de Nyquist

Interf´erence entre symboles (IES)

Cas id´eal - Filtre de Nyquist

Filtre en cosinus sur´elev´e

Filtre en racine de cosinus sur´elev´e

Filtre en cosinus sur´elev´e avec coefficient de retomb´ee

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Interf´erence entre symboles (IES)

Interf´erence entre symboles (IES)

Interf´erence entre symboles :

L’interf´erence entre symboles est un ph´enom`ene qui se

produit si le niveau ´echantillonn´e `a l’instant de d´ecision

ne d´epend pas du seul symbole attendu, mais se trouve

alt´er´e par la superposition d’un ou plusieurs autres sym-

boles voisins.

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Interf´erence entre symboles (IES)

Interf´erence entre symboles (IES)

R´eponse globale du syst`eme, avec filtre adapt´e en r´eception :

(an)

(cid:45)

h(t)

Bruit

(cid:63)

(cid:45)(cid:77)

(cid:45) h(−t)

(cid:45)y(t)

Codage en ligne

Filtrage adapt´e

On a alors,

y(t) =

(cid:88)

n

anr(t − nT ) + bruit

avec la r´eponse globale du syst`eme :

r(t) = h(t) ∗ h(−t)

Symbole

`a d´etecter

IES

bruit

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Interf´erence entre symboles (IES)

Interf´erence entre symboles (IES)

R´eponse en sortie du filtre adapt´e :

y(t) =

(cid:88)

n

anr(t − nT ) + bruit

Echantillonnage aux instants t = kT :

y(kT ) =

(cid:88)

n

anr(kT − nT ) + bk

Symbole

`a d´etecter

IES

bruit

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Interf´erence entre symboles (IES)

Interf´erence entre symboles (IES)

R´eponse en sortie du filtre adapt´e :

y(t) =

(cid:88)

n

anr(t − nT ) + bruit

Echantillonnage aux instants t = kT :

y(kT ) =

(cid:88)

n

anr(kT − nT ) + bk

= ak r(0) +

(cid:88)

n(cid:54)=k

anr ((k − n)T ) + bk

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Interf´erence entre symboles (IES)

Interf´erence entre symboles (IES)

R´eponse en sortie du filtre adapt´e :

y(t) =

(cid:88)

n

anr(t − nT ) + bruit

Echantillonnage aux instants t = kT :

y(kT ) =

(cid:88)

n

anr(kT − nT ) + bk

= ak r(0)

(cid:124) (cid:123)(cid:122) (cid:125)

+

Symbole

`a d´etecter

(cid:88)

n(cid:54)=k

(cid:124)

anr ((k − n)T )

+ bk

(cid:124)(cid:123)(cid:122)(cid:125)

(cid:123)(cid:122)

(cid:125)

IES

bruit

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Interf´erence entre symboles (IES)

Interf´erence entre symboles (IES)

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Interf´erence entre symboles (IES)

Interf´erence entre symboles (IES)

On appelle distortion maximale, le ratio :

|anr(nT )|

|r(0)|

|r(nT )|

(cid:88)

n(cid:54)=0

(cid:88)

n(cid:54)=0

|r(0)|

,

si an ∈ {−1, +1}

Dmax =

=





Dmax ≥ 1 =⇒ IES provoque une erreur de d´etection

Dmax < 1 =⇒ Condition n´ecessaire pour la bonne

r´eception du message

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Interf´erence entre symboles (IES)

Diagramme de l’œil

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Interf´erence entre symboles (IES)

En l’abscence d’IES, l’œil est completement “ouvert” a l’instant de

d´ecision ; tous les trajets passsent par deux points seulement (en

binaire, M points en M-aires).

(a) Ouverture verticale : immunit´e au bruit

(b) Ouverture horizontale : immununit´e au d´ephasage d’horloge

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Cas id´eal - Filtre de Nyquist

1 G´en´eralit´es sur les filtres

2 M´ethode de Fourier

3 M´ethode s´equentielle : Filtres RIF et RII

4 Filtrage adapt´e

5 Effet du bruit et probabilit´e d’erreur

6 Application : Crit`ere et Filtre de Nyquist

Interf´erence entre symboles (IES)

Cas id´eal - Filtre de Nyquist

Filtre en cosinus sur´elev´e

Filtre en racine de cosinus sur´elev´e

Filtre en cosinus sur´elev´e avec coefficient de retomb´ee

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Cas id´eal - Filtre de Nyquist

Cas id´eal - Filtre de Nyquist

C’est le cas Dmax = 0.

Cela revient `a :

(cid:40) r(kT ) = 0 , ∀k (cid:54)= 0

r(0) (cid:54)= 0

Dans ce cas, le filtre, de R.I. r(t), est dit de Nyquist.

Remarque : La condition de Nyquist est trivialement satisfaite

lorsque l’impulsion h(t) est de dur´ee plus petite

que T , et donc |Supp(h)| ≤ 2T

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Cas id´eal - Filtre de Nyquist

Crit`ere de Nyquist fr´equentiel

Le crit`ere de Nyquist implique que :

(cid:88)

n

r(t − nT ) = r(t)

(cid:88)

n

δ(t − nT ) = r(0)δt

soit, dans le domaine fr´equentiel :

(cid:98)r(f ) ∗

1

T

(cid:88)

(cid:16)

δ

f −

n

(cid:17)

n

T

= r(0)

soit,

(cid:16)

(cid:88)

n

T

c’est-`a-dire, (cid:98)r doit ˆetre constant sur une largeur de bande B ≥ 1

2T .

On appelle bande de Nyquist la largeur de bande minimale :

= T r(0)

f −

(cid:98)r

n

(cid:17)

Bmin =

1

2T

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Cas id´eal - Filtre de Nyquist

Choix le plus simple :

R(f ) = T r(0) Rect1/T (f ) =

(cid:26) 1 , si |f | ≤ 1

2T

0 , sinon

On a alors,

r(t) = r(0) sinc

(cid:19)

(cid:18) t

T

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Cas id´eal - Filtre de Nyquist

Inconv´enients des filtres de Nyquist

A l’instant d’´echantillonnage exact t = kT , on a

r(kT ) = sinc(kT ) = 0 , mboxet r(0) = 0

c’est-`a-dire, Dmax = 0.

Mais, si l’instant d’´echantillonnage est nettement moins pr´ecis `a la

r´eception : t = kT + ε, avec ε (cid:54)= 0, alors,

(cid:16)

(cid:12)

(cid:12)

(cid:12)sinc

kπ +

πε

T

(cid:17)(cid:12)

(cid:12)

(cid:12) = +∞

(cid:88)

k

=⇒ Aucune tol´erance vis-`a-vis de l’instant d’´echantillonnage

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Filtre en cosinus sur´elev´e

1 G´en´eralit´es sur les filtres

2 M´ethode de Fourier

3 M´ethode s´equentielle : Filtres RIF et RII

4 Filtrage adapt´e

5 Effet du bruit et probabilit´e d’erreur

6 Application : Crit`ere et Filtre de Nyquist

Interf´erence entre symboles (IES)

Cas id´eal - Filtre de Nyquist

Filtre en cosinus sur´elev´e

Filtre en racine de cosinus sur´elev´e

Filtre en cosinus sur´elev´e avec coefficient de retomb´ee

Publicité

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Filtre en cosinus sur´elev´e

Filtre en cosinus sur´elev´e

Si ak+1 = −ak, alors,

(k + 1)T ,

le cas optimal ´etant :

(cid:88)

n

(cid:88)

n

an r(t − nT ) doit s’annuler entre kT et

anr(t − nT ) = 0 pour t = (cid:0)k + 1

2

(cid:1) T

1`ere solution : r (cid:0)(cid:0)k + 1

2

(cid:1) T (cid:1) = 0 =⇒ r(t) = sinc

(cid:19)

(cid:18) 2πt

T

2`eme solution : On impose

r(t) =





r(0)

2

0

,

si t = ±

T

2

sit = ±T, ±

3

2

T, ±2T, . . . , ±

kT

2

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Filtre en cosinus sur´elev´e

De cette fa¸con, on a, en t = (cid:0)k + 1

2

(cid:1) T





si ak+1 = −ak ,

si ak+1 = ak , ±

r(0)

2

(cid:18) r(0)

2

r(0)

2

r(0)

2

+

= 0

(cid:19)

= ±r(0)

La fonction

(cid:20)

r(t) =

sinc

(cid:19)

(cid:18) 2πt

T

+

1

2

sinc

(cid:18) 2π

T

(cid:18)

t +

T

2

(cid:19)(cid:19)

+

1

2

sinc

(cid:18) 2π

T

(cid:18)

t −

T

2

(cid:19)(cid:19)(cid:21)

v´erifie ce crit`ere, et on a :

r(t) = r(0)

sinc (cid:0) 2πt

T

1 − (cid:0) 2t

T

(cid:1)

(cid:1)2 = r(0) sinc

(cid:19)

(cid:18) πt

T

cosc

(cid:19)

(cid:18) πt

T

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Filtre en cosinus sur´elev´e

avec la fonction cosinus cardinal : cosc(x) =

=⇒ On a alors une d´ecroissance en

1

n2 (et

cos x

1 − x2

1

n2 < +∞)

(cid:88)

n

Dans le domaine fr´equentiel, fonction de trnasfert associ´ee est :

R(f ) =

T r(0)

2

[1 + cos(πf T )]

,

|f | ≤

1

T

“1 + cos” : Filtre en cosinus sur´elev´e

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Filtre en racine de cosinus sur´elev´e

1 G´en´eralit´es sur les filtres

2 M´ethode de Fourier

3 M´ethode s´equentielle : Filtres RIF et RII

4 Filtrage adapt´e

5 Effet du bruit et probabilit´e d’erreur

6 Application : Crit`ere et Filtre de Nyquist

Interf´erence entre symboles (IES)

Cas id´eal - Filtre de Nyquist

Filtre en cosinus sur´elev´e

Filtre en racine de cosinus sur´elev´e

Filtre en cosinus sur´elev´e avec coefficient de retomb´ee

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Filtre en racine de cosinus sur´elev´e

Filtre en racine de cosinus sur´elev´e

Probleme : r(t) ainsi calcul´e est la R.I. globale du systeme,

r(t) = h(t) ∗ h(−t)

−→ D´etermination de h(t) ?

R(f ) =

T r(0)

2

[1 + cos(πf T )] = T r(0) cos2

(cid:19)

(cid:18) πf T

2

=⇒ H(f ) = (cid:112)T r(0) cos

(cid:19)

(cid:18) πf T

2

tel que R(f ) = |H(f )|2. Alors,

H(f ) = (cid:112)T r(0)

ei πf T

2 + ei πf T

2

2

Rect2/T (f )

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Filtre en racine de cosinus sur´elev´e

La transform´ee de Fourier inverse donne alors la R.I.

h(t) =

(cid:112)T r(0)

2

(cid:20)

sinc

(cid:18) 2t

T

(cid:19)

1

2

+ sinc

(cid:19)(cid:21)

(cid:18) 2t

T

+

1

2

soit aussi,

h(t) =

2

π

(cid:112)T r(0) cosc

(cid:19)

(cid:18) 2t

T

−→ Probleme : h(t) ne v´erifie pas le critere de Nyquist...

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Filtre en cosinus sur´elev´e avec coefficient de retomb´ee

1 G´en´eralit´es sur les filtres

2 M´ethode de Fourier

3 M´ethode s´equentielle : Filtres RIF et RII

4 Filtrage adapt´e

5 Effet du bruit et probabilit´e d’erreur

6 Application : Crit`ere et Filtre de Nyquist

Interf´erence entre symboles (IES)

Cas id´eal - Filtre de Nyquist

Filtre en cosinus sur´elev´e

Filtre en racine de cosinus sur´elev´e

Filtre en cosinus sur´elev´e avec coefficient de retomb´ee

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Filtre en cosinus sur´elev´e avec coefficient de retomb´ee

Filtre de Nyquist : R0(f )

Filtre en cosinus sur´elev´e : R1(f )

En pratique, on utilise un compromis entre R0 et R,

Rα , avec 0 ≤ α ≤ 1 , et

Rα = R0 , pour α = 0

Rα = R1 , pour α = 1

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Filtre en cosinus sur´elev´e avec coefficient de retomb´ee

On remarque que pour α = 1, on a

r(t) = sinc

(cid:19)

(cid:18) πt

T

cosc

(cid:19)

(cid:18) πt

T

= rO(t)

=⇒ R(f ) = R0(f ) ∗ C(f )

avec, C(f ) = T.F.

cosc

(cid:18)

(cid:19)(cid:19)

(cid:18) πt

T

On remplace alors C(f ) par

1

α

C

Rα = R0(f ) ∗

πT

2

(cid:19)

,

=

(cid:18) f

α

1

α

(cid:19)

(cid:18) f

α

C

cos (πf T ) Rect1/T (f )

Communication Num´erique

Application : Crit`ere et Filtre de Nyquist

Filtre en cosinus sur´elev´e avec coefficient de retomb´ee

En temporel, on a alors,

rα(t) = r0(t) cosc

(cid:19)

(cid:18) αt

T

= sinc

(cid:19)

(cid:18) πt

T

cosc

(cid:19)

(cid:18) αt

T

et donc le compromis :

• Tol´erance au bruit α = 0

• Tol´erance `a l’instant d’´echantillonnage α = 1