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
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