Total 2
Total 2
Total 2
ELEN017-0
Analyse et conception des systmes de
tlcommunications
Remerciements
Merci J.-M. WAGNER, F. L APIERRE , C. D EMOULIN et V. P IERLOT pour leur aide prcieuse
et pour leur dvouement. Merci galement G.-B. S TAN pour ses commentaires et V. B ROUN pour
la rponse en gain du filtre ADSL.
Notations
Cadre analytique :
N
N = {0, 1, 2, . . .}
Z
Z = {. . . 2, 1, 0, 1, 2, . . .}
R
Ensemble des rels
Constantes :
j
kB
Variables et fonctions :
t
Variable de temps continue
f
Variable frquence
Variable de pulsation = 2 f
T
Priode dun signal
T0
Priode fondamentale dun signal
x
Pas dchantillonnage en x
y
Pas dchantillonnage en y
x(t)
Fonction valeurs continues, dfinie pour tout temps t
h(t)
Rponse impulsionnelle dun systme
H (f)
Transmittance dun systme linaire
f [n]
Fonction chantillonne, fonction discrte
Variables et fonctions de modulation :
m(t)
Signal modulant, normalis 1 : |m(t)| 1
s(t)
Signal modul
ka
Taux de modulation
Indice de modulation
Eb
nergie par bit
D
Dbit binaire [b/s]
R
Rapidit de modulation [baud]
ii
Notations fonctionnelles :
|a|
Valeur absolue
a
Vecteur
kak
Norme de a
a.b
Produit scalaire de
a et b
ab
Produit vectoriel de
a et b
.
a
Divergence de
a
a
Rotationnel de a
2
Laplacien de
Gradient de
X
Complexe conjugu de X
Re (a)
Partie relle de a
Im (a)
Partie imaginaire de a
X (f)
Transforme de F OURIER du signal x(t)
Convolution
iii
iv
Notations
ii
Introduction
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
5
5
5
5
6
7
7
7
7
9
11
11
11
12
12
13
15
16
17
17
18
18
19
19
19
20
20
20
1.4.2
1.4.3
1.4.4
2
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
22
22
23
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
25
25
26
27
28
28
28
29
30
31
31
32
33
34
35
39
39
40
41
41
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
43
43
43
43
43
44
44
44
45
45
45
46
46
46
47
47
47
48
48
48
49
49
50
3.5
Facteur de mrite . . . . . . . . . . . . . . . . . .
3.4.3 Temprature de bruit effective du quadriple . . .
3.4.4 Lien entre facteur de bruit et temprature effective
3.4.5 Cas particulier : quadriple attnuateur rsistif . .
Cascade de quadriples . . . . . . . . . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
50
51
51
51
52
53
55
55
55
55
56
56
56
57
57
58
58
58
59
59
59
60
60
61
61
62
62
62
62
64
66
66
67
67
67
Modulations numriques
4.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.2 Dfinition des modulations numriques . . . . . . . . . . . . . . . . . . . . . . . .
4.2.1 Formulation gnrale . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.2.2 Typologie des modulations . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.3 Modulations linaires classiques . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.3.1 Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.3.2 Calcul de la densit spectrale de puissance . . . . . . . . . . . . . . . . . .
Densit spectrale du puissance du signal modul . . . . . . . . . . . . . . .
Fonction dautocorrlation. . . . . . . . . . . . . . . . . . . . . . .
Densit spectrale de puissance. . . . . . . . . . . . . . . . . . . . .
Densit spectrale de lenveloppe complexe . . . . . . . . . . . . . . . . . .
La densit spectrale en conclusion . . . . . . . . . . . . . . . . . . . . . . .
4.3.3 Modulation damplitude numrique (Amplitude Shift Keying) . . . . . . . .
Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Cas particulier : modulation damplitude numrique deux tats (ASK-2) .
Densit spectrale de puissance de lASK-2 . . . . . . . . . . . . . . . . . .
4.3.4 Modulation de phase numrique (Phase Shift Keying) . . . . . . . . . . . .
Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Cas particulier : modulation PSK-2 ou BPSK . . . . . . . . . . . . . . . . .
Densit spectrale de puissance. . . . . . . . . . . . . . . . . . . . .
4.3.5 Modulation en quadrature de phase (Quadrature Phase Shift Keying) . . . . .
Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Modulateur QPSK . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Dmodulateur QPSK . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Densit spectrale de puissance . . . . . . . . . . . . . . . . . . . . . . . . .
4.4 Modulations linaires dcalage (Offset) . . . . . . . . . . . . . . . . . . . . . . .
4.4.1 Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.4.2 Calcul de la densit spectrale de puissance . . . . . . . . . . . . . . . . . . .
4.4.3 Modulation en quadrature de phase dcalage (Offset Quadrature Phase Shift
Keying) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Modulateur et dmodulateur OQPSK . . . . . . . . . . . . . . . . . . . . .
Densit spectrale de puissance . . . . . . . . . . . . . . . . . . . . . . . . .
4.4.4 Modulation saut de phase minimum (Minimum Shift Keying) . . . . . . .
Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Densit spectrale de puissance . . . . . . . . . . . . . . . . . . . . . . . . .
vii
69
69
70
71
72
72
75
talement de spectre
6.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . .
6.1.1 Techniques dtalement . . . . . . . . . . . . .
6.2 talement direct . . . . . . . . . . . . . . . . . . . . . .
6.2.1 Principes de base de ltalement direct . . . . . .
6.2.2 Gnration des squences pseudo-alatoires . . .
Proprits dune squence de longueur maximale
Choix dune squence de longueur maximale . .
6.2.3 Principe de la transmission en bande de base . .
6.2.4 talement direct par modulation BPSK . . . . .
Bruit . . . . . . . . . . . . . . . . . . . . . . .
Analyse . . . . . . . . . . . . . . . . . . . . . .
6.3 tude des performances . . . . . . . . . . . . . . . . . .
6.3.1 Gain dtalement . . . . . . . . . . . . . . . . .
6.3.2 Probabilit derreur . . . . . . . . . . . . . . . .
6.3.3 Marge dinterfrence . . . . . . . . . . . . . . .
6.3.4 Capacit dun systme talement de spectre . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
77
77
80
81
85
87
87
87
88
88
90
91
93
94
98
98
99
101
101
101
102
103
105
105
105
105
106
106
106
108
108
108
108
111
111
112
112
112
113
113
116
116
116
117
tude du trafic
8.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
8.1.1 Caractrisation du trafic . . . . . . . . . . . . . . . . . . . . . . . . .
Intensit, trafic et charge : quelles mesures ? . . . . . . . . . . . . . . .
Charge. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Charge de rfrence. . . . . . . . . . . . . . . . . . . . . . . .
8.2 Analyse statistique des appels . . . . . . . . . . . . . . . . . . . . . . . . . .
8.2.1 Processus de dnombrement . . . . . . . . . . . . . . . . . . . . . . .
8.2.2 Variable alatoire de P OISSON : dfinition et proprits . . . . . . . . .
8.2.3 Processus de P OISSON . . . . . . . . . . . . . . . . . . . . . . . . . .
8.2.4 Temps entre occurrences . . . . . . . . . . . . . . . . . . . . . . . . .
8.2.5 Temps de communication . . . . . . . . . . . . . . . . . . . . . . . .
Hypothse 1 : temps de communication constant . . . . . . . . . . . .
Hypothse 2 : temps de communication dcroissance exponentielle .
8.3 Lois danalyse du trafic . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
8.3.1 Modle dun systme avec perte (sans mmoire) : statistique dE RLANG
Nombre de tentatives dappel . . . . . . . . . . . . . . . . . . . . . .
Nombre darrts dappel (relchements de ligne) . . . . . . . . . . . .
Charge . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Formule dE RLANG B. . . . . . . . . . . . . . . . . . . . . . .
Interprtation. . . . . . . . . . . . . . . . . . . . . . . . . . . .
Dimensionnement. . . . . . . . . . . . . . . . . . . . . . . . .
Table. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Tentatives et essais reconduits. . . . . . . . . . . . . . . . . . .
8.3.2 Autres modles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. .
. .
. .
. .
. .
. .
. .
. .
. .
. .
. .
. .
. .
. .
B
. .
. .
. .
. .
. .
. .
. .
. .
. .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
119
119
120
120
121
122
123
123
124
126
126
128
128
128
129
129
129
129
129
131
131
132
132
132
134
135
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
137
137
137
138
138
138
139
139
139
140
140
140
141
141
141
142
142
143
143
9.3
9.4
9.5
9.6
Exemples de lignes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
9.3.1 Lignes paires symtriques . . . . . . . . . . . . . . . . . . . . . .
Comportement des cbles tlphoniques basses frquences . . . . .
Comportement des cbles tlphoniques hautes frquences . . . . .
Lignes pour transmissions numriques . . . . . . . . . . . . . . . .
9.3.2 Lignes paires coaxiales . . . . . . . . . . . . . . . . . . . . . . . .
Constitution . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Proprits . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Rseau tlphonique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Fonctions de transfert de puissance diaphoniques et transmission haut dbit
9.5.1 Principe . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
9.5.2 Diaphonie dans le rseau tlphonique . . . . . . . . . . . . . . . . .
Modle des capacits non quilibres . . . . . . . . . . . . . . . . .
Notions de circuit. . . . . . . . . . . . . . . . . . . . . . . .
Calcul des relations entre courants et tensions. . . . . . . . .
Modle des inductances non quilibres . . . . . . . . . . . . . . . .
Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
9.5.3 Calcul du NEXT . . . . . . . . . . . . . . . . . . . . . . . . . . . .
9.5.4 Calcul du FEXT . . . . . . . . . . . . . . . . . . . . . . . . . . . .
9.5.5 Rapport signal sur bruit . . . . . . . . . . . . . . . . . . . . . . . . .
9.5.6 Influence du nombre de perturbateurs . . . . . . . . . . . . . . . . .
Estimation de la capacit de canal . . . . . . . . . . . . . . . . . . . . . . .
9.6.1 Information, incertitude et entropie . . . . . . . . . . . . . . . . . .
Notion dinformation . . . . . . . . . . . . . . . . . . . . . . . . . .
Notion dentropie . . . . . . . . . . . . . . . . . . . . . . . . . . . .
9.6.2 Canal discret sans mmoire . . . . . . . . . . . . . . . . . . . . . .
Modle de canal . . . . . . . . . . . . . . . . . . . . . . . . . . . .
9.6.3 Information mutuelle . . . . . . . . . . . . . . . . . . . . . . . . . .
9.6.4 Capacit de canal . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Estimation de la capacit dun canal binaire symtrique . . . . . . . .
Influence du type de modulation . . . . . . . . . . . . . . . . . . . .
Gnralisation au cas continu et second thorme de S HANNON . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
144
144
144
145
145
145
145
146
146
149
149
150
150
151
151
154
154
155
157
157
158
158
158
158
159
159
160
160
161
161
162
163
.
.
.
.
.
.
.
.
.
.
.
.
.
.
165
165
166
167
167
169
169
169
169
170
171
172
172
174
175
.
.
.
.
.
Glossaire
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
176
176
176
178
178
179
179
181
183
183
183
184
185
186
186
186
.
.
.
.
.
189
189
189
190
190
191
196
xi
xii
Introduction
Ces notes abordent plusieurs questions avances des systmes de tlcommunications ; elles comportent trois parties consacres aux thmes suivants :
1. les aspects de modlisation en vue de lexploitation par simulations numriques.
2. les aspects de systme.
3. les aspects avancs relatifs la couche physique du modle OSI.
Dans la partie relative la modlisation, il sagit danalyser les points cls intervenant lors de lutilisation et lors de la conception de logiciels de tlcommunications. Nous parcourons les thmes de
lestimation spectrale, de la reprsentation des signaux moduls en termes de passe-bas quivalent et
de la manire de considrer un bruit dans un modle. Cette mme partie approfondit la caractrisation
des signaux numriques moduls ainsi que les considrations de transmission de signaux numriques
dans un canal.
Les ressources tant limites, il convient de les partager au mieux. Dans la partie consacre aux
aspects systme, nous tudions les diffrents moyens usuels de partager les ressources (par rpartition
de frquence ou rpartition temporelle) ainsi que le procd plus rcent de partager par talement
de spectre au moyen de codes. Ltude du partage des ressources est doubl dune tude du trafic
dbouchant sur des rgles de dimensionnement applicable au trafic tlphonique.
La dernire partie est consacre ltude dlments de la couche physique. Nous analyserons les
dtails de la transmission sur lignes paires torsades en prsence de diaphonie ainsi que la manire de
traiter les phnomnes dvanouissement et les bilans de puissance dans le contexte de transmissions
mobiles.
Premire partie
Modlisation des systmes de
tlcommunications en vue de lexploitation
par simulations numriques
Chapitre 1
Thorie de lestimation et estimation
spectrale
1.1 Introduction
Si par le pass les dveloppements thoriques primaient, lextraordinaire puissance des ordinateurs fait quune srie de problmes insolubles analytiquement trouvent une solution sous forme de
simulations. On peut songer en particulier tous les phnomnes non-linaires pour lesquels on procdait par linarisation pour obtenir des solutions pratiques. Aujourdhui, les simulations permettent
daborder toute une srie de problmes nouveaux.
(1.1)
(1.2)
Dans le cas de processus alatoires gaussiens stationnaires, les moments des deux premiers ordres suffisent dfinir
compltement la densit de probabilit conjointe et, par consquent, le processus [8].
Cette moyenne est la moyenne statistique de la variable alatoire au temps t. Bien entendu,
elle ne peut sobtenir que par modlisation puisquelle fait intervenir les valeurs de toutes les
ralisations possibles au temps t.
Autocorrlation
XX (t1,t2 ) = E {X (t1)X (t2)}
(1.3)
X (t) = X = constante
XX (t1 ,t2) = XX (t2 t1 ) = XX ( )
(1.4)
(1.5)
On dfinit encore la densit spectrale (de puissance) comme la transforme de F OURIER de la fonction dautocorrlation
Z +
X ( f ) =
(1.6)
XX ( ) e2 j f d
La puissance dun processus stationnaire et ergodique E X 2 (t) vaut alors XX (0).
Lanalyse des systmes de tlcommunications comporte souvent ltude de loccupation spectrale ; cest un critre de comparaison important. On sait par exemple que le passage de lanalogique
au numrique entrane gnralement un accroissement de la bande passante ; en contrepartie, le signal
numrique offre une meilleure rsistance au bruit.
Dans le cas dun signal dterministe plac lentre dun systme linaire invariant en temps,
la dtermination de la bande occupe par le signal est relativement aise. La question est plus dlicate quand il sagit dtudier les signaux numriques. En effet, on est souvent amen mettre des
hypothses sur le fonctionnement du systme de tlcommunications (type de bruit, transmittance
du canal, . . . ) pour arriver une formulation analytique. Dans certains cas, il nest pas possible de
choisir pralablement un jeu dhypothses. On procde alors par modlisation, mesure, validation et
interprtation du rsultat. Cest l tout lenjeu de la thorie de lestimation.
Quant il sagit destimer la densit spectrale de puissance ou la transforme de F OURIER laquelle
elle se ramne sous certaines hypothses, on parle destimation spectrale.
Lestimation spectrale est un domaine important de traitement numrique des signaux et de nombreux ouvrages abordent la question [5, 6, 15]. Il sagit en fait de dterminer la densit spectrale de
puissance, vu limportance physique quelle reprsente. Le contexte de lestimation spectrale tient
principalement trois observations
le processus analyser est de nature alatoire. Son observation rsulte en une squence alatoire qui correspond une trajectoire parmi toutes les trajectoires possibles. Il convient donc de
complter lobservation dune srie dhypothses, typiquement la stationnarit au sens strict ou
au sens large et lergodisme, pour tablir une correspondance entre la trajectoire observe et le
processus sous-jacent ;
lobservation du processus est limite dans le temps. Ds lors, lintgration temporelle queffectue la transforme de F OURIER se limite un intervalle de temps. Cette limitation entrane
un effet bien connu de fentrage, ce qui affecte lallure des rsultats. Plus fondamentalement,
on est en droit de sinterroger sur la validit dune estimation effectue sur base dun intervalle
temporel plus ou moins long dans la mesure o les processus physiques considrs ne sont pas
rellement stationnaires.
lobservation seffectue par chantillonnage temporel. On doit donc tenir compte des effets
induits par lchantillonnage tant au niveau du repli de spectre quau niveau de la ralisation de
limplmentation (ajout de valeurs nulles dans la squence pour augmenter la rsolution, etc).
6
Par la suite, nous allons principalement nous concentrer sur lestimation de paramtres dune squence
chantillonne.
1.2 Dfinitions
1.2.1 Squence stochastique
On part dun processus stochastique X (t) pralablement chantillonn. La squence dchantillons
ainsi obtenus est note X [n]. Lobservation dune trajectoire possible parmi toutes les trajectoires
possibles fournit un signal x[n] ; on parle dune ralisation du processus.
Le processus X [n] est entirement dcrit si, quel que soit le choix des indices n1 , n2 , . . ., nm , on
connat la fonction de densit de probabilit conjointe de X [n1], X [n2 ], . . ., X [nm ]. Ce nest gnralement pas le cas.
On tend sans peine les notions valables pour les processus continus des squences. Ainsi, on
dira dune squence quelle est stationnaire au sens large si
et
E {X [n]} =
(1.7)
(1.8)
Ds lors que la squence est stationnaire au sens strict large, sa moyenne est une constante et sa
fonction dautocorrlation dpend uniquement de lcart temporel entre les variables chantillonnes.
Exemple. Porteuse module par une phase alatoire .
Soit la squence alatoire suivante,
X [n] = Ac cos(2 fc n + )
(1.9)
o Ac et fc sont des constantes et est une variable alatoire uniformment rpartie sur [0, 2 ]. Par
calcul,
E {X [n]} = 0
(1.10)
et
XX [n, n + k] = E {X [n]X [n + k]}
A2c
=
cos(2 fck)
2
(1.11)
(1.12)
ce qui signifie quil sagit bien dun processus stationnaire au sens large.
X (f) =
n=
(1.13)
xs (t) =
n=
x[nTs ] (t nTs )
(1.14)
(1.15)
Z + +
n=
Z +
+
x[nTs ]
(t nTs )e2 j f t dt
n=
(1.16)
(1.17)
n=
X ( f + fs ) =
n=
1
Ts
; la connaissance
(1.18)
n=
Ce phnomne est bien connu : lchantillonnage dans le domaine temporel entrane lapparition
de copies de la transforme au droit des frquences multiples de fs . Il est dusage de dfinir une
frquence normalise ou frquence rduite F par
F=
f
fs
(1.19)
de sorte que F parcourt lintervalle [0, 1[ ou [ 21 , 12 [. Si en plus de lutilisation de la frquence normalise, on adopte lcriture x[n] en lieu et place de x[nTs ], la transforme de F OURIER temps discret
sexprime sous la forme normalise suivante.
Dfinition 2 [Transforme de F OURIER temps discret normalise]
+
X (F) =
x[n]e2 jFn
(1.20)
X (F)e2 jFn dF
(1.21)
X (F)e2 jFn dF
(1.22)
n=
x[n] =
Z 1
ou
1
2
21
Convergence
On montre [5, page 45] que si la squence x[n] est de module sommable, cest--dire si
+
n=
kx[n]k < +
8
(1.23)
n=
kx[n]k2 < +
(1.24)
sans tre de module sommable, alors la srie converge en moyenne quadratique. Il peut ne pas y avoir
convergence uniforme. La fonction kX ( f )k2 est dsigne par le terme de spectre. Dans la littrature,
ce terme est aussi associ la fonction kX ( f )k. Les deux expressions deviennent quivalentes si lon
utilise lchelle logarithmique des dcibels.
Exemple. Transforme de F OURIER dun signal rectangulaire.
Considrons une fonction rectangulaire rectN [n] valant 1 pour n {N, . . . , N} et 0 ailleurs. On
calcule alors
N
X (F) =
x[n]e2 jFn
n=N
2 jFN
= e
(1 + . . . + e2 jF(2N+1) )
(1 e2 jF(2N+1) )
1 e2 jF
sin((2N + 1) F)
=
sin( F)
= e2 jFN
(1.25)
(1.26)
(1.27)
(1.28)
Ce signal est purement rel ; il est reprsent la figure 1.1 pour deux valeurs de N.
La transforme inverse a pour expression
x[n] =
Z 1
X (F)e2 jFn dF
(1.29)
(1.30)
(1.31)
(1.32)
nergie
Tout comme dans le cas dune fonction continue, lnergie du signal est conserve par rapport
sa transforme de F OURIER.
Dfinition 3 Lnergie totale du signal vaut
+
E=
n=
kx[n]k2
(1.33)
n=
kx[n]k =
Z 1
0
kX (F)k2 dF
(1.34)
N=2
5
4
3
2
1
0
1
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
0.6
0.7
0.8
0.9
N=5
12
10
8
6
4
2
0
2
0.1
0.2
0.3
0.4
0.5
10
Squence priodique
Le cas des squences priodiques est particulier car cette priodicit rsulte gnralement dune
opration dchantillonnage ; il sagit alors dune priodicit fictive.
Dfinition 5 Une squence est priodique sil existe une valeur entire N telle que
x[n] = x[n + N]
(1.35)
X (Fk ) =
x[n]e2 jFkn,
Fk =
n=0
k
, k {0, . . . , N 1}
N
(1.36)
n=0
kx[n]k2 =
1 N1
k
kX (Fk )k2 , Fk =
N k=0
N
(1.37)
Comme le signal priodique a une nergie infinie, on utilisera plutt la notion de puissance moyenne
de la squence dfinie par
1 N1
P = kx[n]k2
(1.38)
N n=0
soit encore
k
1 N1
P = 2 kX (Fk )k2 , Fk =
N n=0
N
(1.39)
Prcision
Pour illustrer lutilisation de la transforme de F OURIER discrte dans lobservation de spectres
de signaux, considrons la suite x[n] obtenue par chantillonnage de lexponentielle complexe e2 jF0t
la cadence Fs = 1/Ts (cet exemple est tir de [5]). En posant f0 = F0 /Fs < 1, on a x[n] = e2 j f0 n .
Le fait de rduire la dure dobservation un intervalle N fait apparatre des ondulations dans la
transforme de F OURIER temps discret du signal. Ce phnomne est illustr la figure 1.2 pour
7
.
N = 32 et f0 = 32
35
30
25
20
15
10
0
0.5
0.4
0.3
0.2
0.1
0.1
0.2
0.3
0.4
0.5
7
avec N = 32.
F IG . 1.2 Transforme de F OURIER temps discret de lexponentielle complexe f0 = 32
35
30
25
20
15
10
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
7
32
avec N = 32.
que, si x[n] contient deux sinusodes dont les frquences sont spares de moins de 1/N, les lobes
principaux de chacune delles seront si proches quil sera difficile de les distinguer avec certitude. Et
ceci est dautant plus vrai que les deux amplitudes sont trs diffrentes. Ainsi, si fs = 1/Ts dsigne la
frquence dchantillonnage, la rsolution R exprime en [Hz] est de lordre de grandeur de fs /N qui
est prcisment linverse du temps total danalyse, savoir T = NTs .
Proposition 8 La rsolution en frquence est de lordre de grandeur de linverse du temps total
danalyse.
Le produit R T joue en quelque sorte le rle de facteur de mrite dans lutilisation de la transforme
de F OURIER temps discret pour la recherche de frquences. Typiquement, si on choisit R et T tels
que le produit R T 3, les frquences seront faciles sparer.
Spectre court terme
Comme lanalyse en frquence classique est une opration qui couvre la totalit de laxe des
frquences, cette opration effectue une moyenne sur laxe des temps. Cet effet moyenneur peut aller
jusqu occulter les phnomnes observer. Cela ne signifie pas que linformation soit dtruite. La
transforme de F OURIER est en effet bijective sous les conditions dexistence nonces prcdemment.
Prenons lexemple dune onde continue x(t) constitue de deux portions successives de frquences
f1 et f2 (cf. figure 1.6).
Sa transforme de F OURIER contient linformation concernant lordre dans lequel apparaissent les
sinusodes. Toutefois cette information nest pas explicite, mme si elle se trouve dans linformation
13
25
20
15
10
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
3 dB
14
0.9
7,5
32
avec N = 32.
0.5
0.5
20
40
60
80
100
120
140
160
180
200
0.05
0.1
0.15
0.2
0.25
0.3
0.35
0.4
0.45
0.5
50
40
30
20
10
0
35
30
25
20
3
15
10
2.5
5
2
0
140
120
100
1.5
80
60
40
20
F IG . 1.7 Deux portions de sinusodes et le spectre correspondant dans une fentre danalyse glissante (daprs [5]).
finie, tout comme dans le cas de signaux priodiques, ce qui va nous amener dterminer sa densit
spectrale de puissance.
1.3.1 Estimateur
Soit estimer le paramtre ; cette quantit est non alatoire.
b.
partir des chantillons {x[0], x[1], . . . , x[N 1]}, on calcule une certaine estimation note
Cette estimation est donc une fonction A des chantillons
b = A(x[0], x[1], . . . , x[N 1])
16
(1.40)
30
20
0
10
5
20
40
60
80
100
0.2
0.4
0.6
0.8
0.2
0.4
0.6
0.8
0.2
0.4
0.6
0.8
0.2
0.4
0.6
0.8
30
20
0
10
5
20
40
60
80
100
30
20
0
10
5
20
40
60
80
100
30
20
0
10
5
20
40
60
80
100
MSEb = b2b + b2
17
(1.42)
(1.43)
(1.44)
Sil est souhaitable que lestimateur soit non biais, les moments dordre 2 jouent un rle capital
quand il sagit dexaminer la convergence des estimateurs. En effet, il ne servirait rien davoir un
estimateur non biais mais avec un cart type lev. Comme il faut que biais et cart type soient tous
deux petits, la meilleure mesure reste lerreur quadratique moyenne.
On parle destimateur consistant si son erreur quadratique moyenne tend vers 0 lorsque la dure
dobservation tend vers linfini, soit
lim MSEb = 0
(1.45)
N+
Cette condition, mettre en parallle avec lergodisme, est essentielle quand il sagit dexploiter
les donnes dchantillons prlevs sur une longue priode.
(1.46)
Elle est appele fonction de vraisemblance de ces paramtres lorsquon y insre les valeurs observes
{x[0], x[1], . . . , x[N 1]}.
Lestimateur maximum de vraisemblance des paramtres {1 , 2 , . . . , p } consiste calculer
b1 ,
b2 , . . . ,
b p } qui maximisent la fonction de vraisemblance. On obtient donc les estiles valeurs {
mations en cherchant la solution du systme de p quations p inconnues
bi = i
f
= 0 pour i = 1, 2, p
i
(1.47)
(1.48)
b < 2 ) = 1
p(1 <
(1.49)
b 2 < <
b 1 ) = 1
p(
(1.50)
b qui est la variable alatoire et que soit une variable certaine, mais inconnue,
Bien que ce soit
b
on a tendance interprter cette relation comme la probabilit que appartienne lintervalle [
b 1 ].
2 ,
18
b0 de
b . Le plus petit intervalle
Aprs observation, on aura obtenu une certaine valeur (ralisation)
b0 2 ,
b0 1 ] tel que
[
b0 2 < <
b0 1 ) = 1
p(
(1.51)
ce qui revient
Z b0 1
b0 2
b )d
b = 1
Fb (
(1.52)
est appel intervalle de confiance. Cet intervalle est interprt comme celui dans lequel on risque de
b avec une probabilit 1 ; cependant, une telle interprtation ntant pas
trouver la valeur vraie
correcte, on utilise le vocable confiance pour la quantit 1 . On parle ainsi dintervalle de confiance
90% ou 99%.
bx =
1 N1
x[n]
N n=0
(1.53)
Biais de bx ?
(1.54)
Essayons prsent den dterminer la variance. Soustrayant membre membre les deux quations
prcdentes, il vient
1 N1
bx E {bx } = (x[n] x )
(1.55)
N n=0
et ensuite
2
c
x =
1 N1 N1
E (x[n] x )(x[n ] x )
2
N n=0 n =0
(1.56)
1 N1 N1
Cxx [n n ]
N 2 n=0
n =0
(1.57)
N1
1
(N |i|)Cxx [i]
N 2 i=(N1)
(1.58)
Cette dernire expression sobtient en observant que, dans la double somme 1.57, la diffrence
n n prend N fois la valeur 0, N 1 fois les valeurs 1, etc.
19
Cette expression met en vidence un problme courant dans le domaine de lestimation. Alors
quon tente ici destimer un moment du premier ordre, il faudrait connatre les moments du deuxime
ordre pour dterminer la prcision de lestimateur : plus gnralement, la prcision de lestimateur des
moments dun ordre dtermin dpend des moments dun ordre suprieur. Ceux-ci sont logiquement
inconnus, et il faudra les estimer eux-mmes pour avoir une ide de la prcision de lestimateur.
Revenons au cas prsent. La double somme 1.57 comprend N 2 termes et, dune manire gnrale,
il se pourrait que la variance de lestimateur ne dcroisse pas lorsque la dure N de lenregistrement
augmente. Cependant, si la fonction de covariance Cxx [i] dcrot avec i et, plus prcisment, si
limi+ Cxx [i] = 0, la simple somme apparaissant dans lquation 1.58 crot proportionnellement N
lorsque N +. Dans ce cas, la variance de lestimateur de bx dcrot comme N1 et lestimateur est
2
2 = x.
consistant. Ainsi, si lon savait que x[n] est un bruit blanc dintensit x2 , on aurait c
N
x
On pourrait croire, tort, que la moyenne arithmtique est le seul estimateur de la moyenne.
En fait, dans les cas o lon ne dispose que dun faible nombre dchantillons, il est intressant de
recourir la mdiane. Si pour certaines distributions alatoires, par exemple la gaussienne, mdiane et
moyenne sont confondues, il nempche quen pratique la mdiane est parfois prfre la moyenne
arithmtique.
|x(t)|2 dt < +
(1.59)
Cette hypothse est lgitime en pratique car les signaux observs proviennent de processus physiques
rels et limits dans le temps. On dfinit alors une pseudo fonction dautocorrlation par
xx ( ) =
Z +
x(t)x(t + )dt
(1.60)
dont lexpression est calque sur lexpression correspondante pour un processus stochastique.
supposer que le signal ait des caractristiques stationnaires, on dfinit galement une densit
spectrale de puissance par
x ( f ) =
Z +
xx ( ) e2 j f d
(1.61)
Z +
x(t)e2 j f t dt
20
(1.62)
xx ( ) e2 j f d
Z + Z +
=
x(t)x(t + )dt e2 j f d
Z +
Z +
2 j f
d dt
x(t)
=
x(t + )e
x ( f ) =
Z +
(1.64)
(1.65)
2 j f t
x(t) X ( f )e
dt
= X (f)
(1.63)
Z +
x(t)e2 j f t dt
= X ( f )X ( f )
= kX ( f )k
(1.66)
(1.67)
(1.68)
(1.69)
Ainsi donc, une estimation simple de la densit spectrale de puissance sur base de lobservation
ce qui prsuppose stationnarit et ergodisme quivaut lestimation de la transforme de F OURIER
au carr de lobservation. Ce nest pas la seule faon de faire car on peut tout aussi bien estimer la
fonction dautocorrlation et en prendre la transforme de F OURIER.
Le dveloppement prcdent a montr le lien existant entre une densit spectrale et la transforme
de F OURIER (au carr) dun signal. Il a nanmoins compltement occult la question de lestimation
car il y a plus quune marge entre le traitement dun signal continu et le traitement dune squence. En
fait, les difficults de lestimation rsultent en grande partie des chantillonnages et dune limitation
de la dure dobservation. La raison pour laquelle la littrature propose de nombreuses mthodes
destimation provient des choix qui sont faits pour obtenir une estimation ayant certaines proprits.
Cest donc plus une question de choix quune question de rigueur ou de fiabilit.
Dans les paragraphes suivants, nous prsentons un estimateur de la fonction dautocorrlation et,
ensuite, divers estimateurs de la densit spectrale de puissance. En tout tat de cause, nous considrons une squence x[n] alatoire stationnaire et ergodique que, pour la simplicit, nous supposons
centre. Les estimations sont dduire dun enregistrement exprimental {x[0], x[1], . . . , x[N 1]} de
longueur N.
N n=0
(1.70)
Cette expression nous permet de calculer une estime de la densit spectrale par transforme de
F OURIER
bX (F) =
=
N1
k=(N1)
1 N1 N|k|1
(1.71)
(1.72)
(1.73)
qui nest rien dautre que la transforme de F OURIER au carr de lenregistrement, prolong par des
zros.
Cet estimateur est appel priodigramme. Il est malheureusement biais et
|k| b
E {bX (F)} =
1 N XX [k]e2 jFk
k=(N1)
N1
(1.74)
Le priodigramme est un estimateur asymptotiquement sans biais de la densit spectrale de puissance : pour une dure dobservation N suffisamment grande, bX (F) fluctue autour de la vraie valeur
de X (F). Par contre, lamplitude des fluctuations, qui est donn par la variance de bX (F), ne tend pas
vers 0 lorsque N tend vers linfini. Plus prcisment, on montre que cette variance est de lordre de
grandeur de la valeur estimer. La figure 1.9 montre quatre priodigrammes et une certaine densit
spectrale thorique. On observe que, malgr lallongement de la dure dobservation, lamplitude des
fluctuations ne diminue pas.
N = 100
N = 200
40
40
30
30
20
20
10
10
10
10
20
20
30
0.1
0.2
0.3
0.4
30
0.5
0.1
N = 500
40
30
30
20
20
10
10
10
10
20
20
0
0.1
0.2
0.3
0.3
0.4
0.5
0.4
0.5
N = 1000
40
30
0.2
0.4
30
0.5
0.1
0.2
0.3
2
N1
1
bX (Fk ) =
x[n]e2 jkn/N
N
n=0
(1.75)
On peut aussi sintresser la forme de la fonction fentre que lon applique. En effet,
2
N1
1
bX (F) =
x[n]e2 jFn
N
n=0
2
+
1
2 jFn
=
x[n]wR [n]e
N
n=
(1.76)
(1.77)
(1.78)
Leffet de la fentre nest pas ngligeable car il affecte le spectre du signal observ. Aussi, diffrentes mthodes destimation spectrale tentent-elles dtablir un compromis entre prcision spectrale et effet de fentrage.
La question de lestimation spectrale est trop vaste pour que nous labordions autrement que par
une sensibilisation la difficult intrinsque estimer un spectre. Pour une tude approfondie, nous
renvoyons le lecteur intress des ouvrages spcialiss [5, 6, 15].
23
24
Chapitre 2
Reprsentation des signaux et des systmes
passe-bande
2.1 Introduction
La numrisation des signaux sopre en plusieurs tapes successives. Dans un premier temps, le
signal original est filtr ; il en dcoule un signal spectre limit 1 . Le signal est alors prt pour une
phase de numrisation qui dbute par un chantillonnage une frquence suprieure la frquence
de N YQUIST. Ainsi, la transmission dun signal de parole filtr 4 [kHz] et numris ncessite tout
au plus quelques dizaines de kb/s.
Sur une ligne tlphonique, la transmission seffectue en bande de base. Aprs tout, la ligne tant
physiquement ddie un utilisateur, il ny a aucune raison particulire de partager le canal en modulant le signal. De plus, on a tout intrt conserver un signal en bande de base puisque lattnuation
y est la plus faible.
La situation est tout autre dans lespace. La ressource tant partage, chaque utilisateur se voit
allouer une bande de frquences spcifique. Pour loccuper, un mlangeur transpose le signal de la
bande de base vers la bande de frquences ddie. Aprs transmission, on ramne le signal en bande
de base et on procde la dmodulation.
On voit donc apparatre deux types de considrations frquentielles :
la bande de base. Cest la bande occupe par le signal original.
la bande utile, bande dans laquelle on vient placer le signal initialement en bande de base. On
parle alors de systme passe-bande puisque, du point de vue de lutilisateur, le systme agit par
transparence dans la bande utile mais il annulle les composantes hors bande.
Cette notion de passe-bande est formalise par la dfinition suivante.
Dfinition 9 [Passe-bande] Un signal dterministe g(t) est de type passe-bande sil existe deux
valeurs f0 et W telle que
W
W
f 6 [ f0 , f0 + ], kG ( f )k = 0
(2.1)
2
2
La figure 2.1 reprsente un signal passe-bande. Il est remarquer que le spectre dun signal rel ne
doit nullement tre symtrique par rapport la frquence f0 . Ainsi, si la modulation damplitude
conduit un spectre symtrique autour de la frquence porteuse, la modulation damplitude bande
rsiduelle produit un spectre non symtrique.
1
En toute rigueur, un signal spectre limit est dure infinie. Mais comme les signaux rels sont dure finie, il y a
une contradiction dfinir un signal spectre limit. On utilise nanmoins cette modlisation.
25
kG ( f )k
f0
f0
W
F IG . 2.1 Spectre dun signal passe-bande.
La modlisation dun signal passe-bande peut paratre plus complexe puisque le signal fluctue
une cadence rythme par la plus haute frquence contenue dans le signal. Fort heureusement, on peut
tirer profit du fait que le signal de dpart se situe en bande de base.
Noublions pas que, comme il sagit de reconstituer ce signal en bande de base le signal modul
nest pas une fin en soi, on peut sinterroger sur la possibilit de modliser les signaux et les systmes
des frquences nettement infrieures la frquence porteuse.
La dmodulation dun signal AM porteuse par dtecteur de crte utilise ce principe. En effet,
un choix judicieux des constantes de temps permet de rcuprer lenveloppe du signal. Hormis lobligation de choisir les constantes de temps en fonction de la borne maximale de la bande de base et
en fonction de la frquence porteuse lors de la synthse dun dtecteur, il apparat que le circuit de
dmodulation travaille dans une plage de frquences telle que dfinie par la bande de base ; il ny a
pas de composante haute-frquence.
Avant daborder plus en dtail la phase de modlisation, examinons la question de lchantillonnage des signaux passe-bande.
gs (t) fs
i=
G ( f i fs )
(2.2)
Dans notre cas, le signal est de type passe-bande. Lchantillonnage pourrait amener un chevauchement des copies des spectres si lon ny prend garde. Prenons une frquence dchantillonnage fs
dfinie par 2kfu . Pour viter tout repli de spectre, il faut imprativement que fs > 2W . Ds lors, on
obtient 2kfu > 2W , ce qui implique k < Wfu .
Exemple. Considrons un signal W = 20 [kHz] et fu = 105 [kHz]. Par application du thorme, k <
5, 25. Autrement dit, la frquence dchantillonnage fs = 42 [kHz] convient. La figure 2.2 reprsente
26
le signal chantillonn dans le domaine spectral. On est videmment loin de la frquence de N YQUIST
qui aurait prvalu sil ne sagissait dun signal passe-bande.
111
000
000
111
000
111
kGs ( f )k
111
000
000
111
000
111
W
fu
(2.3)
telle que f0 est une frquence de rfrence contenue dans la bande utile du signal2. Dans cette galit,
g(t) porte le nom de passe-bas quivalent.
Dans un souci de rigueur, il conviendrait plutt de parler de famille de passe-bas quivalents.
lvidence, toute fonction g(t) = (g(t) + jz(t))e2 j f0t j0 respecte lgalit 2.3.
La dtermination dun type de passe-bas quivalent revient donc choisir une fonction particulire
pour z(t). Comme on pourra sen apercevoir le choix nest pas unique. Pour compliquer laffaire, cer
tains auteurs [24] dfinissent le passe-bas quivalent en introduisant un facteur de normalisation 2
dans lexpression 2.3 afin dviter lapparition dun facteur 2 lors de lanalyse de systmes linaires ;
on trouve les deux dfinitions dans la littrature.
2
0
est souvent choisi nul par commodit. Cela ne nuit en rien la gnralit des propos.
27
(2.6)
Dans le cas dune modulation numrique linaire, les signaux m1 (t) et m2 (t) reprsentent une
squence dimpulsions mises en forme. On peut aisment trouver les quivalents passe-bas des
principales modulations numriques.
Modulation numrique damplitude 2n tats
Un signal numrique 2n tats modul en amplitude (2n -ASK) est de la forme
+
s(t) =
k=
(2.7)
o Ak est une variable alatoire pouvant prendre 2n valeurs diffrentes. Le passe-bas quivalent est
tout simplement
+
s(t) =
k=
Ak Rect[0,T ] (t kT )
(2.8)
s(t) =
k=
(2.9)
o k est une variable alatoire pouvant prendre 2n valeurs diffrentes dans lintervalle [0, 2 ]. Le
passe-bas quivalent est
+
s(t) =
k=
Rect[0,T ] (t kT )e jk
(2.10)
Ce choix convient dans bon nombre de cas mais il arrive quun autre choix soit plus judicieux ; cest la facilit de
formulation rsultante qui guidera le choix final.
28
En effet, en prenant c = 0,
+
s(t) = Re
k=
Rect[0,T ] (t kT )e jk e2 j fct
k=
(2.11)
(2.12)
(2.13)
Lexploitation de la symtrie permet de ne considrer que les frquences positives durant les calculs,
condition bien sr de reconstituer les frquences ngatives pour rtablir la nature du signal rel en
fin de calculs.
Soit donc le filtre H ( f ) supprimant les frquences ngatives ; par la mme occasion, il amplifie
les frquences positives dun facteur 2 de manire conserver le mme niveau dnergie. Ce filtre se
caractrise par la transmittance
0 si f < 0
H (f) =
(2.14)
2 si f 0
qui peut galement se mettre sous la forme
H ( f ) = 1 + sgn( f )
tant donn que la transforme de F OURIER de la fonction sgn est donne par
sgn(t)
1
jf
j
t
Dfinition 12 [Signal analytique] La rponse dun tel filtre un signal dentre g(t) est appele
signal analytique. Elle est note ga (t) et vaut
j
ga (t) = g(t) (t) +
(2.15)
t
1
(2.16)
= g(t) + jg(t)
t
29
ga (t) = g(t) + je
g(t)
(2.18)
Le signal analytique est donc compos dune partie relle, qui nest autre que le signal original,
et dune partie imaginaire qui en est la transforme de H ILBERT. Lapparition dune composante
imaginaire ne doit pas nous tonner car elle rsulte de la rupture de symtrie dans le spectre. La
norme du signal analytique est appele enveloppe du signal.
Le signal ge(t) porte aussi le nom de signal en quadrature du signal dentre. Cette appellation
provient des caractristiques du filtre 1t . En effet, la transmittance de ce filtre vaut jsgn( f ). Ds
lors,
Ge( f ) = jsgn( f )G ( f )
(2.19)
Autrement dit, les frquences positives et ngatives subissent toutes un dphasage de 2 . Par la mme
occasion, on remarque que la transforme de H ILBERT occupe la mme bande de frquence que le
signal original.
Proprits de la transforme de H ILBERT
Proprit 14 Lnergie (ou la puissance) dun signal et celle de sa transforme de H ILBERT sont
gales.
Dmonstration
Lnergie de la transforme de H ILBERT dun signal vaut
Z +
Or
Z +
e
2
ke
g(t)k dt =
G ( f )
d f
2
(2.20)
e
2
2
2
G ( f )
= k jsgn( f )k kG ( f )k
= kG ( f )k2
(2.21)
(2.22)
Proprit 15 [36, page 96] Un signal est orthogonal sa transforme de H ILBERT. Pour un signal
dnergie,
Z +
g(t)e
g(t)dt = 0
(2.23)
Z +T
(2.24)
g(t)e
g(t)dt = 0
30
(2.26)
(2.27)
(2.28)
Ce rsultat ne doit pas nous tonner car il dcoule dun dphasage du cosinus par ajout dune phase
2 . Le mme raisonnement conduit aisment au rsultat suivant.
Proprit 16 [Transforme de H ILBERT dun signal modul] Soit un signal g(t) en bande de base,
^ fct) = g(t) sin(2 fct)
g(t) cos(2
(2.29)
(2.31)
(2.32)
(2.33)
Il sagit donc bien dune forme valide de passe-bas quivalent tel que dfinie par la relation 2.3.
31
Dcomposition de R ICE
Il nous reste maintenant dterminer des moyens pratiques pour calculer lenveloppe complexe
dun signal. Pour ce faire, on dcompose lenveloppe complexe en sa partie relle et sa partie imaginaire
eg (t) = gI (t) + jgQ(t)
(2.34)
Nous avons pris le soin de faire rfrence au signal g pour des raisons qui apparatront bientt. gI (t)
sobtient comme suit
gI (t) = Re (eg (t))
= Re ga (t)e2 j f0t
(2.35)
(2.36)
De mme,
(2.37)
(2.38)
(2.39)
(2.40)
Les signaux gI (t) et gQ (t) sont appels composantes de R ICE ou composantes en quadrature. Sur
le plan thorique, on peut les calculer trs simplement partir du schma dcrit la figure 2.3.
g(t)
cos (2 fct)
H ILBERT
+
gI (t)
+ 2
sin (2 fct)
(2.41)
(2.42)
(2.43)
(2.44)
Cette dernire expression traduit un rsultat important : tout signal passe-bande peut sexprimer
comme la somme de deux signaux confins lintrieur de lintervalle de frquences [ W2 , W2 ] et
moduls en quadrature.
La figure 2.4 montre comment reconstituer le signal original partir de ses composantes de R ICE.
32
gI (t)
cos (2 fct)
+
g(t)
+ 2
sin (2 fct)
gQ (t)
F IG . 2.4 Schma de reconstitution dun signal partir de ses composantes de R ICE.
Lenveloppe complexe peut aussi scrire sous la forme dune amplitude instantane et dune
phase instantane
eg (t) = ag (t)e jg(t)
(2.45)
On a alors les trois relations
(2.46)
(2.47)
(2.48)
Ds lors, on trouve un second rsultat important : tout signal passe-bande peut sexprimer sous
la forme dun signal modul en phase et en amplitude. Le passage de g(t) dans un dtecteur de crte
produit ag (t). Dautre part, moyennant la suppression de la modulation rsiduelle, un dtecteur de
phase ou un discriminateur de frquences sera en mesure de dterminer lallure de la modulation
angulaire.
En pratique
Les expressions 2.37 et 2.40 de calcul des composantes de R ICE, pour rigoureuses quelles soient,
ne sont gure pratiques car elles exigent la dtermination pralable de la transforme de H ILBERT.
Un moyen trs facile consiste procder linstar dun dmodulateur damplitude : on multiplie le
signal original par 2 cos(2 f0t) (ou 2 sin(2 f0t)) et on applique un filtre passe-bas, comme illustr
par le schma 2.5.
En effet,
g(t) 2 cos(2 f0t) = 2 [gI (t) cos(2 f0t) gQ(t) sin(2 f0t)]cos(2 f0t)
= 2 gI (t) cos2 (2 f0t) gQ (t) sin(2 f0t) cos(2 f0t)
= gI (t) + gI (t) cos(4 f0t) gQ (t) sin(4 f0t)
(2.49)
(2.50)
(2.51)
gI (t)
2 cos (2 f0t)
g(t)
+ 2
2 sin (2 f0t)
gQ (t)
F IG . 2.5 Schma de calcul pratique des composantes de R ICE.
(2.53)
Filtrer un signal rel g(t) par un filtre de rponse impulsionnelle h(t) peut se faire au choix :
par convolution g(t) h(t) ;
en filtrant le signal analytique ga (t) par un systme ayant ha (t) pour rponse impulsionnelle ;
ou encore
en filtrant lenveloppe complexe eg (t) par eh (t).
=
Comme on peut crire
et
Ds lors,
h( )g(t )d
1
2 j f0t
2 j f0t
eg (t)e
+ eg (t)e
g(t) =
2
1
2 j f0t
2 j f0t
h(t) =
eh (t)e
+ eh (t)e
2
Z
1 2 j f0t +
eh ( )eg (t )d
y(t) =
e
4
Z
1 2 j f0t +
+ e
eh ( )eg (t )d
4
Z
1 2 j f0t +
+ e
eh ( )eg (t )e4 j f0 d
4
Z
1 2 j f0t +
+ e
eh ( )eg (t )e4 j f0 d
4
34
(2.54)
(2.55)
(2.56)
(2.57)
(2.58)
(2.59)
(2.60)
(2.61)
Les troisime et quatrime termes tendent vers 0 lorsque la frquence de rfrence f0 est largement
suprieure aux frquences utilises dans la bande de base. Dans ce cas en effet, le signal eh ( )eg (t
) fluctue peu pendant une priode de lexponentielle et, moyennant lhypothse raisonnable que ce
signal est constant sur toute la dure dune priode, lintgrale se compense demi-cycle par demicycle, pour sannuler. Il reste alors les deux premiers termes
1
y(t) = Re ((eh (t) eg(t))e2 j f0t )
2
(2.62)
1
eh (t) eg(t)
2
(2.63)
(2.64)
fc W2 f fc + W2
2 ja f fc + W2
H1 ( f ) =
(2.65)
2 ja f + fc W2
fc W2 f fc + W2
0
ailleurs
H1 ( f )
j
fc
fc
Eh1 (1)
j
pente = 4 a
W2
H2 ( f )
j
fc
W
2
fc
36
(2.67)
Rt
0 m(t
) dt
(2.68)
= Ac e
(2.69)
Rt
0 m(t ) dt
(2.70)
Es1 ( f ) =
(2.71)
f W2 , W2
ailleurs
(2.72)
m (t) e2 jk f 0 m(t ) dt
(2.74)
es1 (t) = jWaAc 1 +
W
La rponse s1 (t) du filtre rampe est alors donne par
2 j fct
s1 (t) = Re es1 (t) e
Z t
2k f
= WaAc 1 +
m (t) cos 2 fct + 2 k f
m t dt +
W
2
0
(2.75)
(2.76)
Le signal s1 (t) est un signal modul hybride dont lamplitude et la frquence varient avec le signal
modulant m(t). Cependant, si
2k f
<1
(2.77)
m
(t)
W
37
pour tout t, nous pouvons utiliser un dtecteur denveloppe pour retrouver les variations damplitude
de s1 (t) et donc, un terme constant prs, obtenir le signal modulant m(t). La sortie du dtecteur
denveloppe est alors gale
2k f
m (t)
(2.78)
kes1 (t)k = WaAc 1 +
W
(2.79)
Appelons s2 (t) la rponse du filtre rampe complmentaire au signal FM s(t). Alors, en suivant une
procdure similaire celle qui a t faite pour s1 (t), nous pouvons dterminer la sortie du second
dtecteur denveloppe
2k f
m (t)
(2.80)
kes2 (t)k = WaAc 1
W
o es2 (t) est lenveloppe complexe de s2 (t). La diffrence entre les sorties des deux dtecteurs denveloppe est le signal s0 (t) donn par
s0 (t) = kes1 (t)k kes2 (t)k
= 4 k f aAc m (t)
(2.81)
(2.82)
denveloppe
+
signal FM
Dtecteur
H2
signal en
bande de base
denveloppe
sortie
signal FM
Tension de sortie
1,0
2B
0,707
fc 2B
Rponse du filtre
oscillant infrieur
Rponse du filtre
oscillant suprieur
fc fc + 2B
f
Rponse globale
-0,707
-1,0
2B
3B
F IG . 2.9 Rponse frquentielle des deux circuits rsonnants [16, page 179].
39
le spectre du signal FM s(t) nest pas strictement limit lintervalle de frquence [ fc W /2, fc +W /2],
les sorties des circuits rsonnants ne sont pas strictement bande limite et des distorsions
peuvent tre introduites par les filtres passe-bas RC suivant les diodes dans les dtecteurs denveloppe,
la rponse frquentielle des circuits rsonnants nest pas linaire sur toute la bande de frquence
du signal FM.
Cependant, une conception soigne permet de maintenir les distorsions dans des proportions acceptables.
Xa ( f ) = kH ( f )k2 X ( f )
4 X ( f ) si f 0
=
0
si f < 0
(2.84)
(2.85)
Observons que la densit spectrale de puissance du signal analytique peut encore scrire sous la
forme Xa ( f ) = 2H ( f ) X ( f ). Comme la rponse impulsionnelle du filtre H ( f ) est gale
h(t) = (t) +
j
t
(2.86)
il est facile de calculer la fonction dautocorrlation du signal analytique par transforme de F OURIER
inverse
eXX ( ))
Xa Xa ( ) = 2(XX ( ) + j
(2.87)
Les densits spectrales de puissance de lenveloppe complexe et du signal analytique sont alors lies
par
eX ( f ) = Xa ( f + f0 )
(2.89)
On peut encore montrer que lon peut retrouver la fonction dautocorrlation du processus stochastique X (t) partir de celle de lenveloppe complexe
1
Re (Xa Xa ( ))
2
1
Re eX eX ( ) e2 j f0
=
2
XX ( ) =
(2.90)
(2.91)
Il est finalement possible dexprimer la densit spectrale de puissance de X (t) en fonction de celle de
lenveloppe complexe. Il suffit de remarquer que la fonction dautocorrlation de X (t) peut scrire
sous la forme
i
1h
eX eX ( ) e2 j f0 + eX eX ( ) e2 j f0
(2.92)
XX ( ) =
4
et de calculer la transforme de F OURIER de cette dernire expression
X ( f ) =
eX ( f f0 ) + eX ( f f0 )
4
(2.93)
Nanmoins, il est intressant dtudier lenveloppe complexe, non plus en terme de fonction dautocorrlation, mais directement via son expression temporelle.
X (t) tant un processus stochastique, lenveloppe complexe eX (t) est galement un processus stochastique. crit tel quel, le processus stochastique X (t) nest pas stationnaire car sa moyenne dpend du
temps, que lenveloppe complexe soit stationnaire ou non. Comme pour le cas des signaux dterministes, il est ncessaire de stationnariser le signal en introduisant une phase alatoire uniformment
rpartie sur [0, 2 [. Ainsi, on crit
X (t) = Re eX (t) e
j(2 f0t+)
(2.95)
partir de cette expression, il est possible de retrouver la fonction dautocorrlation de X (t). Il suffit
de remarquer que
X (t) =
Il vient
i
1h
eX (t) e j(2 f0t+) + eX (t) e j(2 f0t+)
2
XX (t, t ) = E {X (t) X (t )}
= E {X (t) X (t )}
41
(2.96)
o nous avons utilis le fait que X (t) est rel. En remplaant X (t) par sa valeur, il vient aprs calcul
des esprances
XX (t, t ) =
=
=
=
o
1 n
2 j f0
2 j f0
+ eX (t) eX (t ) e
E eX (t) eX (t ) e
4
n
o
1
2 j f0
E 2 Re eX (t) eX (t ) e
4
o
1 n
Re E eX (t) eX (t ) e2 j f0
2
1
Re eX eX (t, t ) e2 j f0
2
(2.97)
(2.98)
(2.99)
(2.100)
Ds lors, si le processus stochastique X (t) est stationnaire, son enveloppe complexe lest galement
et on peut crire
1
XX ( ) = Re eX eX ( ) e2 j f0
(2.101)
2
Il est important de remarquer que XI (t) et XQ (t) sont galement des processus stochastiques.
42
(2.103)
(2.104)
(2.105)
Chapitre 3
Calcul du bruit dans les systmes de
tlcommunications
3.1 Introduction
La prsence du bruit joue un rle essentiel dans la conception des systmes de tlcommunications. On distingue principalement deux classes de sources de bruit :
les sources externes au systme : bruit atmosphrique, bruit solaire, bruit cosmique, etc,
les sources internes au systme.
Ce ne sont pas tant les sources de bruit et leur nature physique qui nous intressent mais la manire
dont le bruit affecte les signaux utiles.
E ( f ) = 2kB T R
43
(3.1)
R
I(t) =
E(t)
E(t)
R
(a)
G = 1/R
(b)
F IG . 3.1 quivalents de T HVENIN (a) et de N ORTON (b) correspondant au bruit thermique dans
une rsistance.
o kB = 1, 38 1023 [J/K] est la constante de B OLTZMANN.
Lquivalent de N ORTON fournit quant lui
I ( f ) = 2kB T G
(3.2)
Dpendance frquentielle
Ltude plus dtaille du phnomne de bruit thermique fait apparatre une dpendance frquentielle par le biais dun facteur p( f )
E ( f ) = 2kB T Rp( f )
(3.3)
1
hf
(3.4)
h
f
/k
BT 1
kB T e
est le facteur de P LANCK [36]. Dans la gamme de frquences considres, il est raisonnable de supposer que la densit spectrale de puissance est constante.
p( f ) =
E ( f ) = 2kB T Re (Zs ( f ))
(3.5)
(3.6)
(3.7)
44
Zs
E
ZL
et
b
b = EZL
V
Zs + ZL
Ib =
Ps f =
Eb
Zs + ZL
(3.9)
(3.10)
Eb2 Re (ZL )
(3.11)
2 kZs + ZL k2
On sait que la puissance fournie par le diple est maximale ladaptation conjugue. Cette puissance, note Psd , est appele puissance disponible du diple ; elle vaut
Psd =
Eb2
8Re (Zs )
(3.12)
Z T
V (t)I(t)dt
(3.13)
bd ( f ) =
E ( f )
4Re (Zs )
45
(3.14)
Remarquons que cette formulation est aussi bien valable pour les signaux utiles dlivrs par la
source que pour le bruit parasite introduit par cette mme source. Au dnominateur, on retrouve
un facteur 4 et non 8 comme dans lexpression 3.12. Cela sexplique simplement par le fait que la
dfinition de la puissance disponible Psd sappuie sur lexpression de la puissance fournie Ps f , qui elle
fait intervenir les valeurs de crte et donc un facteur 12 supplmentaire.
Exemple. Dans le cas particulier du bruit thermique
bd ( f ) =
E ( f )
2kB T Re (Zs ) kB T
=
=
4Re (Zs )
4Re (Zs )
2
(3.15)
La densit spectrale de bruit disponible dpend donc de la temprature mais pas de la frquence ; cest
donc un bruit blanc de densit spectrale uniforme kB2T .
Rsum
Le tableau 3.1 reprend les principaux rsultats noncs dans cette section.
Signaux sinusodaux
Signaux stochastiques
Bruit thermique
Charge quelconque
2
1
b Ib = Eb Re(ZL )2
Ps f = 2 Re V
Pb f = limT + T1
RT
0
2kZs +ZL k
V (t)I(t)dt
Adaptation conjugue
Eb2
8Re(Zs )
E ( f )
bd ( f ) = 4Re(Z
s)
kB T
bd ( f ) = 2
Psd =
type de T .
bd ( f ) =
46
Temprature de bruit
Si lon a affaire un diple dont la densit spectrale est concentre autour dune position frquentielle centrale, on peut dfinir la bande passante de bruit W . La valeur maximale de T ( f ), note T , est
appele temprature de bruit du diple et la puissance de bruit disponible sur tout laxe des frquences
est donn par
Pbd ( f ) = kB TW
(3.17)
ZS
quadriple
E
entre
linaire
sortie
V
2
1
F IG . 3.3 Schma dun quadriple.
Ce quadriple peut tre un amplificateur ou un attnuateur. Il peut tre purement rsistif ou non.
2
La puissance de la porteuse est souvent note C. Ds lors, le rapport signal bruit est plutt un rapport porteuse
bruit ; on le note CN .
47
Rs
Rin = Rs
Quadriple
linaire
4kB T RsW
bd1 ( f ) = 12 kB T0
F IG . 3.4 Schma dun quadriple bruit.
48
Si lon considre un diple gnrateur port la temprature T0 en entre, on peut crire la densit
spectrale disponible la sortie sous la forme
1
bd2 ( f ) = bd1 ( f ) |T =T0 G( f ) + bdq ( f ) = kB T0 G( f ) + bdq ( f )
2
(3.19)
Ds lors,
F0 ( f ) =
bdq ( f )
bd2 ( f )
= 1+ 1
bd1 ( f ) |T =T0 G( f )
2 kB T0 G( f )
(3.20)
Le facteur de bruit est videmment suprieur ou gal lunit. Il fait lobjet de mesure et caractrise
le quadriple.
Interprtation de la notion de facteur de bruit
Soient in ( f ) et out ( f ) les densits spectrales du signal utile respectivement lentre et la
sortie du quadriple. Le rapport du signal bruit lentre vaut donc
(f)
S
= in
N in bd1 ( f )
(3.21)
S
(f)
= out
N out bd2 ( f )
(3.22)
(3.23)
Ds lors,
S
N in
S
N out
( f ) bd2 ( f )
= in
bd1 ( f ) out ( f )
Mais, comme le signal est amplifi par un gain G( f ), cest--dire que out ( f ) = G( f )in ( f ) , ce
rapport devient
S
bd2 ( f )
N in
=
= F0 ( f )
(3.24)
S
bd1 ( f )G( f )
N out
En dautres termes, le facteur de bruit nest jamais que le quotient du rapport signal bruit en entre
par le rapport signal bruit en sortie. Il reprsente donc la dgradation du rapport signal bruit, aprs
passage dans un quadriple.
Circuits quivalents
partir du facteur de bruit, on peut aussi calculer la puissance de bruit disponible en sortie
(3.25)
Cette formulation permet dexprimer le bruit inject par le quadriple par un bruit quivalent qui
serait plac lentre. Le dessin de la figure 3.5 montre les deux circuits quivalents.
49
bd ( f )
bd1 ( f )
bd2 ( f ) = G( f )bd1 ( f ) + bd ( f )
(F0 ( f ) 1)bd1 ( f )
bd1 ( f )
R + 1
2 kB T0 G( f )F0 ( f )d f
R + 1
2 kB T0 G( f )d f
R +
G( f )F0( f )d f
R +
G( f )d f
(3.27)
Facteur de mrite
Il sagit dune grandeur qui nest pas normalise, contrairement au facteur de bruit.
Si la temprature de bruit du diple tait Ts 6= T0 , on pourrait dfinir un facteur de bruit F( f ), en
reprenant point par point le raisonnement prcdent. En galant la densit spectrale de bruit propre au
quadriple
1
bdq ( f ) = (F0 ( f ) 1) kB T0 G( f )
2
1
= (F( f ) 1) kB Ts G( f )
2
50
(3.28)
(3.29)
Comme le gain nest a priori pas fonction de la temprature, pour toute frquence f ,
(F 1)Ts = (F0 1)T0
(3.30)
do
T0
(F0 1)
(3.31)
Ts
et une relation identique entre F0m et Fm . F, appel facteur de mrite, exprime la dgradation du
rapport signal bruit lorsque la temprature du diple gnrateur est Ts .
F = 1+
bd2 ( f ) =
(3.32)
(3.33)
(3.35)
(3.37)
(3.38)
(3.39)
Ds lors,
F0 = 1 +
(L 1)T0
T0
(3.40)
soit encore
(L 1)T0
=L
(3.41)
T0
pour un attnuateur temprature ambiante T0 . Autrement dit, la facteur de mrite est gal au facteur
dattnuation.
F0 = 1 +
(F02 1)bd1 ( f )
F01
bd1 ( f )
G1
F02
G1 F01 bd1 ( f )
G2
(3.42)
Ce bruit passe travers le second quadriple. la sortie de ce dernier, on a une densit spectrale de
bruit qui est la somme du bruit dentre multipli par le gain et le bruit propre
(3.43)
Ds lors
G1 G2 F01 + G2 (F02 1)
F02 1
= F01 +
G1 G2
G1
Plus gnralement, pour un quadriple n tages,
F0 =
F0 = F01 +
De mme,
n
F02 1 F03 1
F0i 1
+
+ = F01 + n1
G1
G1 G2
i=2 j=1 G j
n
Te2
Te3
Tei
Te = Te1 +
+
+ = Te1 + n1
G1 G1 G2
i=2 j=1 G j
(3.44)
(3.45)
(3.46)
On voit donc quil est prfrable de placer lamplificateur de gain important et de facteur de bruit
faible en tte de chane.
52
Deuxime partie
Aspects systmes des chanes
tlcommunications
53
Chapitre 4
Modulations numriques
4.1 Introduction
linstar de la modulation analogique, la modulation numrique est une opration qui transpose
le spectre dun signal numrique pour lamener autour dune frquence porteuse.
Il existe de nombreux critres guidant le choix dun type de modulation. Il va de soi que limportance relative dun critre dpend de lapplication envisage.
Les critres de comparaison sont classs en trois catgories principales (daprs [3])
1. la rsistance aux distorsions et aux interfrences ; cette classe comporte les critres de
(a) la rsistance au bruit en terme de probabilit derreur, celle-ci tant gnralement une
fonction du rapport nergie bruit Eb /N0 ,
(b) la sensibilit aux interfrences dues des multitrajets,
(c) la sensibilit aux imperfections des filtres qui produit de linterfrence entre les symboles
numriques,
(d) la sensibilit aux non-linarits,
2. loccupation spectrale caractrise par
(a) lefficacit spectrale exprime en (bit/seconde) par H ERTZ [b/s/Hz], qui reprsente le
dbit binaire que lon peut transmettre dans un canal large de 1 [Hz] pour un type de
modulation,
(b) le comportement asymptotique de la densit spectrale de puissance, cest--dire la rapidit
de dcroissance de la courbe de densit spectrale de puissance en fonction de la frquence,
3. la simplicit dimplmentation.
Dans un premier temps, nous aborderons les modulations numriques au moyen dun formalisme
gnral. Nous traiterons ensuite quelques exemples communs de modulations numriques linaires
comme la modulation de phase deux tats BPSK, la modulation de phase en quadrature QPSK ou la
modulation MSK.
j(2 fct+c )
s(t) = Re [m(t)] e
(4.1)
o fc est la frquence porteuse du signal modul et c une phase constante. [m(t)] est une fonction
du signal m(t) et constitue lenveloppe complexe es (t) du signal modul. tant donne la forme de
lquation (4.1), le signal s(t) est rel et passe-bande autour de la frquence fc .
[m(t)] = e j [m(t)]
(4.2)
(4.3)
qui met en vidence la composante en phase I [m(t)] et en quadrature Q [m(t)] du signal modul, et
s(t) = k [m(t)]k cos (2 fct + c + arg [m(t)])
(4.4)
j(2 fc t+c )
dk (t) e j(k 2 fc kT )
(4.5)
k=
o les signaux dk (t) contiennent linformation transmettre et k est une phase constante. Deux types
de modulation linaire seront dtailles :
les modulations classiques, pour lesquelles k = 2 fc kT , et
les modulations dcalage (ou offset), pour lesquelles k = 2 fc kT + k 2 .
es (t) =
dk (t)
(4.7)
Dk gk (t kT )
(4.8)
k=
+
k=
56
Le signal gk (t) est un signal de mise en forme rel (non complexe). Pour la simplicit, nous
choisirons une onde de mise en forme unique gk (t) = g(t), k. Dk est une variable alatoire complexe
qui contient linformation numrique transmettre. Elle prendra gnralement la forme Dk = Ak + jBk
o Ak et Bk sont deux variables alatoires relles.
Lenveloppe complexe sexprime galement par es (t) = sI (t) + jsQ (t) o les signaux rels sI (t) et
sQ (t) reprsentent respectivement la composante en phase et en quadrature
+
sI (t) =
Ak g(t kT )
(4.9)
Bk g(t kT )
(4.10)
k=
+
sQ (t) =
k=
s(t) =
k=
k=
(4.11)
(4.12)
En toute gnralit donc, le signal modul peut tre vu comme la modulation en quadrature de deux
signaux numriques en bande de base (de type NRZ).
o M(t) est un processus stochastique complexe. Remarquons tout dabord que le processus stochastique S(t) nest pas stationnaire au sens large vu car sa moyenne dpend du temps. Il est donc ncessaire de stationnariser le signal. Pour cela, nous ajoutons une phase alatoire uniformment rpartie
sur [0, 2 [
S(t) = Re M(t) e j(2 fct+)
(4.14)
Comme le processus stochastique M(t) et la variable alatoire sont non-corrls, il est ais de
montrer que la moyenne du processus stochastique S(t) est nulle et donc indpendante du temps.
57
Fonction dautocorrlation.
o nous avons utilis le fait que le signal modul est rel. Comme
i
1h
M(t) e j(2 fc t+) + M (t) e j(2 fct+)
S(t) =
2
(4.15)
(4.16)
(4.17)
(4.18)
(4.19)
(4.20)
(4.21)
Densit spectrale de puissance. Ds lors, il vient que, si lenveloppe complexe M(t) est un processus stochastique stationnaire au sens large, le signal modul S(t) est galement stationnaire au sens
large. Il est alors possible de calculer sa densit spectrale de puissance en prenant la transforme de
F OURIER de sa fonction dautocorrlation. tant donn que
i
1h
j2 fc
j2 fc
MM ( ) e
+ MM ( ) e
(4.22)
SS ( ) =
4
il vient finalement
( f f )
M ( f f c ) + M
c
4
o M ( f ) est la densit spectrale de puissance de lenveloppe complexe M(t).
S ( f ) =
(4.23)
M(t) =
k=
Dk g(t kT )
(4.24)
fait, lesprance des deux termes faisant intervenir la variable alatoire est nulle.
58
D ( f T ) =
DD (l) e j2 l f T
(4.26)
l=
Si on ajoute lhypothse de variables alatoires Dk non-corrles2 (CDD (k, k l) = 0), on montre que
"
#
+
1
kG ( f )k2
m
2
2
M ( f ) =
D + kD k
f
(4.27)
T
T
m= T
La densit spectrale de puissance de lenveloppe complexe M(t) est donc relle et symtrique.
La densit spectrale en conclusion
En consquence des relations 4.23 et 4.27, la densit spectrale de puissance du signal modul est
fournie par
M ( f f c ) + M ( f + f c )
S ( f ) =
(4.28)
4
La seule connaissance de lenveloppe complexe et de la frquence porteuse permet donc de dterminer
la densit spectrale dun signal numrique classique.
es (t) =
k=
Ak g(t kT )
(4.29)
(4.30)
Pour arriver interprter le signal modul, on dtermine lenveloppe a(t) et la phase (t) du signal
modul. Ces deux signaux sobtiennent aisment partir de lenveloppe complexe par la relation
es (t) = a(t) e j (t)
(4.31)
tant donn que la variable alatoire Ak peut scrire Ak = kAk k e(1sgn(Ak )) 2 j , nous pouvons dduire
+
a(t) =
(t) =
k=
+
kAk k rect[0,T ] (t kT )
(4.32)
(4.33)
La non-corrlation des symboles successifs garantit un transfert dinformation maximum. Mais en contrepartie de
labsence de redondance, linformation est trs sensible des effets de bruit.
59
Il apparat que lenveloppe du signal varie par paliers, ce qui est un inconvnient si le signal passe
au travers dun canal affectant lamplitude ou dun systme non-linaire, comme par exemple un
amplificateur travaillant saturation. Quant la phase, elle marque des discontinuits par sauts de
. Ces sauts de phase sont galement critiques car il ncessite des filtres de rception large bande
passante sous peine de dtriorer le signal modul.
Cas particulier : modulation damplitude numrique deux tats (ASK-2)
Dans le cas de la modulation ASK-2, la variable alatoire Ak peut prendre deux valeurs : +A ou A.
Pour cette modulation, on voit que lenveloppe du signal est constante. Par contre, la phase subit
toujours des sauts de radians toutes les T secondes.
Il est courant de reprsenter une modulation par un diagramme de constellations ou diagramme
des tats de phase, qui nest rien dautre que le plan complexe de la variable alatoire complexe Dk , dans lequel on reprsente toutes les valeurs que peut prendre Dk . On attribue le signal
g(t) cos(2 fct + c ), correspondant Ak , laxe rel pour faire rfrence la composante en phase
et g(t) sin(2 fct + c ), correspondant Bk , laxe imaginaire pour faire rfrence la composante
en quadrature.
La prsence de g(t) signifie que le diagramme de constellations est indpendant du signal de mise
en forme. La figure 4.1 montre le diagramme de constellations pour la modulation ASK-2.
g(t) sin(2 fct + c )
(A, 0)
(A, 0)
g(t) cos(2 fct + c )
(4.34)
Il en rsulte une densit spectrale de puissance pour lenveloppe complexe donne par
es ( f ) = A2 T sinc2 ( f T )
60
(4.35)
s ( f ) =
A2 T
sinc2 [( f fc ) T ] + sinc2 [( f + fc ) T ]
4
(4.36)
Comme la variable alatoire Ak ne peut prendre que deux valeurs dans le cas de la modulation
ASK-2, le signal modul ne transmet quun seul bit dinformation pendant la dure T . Il vient donc
que la densit spectrale de puissance dun signal ASK-2 est donne par la relation 4.36 pour T = Tb
o Rb = 1/Tb est le dbit binaire, encore appel frquence bit, exprim en [b/s].
s(t) = A
k=
(4.37)
o k est une variable alatoire constante sur lintervalle de temps [kT, (k + 1)T [, pouvant prendre N
valeurs possibles D
2
k = 0 + i , i = 0, ..., N 1
(4.38)
N
Le signal modul a donc une enveloppe constante et subit un saut de phase toutes les T secondes.
Cherchons prsent faire le lien avec les modulations linaires classiques. Pour cela, nous
rcrivons s(t) sous la forme
+
s(t) = A
=
"
k=
+
k=
"
k=
(4.39)
(4.40)
= A
k=
(4.41)
(4.42)
Nous retrouvons donc bien lexpression gnrale de lenveloppe complexe dun signal numrique
linaire classique pour lequel
Dk = A e jk
g(t) = rect[0,T ] (t)
61
(4.43)
(4.44)
a(t) = A
rect[0,T ] (t kT )
(4.45)
k rect[0,T ] (t kT )
(4.46)
k=
+
(t) =
k=
Nous retrouvons bien videmment que lenveloppe du signal modul est constante et que la phase
varie par paliers de manire discontinue.
Cas particulier : modulation PSK-2 ou BPSK
La modulation BPSK est une modulation
j0 2jtats
correspondant k {0, }. La variable alatoire
Dk peut donc prendre les valeurs A e , A e , ce qui correspond {(A, 0) , (A, 0)}. La modulation BPSK correspond ainsi exactement la modulation ASK-2. Son plan de constellations est donc
galement celui de la figure 4.1.
Densit spectrale de puissance. Les modulations ASK-2 et BPSK tant identiques, elles ont mme
densit spectrale de puissance (T = Tb ). Elle est exprime par
s ( f ) =
A2 Tb
sinc2 [( f fc ) Tb ] + sinc2 [( f + fc ) Tb ]
4
(4.47)
3
Dk Ae j 4 , Ae j 4 , Ae j 4 , Ae j 4
(4.48)
Le diagramme de constellations est reprsent la figure 4.2.
Considrons une source binaire fournissant le train dimpulsions suivant
+
I(t) =
k=
Ik (t kTb )
(4.49)
3 On
(4.50)
(4.51)
utilise parfois o
lensemble de ces valeurs dcales de 4 , ce qui conduit lensemble des phases possibles suivant :
A, Ae j 2 , A, Ae j 2 .
62
(+ A2 , + A2 )
01
11
g(t) cos(2 fct + c )
00
10
( A2 , A2 )
(+ A2 , A2 )
(4.52)
T = 2Tb
sI (t)
1
0
+ A
1
0
+1
t
A
sQ (t)
+ A
(t)
3
4
3
4
3
4
3
4
F IG . 4.3 Formation des composantes en phase et en quadrature pour la modulation QPSK (avec
mise en forme par un signal rectangulaire).
63
Nous allons prsent dterminer lenveloppe et la phase instantane pour une impulsion de mise
en forme rectangulaire. Lenveloppe complexe est donne par
es (t) = sI (t) + j sQ(t)
(4.53)
k=
(Ak + j Bk ) rect[0,T ] (t kT )
(4.54)
A +
= (I2k + j I2k+1 ) rect[0,T ] (t kT )
2 k=
(4.55)
Il vient alors
q
s2I (t) + s2Q(t)
q
A +
I 2 + I 2 rect[0,T ] (t kT )
=
2 k= 2k 2k+1
a(t) =
(4.56)
(4.57)
= A
k=
et
(t) =
rect[0,T ] (t kT )
k=
(4.58)
rect[0,T ] (t kT ) tan
I2k+1
I2k
(4.59)
Nous voyons donc que lenveloppe du signal modul est constante et gale A. Le rapport I2k+1 /I2k
peut prendre les valeurs 1, ce qui correspondrait 2 tats seulement. Cependant, la fonction tan1
est dfinie radians prs. En fait, pour dterminer la phase, il faut se rfrer au plan de constellations
(figure 4.2).
La phase varie donc par paliers de manire discontinue. La figure 4.3 montre la valeur instantane
de (t) pour une squence binaire donne. On peut remarquer que la phase peut subir des sauts de
2 ou de toutes les T secondes. Ces diffrentes transitions sont reprsentes en pointills sur
la figure 4.2. Lors des transitions , on voit que lenveloppe instantane du signal peut sannuler
pendant le temps de la transition (passage par zro dans le plan des constellations). En effet, des transitions instantanes sont irralisables physiquement. Si le signal est filtr, il peut perdre sa proprit
denveloppe constante, ce qui est gnant lorsque quil passe dans un systme prsentant des distorsions damplitude. Enfin, des sauts brusques de phase impliquent galement une bande passante plus
importante. La figure 4.4 montre un signal modul QPSK pour une squence binaire donne.
Modulateur QPSK
La figure 4.5 montre le modulateur QPSK. Les diffrents signaux apparaissant sur cette figure
sont galement prsents la figure 4.4.
Dmodulateur QPSK
Le signal modul QPSK correspond la modulation en quadrature de deux signaux numriques
en bande de base. La premire phase de la dmodulation consiste donc retrouver les composantes
en phase et en quadrature du signal modul. Les deux signaux numriques en bande de base obtenus
sont ensuite dmoduls grce au filtre adapt (un pour chaque signal) et recombins pour former la
squence binaire de dpart. Le schma du dmodulateur QPSK est repris la figure 4.6.
64
Modulation QPSK
(a)
1
0
1
(b)
10
12
10
12
10
12
10
12
10
12
10
12
0
1
(c)
1
0
1
(d)
1
0
1
(e)
1
0
1
(f)
1
0
1
F IG . 4.4 Illustration de la modulation QPSK : (a) squence binaire I(t), (b) signal sI (t), (c) signal
sQ (t), (d) sI (t) cos (2 fct), (e) sQ (t) sin (2 fct) et (f) signal modul s(t).
sI (t)
cos (2 fc )t
srie
I(t)
parallle
2
sin (2 fc)t
sQ (t)
65
s(t)
1
2 sI (t)
cos (2 fc )t
parallle
1
2 I(t)
s(t)
2
srie
sin (2 fc)t
1
2 sQ (t)
es (t) =
k=
Dk g(t kT )
(4.60)
(4.61)
(4.62)
es ( f ) = 2A2 Tb sinc2 (2 f Tb )
(4.63)
La densit spectrale de puissance dun signal modul en QPSK sexprime finalement par
A2 Tb
sinc2 [( f fc ) 2Tb ] + sinc2 [( f + fc ) 2Tb ]
s ( f ) =
2
et est identique celle de la modulation BPSK pour un dbit binaire double.
(4.64)
es (t) =
k=
Ak g(t kT ) e jk 2
66
(4.65)
Bien que la variable alatoire Ak soit relle, lenveloppe complexe es (t) est bel et bien un signal
valeurs complexes.
Mettons-nous prsent dans le vas o T = Tb , la dure de limpulsion de mise en forme g(t)
ntant pas encore fixe. Le signal modul prend la forme
+
s(t) = Ak g(t kTb ) cos 2 fct + c + k
(4.66)
2
k=
#
"
+
=
(4.67)
Ak g(t kTb ) cos k 2 cos (2 fct + c)
k=
"
#
+
Ak g(t kTb ) sin k
sin (2 fct + c )
(4.68)
2
k=
En tenant compte du fait que cos (k /2) = 0 pour k impair et que sin (k /2) = 0 pour k pair, les
composantes en phase et en quadrature peuvent sexprimer par
+
(4.69)
sI (t) = Ak g(t kTb ) cos k
2
k=
+
k=
(4.70)
et
A
g(t
kT
)
sin
k
b
k
2
k=
+
sQ (t) =
k=
(4.71)
(4.72)
On voit que, par rapport aux modulations linaires classiques, les squences sI (t) et sQ (t) sont dcales entre elles de la dure dun bit Tb , do le nom de modulation dcalage.
e jk 2 . Afin de rsoudre ce problme, nous allons introduire une nouvelle frquence porteuse fc et une
nouvelle enveloppe complexe v(t) pour laquelle le calcul de la densit spectrale de puissance sera
possible :
s(t) = Re es (t) e j(2 fct+c )
(4.73)
j 2 fc + 4T1 t+c
j2 4Tt
b
b e
(4.74)
= Re es (t) e
j2 4Tt
v(t) = e
fc = fc +
67
1
4Tb
es (t)
(4.76)
(4.77)
Si nous connaissons la densit spectrale de puissance du signal complexe v(t), il est possible de
calculer la densit spectrale de puissance du signal modul par le biais de la formule
1
f f 1
+
c
v
c
v
4Tb
4Tb
v ( f fc ) + v ( f fc )
=
(4.78)
s ( f ) =
4
4
que nous avons dmontre prcdemment. Nous allons prsent montrer que v(t) peut se mettre sous
la forme (4.24) pour laquelle nous sommes mme de calculer la densit spectrale de puissance. En
fait,
+
v(t) =
(4.79)
j 2T (tkTb )
(4.80)
Ak g(t kTb ) e
Ak h(t kTb )
k=
+
k=
+
j2 4Tt
Ak g(t kTb ) e jk 2 e
k=
(4.81)
t
j 2T
h(t) = g(t) e
(4.82)
Le signal h(t) correspond un nouveau signal de mise en forme dont la transforme de F OURIER se
dduit du signal g(t)
1
H (f) = G f +
(4.83)
4Tb
La densit spectrale de puissance du signal v(t) sobtient alors en adaptant la formule (4.27)
"
#
+
kH ( f )k2
m
1
v ( f ) =
A2 + A2
f
Tb
Tb
m= Tb
(4.84)
Nous avons prsent notre disposition tous les outils ncessaires au calcul de la densit spectrale
de puissance de signaux numriques linaires dcalage.
I(t) =
k=
Ik (t kTb )
(4.85)
(4.86)
+
A +
sQ (t) = I2k+1 g (t (2k + 1) Tb ) = A2k+1 (1)k g (t (2k + 1) Tb )
2 k=
k=
o g(t) est une impulsion de mise en forme de dure T = 2Tb , A2k = (1)k I2k A2 et A2k+1 =
(1)k I2k+1 A2 , k. Ces deux squences correspondent respectivement aux bits pairs et impairs de
la squence de dpart. La composante en phase sI (t) est la mme pour la modulation OQPSK et la
modulation QPSK. Par contre, la composante en quadrature SQ (t) est dcale de Tb par rapport la
modulation QPSK. La figure 4.7 illustre la formation de ces squences pour un signal de mise en
forme gal
g(t) = rect[0,T ] (t) = rect[0,2Tb ] (t)
(4.87)
qui correspond modulation OQPSK classique.
Tb
I(t)
+1
T = 2Tb
sI (t)
+ A2
A2
Tb
sQ (t)
(t)
3 3
4
4 4
t
A2
3
4
3 3 3
4
4
4 4
3
4
+ A2
3
4
F IG . 4.7 Formation des composantes en phase et en quadrature pour la modulation OQPSK (avec
une mise en forme par un signal rectangulaire).
Il est maintenant nettement moins vident de calculer lenveloppe et la phase instantane du signal
modul, cela cause du chevauchement entre les impulsions de sI (t) et sQ (t). En toute gnralit,
nous pouvons crire
q
a(t) =
sQ (t)
(t) = arctan
sI (t)
69
(4.88)
(4.89)
pour la phase instantane. En observant la figure 4.7, nous pouvons dduire que lenveloppe du signal
modul est constante, sauf peut-tre aux transitions des signaux sI (t) et sQ (t), cest--dire toutes les Tb
secondes. La phase, quant elle, varie toujours par paliers ; pour sen convaincre, il suffit nouveau
de se rfrer au plan de constellations de la figure 4.2. Deux diffrences essentielles sont cependant
remarquer. Tout dabord, la phase ne change plus de valeur toutes les T = 2Tb secondes comme dans
la modulation QPSK mais bien toutes les Tb secondes. De plus, il ny a plus, chaque transition de la
phase, quun saut de 2 . Le fait de dcaler la squence sI (t) a fait disparatre les transitions de .
En effet, dans la modulation QPSK, les transitions de taient dues aux transitions simultanes
des signaux sI (t) et sQ (t). Dans la modulation OQPSK, ces deux signaux ne varient jamais en mme
temps. Les transitions possibles sont reprsentes en pointills sur la figure 4.8. Celles-ci ne se faisant
plus quhorizontalement ou verticalement, lenveloppe du signal modul ne peut donc plus sannuler,
ce qui est un gros avantage de la modulation OQPSK par rapport la modulation QPSK.
La figure 4.9 montre un signal modul OQPSK pour une squence binaire donne.
g(t) sin(2 fct + c )
(+ A2 , + A2 )
( A2 , + A2 )
01
11
g(t) cos(2 fct + c )
00
10
( A2 , A2 )
(+ A2 , A2 )
es (t) =
k=
Ak g(t kTb ) e jk 2
(4.90)
(4.91)
(4.92)
Modulation OQPSK
(a)
1
0
1
(b)
10
11
12
10
11
12
10
11
12
10
11
12
10
11
12
10
11
12
0
1
(c)
1
0
1
(d)
1
0
1
(e)
1
0
1
(f)
1
0
1
F IG . 4.9 Illustration de la modulation OQPSK : (a) squence binaire I(t), (b) sI (t), (c) sQ (t),
(d) sI (t) cos (2 fct), (e) sQ (t) sin (2 fct) et (f) signal modul s(t).
Nous faisons nouveau lhypothse que les deux valeurs possibles pour Ak sont quiprobables. La
2
moyenne A est donc nulle et la variance est gale A2 = E A2k = A2 .Il vient aprs calcul de la
v ( f ) = 2A Tb sinc
2
f + 4T1b ,
1
f+
2Tb
4Tb
(4.93)
Conformment la relation (4.78), la densit spectrale de puissance dun signal modul en OQPSK
sexprime donc finalement par
s ( f ) =
A2 Tb
sinc2 [( f fc ) 2Tb ] + sinc2 [( f + fc ) 2Tb ]
2
(4.94)
Description
La caractristique qui diffrencie la modulation MSK de la modulation OQPSK est le signal de
mise en forme qui prend maintenant la forme
t
(4.95)
g(t) = rect[0,2Tb ] (t) sin
2Tb
Les squences sI (t) et sQ (t) sont formes exactement de la mme manire que pour la modulation
OQPSK. La composante en phase sI (t) scrit sous la forme
+
(t 2kTb )
sI (t) = A I2k rect[0,2Tb ] (t 2kTb ) sin
(4.96)
2Tb
k=
+
t
k
(4.97)
= A I2k (1) rect[0,2Tb ] (t 2kTb ) sin
2Tb
k=
+
t
(4.98)
= cos
A2k rect[0,2Tb] (t 2kTb)
2Tb 2 k=
o nous avons utilis le fait que sin(a) = cos(a 2 ) et dfini A2k = A I2k (1)k . Le signal sI (t) peut
tre vu comme un signal du type NRZ multipli par une cosinusode. Un dveloppement similaire
donne pour la composante en quadrature de phase
+
t
(4.99)
sQ (t) = sin
A2k+1 rect[0,2Tb] (t 2kTb)
2Tb 2 k=
o A2k+1 = A I2k+1 (1)k . Ces deux expressions vont nous permettre de dterminer lenveloppe instantane et la variation de phase instantane du signal modul. Lenveloppe est donne par
q
a(t) =
s2I (t) + s2Q(t)
(4.100)
s
t
2
2
2
2
=
+ A sin
(4.101)
A cos
2Tb 2
2Tb 2
= A
(4.102)
Lenveloppe du signal modul est donc constante et gale A. La phase instantane peut sexprimer
par
1 sQ (t)
(t) = tan
(4.103)
sI (t)
#
"
+
A
rect
(t
(2k
+
1)T
)
2k+1
b
[0,2T
]
k=
b
(4.104)
= tan1 tan
2Tb 2
A
rect
+
[0,2Tb ] (t 2kTb )
k= 2k
Si maintenant, on observe lvolution de la phase pendant un intervalle de temps Tb dun bit, la fraction
dans largument de la fonction tan1 vaut 1 et la variation de la phase vaut
(t) =
t
2Tb
(4.105)
La phase varie donc linairement de 2 sur la dure de Tb secondes. Contrairement aux modulations
vues jusquici, la modulation MSK prsente des sauts de phase qui se ralisent en douceur. Par
72
01
11
g(t) cos(2 fct + c )
00
10
Tb
/2
t
/2
3 /2
(4.107)
(4.108)
qui reprsente bien une modulation de frquence dont lexcursion est gale f = 2 4T1b =
1
2Tb .
Modulation MSK
(a)
1
0
1
(b)
10
11
12
10
11
12
10
11
12
10
11
12
10
11
12
10
11
12
0
1
(c)
1
0
1
(d)
1
0
1
(e)
1
0
1
(f)
1
0
1
F IG . 4.12 Illustration de la modulation MSK : (a) squence binaire I(t), (b) sI (t), (c) sQ (t),
(d) sI (t) cos (2 fct), (e) sQ (t) sin (2 fct) et (f) signal modul s(t).
es (t) =
k=
Ak g(t kTb ) e jk 2
(4.109)
(4.110)
(4.111)
f + 4T1b ,
i 2
h
1
16 A2 Tb cos 2 f + 4Tb Tb
v ( f ) =
2
2
1 16 f + 1
Tb2
4Tb
74
(4.112)
La densit spectrale de puissance dun signal modul en MSK sexprime donc finalement par
!2
!2
2
4 A Tb
cos [2 ( f + fc ) Tb ]
cos [2 ( f fc ) Tb ]
+
s ( f ) =
(4.113)
1 16 ( f fc )2 T 2
1 16 ( f + fc )2 T 2
b
BPSK
MSK
10
20
OQPSK
dB
30
40
50
60
70
2.5
1.5
0.5
0.5
1.5
2.5
f fc
75
76
Chapitre 5
Modlisation du canal pour transmissions
numriques et interfrence inter-symboles
Dans un canal de communication, une source importante derreurs est la prsence de bruit sur
le canal de transmission ; cette question a dj t largement aborde prcdemment. Nous allons
prsent nous intresser une autre source de distorsion, appele interfrence inter-symboles qui
apparat, pour des communications numriques, lorsque le canal de transmission est dispersif (effet
de la bande passante finie du canal).
Prenons le cas dune transmission par mise en forme NRZ. On sait que la densit spectrale du
signal NRZ occupe tout le spectre, malgr que la dcroissance de la norme en fonction de la frquence
soit significative. Or, il est rare quun systme soit bande passante infinie ou que le canal ait une
transmission de canal idal sur la largeur de bande considre. En consquence, le signal NRZ subit
des distorsions subsquentes au passage dans le canal. Ces distorsions se traduisent inluctablement
par un recouvrement entre symboles successifs.
Pour la facilit des dveloppements, nous nous limitons une transmission binaire en bande de
base.
s(t) =
k=
Ak gE (t kTb )
(5.1)
o les amplitudes Ak (reprsentant linformation binaire transmise linstant t = kTb ) sont lies la
squence {bk } par la relation
+1 si bk = 1
(5.2)
Ak =
1 si bk = 0
Le signal transmis au travers du canal est donc s(t). Celui-ci va tre modifi par convolution avec la
rponse impulsionnelle h(t) du canal de transmission. De plus, un bruit (suppos blanc et gaussien)
vient sajouter au signal dans le canal. Le signal reu, not x(t), est ensuite filtr par un filtre de
rception, de rponse impulsionnelle gR (t). La sortie de ce filtre, note y(t), est chantillonne au
rythme de lmetteur, cest--dire toutes les Tb secondes et est synchronis avec celui-ci (les instants
77
Modulation PAM
Horloge
{Ak }
metteur
Filtre de mise en
forme gE (t)
s(t)
Canal de
transmission h(t)
Bruit blanc
gaussien w(t)
Canal
x(t)
Filtre rcepteur
gR (t)
y(t)
chantillonnage
en ti = iTb
Rcepteur
yi (t)
seuil
Organe de
dcision
1 si y(ti ) >
0 si y(ti ) <
78
dchantillonnage sont dtermins par une horloge qui est extraite du signal y(t)). Finalement la
squence dchantillons ainsi obtenue est utilise pour reconstruire la squence originale {bk } au
moyen dun organe de dcision. Lamplitude de chaque chantillon est compare un seuil .
La sortie du filtre de rception peut scrire sous la forme
+
y(t) =
k=
(5.3)
(5.4)
(5.5)
o est un facteur dchelle constant et p(t) est le rsultat du passage du signal de mise en forme
gE (t) au travers des deux filtres h(t) et gR (t) caractrisant respectivement le canal et le filtre de
rception. Le terme n(t) est le rsidu du bruit additif w(t) ajout par le canal, filtr par le filtre de
rception. Il reprsente donc toujours un bruit blanc gaussien car il est pass au travers dun filtre
linaire. En gnral, w(t) est modlis comme un bruit blanc additif gaussien de moyenne nulle. On
obtient finalement
+
y(t) =
k=
(5.6)
Cette formule montre que le train dimpulsions nest plus rectangulaire. Le facteur dchelle est
choisi de telle sorte que le signal p(t) soit normalis, cest--dire
p(0) = 1
(5.7)
Cette relation justifie lutilisation de comme facteur dchelle pour tenir compte des changements
de lamplitude du signal lors de son passage travers le systme. De plus, cette relation reprsente une
contrainte mathmatique et na aucune signification physique. Avant de poursuivre, il est important de
remarquer que, pour tre rigoureux, il faudrait introduire un dlai t0 dans largument de p(t kTb ) afin
de reprsenter leffet du dlai de la transmission au travers du systme. Cependant, pour simplifier
lexpos, nous choisissons de prendre ce dlai gal zro, et cela sans perdre en gnralit.
Vu quune convolution dans le domaine temporel quivaut un produit dans le domaine frquentiel, nous dduisons de la relation 5.4 que
P( f ) = GE ( f )H ( f )GR ( f )
(5.8)
Le signal y(t) la sortie du filtre de rception est chantillonn aux instants ti = iTb (il est intressant de rappeler que lchantillonnage se fait la fin de la priode, ceci rsultant de la thorie de du
filtre adapt), ce qui fournit les chantillons
y(ti ) =
k=
= Ai +
(5.9)
k =
k 6= i
(5.10)
Le premier terme Ai reprsente la contribution du i-me bit transmis et donc le symbole dintrt.
Le second terme reprsente leffet rsiduel de tous les autres bits transmis sur le i-me bit, quils aient
79
t transmis avant ou aprs le i-me bit. Cet effet rsiduel est appel interfrence inter-symboles.
Limportance de cet effet dpend de ltalement de p(t). On peut dj remarquer que pour minimiser
le terme dinterfrence inter-symboles, il faudra choisir une forme adquate pour la fonction p(t) de
manire liminer ce second terme. Le dernier terme de n(ti ) est le bruit n(t) chantillonn linstant
ti .
En labsence dinterfrence inter-symboles et de bruit, lquation 5.10 se rduit
y(ti ) = Ai
(5.11)
ce qui montre bien que, sous ces conditions idales, le i-me bit transmis est dcod correctement. La
prsence invitable dinterfrence inter-symboles et de bruit dans le systme introduit des erreurs au
niveau de lorgane de dcision. Ds lors, lors la conception du filtre de mise en forme et du filtre de
rception, lobjectif est de minimiser les effets de linterfrence inter-symboles et du bruit, aboutissant
ainsi un taux derreurs aussi faible que possible.
Il existe deux types de transmission :
la transmission o le bruit est prpondrant, cest le cas des communications satellite par exemple, pour lequel on privilgie le bon fonctionnement du filtre adapt,
la transmission o le rapport signal sur bruit est important, comme pour le cas du rseau tlphonique par exemple, et ds lors le fonctionnement du systme est essentiellement limit par
linterfrence inter-symboles, et non par le bruit.
Dans ltude qui suit, nous nous placerons dans le second cas et nous ngligeons ainsi le terme n(t)
et nous nous focalisons sur la question de la rduction, voire de llimination de linterfrence intersymboles.
ps (t) =
m=
p(mTb ) (t mTb )
(5.13)
reprsente alors une version chantillonne du signal p(t). Par le thorme dchantillonnage, il est
facile de calculer la transforme de F OURIER du signal ps (t)
+
Ps ( f ) = f b
m=
P( f m fb )
(5.14)
o fb = 1/Tb est le dbit binaire (ou rythme) exprim en [b/s]. On remarque ds lors une rptition
du spectre autour de tous les multiples entiers de la frquence dchantillonnage. La fonction Ps ( f )
peut encore scrire sous la forme
Ps ( f ) =
Z + +
m=
[p(mTb ) (t mTb )] e2 j f t dt
(5.15)
Z +
= 1
p(0) (t)e2 jt f dt
(5.16)
(5.17)
m=
P( f m fb ) = Tb
(5.18)
Nous pouvons prsent noncer le critre de N YQUIST pour une transmission en bande de base idale
en labsence de bruit.
Proposition 22 [Critre de N YQUIST] La transforme de F OURIER P( f ) de limpulsion p(t) limine totalement linterfrence inter-symboles pour des chantillons pris toutes les Tb secondes si elle
vrifie la condition
+
m=
P( f m fb ) = Tb
(5.19)
Rappelons que la transforme P( f ) dpend du systme dans son ensemble ; cela inclut le filtre de
mise en forme, la transmittance du canal et le filtre de rception.
fb
1
=
2
2Tb
(5.22)
La condition est donc bien satisfaite. Cette quation permet de dire quaucune frquence ncessaire
ne dpasse la moiti du dbit binaire. On en dduit quune fonction p(t) qui satisferait le critre de
N YQUIST est, comme il avait t signal prcdemment, le sinus cardinal :
p(t) =
sin(2 Wt)
= sinc(2Wt)
2 W t
(5.23)
Une telle fonction p(t) (ou P( f )) caractrise le canal idal de N YQUIST. Les figures 5.2.(a) et .5.2.(b)
montre les fonctions p(t) et P( f ). p(t) peut ds lors tre vue comme la rponse impulsionnelle dun
filtre passe bas idal de frquence de coupure correspondant W .
la figure 5.2.(b), on voit galement les diffrents instants dchantillonage. Il est clair que,
si le signal reu y(t) est chantillonn aux instants t = 0, Tb , 2Tb , ..., les impulsions p(t iTb )
(i = 0, 1, 2, . . .) ninterfrent pas entre elles. Cette condition est illustre la figure 5.3 pour la
squence binaire 1011010. On constate ds lors quaux instants dchantillonnage, tous les sinus
cardinal relatif aux symboles perturbateurs sont nuls.
Il existe cependant des raisons pour lesquelles limpulsion en sinus cardinal ne peut tre utilise
en pratique :
1. elle ncessite que P( f ) soit constante sur lintervalle de frquence [W, +W ] et nulle partout
ailleurs. Cela est pratiquement irralisable en raison des transitions abruptes en W (un filtre
passe-bas idal nest pas ralisable),
2. la fonction p(t) dcrot en 1/|t| pour |t| lev ; elle dcrot donc trs lentement, les sinus cardinal
se rpercuteront donc sur des chantillons lointains. Il y a ds lors trs peu de marge derreur
acceptable sur les instants dchantillonnage car si on nchantillonne pas tout fait au bon
moment, tous les sinus cardinaux sont non nuls, ce qui donne une erreur cumule importante.
Pour valuer cet effet derreur dchantillonnage, considrons que lon chantillonne le signal
y(t) en t = ti + t, o t est le dcalage temporel. Pour la clart, on suppose que ti = 0. En
labsence de bruit, nous avons donc
y(t) =
=
Ak p(t kTb )
Ak
k=
+
k=
(5.24)
sin[2 W (t kTb )]
2 W (t kTb )
(5.25)
Sachant que
sin[2 W (t kTb )] = sin(2 W t) cos(2 W kTb )
cos(2 W t) sin(2 W kTb )
= (1)k sin(2 W t)
(5.26)
(5.27)
(5.28)
sin(2 W t)
y(t) = A0 sinc(2W t) +
82
k =
k 6= 0
Ak
(1)k
(2W t k)
(5.29)
2W P( f )
1.0
W=
f
W
1
2Tb
(a)
p(t)
1.0
0.5
Instants dchantillonnage
Tb
intervalles du signal
(b)
83
t
Tb
fb
2
Squence binaire
1.0
Amplitude
0.5
0.0
0.5
1.0
4
2
Temps
84
10
Le premier terme correspond au symbole dsir, tandis que la srie reprsente linterfrence
inter-symboles dues lerreur t. Il est malheureusement possible que cette srie diverge,
provoquant des dcisions errones au rcepteur.
1
,
2W
W f W
(5.30)
Il est possible de dterminer un ensemble de fonctions P( f ) qui vrifient cette dernire condition.
Une forme spciale de la fonction P( f ) souvent utilise en pratique est limpulsion en cosinus
surlev dont le spectre est donn par
1
0 | f | < f1
2W
io
n
h
(| f |W )
1
P( f ) =
(5.31)
f1 | f | 2W f1
4W 1 sin 2W 2 f1
0
| f | 2W f1
= 1
f1
W
(5.32)
o est un paramtre, appel facteur de rolloff , qui indique lexcs en bande passante ncessaire par
rapport au canal idal de N YQUIST. La bande passante du systme BT est maintenant gale
BT = 2W f1 = W (1 + )
(5.33)
Les fonctions P( f ) et p(t) sont respectivement reprsentes aux figures 5.4.(a) et 5.4.(b) pour = 0,
= 0, 5 et 1. Pour = 0, limpulsion correspond au canal idal de N YQUIST.
Limpulsion p(t) est simplement la transforme de F OURIER inverse de P( f )
cos(2 W t)
p(t) = sinc(2Wt)
(5.34)
1 16 2W 2t 2
La fonction p(t) est forme par le produit de deux facteurs : le facteur sinc(2W t) qui correspond
au canal idal et une fonction qui dcrot en 1/|t|2 lorsque |t| est grand. Le premier facteur assure
que la fonction p(t) passe par zro aux instants dchantillonnage ti = it tandis que le second facteur
rduit lextension temporelle de limpulsion de telle sorte que lerreur t est considrablement rduite.
Pour = 1, lextension temporelle de limpulsion est minimale. Lerreur t diminue donc lorsque
augmente. Cependant, la bande passante est dautant plus grande que est grand, pour atteindre sa
valeur maximale 2W soit le double de la bande passante ncessaire pour le canal idal de N YQUIST
, lorsque = 1. Un compromis existe donc entre la bande passante est la diminution de linterfrence
inter-symboles. Typiquement, une valeur de = 0.3 convient.
La synthse dun systme de tlcommunications numrique est rendue difficile par le fait que
cest la fonction p(t) qui doit possder une rponse impulsionnelle en forme de cosinus surlev. Or,
p(t) rsulte de la mise en cascade de trois systmes (gE (t), h(t) et gR (t)). La synthse de ces systmes
doit donc tre telle que lensemble doit tre caractris par un cosinus surlev.
85
2W P( f )
=0
= 0.5
=1
0.8
0.6
0.4
0.2
f
W
-2
-1.5
-1
-0.5
0.5
1.5
(a)
p(t)
1.0
=1
0.5
= 0.5
=0
(b)
86
t
Tb
Chapitre 6
talement de spectre
6.1 Introduction
Loccupation spectrale est un souci rcurrent pour la mise au point de techniques de modulation.
Lors dune transmission en bande de base, on cherche donner une forme adquate la densit
spectrale en fonction de lapplication vise. Pour des communications en espace libre, on tente de
trouver le compromis entre loccupation spectrale et la rsistance au bruit.
Un systme talement de spectre se caractrise par le fait que la bande passante B est largement
suprieure au dbit dinformation Rb exprim en [b/s]. Cet accroissement de bande saccompagne,
comme nous le verrons par la suite, dune meilleure rsistance au bruit.
Historiquement, ce sont les militaires qui ont dvelopp en premier les techniques dtalement de
spectre. Il existe nanmoins aujourdhui des systmes de transmission numrique commerciaux qui en
utilisent le principe. Le tableau 6.1 fournit une comparaison des diffrentes fonctions de ltalement
de spectre.
Objectifs
Lutte contre le brouillage
Accs multiple
Dtection difficile
Protection des donnes
Appel slectif
Identification
Navigation
Protection contre les multi-trajets
Faible densit de flux mis
Domaine militaire
Application commerciale
ltalement direct ou Direct Spreading (DS). Cette technique porte aussi le nom Code-Division
Multiplexing (CDM) car elle consiste taler le spectre en multipliant le signal utile par un
signal numrique dont la cadence bit est nettement suprieure appel code.
Dans tout systme de transmission, quil soit partag ou non, il est ncessaire de crer le signal
transmettre mais aussi de dfinir des moyens et conventions pour interprter correctement le signal
la rception. Il faut donc dtailler un moyen daccder au message utile aprs talement. Ces techniques sont reprises gnralement sous le vocable de Code-Division Multiple Access (CDMA), qui
est souvent utilis pour dsigner la technique dtalement direct elle-mme.
Dans la suite, nous analysons la technique de ltalement ; la question de laccs au signal nest
pas aborde dans ce chapitre.
(6.1)
Ce choix permet, entre autres, dassurer une cohrence entre la synchronisation des informations
portes par les deux signaux.
Considrons la figure 6.1. Le premier signal reprsente le signal de donnes binaires transmettre
(ici, deux bits) ainsi que loccupation spectrale associe (idalise). Le deuxime signal reprsente un
exemple de squence dtalement. On constate que son occupation spectrale est dilate par N.
Dans un systme talement direct, le taux de transmission de symboles (1/Tb) dun signal bipolaire binaire est modifi en multipliant cette onde par un code bipolaire binaire pseudo-alatoire priodique dont la dure Tc est nettement plus courte que celle de dpart ; on obtient alors un signal tel
que celui reprsent la troisime ligne de la figure 6.1. Son spectre est approximativement aussi
large que celui de la squence dtalement.
La figure 6.2 donne une reprsentation plus fidle des densits spectrales de puissance du signal
de donnes et de la squence dtalement. Elle met aussi en lumire le thorme de conservation de
la puissance (conservation de laire sous la courbe).
Un raisonnement simple, par labsurde, permet de se convaincre de la conservation de puissance.
En effet, le signal tal rsulte de la multiplication de la squence initiale par le code dtalement.
Remarquons quune deuxime multiplication du signal tal par le code dtalement fournit le signal
de dpart. Ds lors, sil y avait perte ou gain de puissance durant ltalement, lapplication successive
de deux phases dtalement ne fournirait pas le signal de dpart.
Les squences dtalement ont pour but de faire apparatre, pour lutilisateur courant, les autres
utilisateurs comme du bruit. Ces squences doivent ds lors avoir les caractristiques les plus proches
possible de celles dun bruit blanc. Nanmoins, pour que le rcepteur puisse retrouver le signal utile, il
doit aussi tre capable de rgnrer la squence ayant servi ltalement. Il ne peut donc sagir dune
88
Tb
(un bit)
W 0 W
Squence dtalement
F
t
NW
NW
NW
NW
Tc
Signal tal
F
t
-1
T1
1
T
t
+1
Symbole
tal
-1
Tc
f
T1c
1
Tc
Priode du chip
Tc
DOMAINE FREQUENTIEL
DENSITE SPECTRALE DE PUISSANCE
DOMAINE TEMPOREL
vritable squence alatoire mais bien dune squence dterministe ayant des proprits aussi proches
que possibles de celles dune squence alatoire ; on parle donc de squences pseudo-alatoires. De
plus, pour des raisons pratiques de ralisation, les squences sont priodiques.
Pour construire des gnrateurs de squences pseudo-alatoires, on peut partir de ltude des proprits des squences alatoires et imposer algbriquement des proprits similaires aux gnrateurs.
Cest la dmarche que nous suivrons ici. Une remarque simpose nanmoins : il est possible de construire plusieurs types de gnrateurs de squences alatoires. Ltude mene dans ce document ne se
veut donc pas exhaustive. Le lecteur est invit consulter des ouvrages spcialiss pour une typologie
des gnrateurs de squences pseudo-alatoires.
Flip-flop
....
Horloge
Squence
de sortie
Additionneur
modulo-2
Flip-flop
s0
2 s
2
s1
3 s
3
Squence
de sortie
Horloge
F IG . 6.4 Exemple de registre dcalage linaire.
Proprits dune squence de longueur maximale
Considrons des symboles 0 et 1 respectivement reprsents par les niveaux 1 et +1. Les
squences de longueur maximale ont des proprits proches de celles des squences binaires purement alatoires1 , savoir une moyenne nulle et une fonction dauto-corrlation en forme de delta de
D IRAC.
tablissons ces proprits. La priode dune squence de longueur maximale est dfinie par
N = 2m 1
(6.2)
Appelons c(t) londe rsultant de la squence de longueur maximale, comme illustr la figure 6.5(a)
pour N = 7.
La priode de londe c(t) est gale
Tb = NTc
(6.3)
o Tc , appele priode chip, est la dure dun bit de la squence de longueur maximale.
Deux proprits des squences pseudo-alatoires sont immdiates :
1. [Esprance] Pour une squence longueur maximale, tous les tats internes des registres sont
possibles lexception de ltat nul (il y a 7 tats internes possibles dans notre exemple). Si
lon retire ltat avec tous les bits mis 1, pour tout tat interne il y a un dual (les 0 sont 1
et rciproquement). Un tat et son dual fournissent en sortie des bits de signe oppos et donc
ces derniers se compensent. Pour le calcul de lesprance, ll ne reste alors qu considrer ltat
avec tous les bits 1, qui lui fournit 1 en sortie. Ds lors, lesprance vaut
1
(6.4)
N
Contrairement une squence alatoire, lesprance nest pas nulle mais elle diminue avec la
taille du code dtalement.
2. [Priodicit de la fonction dauto-corrlation] La fonction dautocorrlation dune squence
de longueur maximale est priodique.
E {c(t)} =
Par ltablissement de la fonction dautocorrlation, nous partons de la dfinition de la fonction dautocorrlation dun signal priodique c(t) de priode Tb , qui est donne par
1
cc ( ) =
Tb
Z +Tb /2
Tb /2
c(t)c(t + )d
(6.5)
Lutilisation dune squence rellement alatoire ne permettrait pas de restituer le signal utile au rcepteur prcisment en raison de laspect alatoire !
91
o appartient lintervalle [Tb /2, +Tb /2]. En appliquant cette formule la squence de longueur
maximale c(t), on obtient
1 N+1
NTc | | , | | Tc
(6.6)
cc ( ) =
1
N
ailleurs
Ce rsultat est montr la figure 6.5(b) pour m = 3 et N = 7.
Squence binaire
00111010011101
+1
-1
t
Tc
NTc
(a)
cc ( )
1.0
Tc
NTc
Tc
0
NTc
N1
(b)
c ( f )
T2c
T1c
0
1
NTc
1
Tc
2
Tc
(c)
F IG . 6.5 Proprits dune squence de longueur maximale.
92
tant donn quune priodicit dans le domaine temporel correspond un chantillonnage uniforme dans le domaine frquentiel, et la forme de la fonction dautocorrlation de c(t), on peut dterminer la densit spectrale de puissance de c(t)
+
sin( Nn ) 2
1
1+N
n
(6.7)
c ( f ) = 2 ( f ) + 2
f
N
N n=,
Nn
NTc
n6=0
qui est reprsente la figure 6.5(c) pour m = 3 et N = 7. Au vu de ces rsultats, il est maintenant
possible de comparer une squence de longueur maximale et une squence rellement alatoire :
sur une priode de la squence de longueur maximale, la fonction dautocorrlation de c(t) est
similaire celle de la squence alatoire. En effet, la fonction dautocorrlation dune squence
alatoire vaut
(
1 |Tc| | | Tc
(6.8)
0
| | > Tc
Il apparat donc que les densits spectrales des deux squences ont une mme enveloppe,
savoir sinc2 ( f T ).
la diffrence fondamentale rside dans le fait que le spectre de la squence alatoire est continu
tandis que celui de la squence de longueur maximale est compos de raies espaces de 1/NTc .
Lorsque lon augmente la valeur de N, cest--dire le longueur de la squence de longueur maximale,
la squence pseudo-alatoire se rapproche plus dune squence alatoire pure. Cependant, le prix
payer avec une grande valeur de N, est une augmentation de la bande passante ncessaire. Un
compromis doit ds lors tre tabli.
Choix dune squence de longueur maximale
La question qui se pose prsent consiste savoir comment choisir la fonction logique du registre
dcalage pour obtenir une priode dsire N. La littrature fournit des tables avec les connections
raliser entre les lments du circuit pour diffrentes valeurs de m. Une telle table est montre la
figure 6.6.
Registre dcalage
de longueur m
2
3
4
5
6
7
8
F IG . 6.6 Table permettant de construire le registre dcalage adquat [16, page 584].
Plus m augmente, plus il y a de possibilits de circuits.
Exemple. Considrons une squence de longueur maximale ncessitant lutilisation dun registre de
longueur m = 5. Pour cette valeur de m, nous choisissons dans la table une boucle rtroaction du
93
type [5, 2]. La configuration correspondante est montre la figure 6.7(a). En supposant un tat initial
10000, la figure 6.2 montre lvolution du registre pour une priode de la squence de longueur maximale. On remarque quen fin de priode, le registre est revenu son tat initial. On aurait galement
pu choisir la configuration [5, 4, 2, 1] comme le montre la figure 6.7(b) et on aurait une autre squence
pseudo-alatoire de mme longueur.
Additionneur
modulo-2
(a)
Flip-flop
Squence
de sortie
Horloge
(b)
Flip-flop
Horloge
Squence
de sortie
Z +
C ( )B( f ) d
(6.10)
(6.11)
Conformment lusage repris dans la littrature scienfitique, le terme modulation est utilis tantt pour dcrire
ltalement, tantt pour dsigner lopration qui fait intervenir une porteuse AC cos(2 fct).
94
Bit de rtroaction
0
1
0
1
1
1
0
1
1
0
0
0
1
1
1
1
1
0
0
1
1
0
1
0
0
1
0
0
0
0
1
1
0
1
0
1
1
1
0
1
1
0
0
0
1
1
1
1
1
0
0
1
1
0
1
0
0
1
0
0
0
0
1
Registres
0 0 0
1 0 0
0 1 0
1 0 1
0 1 0
1 0 1
1 1 0
1 1 1
0 1 1
1 0 1
1 1 0
0 1 1
0 0 1
0 0 0
1 0 0
1 1 0
1 1 1
1 1 1
1 1 1
0 1 1
0 0 1
1 0 0
1 1 0
0 1 1
1 0 1
0 1 0
0 0 1
1 0 0
0 1 0
0 0 1
0 0 0
0 0 0
0
0
0
0
1
0
1
0
1
1
1
0
1
1
0
0
0
1
1
1
1
1
0
0
1
1
0
1
0
0
1
0
Bit de sortie
0
0
0
0
1
0
1
0
1
1
1
0
1
1
0
0
0
1
1
1
1
1
0
0
1
1
0
1
0
0
1
TAB . 6.2 volution du registre de la figure 6.7(a) au cours dune priode de la squence [16, page
585].
95
{bk }
m(t)
b(t)
NRZ
m(t)
r(t)
c(t)
{ck }
NRZ
i(t)
(b)
(a)
z(t)
r(t)
X
R Tb
0
dt
Organe de
dsision
1 si v > 0
0 si v < 0
Seuil = 0
c(t)
(c)
F IG . 6.8 Modle dune transmission spectre tal en bande de base.
Or, on effectue le produit des deux courbes et on intgre sur toutes les valeurs diffrentes de zro. On
obtient ds lors un signal large bande. De plus, la figure 6.9 montre trs clairement que la priode du
signal tal est de Tc et donc la bande de base du signal tal vaut 2T1c , ce qui reprsente un signal
large bande.
On dit que la squence c(t) joue le rle de code dtalement. La formation du signal m(t) est
montre la figure 6.9.
Pour une transmission en bande de base (sans utilisation dune porteuse), on transmet le signal tal
m(t) sur le canal.
Le rcepteur reoit le signal tal m(t) bruit par des interfrences caractrises par le signal i(t)
(cf. figure 6.8(b))
r(t) = m(t) + i(t)
= c(t)b(t) + i(t)
(6.12)
(6.13)
Le signal dinterfrence se compose de signaux mis par dautres utilisateurs dans la mme bande
de frquences ainsi que du bruit.
Pour recouvrer le signal original et la squence binaire sous-jacente, le signal reu r(t) est appliqu
lentre dun dmodulateur qui consiste en une multiplication suivi dun intgrateur et dun organe
de dcision, comme lillustre la figure 6.8(c). Au niveau du mlangeur, le signal reu r(t) est multipli
par une rplique exacte de la squence pseudo-alatoire utilise au rcepteur.
Nous faisons lhypothse que le rcepteur travaille en synchronisme parfait avec lmetteur. La
sortie de ltage de multiplication est alors donne par
z(t) = c(t)r(t)
= c2 (t)b(t) + c(t)i(t)
(6.14)
(6.15)
Cette dernire quation montre que le signal b(t) est multipli deux fois par la squence pseudoalatoire c(t), tandis que le signal i(t) nest multipli quune seule fois. Vu que c(t) vaut +1 ou 1,
96
0
-1
Tb
(a)
Code dtalement c(t)
+1
0
-1
Tc
NTc
(b)
0
-1
(c)
F IG . 6.9 Formation du signal m(t) par talement du spectre de b(t).
97
le signal c2 (t) est gal 1 pour tout instant t. Donc, le signal z(t) se rduit
z(t) = b(t) + c(t)i(t)
(6.16)
Nous voyons donc que le signal b(t), qui contient linformation utile, se retrouve la sortie de ltage
de multiplication. Si le terme additif c(t)i(t) semble gnant, il faut bien voir que la multiplication de
i(t) par le code dtalement c(t) a pour effet dtaler le spectre de i(t). Ds lors, le signal c(t)i(t) est
large bande tandis que le signal b(t) est bande troite. Ce signal large bande contient notamment
les signaux tals des autres utilisateurs, ceux-ci tant considrs comme du bruit.
En appliquant le signal z(t) lentre dun filtre passe-bas, en loccurrence lintgrateur, la majeure partie de lnergie de c(t)i(t) est filtre. Leffet de linterfrence i(t) est donc fortement rduit
la sortie du rcepteur. En ralit, partir du signal z(t), tout se passe comme pour la dmodulation dun signal numrique par le filtre adapt suivi de lorgane de dcision. Le signal b(t) est le
signal numrique dmoduler et c(t)i(t) peut tre vu comme lapproximation dun bruit blanc additif gaussien. Plus la squence pseudo-alatoire est longue, plus cette approximation est vrifie.
Cependant, le prix payer pour cette augmentation de la protection contre les interfrences est une
augmentation de la bande passante requise, de la complexit du systme, ainsi que des dlais.
En rsum, lutilisation dun code dtalement dans lmetteur produit un signal large bande transmis sur le canal. Celui-ci apparat comme du bruit pour un rcepteur qui na pas connaissance du code
dtalement initial.
Squence de
donnes binaires
{bk }
Encodeur
NRZ
b(t)
m(t)
Modulateur
B-PSK
x(t)
c(t)
Gnrateur
de code PN
Porteuse
(a)
Dtecteur cohrent
Signal reu
y(t)
Modulateur
produit
Filtre
passe-bas
R Tb
0
dt
Organe de
dsision
1 si v > 0
0 si v < 0
Gnrateur
local de code
PN
Porteuse locale
(b)
Analyse
La figure 6.10 reprsente un modulateur o ltalement de spectre se fait avant la modulation de
phase numrique. Pour une analyse plus dtaille, il est intressant dintervertir conceptuellement
ces deux tapes de modulation comme le montre le schma de la figure 6.12. Cette permutation est
licite car les oprations dtalement de spectre et de modulation BPSK sont des multiplications, et
donc commutatives. Il en va de mme pour les oprations de dmodulation. Le fait de permuter les
deux oprations permet de modifier la rponse impulsionnelle du filtre adapt qui, pour rappel, est
donne par
h(t) = g(Tb t)
(6.17)
Combinant les deux oprations, g(t) devient cos(2 fct). La rponse impulsionnelle du filtre adapt
sera donc un cosinus retourn.
Les interfrences sont modlises par le signal j(t) qui est de type passe-bande. La sortie du canal
99
m(t)
+1
t
-1
Tb
(a)
Porteuse
Ac
t
Ac
(b)
s(t)
Ac
t
Ac
(c)
F IG . 6.11 Modulation BPSK du signal.
Emetteur
Signal de
donnes b(t)
Modulation
binaire PSK
Rcepteur
Canal
s(t)
x(t)
c(t)
Porteuse
y(t)
u(t)
Dtecteur
cohrent
c(t)
j(t)
Gnrateur
de code PN
Porteuse locale
Gnrateur
local de
code PN
100
Estimation
de b(t)
(6.18)
(6.19)
(6.20)
(6.21)
(6.22)
Cette dernire quation montre que le signal lentre du dmodulateur cohrent BPSK est un signal
modul en BPSK augment de linterfrence c(t) j(t). La dmodulation peut alors se faire, comme
pour le cas de la transmission en bande de base, en utilisant un filtre adapt suivi dun organe de
dcision, le terme dinterfrence tant considr comme un bruit blanc additif gaussien.
(6.24)
Ce gain mesure limpact sur le rapport signal sur bruit obtenu par usage de la technique dtalement.
On constate que le gain dtalement gale le facteur dtalement N. Plus la squence dtalement est
longue, plus le gain dtalement est important.
o Eb est lnergie par bit du message b(t), Tc la dure dun chip et I la puissance moyenne du signal
dinterfrence
Z
1 Tb 2
I=
j (t) dt
(6.26)
Tb 0
En fait, si on compare cette probabilit derreur celle obtenue pour la dmodulation BPSK, tout
se passe comme si le signal BPSK tait plong dans un bruit blanc gaussien de densit spectrale de
puissance gale ITc .
(6.27)
Eb
N0 req
I
S
B
Rb
(6.28)
(6.29)
Eb
N0
Exemple. Soit transmettre un signal vocal numrique une cadence bit de 9600 [Hz]. Prenons ensuite une largeur de canal de B = 1, 2288 [MHz] telle quutilise dans le systme amricain talement
de spectre IS-95. Pour un rapport nergie bruit requis de 6 [dB], ce qui correspond une certaine
probabilit derreur, on calcule une marge de
1, 2288 106
10 log
6 = 15, 1 [dB]
(6.30)
9600
En conclusion, les bits dinformation la sortie du rcepteur peuvent tre dtects de manire
fiable mme lorsque le bruit ou les interfrences lentre du rcepteur sont majors de 15, 1 [dB].
102
Cette marge disponible va donc permettre de rajouter des utilisateurs sans hypothquer la qualit de
la communication.
Ltalement a toutefois ses limites car il faut bien borner la valeur de N. Pour accrotre davantage
encore le nombre dutilisateurs, il est ncessaire de recourir une technique dtalement de spectre
avec saut de frquence.
(6.32)
Soit M le nombre de mobiles traits par la cellule. En supposant que tous les metteurs travaillent
mme niveau de puissance et en ngligeant le bruit thermique, les autres mobiles produisent une
puissance dinterfrence
I = C(M 1)
(6.33)
ce qui conduit
C
1
=
I
M 1
(6.34)
1
1
B
B
+1
Rb Eb /N0
Rb Eb /N0
(6.35)
On voit donc que la capacit de la cellule est bien proportionnelle au gain dtalement RBb . M est
donc approximativement gal la marge dinterfrence. Dans lexemple du paragraphe prcdent, on
observe un nombre dutilisateurs gal M = 101,51 32. On voit ds lors que, plus on veut mettre
des utilisateurs sur la mme bande de frquences pour une mme probabilit derreur, plus le gain
dtalement doit tre important et donc plus B doit tre lev par rapport Rb , et donc, plus Tc doit
tre petit par rapport Tb . Une interprtation intressante est montre la figure 6.13 ; la puissance de
lutilisateur ramene dans une bande large de W [Hz] dpasse largement la puissance dinterfrence
des autres utilisateurs.
103
Puissance = C
11
00
1
0
1
0
1
0
111
000
N
00
11
00
11
(M 1)utilisateurs
00
111
000
0101010101010101010111
00
0111
1
0
1
0
1
0
1
0
1
0
1
0
1
0
1
0
1
0
1
0
1
0
1
0
00 0110
1010 0111
W
1010
10 B
marge
Puissance = C
Surface utile
Surface gographiquetotale
(6.38)
B
1 1
GFe
Rb Eb /N0
Pour le GSM, Fe 1.
104
(6.39)
Chapitre 7
Multiplexage et accs au multiplex
7.1 Introduction
La paire torsade qui relie un abonn au central tlphonique est habituellement ddie un
utilisateur unique1. Cet utilisateur peut ds lors, sa guise et toute heure, utiliser la totalit du
spectre frquentiel ; par ailleurs, la technologie ADSL en tire profit. Cette souplesse a un prix car, la
ligne restant inactive pendant une trs longue priode, lutilisation moyenne de quelques dizaines de
minutes par jour doit suffire rentabiliser linvestissement consenti par loprateur.
7.1.1 Multiplexage
Dans lair ou sur un rseau de tldistribution, le partage du spectre et le partage temporel des
ressources font partie intgrante du dimensionnement du rseau. Ces partages rsultent de procds
de multiplexage qui consistent combiner plusieurs signaux pour les transmettre sur un mme
support.
Historiquement, on en distingue principalement deux :
le multiplexage en frquences (Frequency Division Multiplexing - FDM). Cette technique de
multiplexage alloue une bande de frquences spcifique chaque signal.
le multiplexage temporel (Time Division Multiplexing - TDM). Il consiste rguler les moments
doccupation du canal pour chaque signal.
La technique plus rcente du multiplexage par talement de spectre a ouvert de nouvelles voies pour
le partage de ressources. Cette technique ralise la fois un partage frquentiel et temporel.
Pour tre complet, signalons quil existe des techniques de multiplexage propres certains supports. Ainsi, dans le domaine de la transmission par fibre optique, on a dvelopp des techniques
de multiplexage par longueurs donde (Wave Division Multiplexing - WDM) avec quelques variantes
particulirement adaptes la transmission trs haut dbit (Dense Wave Division Multiplexing DWDM) .
Dans de trs rares cas, il existe des concentrateurs de lignes qui combinent des signaux entre les abonns et le central
tlphonique.
105
des stratgies pour disposer du canal. Cest par le biais de stratgies quon espre atteindre au
haut de performance (haut dbit, faible dlai, faible taux de congestion, . . .). La question des
stratgies ne sera pas aborde ici.
f1
f1
f1
X2 ( f )
f2
f2
f2
X3 ( f )
f
f3
f3
f3
Signal multiplex
f3
f2 f1
f1
f2
f3
Signal multiplex
f3
f2 f1
f1
f1
f2
f3
f2
f1
f3
f2
f3
X1 ( f )
X2 ( f )
X3 ( f )
107
La borne infrieure de la bande passante est fournie par la somme des bandes des messages constituants. Ainsi, soit W la bande totale et Wi les bandes individuelles,
W > Wi
(7.1)
C
Canal physique
Utilisateur 2
Non occup
A,B,C : porteuses
B
A
Temps
F IG . 7.3 Partage des ressources par multiplexage en frquence : le FDMA (Frequency Division
Multiple Access).
(7.2)
o chaque signal est modul par une modulation bande latrale rsiduelle. Les groupes sont ensuite
rassembls pour former des super-groupes et ainsi de suite. Le second chelon de la hirarchie est
reprsent la figure 7.5.
108
Frquence porteuse
des voix (en kHz)
Frquence porteuse
des groupes (en kHz)
108 kHz
108
104
100
96
92
88
84
80
76
72
68
64
4 kHz
0
612
564
516
468
420
12
11
10
9
8
7
6
5
4
3
2
1
5
4
3
2
1
552 kHz
504
456
408
360
312
Supergroupe
de 5 groupes
60
1 voix
Groupe basique
de 12 voix
groupes basiques
frquences porteuses
612
564
516
468
420
60
411.920
108
groupe No
312
2
360
3
408
4
456
5
504
552
109
kHz
un mme canal (appel voie haute vitesse HV), un ensemble dinformations provenant de diffrents
canaux faibles dbits (appels voies basses vitesses BV) lorsque celles-ci doivent communiquer
simultanment dun mme point de dpart un mme point darrive. Le principe est illustr la
figure 7.6.
1
2
M
U
L
T
I
P
L
E
X
Voies BV
n
1
2
M
U
L
T
I
P
L
E
X
Voie HV
"mots" contenant
linformation de chaque voie
F IG . 7.7 Structure temporelle de multiplexage par rpartition temporelle des ressources.
Le multiplexage peut tre effectu soit aprs lopration dchantillonnage de chaque voie et avant
le codage, soit aprs celui-ci.
Lors du dmultiplexage des voies, cest--dire de la r-affectation des donnes en fonction de leur
source dorigine, il faut que lon puisse distinguer les diffrents chantillons. Pour cela, on intercale
dans la succession de donnes un mot de repre, dont loccurrence est gale la priode dchantillonnage (figure 7.8).
mots de repre
1
Te
F IG . 7.8 Mot de repre.
Lensemble constitu du mot de repre et de linformation de mme rang pour toutes les voies est
appel trame. Sa dure est gale Te .
110
synchronisation
utilisateur 1
source 2
canal
utilisateur 2
..
.
..
.
signal en bande de base
source N
s2
s1
s1
s2
sN
.....
sN
s1
s1
s2
utilisateur N
..... sN
.....
s2
.....
t
(7.3)
Si le signal final est en bande de base et que cette bande vaut W , il faut imprativement 2W chantillons au moins pour le caractriser. Pour lintervalle T , on a donc
ns = 2W T = 2Wi T
(7.4)
W = Wi
(7.5)
Ds lors,
Canal physique
Frquence
Utilisateur 1
Utilisateur 2
Non occup
1,2,3,4 : timeslots
1 2
1 2
3 4
Temps
F IG . 7.10 Partage des ressources par multiplexage en temps : le TDMA (Time Division Multiple
Access).
......
16
......
30
31
......
16
......
30
31
3, 9[ s]
1 2 3 4 5 6 7 8
488[ns]
Rle des IT
Les IT de 1 15 et 17 31 sont affects aux voies de parole tlphonique ou la transmission
numrique (donnes, fax, . . .). Ils ont chacun un dbit de 64 [kb/s].
LIT0 contient le mot de repre voqu prcdemment, appel mot de verrouillage de trame,
et not VT. Il ne joue pas le mme rle selon quil se situe dans une trame paire ou impaire.
Le verrouillage de trame permet de synchroniser les quipements de rception sur le cycle
des donnes mises la mise en marche du systme, puis dune manire priodique afin de
vrifier le synchronisme du fonctionnement gnral. Dans ce cas prcis, on parle de verrouillage
regroup, cest--dire dont les informations sont groupes dans un mme IT en tte de trame.
Il existe un autre type de verrouillage dit rparti si le mot possde des informations localises
dans plusieurs IT diffrents, ou bien sur plusieurs trames diffrentes.
LIT16 est rserv la signalisation des 30 voies. Celle-ci a pour but de transmettre les informations relatives ou pas aux donnes (maintenance des systmes, tests, commandes, gestion du
rseau...), pour la commande des oprations de commutation. Elle est ici de type hors octet
regroup, cest--dire regroup dans lIT 16. Chaque IT16 a pour rle la signalisation simultane de deux voies. Il faut donc 15 IT16, soit 15 trames, pour signaliser toutes les voies de
donnes. Lensemble de ces 15 trames, auxquelles on ajoute une 16-ime trame appele trame
0, constitue une multitrame. Cette trame 0 comprend dans lIT16 un mot de repre, appel mot
de verrouillage multitrame (VMT) ncessaire lors du multiplexage la synchronisation des
trames. Il occupe les quatre premiers bits de cet IT (figure 7.13).
b j (t)=0
Z Tb
bi (t )c j (t)ci (t )dt
Z Tb
0
c j (t)ci(t )dt
(7.6)
(7.7)
(7.8)
1 Tb
ji ( ) =
(7.9)
c j (t)ci(t )dt
Tb 0
Cette moyenne temporelle est appele fonction dintercorrlation partielle. Pour avoir une interfrence nulle, il faudrait que cette fonction soit nulle pour toute valeur . En gnral, les squences
113
MULTITRAME
125[ s]
T 0 T 1 T 2 T 3 T 4 T 5 T 6 T 7 T 8 T 9 T 10 T 11 T 12 T 13 T 14 T 15
125[ s]
IT IT IT IT
0 1 2 3
IT
16
IT
31
3, 9[ s]
Mot VT
X 0 0 1 1 0 1 1
0 0 0 0 X X X
Mot altern :
1 trame sur 2
X 1 X X X X X
488[ns]
IT 1 15
IT 17 31
P 1 2 3 4 5 6 7
IT 16 de la T 0
30 voies tlphoniques
Mot VMT
IT 16 trame T 1
a b c d a b c d
Signalisation
Signalisation
Voie 1
Voie 16
X : bits libres
: bits dalarmes
IT 16 trame T 2
a b c d a b c d
Signalisation
Signalisation
Voie 2
Voie 17
IT 16 trame T 15
a b c d a b c d
Signalisation
Signalisation
Voie 15
Voie 30
Code
Canal physique
Utilisateur 1
Utilisateur 2
Non occup
Frquence
Temps
F IG . 7.14 Multiplexage de ressources par rpartition de code (Code Division Multiple Access).
114
utilisateur j
b j (t) = 0
utilisateur i
bi (t )
ci (t )
Transmetteur du
i me utilisateur
c j (t)
Transmetteur du
j me utilisateur
cos (2 fct)
cos (2 fct)
cos (2 fct)
cos (2 fct)
canal
ci (t )
R Tb
0
Rcepteur du
i me utilisateur
c j (t)
R Tb
dt
vj
vi
organe
dcision
soit 1 si
vi > 0
dt
seuil 0 seuil 0
soit 0 si
vi < 0
Rcepteur du
j me utilisateur
organe
dcision
soit 1 si
vj > 0
soit 0 si
vj < 0
115
longueur maximale nont pas de bonnes proprits dintercorrlation partielle comme le montre la
figure 7.16 ; la fonction dintercorrlation est loin dtre nulle pour toute valeur de .
20
15
10
5
Ri j ( ) 0
-5
-10
-15
-20
-60
-40
-20
20
40
60
Delai
7.5 Discussion
7.5.1 Combinaison de techniques de multiplexage
Les techniques de multiplexage ne sont pas ncessairement utilises isolment. En effet, dans
le cas pratique du GSM, loprateur alloue une srie de canaux frquentiels. Ces canaux sont ensuite partags entre utilisateurs par multiplexage rpartition temporelle. On parle de technique
FD/TDMA. Un tel schma est illustr la figure 7.17. De plus, il y a sparation des bandes de
frquence suivant que le signal va de la station de base vers lutilisateur ou linverse. Cette technique
est parfois appele FDD/TDD (le second D signifiant Duplex).
Dans le cas du standard amricain IS-95, le multiplexage temporel est remplac par un multiplexage par spectre tal ; il sagit alors de FD/CDMA.
116
Frquence
1 2
1 2
3 4
Utilisateur 1
Utilisateur 2
Non occup
1,2,3,4 : timeslots
A,B,C : porteuses
Temps
F IG . 7.17 Partage de ressources par multiplexage en temps et en frquence, combinaison du TDMA
et du FDMA.
DCS-1800
TDMA/FDMA
IS-54
TDMA/FDMA
IS-95 (DS)
CDMA/FDMA
1710-1785
1805-1880
(Europe)
869-894
824-849
(USA)
869-894
824-849
(USA)
200
200
GMSK
30
30
/4 DQPSK
1250
1250
BPSK/QPSK
1W/125mW
RPE-LTP
13
600mW/200mW
VSELP
7,95
600mW
QCELP
8 (var.)
270,833
270,833
4,615
48,6
48,6
40
20
117
118
Chapitre 8
tude du trafic
8.1 Introduction
La difficult majeure de lexploitation dun rseau consiste en son dimensionnement. En effet, tout
est question de compromis entre le risque de blocage d une occupation de toutes les ressources et
le cot quentrane un sur-dimensionnement dun rseau.
Le rseau tlphonique est critique du point de vue de la gestion du trafic, en raison de sa structure
de rseau commut (cf. illustration la figure 8.1). En effet, une communication occupe une voie de
bout en bout pour la totalit de la dure dun appel. Il faut donc viter quun commutateur nait plus
aucun circuit disponible lors de ltablissement dun nouvel appel.
Commutateur local
Commutateur
intermdiaire
Commutateur local
Tlphone
fils de cuivre
Tlphone
fibre optique
lien arien
fils de cuivre
Commutateur
Commutateur
pour les appels prioritaires pour lesquels le systme libre les circuits ncessaires.
120
temps
Nombre de communications
5
4
3
2
1
0
temps
T
N
N
i=1 0 1i (t)dt
0 i=1 1i (t)dt
I=
=
(8.1)
T
T
Lintensit de trafic reprsente donc loccupation moyenne du lien durant une certaine priode. Bien
que lintensit soit adimensionnelle, comme lexpression dun rapport entre deux temps, on parle
dE RLANG, not [E], du nom du thoricien Danois, pre de la thorie. Remarquons que la capacit
maximale dun circuit est de 1 E RLANG. Dans un systme perte, la capacit nexcde jamais le
nombre de circuits.
En pratique, on utilise deux paramtres importants pour caractriser le trafic :
N
i=1
RT
0
1i (t)dt
#T
(8.2)
Dfinition 26 [Charge de trafic coul] On dfinit la charge coule dun faisceau, exprime en
E RLANG, comme le produit du taux dappels entrants par la dure moyenne dappel. On la note Ae .
Ae = e t m
2 Il
121
(8.3)
Il est important de remarquer que lon exclut, dans la charge, la charge due la signalisation. Dautre
part, la charge, qui est une combinaison de deux paramtres moyens, ne fournit aucune information
quant une corrlation ventuelle entre le taux dappel et la dure.
Exemple : lanalyse des appels dune socit rvle un taux dappel de 40 appels par heure et une
dure moyenne par appel de 5 minutes. La charge vaut donc
40
200
5 =
= 3, 3 [E]
(8.4)
60
60
ce qui signifie quil y a en moyenne 3,3 lignes occupes si la distribution est uniforme.
Charge de rfrence. La charge maximale disponible sur un faisceau N liens vaut thoriquement
N. Elle correspond une occupation permanente des canaux. En pratique, comme les appels ont lieu
de manire alatoire, il arrive que le faisceau soit congestionn, cest--dire que les N liens soient
occups. La charge pratique est donc infrieure N.
Lors de ltude du trafic, il faut distinguer deux types de charge :
1. la charge offerte, cest la charge qui serait transporte par le rseau sil pouvait honorer toutes
les demandes de connexion.
2. la charge coule ; la charge rellement mesure dans le rseau.
La charge offerte est utilise dans la dfinition de la charge A.
Dfinition 27 La charge offerte A vaut
A = tm
(8.5)
D8
D1
3
2
1
t
T
F IG . 8.4 Ralisations dun processus de comptage.
Plutt que cette courbe, cest le nombre doccurrences pendant tout intervalle [t1,t2] qui nous
intresse ; aprs tout, on devra bien se choisir une dure typique pour lanalyse. Imaginons que lon
puisse dcouper laxe temporel en intervalles de largeur T tellement troits quils contiennent tout
au plus une occurrence. Chacun de ces intervalles, rfrencs par ]nT, (n + 1)T ], est le lieu dune
variable alatoire binaire Dn . Ds lors, Dn = 1 implique quune tentative dappel a t mene pendant
lintervalle de temps ]nT, (n + 1)T], Dn = 0 mentionne quaucune tentative dappel na t mene
pendant cet intervalle. supposer que ces variables alatoires soient indpendantes, la rsultante nest
autre quune loi binomiale.
Par gnralisation, pour toute constante positive > 0, il est possible de choisir une valeur pour
T telle que T < 1 (pour pouvoir dfinir une probabilit). Si T reprsente la probabilit
de russite3 de chacune des variables D1 , D2 , . . ., considrant un intervalle de temps T = mT , on
obtient les probabilits suivantes :
3 Et
123
1. probabilit davoir n tentatives dappels pendant la dure T et donc sur les m intervalles de
temps T (n russites pour la loi binomiale) :
T n
n
(8.8)
( T ) =
m
2. probabilit davoir m n intervalles de temps T sans tentative dappel (m n checs)
T mn
mn
(1 T )
= 1
m
(8.9)
m!
.
n!(mn)!
(8.15)
(8.16)
124
et donc
Ds lors,
m
1 mT
T m
n = lim 1
lim
= e T
m+
m+
T
m
1 m
lim fDm (n) =
m+
( T )n T
n! e
n = 0, 1, . . . , m
n 6= 0, 1, . . . , m
(8.17)
(8.18)
qui nest autre que la variable alatoire de P OISSON de paramtre = T . tant donn que lon
fait tendre m vers linfini, lindice de la variable alatoire Dm na plus de raison dtre ; nous crirons
donc D.
Dfinition 28 [Variable alatoire de P OISSON] On appelle variable alatoire de P OISSON une variable alatoire D discrte dont la densit de probabilit est
( T )n
T n = 0, 1, . . .
n! e
(8.19)
fD (n) =
0
n 6= 0, 1, . . .
Elle reprsente donc ici la probabilit davoir n tentatives dappels pendant la dure T .
La loi de P OISSON fait intervenir un taux moyen doccurrences et un intervalle de temps dobservation T ; ces dnominations sont justifies par le thorme suivant (on dfinit = T ) :
Thorme 29 [35, page 63] Soit effectuer n tirages. Pour chaque tirage, la probabilit de lvnement pertinent est dfinie gale n , o > 0 est une constante et n > . Soit une variable alatoire
Sn qui dnombre les occurrences dvnements pertinents pour n tirages effectus. Lorsque n +,
la densit de probabilit de Sn converge vers une loi de P OISSON de paramtre .
Dmonstration.
Elle dcoule des dveloppements prcdents.
Thorme 30 Lesprance et la variance dune variable alatoire de P OISSON valent toutes deux
D =
D2 =
Dmonstration.
Dmontrons la relation 8.20.
D = i fD (i) = i
i=0
i=0
(8.20)
(8.21)
i
e
i!
(8.22)
D =
j+1
j
e =
j=0 j!
| {z }
+
( j + 1) ( j + 1)! e =
j=0
(8.23)
Le fait que la somme soit gale 1 dcoule du fait que la fonction de rpartition de la loi de P OISSON
vaut 1 linfini (la somme des probabilits dfinies par la variable alatoire de P OISSON vaut 1). La
relation 8.21 se dmontre dune manire identique.
Ainsi donc, = T est gal la valeur moyenne de la loi de P OISSON, do le fait que = T est
la valeur moyenne du nombre doccurrences par unit de temps ; on parle parfois dintensit pour .
125
Question 31 Soit un taux moyen dappels de 4 par minute. Quelle est la probabilit P5 davoir 5
appels ou plus pendant une minute ?
Rponse
Par la variable de P OISSON,
+
i=5
i=0
P5 = fD (i) = 1 fD (i)
(8.24)
41 42 43 44 4
e = 0, 37
P5 = 1 1 + + + +
1! 2! 3! 4!
(8.25)
Comme = T = 4 1 = 4,
(8.27)
Pour toute valeur x > 0, Xn > x si et seulement sil ny a aucune occurrence durant lintervalle de
temps ]tn1 , tn1 + x]. Par dfinition du processus de P OISSON, la valeur des variables X1, . . . , Xn1
126
D(t)
occurrence 4
5
4
occurrence 1
3
2
1
t
X1
X2
X3
X4
X5
(8.28)
Or, cette probabilit reprsente la probabilit de ralisation dune variable alatoire de P OISSON dont
le paramtre n est nul, donc la probabilit de ne pas avoir de tentative dappel pendant lintervalle de
temps ]tn1 , tn1 + x]. Ds lors,
(8.29)
p(D(tn1 + x) D(tn1 ) = 0) = e x
Xn est donc bien indpendante des autres Xi et sa fonction de rpartition vaut
1 e x x > 0
FXn (x) = p(Xn x) = 1 p(Xn > x) =
0
x0
(8.30)
(8.31)
Question 35 Soit un central tlphonique desservant 10.000 utilisateurs. Supposons que chaque
utilisateur gnre une tentative dappel par heure. Dterminez la cadence (en nombre de fois par
seconde) laquelle deux appels arrivent moins de 0,1 [s] entre eux.
Rponse
Le taux moyen dappels vaut
10.000
= 2, 78 [appel/s]
3600
Selon lquation 8.31, la probabilit de ne pas avoir de tentative pendant 100 [ms] vaut
Ce qui revient dire que 24,3% des tentatives dappel sont cartes de moins de 100 [ms]. La cadence
sobtient alors par
(1 P0 (0, 278)) = 0, 67 [ f ois/s]
(8.32)
8.3.1 Modle dun systme avec perte (sans mmoire) : statistique dE RLANG
B
La statistique E RLANG B va nous permettre de dterminer la probabilit de tentative dappels lors
dune congestion du rseau (probabilit que les N lignes du rseau soient occupes) en fonction de la
charge et du nombre de lignes caractrisant le tronon considr.
Pour analyser la charge dun rseau, il faut prendre en compte le nombre de tentatives dappel, le
nombre de connexions en cours et le nombre darrts dappel.
Nombre de tentatives dappel
Le nombre de tentatives dappel peut tre vu comme un processus de dnombrement. Pendant
lintervalle de temps t (intervalle de temps dobservation), le nombre de tentatives dappel NA est
donc une variable alatoire de P OISSON telle que
p(NA = n) =
( t)n t
e
, n = 0, 1, . . .
n!
(8.35)
( t)n t
e
, n = 0, 1, . . .
n!
(8.36)
Charge
tout moment, des tentatives dappel et des arrts peuvent avoir lieu sur la ligne ; la charge
fluctue donc entre 0 et N appels en cours. Pour un intervalle de temps trs court dt, la probabilit
davoir exactement une tentative dappel (n = 1) vaut dte dt dt alors que la probabilit davoir
un appel qui se termine vaut dte dt dt.
Ds lors, supposer que k lignes sur un total de N lignes soient occupes linstant t, on calcule
respectivement trois probabilits sur un intervalle de temps dt
P1 la probabilit dune tentative dappel,
P2 la probabilit dun arrt et
P3 la probabilit dun statu quo en matire doccupation de lignes.
Elles valent respectivement
P1 = dt
(8.37)
P2 = Ck1 ( dt)1(1 dt)k1 k dt
(8.38)
car on tient compte du fait que nimporte quelle ligne parmi les k lignes peut se librer pendant
lintervalle de temps dt. Finalement, lvnement de statu quo se produit quand il ny a pas de tentative
dappel et quaucune des k lignes occupes nest sujet un relchement. Ainsi, la probabilit P3 est
donne par
P3 = (1 dt)(1 k dt)
= 1 dt k dt + k (dt)2
1 dt k dt
129
(8.39)
(8.40)
(8.41)
Dfinissons alors p(k;t + dt) comme la probabilit quil y ait k lignes occupes linstant t + dt. Il
est alors possible de dterminer p(k;t + dt) en utilisant les probabilits P1 , P2 et P3
p(k;t + dt) = P3 p(k;t) + P1 p(k 1;t) + P2 p(k + 1;t)
(1 dt k dt) p(k;t) + dt p(k 1;t) + (k + 1) dt p(k + 1;t)
(8.42)
Le premier terme correspond au fait quil y avait dj k lignes occupes linstant t mais quaucune tentative dappel ni darrt nait eu lieu. Le deuxime terme correspond au fait quil y avait
k 1 lignes occupes linstant t mais quune tentative dappel (ds lors russie) ait eu lieu. Enfin,
le dernier terme correspond au fait quil y avait k + 1 lignes occupes linstant t et quune ligne se
soit libre pendant lintervalle de temps dt. Il existe cependant deux cas particuliers correspondant
k = 0 et k = N
p(0;t + dt) = (1 dt) p(0;t) + dt p(1;t)
(8.43)
reprsentant le cas o aucune ligne nest occupe (il ny a donc pas de possibilit de relchement
dune ligne) et
(8.44)
p(N;t + dt) = (1 dt N dt) p(N;t) + dt p(N 1;t)
pour le cas o toutes les lignes sont occupes (une nouvelle tentative dappel ne peut donc tre
gnre).
En rgime, on peut faire lhypothse que les probabilits ne sont pas fonction du temps et crire
p(k;t + dt) = p(k;t) = Pk ,
k = 0, 1, 2, ..., N
(8.45)
0 = [( k ) Pk + Pk1 + (k + 1) Pk+1 ] dt
(8.46)
(8.47)
(8.48)
P0 = P1 , k = 0
( + N ) PN = PN1 , k = N
(8.49)
(8.50)
(8.51)
On peut montrer [22] que lexpression de Pk vrifiant toutes ces conditions est donne par
Pk =
( / )k
k!
N ( / )i
i=0 i!
(8.52)
Cette formule reprsente ainsi la probabilit davoir k lignes occupes. Elle est valable k [0, N].
Formule dE RLANG B. Ltat qui rsulte dune occupation de toutes les lignes est appel congestion. Si un appel est rejet en raison dune occupation des N lignes, la probabilit de cet vnement
de blocage est celle de PN (k = N)
B = PN =
( / )N
N!
N ( / )i
i=0 i!
(8.53)
E{k} =
( / )k
k!
( / )i
i=0 i!
N1
=
k=0 N
(1 B) = A(1 B)
(8.55)
o A est la charge offerte du faisceau et A(1 B) la charge relle. Dans la mesure o la probabilit de
blocage est petite, loccupation moyenne des N lignes est gale .
La probabilit de blocage sexprime donc par lexpression
B = PN =
AN
N!
N Ak
k=0 k!
, A=
(8.56)
10
10
10
10
10
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
A/N
F IG . 8.6 Loi de probabilit E RLANG B (daprs [22, page 271]).
131
N
A
A
N
10
4,5
0,45
20
12,0
0,60
30
20,3
0,68
40
29,0
0,73
50
37,9
0,76
60
46,9
0,78
70
56,1
0,80
80
65,4
0,82
90
74,7
0,83
100
84,1
0,84
A
1B
(8.57)
La table de la figure 8.2 peut tre utilise pour trouver A . Il est important de noter quen principe le
rseau doit tre dimensionn pour cette charge A et non pour la charge A. Alors A peut tre obtenue
par A = (1 B)A . Par exemple, pour B = 0, 02 et N = 30, la table donne A = 21, 9 qui doit tre
vu comme la valeur de A . La valeur de A est alors gale (1 0, 02) 21, 9 = 21, 5. Si on dsire
offrir chaque utilisateur une charge moyenne A0 = 0, 03, on doit limiter le nombre dutilisateurs
M = A/A0 = 715, contre 730 lorsque quaucune tentative ayant chou nest reconduite.
132
B
N
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
0,01
0,01
0,153
0,46
0,87
1,4
1,9
2,5
3,1
3,8
4,5
5,2
5,9
6,6
7,4
8,1
8,9
9,7
10,4
11,2
12,0
12,8
13,7
14,5
15,3
16,1
17,0
17,8
18,6
19,5
20,3
0,005
0,005
0,105
0,35
0,7
1,1
1,6
2,2
2,7
3,3
4,0
4,6
5,3
6,0
6,7
7,4
8,1
8,8
9,6
10,3
11,1
11,9
12,6
13,4
14,2
15,0
15,8
16,6
17,4
18,2
19,0
0,003
0,003
0,081
0,29
0,6
1,0
1,4
1,9
2,5
3,1
3,6
4,3
4,9
5,6
6,2
6,9
7,6
8,3
9,0
9,8
10,5
11,2
12,0
12,7
13,5
14,3
15,1
15,8
16,6
17,4
18,2
0,001
0,001
0,046
0,19
0,44
0,8
1,1
1,6
2,1
2,6
3,1
3,7
4,2
4,8
5,4
6,1
6,7
7,4
8,0
8,7
9,4
10,1
10,8
11,5
12,2
13,0
13,7
14,4
15,2
15,9
16,7
N
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
0,01
21,2
22,0
22,9
23,8
24,6
25,5
26,4
27,3
28,1
29,0
29,9
30,8
31,7
32,5
33,4
34,3
35,2
36,1
37,0
37,9
38,8
39,7
40,6
41,5
42,4
43,3
44,2
45,1
46,0
47,0
0,005
19,9
20,7
21,5
22,3
23,2
24,0
24,8
25,7
26,5
27,4
28,2
29,1
29,9
30,8
31,7
32,5
33,4
34,2
35,1
36,0
36,9
37,7
38,6
39,5
40,4
41,2
42,1
43,0
43,9
44,8
133
B
0,003
19,0
19,8
20,6
21,4
22,2
23,1
23,9
24,7
25,5
26,3
27,2
28,0
28,8
29,7
30,5
31,4
32,2
33,1
33,9
34,8
35,6
36,5
37,3
38,2
39,0
39,9
40,8
41,6
42,5
43,4
0,001
17,4
18,2
19,0
19,7
20,5
21,3
22,1
22,9
23,7
24,4
25,2
26,0
26,8
27,6
28,4
29,3
30,1
30,9
31,7
32,5
33,3
34,2
35,0
35,8
36,6
37,5
38,3
39,1
40,0
40,8
Reconduite des
appels chous
E RLANG C
M OLINA
E NGEST
B ERNOULLI
10
E RLANG C
M OLINA
10
B
E RLANG B
2
10
Reconduite des
tentatives dappel
10
10
15
20
25
30
35
40
45
50
A [E RLANG]
F IG . 8.7 Comparaison des probabilits de blocage (daprs [22, page 332]).
134
Troisime partie
lments de la couche physique
135
Chapitre 9
Transmission sur ligne paires symtriques
et rseau tlphonique
9.1 Introduction
La transmission peut se raliser sur diffrents supports (cble coaxial, paire torsade, fibre optique,
. . .) ou dans lair. Si les lois de M AXWELL rgissent lensemble des phnomnes de propagation lectromagntiques, les proprits des matriaux influencent considrablement les performances, rendant
les supports non interchangeables bien des gards.
La figure 9.1 fournit les ordres de grandeur des performances distance/dbit que lon peut obtenir
pour diffrents supports de transmission physiques. Trs clairement, cest lutilisation de la fibre
monomode qui reprsente le meilleur choix actuel.
Distance maximale [km]
1000
Fibre monomode
100
Cble
10
Paire torsade
1
1
10
100
1000
10000
Dbit [Mb/s]
F IG . 9.1 Comparaison de divers supports de transmission (daprs [34]).
9.1.1 Limitations
Plusieurs phnomnes affectent la transmission dun signal numrique par une onde lectromagntique. Citons en quelques-uns :
Attnuation Lattnuation est un phnomne quivalent une perte dnergie du signal se
propageant. On lexprime gnralement en [dB] par kilomtre.
137
Distorsion Un canal qui fournit en sortie la version du signal original ainsi que des versions attnues
dcales dans le temps, introduit un effet de distorsion damplitude (lamplitude de la transmittance du canal nest pas constante en fonction de la frquence) et de phase (la phase nest pas
linaire en fonction de la frquence).
Dispersion La dispersion est le phnomne lorigine dun talement de londe, ce qui dans le cas
de communications numriques, se rpercute par une confusion entre symboles successifs (le
canal possde une bande passante finie).
Bruit Il est totalement impossible dviter un bruit dans un systme de communications. En fait, les
quipements gnrent eux-mmes un bruit. cela vient se rajouter le bruit du canal. Aujourdhui, on prfre utiliser des techniques de transmission numriques pour lutter contre les effets
dun bruit additif de canal. On peut ainsi esprer un taux derreur de lordre de 1012 pour une
transmission sur fibre optique et 107 sur une ligne paire de cuivre.
(9.1)
(9.2)
Exemple. Une taille de paquet de N = 105 bits et un taux derreur par bit Pe = 107 conduisent
PP 102 .
9.2.1 Dfinitions
Dfinition 37 On appelle ligne tout support physique de transmission constitu dun milieu matriel
fini.
Les deux formes de lignes les plus courantes sont la paire torsade constitue de deux conducteurs
arrangs en hlice et la paire coaxiale forme de deux conducteurs concentriques spars par un
isolant.
Dfinition 38 On appelle cble de transmission, tout support physique constitu dun ensemble de
lignes.
138
l
S
(9.3)
ln(s/r)
C=
(9.5)
s
r2
r1
Coaxe
r
Paire de fils
139
Inductance
Quant la valeur de linductance, elle sobtient en tirant profit de la proprit gnrale suivante
LC =
(9.6)
Dans ce qui suit, la ligne est considre comme idale, cest--dire comme ayant les mmes proprits lectriques sur toute sa longueur.
9.2.3 Propagation
La propagation dun onde lectromagntique le long dune ligne de transmission prend la forme
dune onde lectromagntique transverse (TEM, Transverse Electromagnetic Mode). Londe se
propage principalement dans le dilectrique qui spare les deux conducteurs dune ligne de transmission. La figure 9.4 montre les formes des champs dans une srie de structures lectriques deux
conducteurs.
coaxe
2 plans
bifilaire
Lignes de champ magntique
Lignes de champ lectrique
F IG . 9.4 Configuration du champ lectromagntique en mode TEM pour quelques types de lignes.
Rdz
Gdz
Cdz
140
Paramtres primaires
R, L, C et G sont appels paramtres primaires de la ligne avec
R = rsistance linique lmentaire, reprsentant la rsistance de la ligne par unit de longueur
[/m]. Elle dpend en particulier de la section et de la nature du conducteur,
L = inductance linique [H/m], modlisant la prsence dun flux variable autour et entre les
structures conductrices,
C = capacit linique [F/m], caractrisant la capacit du dilectrique constituant la ligne,
G = admittance linique [1 /m], reprsentant les pertes dilectriques et les dfauts disolation
de la ligne. Elle dpend de la nature des isolants.
quations des tlgraphistes
En mettant bout bout des segments de ligne infinitsimaux et sur base du schma de la figure 9.6,
on obtient aisment le systme dquations suivantes, dites quations des tlgraphistes,
V
I
= RI L
z
t
I
V
= GV C
z
t
Ldz
I(z)
V (z)
(9.7)
(9.8)
Rdz
Gdz
Cdz
V (z + dz)
= V (z) dV
Charge
Gnrateur
2V
V
V 2
=
RGV
+
(RC
+
LG)
+
LC
z2
t
t2
(9.9)
2V
V 2
=
LC
(9.10)
z2
t2
ce qui correspond une quation dondes bien connue dont la solution est une combinaison linaire
de signaux sinusodaux
V (z,t) = (A cos kz + B sin kz)(C cos 2 f t + D sin 2 f t)
o A, B, C et D sont des constantes dont les valeurs dpendent des conditions initiales.
141
(9.11)
2V
= (R + jL )(G + jC )V (z) = 2V (z)
z2
(9.12)
(9.13)
Londe est donc constitue dune onde incidente (Vi e z ) et dune onde rflchie (Vr e z ). De plus,
on constate que les deux ondes subissent une attnuation e z lie au facteur . On voit tout de suite
que lattnuation crot avec la longueur de la ligne.
La prsence dune attnuation ne signifie pas que toute transmission soit impossible mais bien que
le signal est attnu ds quil y a des pertes dans le conducteur, impliquant une longueur maximale de
la ligne. Lanalyse en dtail de la question montre que lattnuation dpend de la frquence. En fait,
elle augmente avec la frquence. Il est ds lors plus intressant dutiliser les basses frquences pour
la transmission. Nanmoins, rien nempche dutiliser les zones dattnuation plus importantes. Cest
le mode de fonctionnement des modems haut dbit ADSL dont le spectre dutilisation est montr
la figure 9.7.
Tlphone : 0 4kHz
f
en remonte : 25 100kHz
F IG . 9.7 Spectre dun signal ADSL.
Paramtres secondaires
Les paramtres primaires ne modlisent la ligne que dune manire microscopique. On leur prfre
souvent les paramtres dits secondaires suivants pour dterminer les proprits macroscopiques du
support :
impdance caractristique Zc
Cest une donne complexe qui reprsente la valeur de limpdance connecter en bout de
ligne de manire obtenir une impdance dentre de la ligne gale cette impdance ; on la
nomme Zc . En pratique, on doit absolument tenir compte de la valeur de limpdance lors du
raccordement de lignes ou dquipements un rseau. En effet, si deux lignes nont pas la mme
impdance, le droit du raccord est le lieu de rflexions parasites qui diminuent considrablement
les performances de la transmission. En effet, le coefficient de rflexion est donn par
ZL Zc
ZL + Zc
(9.14)
o ZL est limpdance de charge. On voit donc que le coefficient de rflexion sannulle lorsque
limpdance de charge dune ligne (gale limpdance dentre de la ligne raccorde cette
142
ligne) est gal son impdance caractristique. Il convient de remarquer que ladaptation de
ligne nest pas quivalente ladaptation conjugue ; la premire conduit une absence de
rflexion, la seconde un transfert maximum de puissance dans la charge. Ces deux types
dadaptation sont nanmoins compatibles lorsque la charge est purement relle.
coefficient de propagation (cf. supra)
Par dfinition : = + j o
= affaiblissement linique en Nper/mtre [N p/m]1
= dphasage linique (en [rad/m])
Le facteur dattnuation reprsente les pertes subies par le signal lectrique lors de la propagation
le long de la ligne. Il se mesure en injectant un signal lune des extrmits de la ligne et en mesurant
le signal reu lautre extrmit. est li la longueur donde et la vitesse de propagation v de
londe lectromagntique dans le support par
2 f
2
= =
v
v
(9.15)
(9.16)
(9.17)
Ces quations sont gnrales et valables sur tout type de ligne. Toutefois, certaines simplifications sont
possibles en considrant un caractre plutt inductif ou pas de ligne, une frquence dutilisation leve
ou non. Dautre part, on peut raisonnablement admettre que ladmittance linique est ngligeable,
autrement dit G = 0, en prsence dun isolant entre les conducteurs. Et donc,
s
R + j L
(9.18)
Zc
jC
p
(R + j L) jC
(9.19)
j L
p
p
R
(9.22)
jC j L 1 +
j2 L
r
1
C
(9.23)
R
+ j LC
2
L
1 Le
143
(9.24)
=k
o Lo dsigne la longueur de la ligne.
Exprimons prsent lattnuation en dcibels
f Lo
h i
h i 20 ln ek f Lo
A [dB] = 20 log ek f Lo =
ln 10
(9.25)
(9.26)
f Lo
(9.27)
Ds lors, passer de Lo 2Lo double les pertes. Cette quation est importante pour ltude du transfert
de puissance entre paires voisines.
et
Zc =
R
1
C
f
(9.29)
Exemple. Ordre de grandeur des paramtres primaires pour une paire de diamtre 0,4 [mm] : R =
290 [/km], C = 50 [nF/km], = 0, 2 [N p/km] 800 [Hz].
Ces lignes sont utilises pour des liaisons faible distance, par exemple entre abonns et centraux
tlphoniques.
144
R [/km]
52,3
98,7
L [ H/km]
766
67,5
|Z|c []
188
156
[mN p/km]
151
363
Bande passante
4
5
20 [MHz]
100 [MHz]
16 [MHz]
Exemples dutilisation
Distribution tlphonique (voix)
Voix numrique, rseaux locaux Ethernet 10 [Mb/s] et Any
Lan
Rseaux Token Ring
Rseaux locaux Ethernet 10 et 100 [Mb/s], Token Ring et
Any Lan
Lme centrale peut tre compose dun ou de plusieurs fils de cuivre ou dacier. Lisolant sert
sparer lme du conducteur extrieur. Les proprits du cble sont principalement lies la nature de
cet isolant. Le conducteur extrieur est ralis partir de tresses, dun tube ou dune feuille de cuivre
ou daluminium. La gaine extrieure, gnralement en PVC, tflon ou polythylne (PE), sert de protection mcanique et chimique. Elle peut elle-mme tre incluse dans une enveloppe supplmentaire
(en plomb, fibre de verre ou mtallique) pour des conditions extrmes dutilisation.
Proprits
La figure reprsente une vue de face dun cble coaxial.
gaine
d
D
dilectrique
(r )
me
F IG . 9.8 Vue de face dun cble coaxial.
Le rapport D/d des diamtres des conducteurs est directement li lattnuation . Certaines
tudes ont montr quune attnuation minimale par unit de longueur est obtenue pour D/d = 3, 6.
Ainsi, la plupart des cbles commercialiss respectent ce rapport. Pour des conducteurs en cuivre :
0, 692 103 r f d1 + D1
[N p/km] =
log(D/d)
avec d et D en [mm], f en [Hz].
Le tableau 9.3 donne deux exemples des proprits de diffrents cbles coaxiaux normaliss par
lITU (cf. [13, page 66]).
Type
RG58W
RG35BU
d [mm]
0,9
2,65
D [mm]
4,95
17,27
Zc []
50
75
[dB/100m] 3 [GHz]
140
37
signal analogique est converti sous forme numrique et achemin tel quel jusquau dernier central
tlphonique. Le signal numrique est ensuite interpol et transmis sur la paire de labonn destinataire. Ce schma est illustr la figure 9.9.
Analogique
Analogique
Rseau numrique
Analogique
Modem
Modem
Rseau numrique
F IG . 9.10 Transmission par modem dans la bande [300 Hz, 3400 Hz].
Fort heureusement, la transmission par le biais de modems dans la bande vocale nest pas le
seul moyen de transmettre des informations numriques. Il est par exemple possible de recourir au
Rseau Numrique Intgration de Services (RNIS). Il sagit dun prolongement de laccs numrique
jusqu labonn.
Comme reprsent la figure 9.11, laccs au RNIS se prsente sous la forme dun bus offrant
plusieurs canaux en parallle. Pour un accs de base, labonn dispose de 2 canaux 64 [kb/s]. Ces
2
linverse, la commutation dinformations suivant le protocole IP se fait par paquets. Il ny a pas de chemin prtabli
pour toute la dure de lchange.
147
canaux sont compatibles avec les canaux 64 [kb/s] utiliss entre centraux. Ds lors, il est tout fait
possible davoir une communication numrique lorigine et une terminaison analogique.
Numrique
Analogique
Numrique
NT1
Rseau numrique
F IG . 9.11 RNIS.
Lusage dune bande de frquences limite la bande [300 Hz, 3400 Hz] na de sens que sil sagit
dtablir une communication entre deux points distants. Bien entendu, rien ninterdit dutiliser une
bande plus large entre un abonn et le central tlphonique car le support est spcifique un abonn.
Le RNIS utilise cette astuce pour transmettre des signaux numriques dans une bande de frquences,
excdant la bande vocale, entre labonn et le central.
La transmission numrique haut dbit dans le rseau tlphonique par ADSL pousse le principe
plus loin : la transmission de linformation numrique se fait hors bande vocale de manire garantir
la coexistence avec le signal vocal analogique usuel. Cette coexistence nest effective que sur la
paire ddie car au franchissement du central tlphonique, les signaux sont spars par filtrage et
injects dans des rseaux de transmission spcifiques. Cette sparation, effectue au niveau du central
tlphonique par un splitter, est illustre la figure 9.12.
Modem
Splitter
ATM / IP
Rseau numrique
F IG . 9.12 Configuration dune connexion ADSL.
Du ct de labonn, on place un filtre qui spare les bandes de frquences ; ce filtre nest en
principe pas indispensable puisque les signaux vocaux et de donnes occupent des bandes distinctes.
Ce filtre ainsi que sa courbe de gain3 (tension de sortie sur tension dentre, exprim en [dB]) sont
prsents la figure 9.13.
3 Le
148
20
10
0
Gain [dB]
10
20
30
40
50
60
70
10
100
1000
Frquence [Hz]
10000
100000
F IG . 9.13 Filtre de sparation pour les signaux dune ligne convertie lADSL et sa courbe de gain.
Paire 1
Paire 2
Paire 3
Z
Z
Ligne
perturbatrice
Ligne
utile
Puissance Z
du signal
Puissance
du signal
Paradiaphonie
Tldiaphonie
Cble multi-paires
Modem
Paire perturbe
FEXT
NEXT
Paire perturbatrice
Commutateur
de donnes
Modem
F IG . 9.16 Diaphonie dans le rseau tlphonique.
En raison de la structure du rseau, les caractristiques locales ne sont pas invariantes dans tout
le rseau. Lapproche analytique consiste simplifier la situation pour aboutir des conclusions de
principe. On cherchera donc tablir une fonction de transfert de puissance entre lignes de rfrence.
On distingue deux types de modle pour expliquer les phnomnes de diaphonie :
le modle des capacits non quilibres et
le modle des inductances non quilibres.
Nous allons tudier ces deux modles.
Modle des capacits non quilibres
Ce modle considre deux paires voisines, places au dessus dun cran. Comme le montre la
figure 9.17, les quatre fils sont caractriss par une impdance propre, une admittance avec les autres
conducteurs et une capacit une masse commune. On supposera raisonnablement que la conductance
entre fils est nulle parce que les conducteurs sont enrobs dun isolant.
150
Z1 x
V1 (x)
I1 (x)
I1 (x + x)
Y1G x
V1 (x + x)
Y12 x
V2 (x)
Y2G x
I2 (x)
I2 (x + x)
V2 (x + x)
Z2 x
Y13 x
Y14 x
Y23 x Y24 x
Z3 x
V3 (x)
I3 (x)
I3 (x + x)
Y3G x
V3 (x + x)
Y34 x
V4 (x)
Y4G x
I4 (x)
I4 (x + x)
V4 (x + x)
Z4 x
F IG . 9.17 Modle capacitif de deux paires.
Notions de circuit. Pour un systme constitu de N + 1 conducteurs, on peut dfinir N tensions et
N courants ; ces valeurs correspondent N circuits. Lutilisation usuelle dun cble N conducteurs
revient considrer N/2 paires, elles-mmes organises en quarte. Ces paires servent transmettre
un signal en mode diffrentiel, cest--dire de type Vi Vi+1 . Pour arriver N circuits au total, il reste
un solde thorique de N/2 circuits. Mais ces derniers sont dlaisss en pratique en raison des effets
diaphoniques inacceptables quils induisent.
Calcul des relations entre courants et tensions. Pour la facilit, on crit les capacitances comme
des admittances : Y = jC (on admet que G = 0). Pour les calculs, nous utiliserons des tensions
rfrences par rapport la masse.
Lapplication de la loi des mailles donne les quations suivantes
V1 (x + x)
V2 (x + x)
V3 (x + x)
V4 (x + x)
=
=
=
=
V1 (x) I1 (x)Z1 x
V2 (x) I2 (x)Z2 x
V3 (x) I3 (x)Z3 x
V4 (x) I4 (x)Z4 x
(9.30)
(9.31)
(9.32)
(9.33)
Dans toutes ces quations, on peut passer la limite x 0 aprs avoir divis tous les membres par
x. Se rappelant que la drive est dfinie par la formule suivante
d
f (x + x) f (x)
f (x) = lim
dx
x0
x
Les 8 quations sexpriment alors sous la forme matricielle suivante
V1
0
0
0
0 Z1
0
0
0
V2 0
0
0
0
0 Z2
0
0
V3 0
0
0
0
0
0 Z3
0
d
V
0
0
0
0
0
0
0
Z
4
4 =
A1 Y12 Y13 Y14
I
0
0
0
0
dx
1
I2 Y12 A2 Y23 Y24
0
0
0
0
I3 Y13 Y23 A3 Y34
0
0
0
0
I4
Y14 Y24 Y34 A4
0
0
0
0
o
A1
A2
A3
A4
=
=
=
=
(9.34)
V1
V2
V3
V4
I1
I2
I3
I4
= A
S
(9.35)
(9.36)
(9.37)
(9.38)
(9.39)
Les tensions utiles sexprimant comme une diffrence de tensions, on procde aux changements de
variables suivants
V1M = V1 V2
V2M = V3 V4
V1 +V2
V1L =
2
V3 +V4
V2L =
2
I1 I2
I1M =
2
I3 I4
I2M =
2
I1L = I1 + I2
I2L = I3 + I4
Ces quations peuvent galement tre regroupes sous forme matricielle
V1M
V1L
V2M
V2L
= T
S =T
S
I1M
I1L
I2M
I2L
152
(9.40)
(9.41)
(9.42)
(9.43)
(9.44)
(9.45)
(9.46)
(9.47)
(9.48)
1
1 0 0 0 0 0 0
2
1 1 0 0 0 0 0 0
2
0 0 1 1 0 0 0 0
1
0 0 2 1 0 0 0 0
T =
0 0 0 0 1 12 0 0
0 0 0 0 1 12 0 0
0 0 0 0 0 0 1 21
0 0 0 0 0 0 1 21
(9.49)
Ds lors,
S = T 1 S
(9.50)
d
S = AS
dx
d 1
T S = T 1 A S
dx
S = T 1 AT S
dx
(9.51)
(9.52)
(9.53)
La consquence du changement de variables est le dcouplage des quations, cest--dire quil est
possible dexprimer les courants uniquement en fonction des tensions. Ds lors, en remplaant Y par
jC, on obtient
I1M
d
I1L = j
dx I2M
4
I2L
a11
a21
a31
a41
a22
a23
a24
a33
a34
a44
=
=
=
=
=
=
=
=
=
=
a11
a21
a31
a41
a12
a22
a32
a42
a13
a23
a33
a43
a14
V1M
V1L
a24
a34 V2M
a44
V2L
(9.54)
(9.55)
(9.56)
(9.57)
(9.58)
(9.59)
(9.60)
(9.61)
(9.62)
(9.63)
(9.64)
Le paramtre a31 de ces quations dtermine le couplage entre la tension sur une paire et le courant
induit dans lautre paire. Ce terme serait nul si toutes les capacitances taient gales par couple. Cest
153
ce dsquilibre de capacits qui est lorigine dun effet diaphonique ; il sexprime par 4
d
j
I2M = CM1 M2V1M
dx
4
(9.65)
o CM1 M2 est gal a31 (ou a13 ). Cette quation dtermine donc le courant lvolution du courant
perturbateur dans la paire perturbatrice d une tension dans la paire perturbatrice. Les termes a12 et
a21 sont eux lis au mode commun.
Modle des inductances non quilibres
Le modles des inductances non quilibres est similaire au prcdent. Les circuits sont reprsents par une srie deffets inductifs tels que repris la figure 9.18.
Zc
I1m
V1m
Zc
x M3
x M1
x M2
Zc
I2m
x M4
Zc
V2m
(9.67)
I2M =
CM1 M2 V1M
dx
4Zc2
4
4
Lquation 9.65 nglige le terme de diaphonie d la tension V1L de mode commun car, lors de linstallation des
lignes, loprateur sarrange pour avoir un rapport de mode diffrentiel mode commun de 35 55 [dB].
154
o QM1 M2
d
I2M = j QM1 M2V1M
dx
tient compte des effets capacitifs et inductifs.
(9.68)
Modem
Paire perturbe
NEXT
Paire perturbatrice
0
Modem
x1
V0 ( f )e ( f )x1
V0 ( f )
Commutateur
de donnes
(9.69)
Le courant diaphonique induit sur la paire perturbe sobtient en exploitant la relation 9.68
d
I2 ( f , x1 ) = j QM1 M2 V1 ( f , x1 ) = j QM1 M2 V0 ( f )e ( f )x1
dx
(9.70)
Pour que ce signal perturbe la seconde paire, il faut quil parvienne son extrmit. lorigine, le
courant induit vaut
d
d
I2 ( f , 0) = I2 ( f , x1 )e ( f )x1 = j QM1 M2 V0 ( f )e2 ( f )x1
dx
dx
(9.71)
Pour obtenir la totalit du signal perturbateur sur la seconde paire, il suffit dintgrer sur lensemble
des positions x1 , cest--dire sur toute la longueur de la zone dinterfrence L. Ds lors,
I2 ( f ) =
Z L
0
(9.72)
Cette expression reprsente donc le courant gnr dans la paire perturbe d une tension dans
la paire perturbatrice, ramene lentre du cble. En terme de puissance, le signal perturbateur
sexprime par
P2 ( f ) = V2 ( f )I2 = ZL I2 ( f )I2( f )
(9.73)
155
o ZL est limpdance de la charge qui termine la ligne. La rsolution de lintgrale sous forme
analytique nest malheureusement pas possible car la quantit QM1 M2 nest pas connue ; on peut donc
la traiter comme une variable alatoire et ds lors rechercher lexpression dune puissance moyenne.
En supposant que la charge ZL est adapte limpdance de ligne Zc , il y a donc uniquement un transit
de puissance active. Ainsi, la puissance moyenne est donne par lexpression suivante
Z L
Z L
2
2 ( f )x
2 ( f )y
dx
j Q M1 M2 (y)e
E {P2 ( f )} = E RLV0 ( f )
j QM1 M2 (x)e
dy
0
0
Z L Z L
2 2
2 ( f )x 2 ( f )y
dxdy
(9.74)
e
QM1 M2 (x)QM1 M2 (y)e
= RL V0 ( f )E
0
= RL 2V02 ( f )
Z LZ L
0
(9.75)
Une hypothse raliste consiste considrer que QM1 M2 (x) et QM1 M2 (y) sont deux variables non corrles, ce qui implique que la variable alatoire QM1 M2 (x)QM1 M2 (y) naura une valeur non nulle que
si x = y do
E QM1 M2 (x)QM1 M2 (y) = k (x y)
(9.76)
Cette hypothse est un artifice de calcul commode mais son interprtation physique est malaise.
Notons nanmoins que QM1 M2 (x) peut prendre des valeurs ngatives ; lesprance peut donc sannuler.
Grce cette hypothse,
E {P2 ( f )} = RL 2V02 ( f )
=
Z LZ L
RL 2V02 ( f )k
Z L
0
k (x y)e2 ( f )x e2
e2 ( f )x e2
( f )x
( f )y
dxdy
dx
(9.77)
Z L
e4 ( f )x dx
#L
"
4 ( f )x
e
= RL 2V02 ( f )k
4 ( f )
0
RL 2V02 ( f )k 4 ( f )L
1
=
e
4 ( f )
=
RL 2V02 ( f )k
(9.78)
(9.79)
(9.80)
RL 2V02 ( f )k
4 ( f )
haute frquence (ce qui est le cas pour lADSL), on sait que est proportionnel
HNEXT ( f ) =
3
RL 2 k
k f 2
= = KNEXT f 2
4 ( f ) 0 f
(9.81)
f . Do
(9.82)
Modem
j QM1 M2 V0 ( f )e ( f )L
Paire perturbe
FEXT
Paire perturbatrice
0
Modem
x1
V0 ( f )e ( f )x1
V0 ( f )
Commutateur
de donnes
(9.83)
(9.84)
o E ( f ) kHc ( f )k2 reprsente la densit spectrale de puissance du signal utile la sortie de la paire
et P ( f ) HFEXT ( f ), la densit spectrale de puissance due au FEXT lentre du rcepteur.
paire perturbe
PR
PE
FEXT
PP
paire perturbatrice
E ( f ) kHc ( f )k2
(9.85)
(9.86)
(9.87)
(9.88)
Une formule identique peut tre drive dans le cas du NEXT, ce qui donne deux expressions
du rapport signal sur bruit. On peut constater que lexpression du canal nintervient pas dans ces
formules.
(9.89)
i=1
Si, en plus, on fait lhypothse que tous les signaux perturbateurs ont la mme densit spectrale de
puissance,
(9.90)
Cette formule survalue la puissance perturbatrice car les paires trs loignes de la paire perturbe
influencent peu cette dernire. Le facteur considrer est ds lors infrieur N. Il est donc ncessaire
de mener une tude statistique. U NGER a propos la formule empirique suivante
(9.91)
Supposons que la source mette un symbole. Si une des deux probabilits est unitaire, il ny a
pas de surprise la sortie et donc aucune information (on connat le symbole qui va tre mis). Dans
le cas o lmission des symboles est quiprobable (p0 = p1 = 0, 5), une grande incertitude subsiste
sur le symbole mis et la rception du symbole mis apporte beaucoup dinformation. On peut ainsi
caractriser linformation par la probabilit dmission des symboles.
Dfinition 39 [Information dun symbole] Ds lors, le gain en information apport par lobservation de lvnement S = sk , de probabilit pk est donn par linformation de symbole telle que dfinie
par
1
= log2 pk
(9.92)
i(sk ) = log2
pk
Linformation est nulle quand pk vaut 1 ce que lon peut interprter comme le fait que la ralisation dun vnement certain napporte aucune information et elle augmente mesure que pk
diminue. Cette grandeur sexprime en bit dinformation. Il convient de ne pas confondre le bit et
le bit dinformation qui reprsentent respectivement les mesures du dbit et du dbit dinformation.
Ainsi, il est tout fait possible davoir un nombre de bits dinformation qui ne soit pas un nombre
entier. En pratique, on ne retrouve gure cette nuance puisque le terme bit est utilis pour les deux
notions. Il faudra donc recourir au contexte pour rtablir la distinction.
Notion dentropie
Soit une source S comprenant K symboles. On cherche caractriser linformation moyenne
fournie par chaque symbole mis.
Dfinition 40 Cette grandeur est appele entropie de la source et note H(S). Ds lors,
K1
H(S) =
pk log2
k=0
1
pk
(9.93)
Dans le cas dune source binaire, on peut considrer deux cas particuliers :
1. le cas o tous les symboles sont quiprobables. On montre alors que H(S) = 1, ce qui implique
que linformation apporte lors de lobservation du symbole vaut un bit.
2. le cas o une probabilit pk vaut 1. Alors H(S) = 0 ; linformation apporte lors de lobservation
du symbole est nulle car on sait a priori quel symbole va tre mis par la source.
x0
y0
y1
x1
Canal : p(yk |x j )
yK1
xJ1
H(X |Y = yk ) =
j=0
1
p(x j |yk )
(9.94)
Cette formule nous donne ds lors lincertitude subsistant sur la valeur de lentre connaissant la
valeur yk de la sortie. Il serait intressant davoir une moyenne de cette information pour toutes les
valeurs de Y et donc pour toutes les sorties possibles, cest lentropie conditionnelle H(X |Y ), donne
par
K1
H(X |Y ) =
=
H(X |Y = yk )p(yk )
k=0
K1 J1
k=0 j=0
K1 J1
(9.95)
p(x j , yk ) log2
k=0 j=0
1
p(x j |yk )
1
p(x j |yk )
(9.96)
(9.97)
On dfinit alors linformation mutuelle moyenne I(X ;Y ) comme la quantit dinformation, et donc
dentropie, que lon retrouve en sortie du canal (et donc qui na pas t affect par le bruit). Ds lors,
I(X ;Y ) = H(X ) H(X |Y)
(9.98)
2. le cas dun canal fortement bruit au point que tous les symboles soient quiprobables. Alors
H(X |Y ) = H(X ), car lanalyse de la sortie ne permet pas dobtenir la moindre information sur
lentre. Ainsi donc, I(X ;Y ) est nulle ; le canal ne transmet aucune information.
On peut montrer que linformation mutuelle moyenne possde notamment la proprit dtre
symtrique, ds lors I(X ;Y ) = I(Y ; X ), on en dduit que
I(X ;Y ) = H(X ) H(X |Y)
I(Y ; X ) = H(Y ) H(Y |X )
(9.99)
(9.100)
(9.101)
Cette dfinition reprsente donc la meilleure utilisation que lon peut faire du canal et elle est indpendante de la distribution de probabilits des symboles de lalphabet dentre. Il sagit videmment
dune capacit thorique que lon pourra utiliser comme rfrence pour la mise au point de technique
de transmission.
Remarquons enfin que la relation 9.101 dfinit une capacit de canal par symbole. Si le canal est
capable de traiter s symboles par seconde, la capacit par seconde est donne par
C = sCs
(9.102)
Y
p
x0 = 0
p(x1 ) = 1
x1 = 1
y0 = 0
y1 = 1
(9.103)
(9.104)
(9.105)
Calculons lentropie conditionnelle H(Y |X ). Pour ce faire, appliquons la formule 9.96 adapte.
Pour la paire entre-sortie (x0 , y0 ), on obtient
H(Y |X )(x0,y0 ) = p(x)0 p(y0 |x0 ) log2 (p(y0 |x0 ))
(9.106)
(9.107)
En effectuant de mme pour toutes les combinaisons entre-sortie ((x1 , y0 ), (x0, y1 ) et (x1 , y1 )) et
en sommant les diffrentes contributions, on trouve
1
H(Y |X )(xi,y j )
(9.108)
= p log2 p q log2 q
(9.109)
H(Y |X ) =
i=0 j=0
On constate donc que H(Y |X ) ne dpend que des caractristiques du canal, on parle dentropie de
canal. Ds lors,
I(X ;Y ) = H(Y ) + p log2 p + q log2 q
(9.110)
Pour maximiser cette quantit, on ne peut plus jouer que sur H(Y ). Le maximum que lon peut
obtenir avec un canal binaire est une valeur unitaire pour H(Y ). La capacit maximale de canal est
alors
Cs = 1 + p log2 p + q log2 q
(9.111)
Sur le schma de la figure 9.23, q reprsente la probabilit de prendre une dcision errone ; cest
la probabilit derreur. Dans le cas dune modulation deux tats, la probabilit derreur est lie au
rapport signal bruit par la relation
r
1
Eb
q = pe = er f c
(9.112)
2
N0
Lvolution de la capacit de canal en fonction du rapport signal sur bruit est reprsente la
figure 9.24. On constate que, lorsque le rapport signal sur bruit tend vers linfini (cest--dire que la
probabilit derreur tend vers 0), la capacit maximale du canal tend vers 1 bit.
Influence du type de modulation
Prenons le cas dune modulation quatre tats (4-AM par exemple). La matrice de canal se complique car il faut calculer plusieurs probabilits derreur pour un symbole dentre donn.
Cependant, lexpression de la capacit de canal reste de la forme
Cs = H(Y ) + Hc
(9.113)
o la valeur maximale de H(Y ) est de deux, la consquence dune modulation quatre tats tant
que lon a deux bits par symbole et o Hc est un terme associ aux transitions de canal. Il est alors
intuitivement acceptable que la capacit de canal tende vers deux quand le rapport signal bruit tend
vers linfini. Le raisonnement peut tre rpt pour des modulations 16 voire 64 tats. Les rsultats
sont rsums la figure 9.24.
162
Cs (bit/symbole)
16-AM
p = 105
8-AM
p = 105
4-AM
p = 105
2-AM
p = 105
10
20
30
Eb
N0 [dB]
X2 .
Thorme 42 [S HANNON] Des calculs montrent que la capacit dun canal blanc additif gaussien
de moyenne nulle et de variance N2 sexprime par
1
X2
Cs = log2 1 + 2
2
N
o
X2
N2
(9.115)
Cette formule porte le nom de second thorme de S HANNON. On constate que, dans le cas dun
bruit important (N2 X2 ), le canal ne parvient pas transmettre de linformation. La figure 9.25
compare la limite thorique de S HANNON avec les capacits maximales atteintes par les diffrentes
techniques de modulation. Il est intressant de noter que la limite thorique de S HANNON correspond
au cas dune modulation comportant un nombre infini dtats.
163
Cs (bit/symbole)
16-AM
p = 105
limite de S HANNON
8-AM
3
p = 105
4-AM
p = 105
2-AM
p = 105
10
20
30
164
Eb
N0 [dB]
Chapitre 10
Ingnierie des radiocommunications mobiles
terrestres
Ce chapitre aborde ltude de la partie radio des rseaux de radiocommunications mobiles terrestres. Les aspects systmes ne seront pas couverts.
10.1 Introduction
La figure 10.1 reprsente les lments qui interviennent dans le calcul du bilan de puissance
dune liaison sans fil. On distingue les pertes dmission (LE ) et de rception (LR ), dues aux lignes
de transmission reliant les quipements aux antennes ainsi qu certains quipements intermdiaires
(duplexeurs, connecteurs, . . .), les gains dantenne et laffaiblissement dans lespace (Ae ).
GE
Ae
GR
LE
LR
metteur
Rcepteur
F IG . 10.1 lments intervenant dans le calcul du budget de puissance dune liaison sans fil.
Ds lors, la puissance de rception PR , exprime en [dB], vaut
PR = PE LE + GE Ae + GR LR
(10.1)
fait, cette situation est rare dans le contexte dune radiocommunication mobile.
165
Eb
N0
Eb
N0
(10.2)
Eb
C=
+N
N0 seuil
(10.3)
Comme en GSM, la largeur de modulation est de 271 [kHz], le bruit thermique est gal 1, 1
1012 [mW ], soit 120 [dBm] une temprature de 290 [K].
Ce rapport est calcul au droit de la sortie de ltage dentre (typiquement aprs dmodulation).
Comme lamplificateur dentre amplifie le bruit dun facteur typique de 10 [dB], le bruit de fond est
born par 110 [dBm]. Il en dcoule que la sensibilit du mobile vaut 102 [dBm]. Le tableau 10.1
reprend les quelques valeurs de sensibilit pour diffrents types de rcepteurs.
2 En
166
Type de rcepteur
Station de base
Portable 8 [W ]
Portable 2 [W ] (GSM 900)
GSM bi-bande
Sensibilit en [dBm]
-104
-104
-102
-102
167
P1 [dB]
35 dB / dcade
10
100
1000
10000
d [km]
10
100
1000
10000
d [km]
F IG . 10.3 volutions de la puissance en fonction de la distance pour une mme configuration apparente.
168
Dautre part, en pratique, il suffit de parcourir un autre rayon du cercle centr sur lantenne dmission pour observer une autre courbe daffaiblissement, la configuration de lenvironnement tant diffrente.
10.2.2 Modle
Dans les cas les plus courants, le rcepteur et lmetteur ne sont pas en visibilit directe. On
approxime frquemment laffaiblissement du canal par un modle 3 tages comprenant :
1. une attnuation mdiane due la distance,
2. un terme alatoire prenant en compte les effets de masque, dus la prsence dobstacles, et
3. un autre terme alatoire dcrivant les vanouissements.
10.2.4 Masquage
Base physique
Le masquage provient de la prsence dobstacles. En premire approximation, on peut considrer quun obstacle absorbe un certain pourcentage de la puissance. Ainsi, la puissance transmise ou
diffracte ne reprsente quun certain pourcentage de la puissance incidente. De mme, la prsence
dune vitre se traduit par une certaine diminution de puissance.
Pour modliser les effets de masquage (shadowing en anglais), aussi appel vanouissement lent,
nous partons de lhypothse suivante. Supposons que toutes les contributions de masquage laffaiblissement A, en terme de puissance, soient dues des facteurs multiplicatifs A1 , A2 , . . . , AN , tous
infrieurs 1, tels que
A = A1 A2 . . . AN
(10.4)
169
(10.5)
(10.6)
(10.7)
Si toutes les contributions sont des variables alatoires de mmes esprance et variance, L est une
variable alatoire normale par application du thorme de la limite centrale :
L [dB] = L50% [dB] + s [dB] N(0, 1) = N(L50% , s2 ) [dB]
(10.8)
(10.9)
(10.10)
V = 10s [dB]X/10
(10.11)
o Lo = 10L50% [dB]/10 et
est une loi log-normale.
Densit de probabilit dune loi log-normale. Pour calculer la densit de probabilit de la lognormale, on part de lexpression de la gaussienne. Soit X une variable alatoire gaussienne centre,
de variance unitaire. Sa densit de probabilit vaut, pour x [, +],
x2
1
fX (x) = e 2
2
(10.12)
s X
10
= e s X
(10.13)
ln(10)
= 0, 23
10
La fonction de rpartition de V vaut [22, page 152]
(10.14)
Z ln v/ s
1
fX (x)dx =
2
Z ln v/ s
ln v
)
s
x2
e 2 dx
(10.15)
(10.16)
Cette dfinition, courante dans la littrature scientifique, peut prter confusion. En effet, en toute rigueur A
reprsente linverse de laffaiblissement. Par la suite, nous adopterons nanmoins la dfinition 10.4.
170
fV (v) =
2
1 x2
e
fX (x)
= 2 x
| v/ x| x= ln v
s e s
(ln v)2
2 2
2 s
1 1
s v 2 e
(10.17)
x= lnv
si v 0
si v < 0
(10.18)
s = 1 [dB]
s = 2 [dB]
s = 6 [dB]
0.8
0.6
fV (v)
0.4
0.2
0.5
1.5
2.5
v
F IG . 10.4 Densit de probabilit dune loi log-normale.
La valeur mdiane de V , v50% , est tire de la valeur mdiane de X par v50% = ex50% = 100 = 1
quelle que soit la valeur de s .
En pratique, tous les termes ne contribuent pas de la mme manire laffaiblissement. Dautre
part, on peut montrer que les termes de diffraction nobissent pas une simple loi additive, en
raison dune dpendance entre certains termes. Nanmoins, diverses mesures montrent que la loi lognormale est gnralement suffisante pour modliser le masquage.
Impact sur la zone de couverture en bordure dune cellule
En prsence dun effet de masquage, laffaiblissement L est une variable alatoire normale
N(L50% , s2 ) (cf. quation 10.8) :
L [dB] = N(L50% , s2 ) [dB]
(10.19)
(10.20)
L50% est la valeur mdiane5 de laffaiblissement de parcours telle que fournie par les modles
empiriques rsultant de mesures. En effet, ces mesures fournissent une valeur mdiane qui,
outre laffaiblissement de parcours, englobe un effet moyen de masquage. Cest la raison pour
laquelle on classe les modles empiriques en fonction de la densit de lenvironnement urbain.
Ls , la variable alatoire de masquage ; cest--dire une variable alatoire gaussienne de moyenne
nulle et dcart-type s .
La densit de probabilit de Ls est donc
f (ls ) =
s 2
ls2
2s2
(10.21)
Pour des raisons defficacit, il est ncessaire de garantir un niveau de puissance suprieur la
valeur mdiane ; lingnieur de rseau na gure le choix que dajouter une marge supplmentaire de
puissance lmission afin daugmenter la couverture. lajout de cette marge de ms [dB] correspond
une probabilit de couverture suprieure celle de la puissance mdiane pour une mme sensibilit
de rcepteur. Cette probabilit se calcule par
p(ls < ms ) =
Z ms
s 2
ls2
2s2
dls
(10.22)
2
ls
ms
1
1
2
e s2 dls
+
=
2
0 s 2
1 1
ms
=
+ erf
2 2
2s
Z x
0
et dt
(10.23)
(10.24)
(10.25)
Leffet de masque reprsente une variation sur les conditions de propagation, tantt favorable (par
exemple visibilit directe), tantt dfavorable (par exemple lors de la prsence dun obstacle important
entre lmetteur et le rcepteur). En environnement urbain, lcart type de la loi a une valeur typique
de 6 [dB] (daprs [20]).
La figure 10.5 montre la probabilit de couverture (en %) en fonction de la marge additionnelle
de puissance.
10.2.5 vanouissement
Pour obtenir la statistique qui dcrit le phnomne dvanouissement, il faut distinguer deux cas,
suivant quil y ait ou non visibilit directe ou du moins, un trajet nettement moins attnu que les
autres entre lantenne de la station de base et lantenne du mobile.
vanouissement en labsence dun trajet direct
Considrons une srie de trajets dus la diffraction et la rflexion (cf. figure 10.6). Dans un
premier temps, nous considrons le cas dun metteur et dun rcepteur qui ne sont pas en ligne de
vue.
5 La
valeur mdiane dune variable alatoire normale est thoriquement gale sa moyenne. Pour les mesures, le
recours la mdiane est prfrable si lon ne dispose pas dun large ensemble de ralisations car elle limine limpact des
valeurs aberrantes.
172
10
12
14
Marge en [dB]
F IG . 10.5 Pourcentage de couverture tel que dfini par un seuil de puissance.
A cos(2 fot)
x(t)
173
Lmetteur envoie un signal cosinusodal A cos(2 fot) avec une amplitude et une frquence fo
constantes. Le signal reu au droit du rcepteur X (t) est compos dune srie dondes dphases et
attnues
X (t) = Ci cos(2 fot + i )
(10.26)
i
(10.27)
(10.28)
(10.29)
(10.30)
Le canal ntant pas invariant dans le temps (par exemple en raison de dplacements), les composantes en phase et en quadrature sont des fonctions temporelles. Dautre part, X (t) est la ralisation
dun processus stochastique qui est la somme de termes indpendants. Ds lors, en vertu du thorme
de la limite centrale, XI (t) et XQ (t) reprsentent des variables alatoires gaussiennes.
Si lon veut dterminer lintensit de champ lectrique, il faut analyser lvolution de lamplitude
du signal. Cette amplitude est dfinie par6
q
R(t) = XI2 (t) + XQ2 (t), R(t) 0
(10.31)
(t) = tan1
et donc
XQ (t)
, (t) [0, 2 [
XI (t)
(10.32)
(10.33)
Pour exprimer les densits de probabilit de ces variables alatoires partir de XI (t) et XQ (t),
supposons tout dabord que les composantes en phase et en quadrature soient des variables gaussiennes centres, de mme variance X2 et non corrles. Les variables sont centres car elles rsultent
de la somme de variables alatoires de type Ci cos i . Or, en supposant que les variables Ci et cos i
soient indpendantes, E {Ci cos i } = E {C}i E {cos i } = 0 si lon considre un argument uniformment distribu sur lintervalle [0, 2 ]. Physiquement, une moyenne correspond des contributions
damplitude ngative, plus prcisment des termes en dcalage de phase de 900 2700.
Densit de probabilit. Sagissant de gaussiennes, la non corrlation entrane lindpendance [7,
page 162]. La densit de probabilit conjointe est ds lors gale
1
e
fXI XQ ( , ) =
2X2
2 + 2
2
2X
(10.34)
Le changement de variables
xI = r cos
xQ = r sin
6
(10.35)
(10.36)
Les rsultats suivants rsultent de la dcomposition en composantes de R ICE et de la reprsentation par un passe-bas
quivalent dun signal passe-bande.
174
permet de trouver la densit de probabilit conjointe fR, (r, ). En effet, on sait quen toute gnralit,
un changement de variables fait intervenir le Jacobien7
XI XQ
(10.37)
fR,(r, ) = XR XRQ fXI XQ (r cos , r sin )
I
Ds lors,
cos
sin
fR, (r, ) =
r sin r cos
fX X (r cos , r sin )
I Q
(10.38)
r
2
2X
=
e
2
2
X
2
r 2r 2
e X , r 0, [0, 2 [
=
2X2
0
r<0
(10.39)
(10.40)
Densit de probabilit marginale de lenveloppe. La densit de probabilit marginale des amplitudes sobtient en intgrant la densit conjointe sur . Un calcul simple conduit
2
r 2r 2
e X, r0
fR (r) =
(10.41)
2
X
0
r<0
Cest la densit de probabilit de R AYLEIGH. Elle est illustre, sous forme normalise, la figure 10.7.
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
175
pour n = 2k + 1
pour n = 2k
(10.42)
En particulier,
R2
/2 X 1, 253 X
2
= 2
X 0, 429 X2
2
E {r} =
r2
r
2X
dr
=
Rrms =
r2
e
2X 1, 41 X
2X2
0
(10.43)
(10.44)
(10.45)
On parle gnralement dvanouissement lorsque lamplitude du signal passe en dessous de Rrms [22].
Pour obtenir la densit marginale f ( ), on intgre
1
, [0, 2 [
2
ce qui signifie que est une variable alatoire uniforme et indpendante de R, do
f ( ) =
(10.46)
fR,(r, ) = fR (r) f ( )
(10.47)
(10.48)
Ce cas est videmment plus favorable que le cas prcdent puisque le rcepteur reoit galement le
signal en ligne de vue. Ltude mene ci-aprs vise mesurer limportance de lajout du signal direct.
Il ne faudrait pas pour autant en dduire des rgles de dimensionnement car on doit tenir compte des
conditions les plus dfavorables pour ltablissement du bilan de puissance, cest--dire en labsence
dun trajet direct.
Pour un signal Z(t) bande troite, il est possible de recourir la dcomposition de R ICE
Z(t) = ZI (t) cos(2 fot) ZQ (t) sin(2 fot)
(10.49)
(10.50)
(10.51)
(10.52)
(10.53)
60
Nombre doccurrences
50
40
30
20
10
0
85
80
75
70
65
60
Niveau du signal [dBm]
55
50
45
177
et
(t) = tan1
ZQ (t)
, (t) [0, 2 [
ZI (t)
(10.54)
Ds lors,
ZI (t) = R(t) cos(t)
ZQ (t) = R(t) sin(t)
(10.55)
(10.56)
Pour une valeur donne de , les variables alatoires gaussiennes XI (t) et XQ (t) tant indpendantes,
elles conservent ce caractre. Il en rsulte que
ZI (t) = N(A cos , X2 )
ZQ (t) = N(A sin , X2 )
(10.57)
(10.58)
( A cos o )2 +( A sin o )2
2
2X
(10.59)
Densit de probabilit.
r
=
e
2
2X
r
e
=
2X2
(10.61)
(10.62)
Z 2
0
fR,(r, |o )d
2
r +A2
r
2X
e
=
2X2
Z 2
(10.63)
rA cos(o )
2
X
2 2
r +A2
Ar
r
2X
e
I
2
=
0
2X2
X2
(10.64)
(10.65)
x2n
2n 2
n=0 2 (n!)
(
2
x4
1 + x4 + 64
, x < 2, 5
e1,1x
2, 5 < x < 10
1,9x ,
I0 (x) =
(10.66)
(10.67)
12
10
Io(x)
0.5
1.5
2
x
2.5
3.5
2 2
r r 2+A
2
X
I0 Ar2 , r 0
e
fR (r) =
(10.68)
X2
X
0
r<0
Il sagit de la densit de probabilit de R ICE ; elle est illustre (sous forme normalise) la figure 10.10. Lorsque A = 0, il sagit de la loi de R AYLEIGH.
0.7
a=0
0.6
a=1
0.5
a=2
a=4
a=8
0.4
0.3
0.2
0.1
10
180
A
X ).
Prenons le cas dune source qui met un champ lectrique sinusodal de frquence fc . En labsence
de mouvement, le champ reu au droit du mobile vaut, suivant la thorie des ondes [26],
E(t) = E0 e j(ct z)
(10.69)
= E0 e
= E0 e
(10.70)
j(c t z V t cos )
j2 ( fc V
c
(10.71)
cos )t j z
(10.72)
V
cos = fm cos
c
(10.73)
appele frquence D OPPLER. Comme V et fluctuent au cours du temps, on assiste une distorsion
de frquence variable et donc de phase pour lanalyse de laquelle il faudra galement recourir des
outils statistiques. La modification de la frquence porteuse dans une place de frquences porte le nom
dtalement D OPPLER. La figure 10.12 montre les spectres lmission et la rception rsultant de
leffet D OPPLER.
Spectre lmission
fc
Spectre la rception
fc
(10.74)
En associant un secteur angulaire une contribution spcifique de ltalement D OPPLER, la puissance correspondante en termes de frquence vaut
PR ( f ) = D ( f ) d f
(10.75)
On remarquera aisment que deux secteurs angulaires dorientation oppose contribuent une mme
frquence. Ds lors,
D ( f ) d f = G( ) f ( ) d + G( ) f ( ) d
(10.76)
Ce qui mne
|D ( f )| =
Par la relation 10.73,
G( ) f ( ) + G( ) f ( )
df
d
(10.77)
d f
= fm | sin |
d
p
= fm 1 cos2
s
2
f
= fm 1
fm
(10.78)
(10.79)
(10.80)
G( )
2r
fm
G( )
2
f
fm
2 =
G( )
r
2
fm 1 ffm
pour | f | fm
(10.81)
et 0 pour | f | > fm . Cette forme de spectre porte le nom de spectre D OPPLER classique ; elle est
illustre la figure 10.13.
0.02
0.015
|D ( f )| 0.01
0.005
0
-100
-50
50
100
f [Hz]
F IG . 10.13 Spectre D OPPLER classique ( fc = 900 [MHz], V = 100 [km/h] et G = 1).
Bien que leur obtention soit mathmatiquement complexe, il existe des expressions des densits
spectrales de la phase du signal, de la frquence instantane, de lenveloppe et dautres quantits
encore. En particulier, lanalyse de lenveloppe montre que sa fonction dautocorrlation se resserre
lorsque la vitesse augmente, ce qui signifie que les enveloppes prises deux temps voisins sont moins
corrles et donc que lvanouissement marque des sauts plus brusques [1].
182
exposant
3,2
3,5
3,8
d1
hb
d2
hm
s
2
2
hb + hm
hb hm
d2 d1 = d
+1
+ 1
d
d
184
(10.86)
2hb hm
d
(10.87)
(10.88)
(10.89)
(10.90)
(10.91)
PE GE GR
4 d
2
h h2
PE GE GR b 4m
d
2
4 hb hm
d
2
(10.92)
(10.93)
Ou encore, en dcibels,
PR [dB] = 10 log10 (PE GE GR ) + 20 log10 hb + 20 log10 hm 40 log10 d
(10.94)
Cette formule montre que laffaiblissement en fonction de la distance est deux fois plus rapide
quen espace libre et quil ne dpend plus de la frquence, si ce nest au travers des gains dantenne.
En pratique, il apparat que le modle simple que nous avons choisi nest pas trs reprsentatif des
situations relles. Lanalyse montre nanmoins certaines tendances que lon retrouve dans certains
modles empiriques.
(10.95)
avec
f la frquence, d la distance, hb , hm , des hauteurs ; ces grandeurs sont exprimes respectivement
en [MHz], [km] et [m].
185
a(hm ) = (1, 1 log( f )0, 7)hm (1, 56 log( f )0, 8) pour une ville de taille moyenne ; ce facteur
de correction dpend de la hauteur de lantenne du mobile mais galement du type denvironnement.
Cm = 0 [dB] pour les villes de taille moyenne et les banlieues, et Cm = 3 [dB] pour les grands
centres mtropolitains.
Nous renvoyons le lecteur des ouvrages spcialiss [20, 22, 33] pour les conditions dapplicabilit
de la formule et la description dautre modles.
Rayon transmis
Rayons diffracts
Rayon rflchi
Rcepteur
Faade
metteur
187
188
Annexe A
Quelques rsultats
A.1 Gain et attnuation
Dfinition 43 [Gain et attnuation dun systme linaire] Pour un systme linaire, le gain en puissance est le rapport entre la puissance de sortie et la puissance en entre, lattnuation tant le rapport
inverse.
On dfinit parfois, par le biais de la rponse impulsionnelle, le carr du rapport des tensions ou
courants de sortie et dentre en [dB]. Nanmoins, kH ( f )k2 nest pas toujours gal au gain en puissance du circuit.
(A.1)
Lexpression temporelle est la partie relle du phaseur pralablement multipli par lexponentielle
imaginaire e j t . En effet,
b j t )
x(t) = Re(Xe
= Re(X e j( t + ))
= X cos( t )
(A.2)
(A.3)
(A.4)
(Xb e j t )
= j Xbe j t
t
et
A.2.1
Impdance complexe
b j t
b j t dt = X e
Xe
j
(A.5)
(A.6)
Soit un circuit formant un diple lectrocintique que lon supposera passif, cest--dire quil ne
comporte pas de gnrateur. Supposons tous ses lments linaires ; toute excitation sinusodale
v(t) = V cos( t) entre ses bornes, correspond un courant i(t) = I cos( t ).
En utilisant la notation phasorielle, limpdance complexe est dfinie par la relation
Vb = ZbIb
(A.7)
Zb est donc de la forme Ze . Son module Z, simplement appel impdance, est gal au rapport des
amplitudes V et I ; il se mesure en ohms [].
A.2.2
Soient v(t) = V cos( t) la tension aux bornes du diple tudi et i(t) = I cos( t ) lintensit
qui le parcourt. Le diple reoit la puissance lectrocintique instantane
(A.8)
v(t)i(t) = V I cos( t) cos( t )
VI
(A.9)
=
(cos + cos(2 t ))
2
Cette puissance varie au cours du temps de faon sinusodale la frquence angulaire 2 , autour
dune valeur Pm , appele puissance moyenne
VI
cos
(A.10)
2
Pm est exactement la moyenne de la puissance reue par le circuit lorsquelle est calcule sur un
nombre entier de priodes.
En notation phasorielle,
b j t
b e j t Re Ie
(A.11)
Pb = Re V
!
!
b e j t + Vb e j t
b j t + Ib e j t
V
Ie
=
(A.12)
2
2
!
!
Vb Ib
Vb Ibe2 j t
= Re
+ Re
(A.13)
2
2
Pm =
Pour un circuit passif, dimpdance complexe Zb = R + jS, la puissance moyenne vaut alors
!
RI 2
(R + jS) kIk2
=
(A.14)
Pm = Re
2
2
(A.15)
On dit alors que limpdance du circuit dutilisation est adapte celle du gnrateur. La puissance
fournie est maximale et vaut
V2
(A.16)
8R
Notons quune puissance gale est dissipe dans le gnrateur, ce qui diminue pour de nombreuses
applications lintrt de cette adaptation.
Pd =
191
192
Index
A
accs
multiple, 105
adaptation dimpdance, 191
ASK, 59
attnuation, 137
B
bande
de base, 25
de garde, 106
utile, 25
bande vocale, 147
biais, 18
bit, 159
dinformation, 159
bruit, 98, 138
bande passante, 47
facteur, 48, 50
temprature, 46, 47, 51
thermique, 43
C
cble, 138
canal, 159
capacit, 103, 161
CDM, 88
CDMA, 88, 105
cellule, 183
charge, 121
coule, 121, 122
dun rseau, 120
de rfrence, 122
offerte, 122
chip, 88
circuit, 120, 147
coaxe, 138
code, 88
de Gold, 116
dtalement, 96
composante
en phase, 28
en quadrature, 28
congestion, 119, 131
consistant, 18
constellation, 60
convergence, 9
cosinus surlev, 85
couverture, 172
critre
de Nyquist, 80
D
dcibel, 143
dnombrement, 123
dtecteur denveloppe, 35
densit
Rayleigh, 175
Rice, 179
DFT, 11
diagramme des tats de phase, 60
diaphonie, 149
dimensionnement, 132
discriminateur, 35
dispersion, 138
distorsion, 138
Doppler, 166, 179
spectre, 182
DS, 88
dtFT, 8
DWDM, 105
E
chantillonnage, 26
effet de fentrage, 23
effet de masque, 169
effet Doppler, 166, 179
efficacit spectrale, 55
nergie, 11
entropie
conditionnelle, 160
entropie de source, 159
enveloppe, 30
enveloppe complexe, 31
193
modulation dcalage, 67
MSK, 72
multiplexage, 105
en frquences, 105, 106
par rpartition de codes, 113
temporel, 105, 108
N
Nper, 143
Nyquist
critre, 80
O
Okumura-Hata, 185
OQPSK, 69
P
priodigramme, 22
paire coaxiale, 138
paire torsade, 138
paquet, 147
paquets, 120
paradiaphonie, 149
paramtres
primaires, 141
Parseval, 11
passe-bande, 25
reprsentation, 31
passe-bas
quivalent, 27
path loss, 169
PCM, 112
phaseur, 189
pico-cellule, 183
Poisson
esprance, 125
processus, 126
variable alatoire, 124, 125
variance, 125
prcision, 12
pseudo-alatoire, 90
PSK, 61
PSK-4, 62
puissance
disponible, 45
Q
QPSK, 62
dmodulation, 66
quarte, 149, 151
R
rsistivit, 139
rsolution, 12, 13
Rayleigh, 175
Rice, 179
dcomposition, 32
rolloff, 85
S
squence
de longueur maximale, 90
pseudo-alatoire, 91
squence pseudo-alatoire, 90
sensibilit, 166
sensibilit du mobile, 166
shadowing, 169
Shannon
thorme, 163
signal analytique, 29
signalisation, 113, 120
source, 159
spectre, 9
damplitude, 9
Doppler, 182
splitter, 148
STP, 145
T
tldiaphonie, 149
TDD, 116
TDMA, 105
tlgraphistes
quations, 141
temprature
de bruit, 47, 51
de bruit ponctuelle, 46
thorme de Shannon, 163
trafic
coul, 121
tude, 119
intensit, 121
trame, 110
transforme de Hilbert, 30
transmission sur ligne, 137
treillis de phase, 73
twisted pair, 145
W
WDM, 105
195
Glossaire
ADSL Asymmetric Digital Subscriber Line. Technologie permettant la transmission de signaux
numriques haut dbit (jusqu plusieurs Mgabits par seconde) sur une paire torsade. Le
dbit est asymtrique ; il est plus important du central vers labonn quen retour. Le systme
ADSL est compatible avec le signal tlphonique. . . . . . . . . . . . . . . . . . . . . 105, 142, 148
BPSK Binary Phase Shift Keying. Modulation numrique de phase 2 tats. Il sagit dun cas particulier des modulations numriques de phase plusieurs tats (PSK). . . . . . . . . . . . . . . . . 55
CDMA Code Division Multiple Access. Technologie de transmission numrique permettant la transmission de plusieurs flux simultans par rpartition de code. Cette technologie permet une utilisation permanente de la totalit de la bande de frquences alloue lensemble des utilisateurs.
La technologie prvoit un mcanisme daccs aux ressources. . . . . . . . . . . . . . . . . . . . 105
ellule En radiocommunications, zone gographique lmentaire dun rseau radiocellulaire
laquelle on affecte un ensemble de frquences non rutilisables dans les zones contigus. Cest
galement le nom donn un paquet ATM qui a une taille de 53 bytes dont 48 sont destines
recevoir les donnes dun utilisateur. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 183
DFT
Discrete F OURIER Transform. La transforme de F OURIER discrte sobtient par chantillonnage de la transforme continue. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
Duplex Terme utilis en tlcommunications pour dsigner une communication bidirectionnelle.
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 116
FDM
GSM
Global System for Mobile Communications. Standard de tlphonie mobile adopt en Europe,
en Asie et en Australie. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 116
IS-95 Norme amricaine de rseau cellulaire (dit de seconde gnration ou 2G) base sur la mthode daccs CDMA. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 102
LOS
Line Of Sight. Terme qui dsigne quun metteur et un rcepteur sont en ligne de vue, cest-dire que, parmi tous les trajets reliant lmetteur et le rcepteur, un trajet est nettement moins
attnus que les autres. La visibilit en ligne de vue nexclut pas les trajets multiples, contrairement la visibilit en espace libre. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 165
modem modulateur dmodulateur. Appareil transmettant des signaux numriques sur le rseau tlphonique analogique. Offre les fonctions de numrotation, de connexion, et ventuellement de
compression et de correction derreur. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 147
modulation Technique consistant modifier lamplitude, la phase ou la frquence dune onde
sinusodale, appele porteuse, au moyen dun signal transmettre, appel signal modulant. Grce
la modulation, on peut translater le contenu frquentiel dun signal modulant ; ce procd permet
de partager le spectre de frquences entre plusieurs utilisateurs. Pour retrouver le signal modulant
original, il faut procder une dmodulation. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
196
MSK
Minimum Shift Keying. Technique de modulation numrique consistant effectuer une fonction
XOR entre 2 bits successifs pralablement une modulation de frquence 2 tats. . . . . . . 55
multiplexage Terme technique utilis en tlcommunications pour dsigner un procd qui consiste partager des ressources entre plusieurs utilisateurs. Dans un autre contexte et lorsquon
parle de signaux, il sagit dune manire de combiner plusieurs signaux. . . . . . . . . . . . . 105
NEXT Near-End Cross(X) Talk. Mesure de la paradiaphonie dun cble de cuivre. . . . . . . . . 150
NRZ Non Return to Zero. Mthode de reprsentation de signaux numriques. Cette mthode
reprsente les symboles numriques par des niveaux de tension constants pendant toute la dure
du symbole. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55, 57
OSI
PSK
Phase Shift Keying. Modulation numrique de phase plusieurs tats. Technique de modulation consistant slectionner des chantillons dune porteuse amplitude constante mais avec
plusieurs tats de phase possibles. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
QPSK Quadrature Phase Shift Keying. Technique de modulation numrique consistant appliquer
une modulation des tats de phase tant sur une porteuse en cosinus que sur la porteuse en sinus. 55
RNIS Rseau Numrique Intgration de Services. Dsigne le rseau tlphonique numrique. Au
niveau du rseau, les signaux numriques utiles sont transmis des multiples de 64[kb/s]. 147
TDM
Time Division Multiplexing. Mcanisme de rpartition de ressources par multiplexage temporel. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105,
108
TDMA Time Division Multiple Access. Technique de rpartition de ressources par multiplexage temporel. Cette technique prvoit un mcanisme daccs. . . . . . . . . . . . . . . . . . . . . . . . . 105
UTP
Unshielded Twisted Pair. Dsigne un cble non blind compos de 4 paires en cuivre. La norme
EIA-TIA 586 dfinit 5 catgories de cbles de ce type. . . . . . . . . . . . . . . . . . . . . . . . . 145
197
198
Bibliographie
[1] J. Anderson. Digital transmission engineering. IEEE Press, Prentice Hall, 1999. 182
[2] J. Bellamy. Digital telephony. John Wiley & Sons, third edition, 2000.
[3] J. Bic, D. Duponteuil, et J. Imbeaux. Elements of digital communications. John Wiley & Sons,
1991. 55
[4] J. Bingham. ADSL, VDSL, and multicarrier modulation. Wiley, 2000.
[5] G. Blanchet et M. Charbit. Traitement numrique du signal : simulation sous matlab. Herms,
1998. 6, 9, 12, 15, 16, 23
[6] J. Candy. Signal processing : The modern approach. McGraw-Hill, 1988. 6, 23
[7] M. Capinski and E. Kopp. Measure, integral and probability. Springer, 1999. 174
[8] M. Charbit. lments de thorie du signal : les signaux alatoires. Ellipses, 1990. 5
[9] W. Chen. DSL Simulation techniques and standards development for digital subscriber line
systems. Macmillan Technical Publishing, 1998.
[10] P. Clerc et P. Xavier. Principes fondamentaux des tlcommunications. Ellipses, 1998. 114
[11] P. Delogne. Estimation spectrale de squences alatoires. Chapitre 7 du cours de traitement
numrique du signal, Universit catholique de Louvain, 40 pages, Mai 1994.
[12] C. Demoulin and M. Van Droogenbroeck. Principes de base du fonctionnement du rseau gsm.
Revue de lAIM, (4) :318, 2004. 87
[13] P. Fraisse, R. Protire, et D. Marty-Dessus. Tlcommunications 1 : transmission de linformation. ellipses, 1999. 146
[14] D. Ginsburg. Implementing ADSL. Addison-Wesley, 1999.
[15] M. Hayes. Statistical digital signal processing and modelling. John Wiley & Sons, 1996. 5, 6,
23
[16] S. Haykin. Communication system. John Wiley & Sons, 1994. 39, 93, 94, 95, 101, 116
[17] S. Keshav. An engineering approach to computer networking : ATM networks, the Internet, and
the telephone network. Addison-Wesley, 1997. 116
[18] D. Kingdom. Frequency-division multiplexing. In J. Flood and P. Cochrane, editors, Transmission systems, volume 27 of IEE Telecommunications, chapter 5, pages 110149. Peter Peregrinus
Ltd., 1991. 109
c
[19] X. Lagrange. Les r
seaux
radiomobiles. Herms, 2000. 183
[20] X. Lagrange, P. Godlewski, et S. Tabbane. Rseaux GSM-DCS. Herms, troisime edition, 1997.
167, 172, 186
[21] F. Lapierre. Estimation des seuils de design pour la couverture Indoor. Travail de fin dtudes,
Universit de Lige, Juin 2000. 177
199
[22] J. Lee and L. Miller. CDMA systems engineering handbook. Artech House, 1998. 87, 131, 132,
134, 170, 176, 186
[23] W. Lee. Mobile communications design fundamentals. John Wiley & Sons, second edition,
1993.
[24] H. Meyr, M. Moeneclaey, and S. Fechtel. Digital communication receivers : synchronisation,
channel estimation and signal processing. John Wiley & Sons, 1998. 27
[25] A. Papoulis. Probability, random variables, and stochastic processes. McGraw-Hill, 1984. 175
[26] C. Paul and S. Nasar. Introduction to electromagnetic fields. McGraw-Hill, 1987. 181, 184
[27] D. Rauschmayer. ADSL/VDSL principles : a practical and precise study of asymetric digital
subscriber lines and very high speed digital subscriber lines. Macmillan Technical Publishing,
1999.
[28] T. Saaty. Elements of queueing theory with applications. Dover, 1961.
[29] S. Saunders. Antennas and propagation for wireless communication systems. John Wiley &
Sons, 1999.
[30] C. Servin. Tlcoms 2 : de lingnierie aux services. Interditions, 1998. 122, 145
[31] M. Simon and M.-S. Alouini. Digital communication over fading channels : a unified approach
to performance analysis. Wiley, 2000. 179
[32] T. Starr, J. Cioffi, and P. Silverman. Understanding digital subscriber line technology. Prentice
Hall, 1999.
[33] R. Steele and L. Hanzo. Mobile radio communications : second and third generation cellular
and WATM systems. John Wiley & Sons, second edition, 1999. 186
[34] J. Walrand and P. Varaiya. Communication networks : a first course. McGraw-Hill, 1998. 137
[35] R. Yates and D. Goodman. Probability and stochastic processes : a friendly introduction for
electrical and computer engineers. John Wiley & Sons, 1999. 125
[36] R. Ziemer and W. Tranter. Principles of communications : systems, modulation and noise. John
Wiley & Sons, fourth edition, 1995. 30, 44
200