04.01.2015 Views

Traitement du Signal Bruit de la Numérisation - PRIMA

Traitement du Signal Bruit de la Numérisation - PRIMA

Traitement du Signal Bruit de la Numérisation - PRIMA

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

<strong>Traitement</strong> <strong>du</strong> <strong>Signal</strong><br />

James L. Crowley<br />

Deuxième Année ENSIMAG première Semestre 2002/2003<br />

Séance 5 : 21 octobre 2002<br />

<strong>Bruit</strong> <strong>de</strong> <strong>la</strong> Numérisation<br />

Formule <strong>du</strong> Jour : .......................................................................1<br />

<strong>Bruit</strong> <strong>de</strong> Echantillonage et <strong>de</strong> Numerisation ..........................2<br />

<strong>Bruit</strong> d'Echantillonage........................................................... 3<br />

Quanitification:..................................................................... 5<br />

Métho<strong>de</strong>s <strong>de</strong> conversion A/N...................................................9<br />

Métho<strong>de</strong>s Indirect : VCO....................................................... 9<br />

Convertiseeur parallèle. (Ang<strong>la</strong>is “F<strong>la</strong>sh Converteur”).............. 9<br />

Approximation sucessive ....................................................... 10<br />

Conversion Numérique-analogique......................................... 11<br />

Formule <strong>du</strong> Jour :<br />

La qualité d'un signal est souvent représentée par le "rapport signal/bruit"<br />

(SNR en ang<strong>la</strong>is).<br />

pour x(t) = s(t) + b(t).<br />

Le rapport signal sur bruit est défini par<br />

ξ= W s<br />

W n<br />

où W s est l'énergie <strong>du</strong> signal s(t) et<br />

W b est l'énergie <strong>du</strong> bruit b(t)<br />

Le SNR est souvent représenté avec une échelle logarithmique appelée<br />

décibels et noté dB.<br />

ξ dB = 10 log 10 ξ<br />

Un facteur <strong>de</strong> 3 dB est équivalent à un facteur <strong>de</strong> 2 en puissance


<strong>Bruit</strong> <strong>de</strong> <strong>la</strong> numérisation Séance 5<br />

<strong>Bruit</strong> <strong>de</strong> Echantillonage et <strong>de</strong> Numerisation<br />

Soit<br />

echantillonage : s(n) = Echant{s(t)} (paramètre : T e )<br />

quantification : s q (n) = Quant{s(n∆T)} (paramètres : V, q)<br />

Conversions Analogique : s^(t) = Analog{s q (n)} (paramtres :T a ,V a q a }<br />

Cas Continu :<br />

s(t)<br />

s(nT e )<br />

s q (nT e )<br />

s(t) ^<br />

Echant{s(t)}<br />

Quant{s(n)}<br />

Analog{s q (n)} –<br />

T e V, q V a , q a<br />

e(t)<br />

En signal continu, le bruit <strong>de</strong> numérisation peut être calculé comme :<br />

e(t) = s(t) – s^(t)<br />

Son energie est W T = ∫<br />

0<br />

Te(t) 2 dt<br />

Pour le cas continue, l'etu<strong>de</strong> <strong>de</strong>s <strong>de</strong> <strong>la</strong> bruit est fait avec <strong>la</strong> theorie <strong>de</strong> <strong>la</strong> probabilité<br />

en utilisant distribution continu <strong>de</strong>s probabilité est les intégrales.<br />

5-2


<strong>Bruit</strong> <strong>de</strong> <strong>la</strong> numérisation Séance 5<br />

<strong>Bruit</strong> d'Echantillonage<br />

Soit f e = 1 T e<br />

= 1.<br />

X e (f) = X(f) * δ fe(f) =<br />

∞<br />

∑ X(f – kfe )<br />

k=–∞<br />

Le "bruit" d'echantillonage est<br />

W b 2<br />

= 2<br />

⌠<br />

⎮<br />

⎮<br />

⌡<br />

-f e /2<br />

f e /2<br />

∫ X e (f) 2 df<br />

-f e /2<br />

f e /2<br />

∞<br />

∑<br />

k=1<br />

X(f – kfe ) 2 df<br />

Mais ceci est équivalent a tout l'energie au <strong>de</strong><strong>la</strong> <strong>de</strong> f e<br />

2<br />

∞<br />

∞<br />

W b = 2 ∫ X e (f) 2 df = 2 ∫ X(f) 2 df<br />

f e /2<br />

f e /2<br />

Exemple : Soit x(t) avec Tranformée X(f) = 1 pour | f | ≤ F max<br />

Si l'echantillonage est fait à T e = 1 f e<br />

tout energie entre f e<br />

2 et F max <strong>de</strong>vient bruit.<br />

Après échantillonage, l'energie <strong>du</strong> signal : W s = 2 ∫<br />

0<br />

f e /2<br />

X(f) 2 df<br />

Energie <strong>du</strong> bruit<br />

∞<br />

: W b = 2 ∫ X e (f) 2 df= 2<br />

f e /2<br />

F max<br />

∫ X(f) 2 df<br />

f e /2<br />

5-3


<strong>Bruit</strong> <strong>de</strong> <strong>la</strong> numérisation Séance 5<br />

–F max<br />

X(f)<br />

–F e /2<br />

F e /2<br />

F max<br />

f<br />

Rapport ξ= W s<br />

W n<br />

=<br />

f e /2<br />

∫ X(f) 2 df<br />

0<br />

F max<br />

∫X(f) 2 df<br />

f e /2<br />

L'information d'un signal, x(t), est donc assimilé a sa variance,<br />

σ x<br />

2 = E{(x(t)–µ x ) 2 }<br />

Dans ce qui suit, nous suposons que a moyenne <strong>du</strong> signal, µ x = µ s = µ n = 0.<br />

Donc σ x<br />

2 = W x<br />

La variance d'un signal est une caractéristique liée à <strong>la</strong> quantité <strong>de</strong> l'information<br />

représentée.<br />

Le rapport signal sur bruit est défini par<br />

ξ q =<br />

W s<br />

W<br />

= σ s 2<br />

n σ n<br />

2<br />

où σ s 2 est <strong>la</strong> variance <strong>du</strong> signal s(t) et<br />

σ n<br />

2 est <strong>la</strong> variance <strong>du</strong> bruit n(t)<br />

5-4


<strong>Bruit</strong> <strong>de</strong> <strong>la</strong> numérisation Séance 5<br />

Quanitification:<br />

Soit un signal x e (nT e ) <strong>de</strong> p<strong>la</strong>ge utile (le range dynamique) V.<br />

Pour le suivant, nous <strong>la</strong>isserons T e = 1.<br />

Le signal est décomposée en 2 B = V q<br />

intervalles <strong>de</strong> <strong>la</strong>rgeur q, par<br />

x q (n) = Max<br />

k {(x e(nT e ) – k·q + q 2 < 0}<br />

3q<br />

2q<br />

x q<br />

q<br />

s<br />

Pour une quantifiaction uniforme, <strong>la</strong> <strong>de</strong>nsité <strong>de</strong> probabilité est uniforme p(e q ) = 1 q<br />

Pour les valeurs quantifiées par arrondi, µ q = 0.<br />

–q<br />

–2q<br />

–3q<br />

e q = s e – s q<br />

q<br />

2<br />

–q<br />

2<br />

s e<br />

1<br />

q<br />

–q<br />

2<br />

p(e q )<br />

q<br />

2<br />

e q<br />

µ = 0 La variance reste σ 2 =<br />

1<br />

12 q2 5-5


<strong>Bruit</strong> <strong>de</strong> <strong>la</strong> numérisation Séance 5<br />

Soit un signal s(n), L'effet <strong>de</strong> <strong>la</strong> quantification est mesuré par :<br />

q<br />

s(n)<br />

Re-Quant{s(n)}<br />

+<br />

+ e(n)<br />

–<br />

L'erreur est : e(n) = s(n) – s q (n)<br />

N–1<br />

Son energie est W N = ∑ e(n) 2 dt<br />

n=0<br />

Son valeur moyenne est Sa moyenne est<br />

q/2<br />

µ = E{p(e)} = ∫ p(e) · e <strong>de</strong> =<br />

-q/2<br />

q/2<br />

⌠<br />

= ⎮<br />

⌡<br />

-q/2<br />

1<br />

q . · e <strong>de</strong> = 1<br />

q . 1 2 e2 ⏐ q/2<br />

-q/2 = 1 q . 1 2 (q2 4 – q2<br />

4 ) q2 = 0<br />

Sa variance est σ 2 =<br />

1<br />

12 q2<br />

σ 2 = E{p(e)} =<br />

=<br />

=<br />

q/2<br />

⌠<br />

⎮<br />

⌡<br />

-q/2<br />

1<br />

q/2<br />

∫ p(e) · (e–µ) 2 <strong>de</strong>)<br />

-q/2<br />

1<br />

q . · (e – q 2 )2 <strong>de</strong>) =<br />

3q [( q 2 – q 2 )3 – ( – q 2 – q 2 )3 ]<br />

1<br />

q · 1 3 (e – q 2 )3 ⏐ q/2<br />

-q/2<br />

=<br />

1<br />

3q<br />

q3<br />

[(<br />

8 ) + ( q3<br />

8 ) ] =<br />

1<br />

3q (q3 4 ) =<br />

1<br />

12 q2 5-6


<strong>Bruit</strong> <strong>de</strong> <strong>la</strong> numérisation Séance 5<br />

Le rapport signal sur bruit <strong>de</strong> quantification est défini par<br />

ξ q = W s<br />

W n<br />

N 2s<br />

L'energie W s (N 1 , N 2 ) = ∑ 2(n).<br />

n=N1<br />

Si on suppose que E{s} = 0, <strong>la</strong> variance vaut l'energie moyenne.<br />

σ s<br />

2 =<br />

Si E{n q } = 0,<br />

1<br />

N 1 –N 2<br />

W s =<br />

1<br />

N 1 –N 2<br />

N 2s(n)<br />

∑ 2.<br />

n=N1<br />

σ q<br />

2 =<br />

1<br />

N 1 –N 2<br />

W q =<br />

Dans ce cas, ξ q = W s<br />

W n<br />

=<br />

1<br />

N 1 –N 2<br />

σ s<br />

2<br />

σ q<br />

2 .<br />

N 2eq<br />

∑ (n) 2 .<br />

n=N1<br />

Donc, le rapport <strong>Signal</strong>-<strong>Bruit</strong> peut être estimé par les <strong>de</strong>nsité <strong>de</strong> probabilité.<br />

Par <strong>la</strong> formule précè<strong>de</strong>nt :<br />

ξ q = σ s 2<br />

q 2<br />

12<br />

12 σ s 2<br />

q 2<br />

et en décibels :<br />

ξ qdB = 10 log 10 ξ q = 10 log 10 (12 σ x 2<br />

q 2 ) = 2 . 10 log 10 ( σ x<br />

q ) + 10 log 10(12) dB<br />

≈ 20 log 10<br />

σ x<br />

q<br />

+ 10.8 dB<br />

Si <strong>la</strong> p<strong>la</strong>ge utile V <strong>du</strong> signal est décomposée en 2 B intervalles <strong>de</strong> <strong>la</strong>rgeur q,<br />

2 B = V q ⇔ q = V 2 B ⇔ log 10{ σ x<br />

V/2 B} = 10 log 10{2 B } + 20 log 10 { σ x<br />

V } 5-7


<strong>Bruit</strong> <strong>de</strong> <strong>la</strong> numérisation Séance 5<br />

et donc <strong>la</strong> ré<strong>du</strong>ction en rapport <strong>Signal</strong>-bruit entrainée par <strong>la</strong> numérisation est<br />

ξ qdB =<br />

10.8 + B 10 log 10 (2) – 20 log 10<br />

V<br />

σ x<br />

dB<br />

ξ qdB =<br />

6B + 10.8 – 20 log 10<br />

V<br />

σ x<br />

dB<br />

Ainsi, pour un convertisseur analogique-numérique, où B représente le nombre <strong>de</strong><br />

bits <strong>de</strong>s valeurs <strong>de</strong> sortie, le rapport signal sur bruit <strong>de</strong> une quantification mesuré en<br />

décibels varie linéairement avec B et augmente <strong>de</strong> 6 <strong>de</strong>cibel avec chaque bit<br />

supplémentaire.<br />

5-8


<strong>Bruit</strong> <strong>de</strong> <strong>la</strong> numérisation Séance 5<br />

Métho<strong>de</strong>s <strong>de</strong> conversion A/N<br />

Métho<strong>de</strong>s Indirect : VCO<br />

La valeur <strong>de</strong> <strong>la</strong> tension <strong>du</strong> signal d’entrée est convertie en une fréquence par une<br />

osci<strong>la</strong>teur commandé en tension. en ang<strong>la</strong>is : Voltage Controlled Oscil<strong>la</strong>teur (VCO).<br />

La fréquence est proportionnelle à <strong>la</strong> tension.<br />

Ensuite on compte le nombre <strong>de</strong> passage à zéro pendant ∆T.<br />

Convertiseeur parallèle. (Ang<strong>la</strong>is “F<strong>la</strong>sh Converteur”).<br />

La tension x(n∆T) d’entré est comparé à 2 B – 1 valeurs <strong>du</strong> type : k U 0<br />

2 B<br />

dé<strong>du</strong>ite d’une tension <strong>de</strong> référence U 0 .<br />

Le résultat est tra<strong>du</strong>ite en mot binaire par un déco<strong>de</strong>ur logique.<br />

U 0<br />

(tension <strong>de</strong> référence)<br />

x(t)<br />

R<br />

+<br />

–<br />

2 n –1 comparateurs<br />

R<br />

•<br />

•<br />

•<br />

R<br />

+<br />

–<br />

déco<strong>de</strong>ur<br />

Logique<br />

Sortie<br />

numérique<br />

(codée<br />

en<br />

binaire)<br />

+<br />

–<br />

R<br />

+<br />

–<br />

R<br />

Utiliser pour <strong>la</strong> vidéo, radar, etc.<br />

5-9


<strong>Bruit</strong> <strong>de</strong> <strong>la</strong> numérisation Séance 5<br />

Avantage : Rapi<strong>de</strong> est simple.<br />

Inconvenance : Manque <strong>de</strong> précision.<br />

Approximation sucessive<br />

La tension d’entrée est comparé successivement à une succesion <strong>de</strong> 2 B – 1 valeurs<br />

<strong>de</strong> référence pondérées k U 0<br />

2 N.<br />

C’est un système bouclé qui inclut un convertisseur numérique-analogique :<br />

x(t)<br />

+<br />

Σ<br />

–<br />

+<br />

–<br />

Horloge et logique<br />

<strong>de</strong> comman<strong>de</strong><br />

interne<br />

• • •<br />

} Sortie<br />

numérique<br />

Convertisseur<br />

N/A<br />

Tension <strong>de</strong><br />

Référence U 0<br />

5-10


<strong>Bruit</strong> <strong>de</strong> <strong>la</strong> numérisation Séance 5<br />

Conversion Numérique-analogique<br />

Un convertisseur numérique-analogique est un dispositif pro<strong>du</strong>isant une gran<strong>de</strong>ur<br />

<strong>de</strong> sortie y qui possè<strong>de</strong> 2 N valeurs distinctes. Il est à loi uniforme si ces valeurs<br />

sont régulièrement réparties sur une p<strong>la</strong>ge <strong>de</strong> valeurs al<strong>la</strong>nt <strong>de</strong> zéro à 2 N · q selon <strong>la</strong><br />

loi :<br />

y = q (d 1 2 B-1 + d 2 2 B-2 + ... + d B 2 0 )<br />

où q est le pas <strong>de</strong> quantification.<br />

Des interrupteurs commandés par les variables binaires d k (d k = 0 interrupteur<br />

ouvert, d k = 1 interrupteur fermé) contrôlent le passage <strong>de</strong> courants pondérés I 0 / 2 N<br />

provenant <strong>de</strong> source <strong>de</strong> cournat dépendant d’une référence I 0 vers un point <strong>de</strong><br />

sommation.<br />

I 0 /2 d 1<br />

•<br />

•<br />

•<br />

I 0 /4 d 2<br />

I = d 1 I 0 /2 + d 2 I 0 /4 + ... + d B I 0 /2 B<br />

•<br />

•<br />

I •<br />

0 /2 B d B<br />

Dans <strong>la</strong> pratique, tous les interrupteurs ne réagissent pas exactement au même<br />

instant. Il ne résulte <strong>de</strong>s parasites <strong>de</strong> commutations (“Glitches” en ang<strong>la</strong>is) qui<br />

doivent être éliminés par filtrage.<br />

5-11

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!