TP En206 2011 2012 PDF
TP En206 2011 2012 PDF
TP En206 2011 2012 PDF
pratiques de
Communications
Numriques
ENSEIRB-MATMECA E2
Cette partie commence par ltude de la modulation de phase BPSK (Binary Phase Shift Keying). Puis, ltude
est gnralise aux autres modulations de phase, aux modulations damplitude et aux modulations
damplitude en quadrature, et enfin aux modulations de frquence. La dernire section est consacre ltude
de lefficacit spectrale.
1
Le s de Fs fait rfrence sampling (chantillonnage).
2
BITS SIGNAL PASSE-BAS
(t )=k (t kTb ) sl (t )=Ake g(t kT )=A(t )e =i(t )+ jq(t )
jk j (t )
k k
g(t)
SYMBOLES a(t )= A e (t kT )
j k
k
k
SIGNAL PASSE-BANDE
[ ] ( ) (
s(t )=A(t )cos 2fct + (t ) =i (t )cos 2fct q(t )sin 2fct )
Figure 2 : gnration de modulations numriques - schma dtaill
Conseil de programmation : prvoir une partie initialisation de variables au dbut du code MATLAB. Lobjectif
est de programmer de faon pouvoir passer facilement aux modulations M-PSK. Par exemple, utiliser une
variable M initialise 2 au lieu dutiliser directement 2 dans le code.
NB : les symboles gnrs par la fonction pskmod() sont complexes, mme pour une BPSK. Cela permet
notamment de gnrer des symboles binaires de la forme exp(j) et exp[j(+)].
Mise en forme des symboles BPSK : gnration des signaux a(t) et sl(t)
Gnrer un peigne de Dirac a(t) de dbit D, dbit symbole, dont les amplitudes correspondent aux symboles
BPSK.
Rappel : D=Db/log2(M).
Indication : gnrer un vecteur nul de dure NT o N reprsente le nombre de symboles, avec N=Nb/log2(M), et
T dsigne la priode symbole (T=1/D), avec linstruction zeros(). Puis, placer, intervalles de temps
rguliers, les valeurs des symboles BPSK.
Remarque sur lhypothse de signaux bande troite : lhypothse de signaux bande troite (fc trs grande
devant D) permet de faire un certain nombre de simplifications dans les calculs sur les signaux numriques. Ici,
cette hypothse nest pas vrifie. En effet, pour visualiser les signaux et leurs densits spectrales de puissance
sans devoir systmatiquement faire un zoom avant et pour rduire le temps de simulation, cette hypothse a
t abandonne. Cependant, pour que les rsultats tablis en cours restent valables, cette hypothse a t
remplace par une autre hypothse, savoir que fc=kD, o k est un entier non nul.
3
Visualisation des rsultats
Reprsentations en temps
Diviser une figure en trois parties et reprsenter les signaux suivants :
Un exemple est fourni Figure 3. Laxe temporel sera exprim en multiples de la priode symbole T. Pour la
reprsentation du peigne de Dirac des symboles, seuls les chantillons non nuls seront reprsents. Expliquer
pourquoi le signal passe-bande na pas la forme exacte dun signal sinusodal.
1
Band-Pass Signal
0
-1
0 2 4 6 8 10
[t/T]
1
Re[Low Pass Signal]
0 Re[Symbols]
-1
0 2 4 6 8 10
[t/T]
1
Im[Low Pass Signal]
0 Im[Symbols]
-1
0 2 4 6 8 10
[t/T]
Figure 3 : reprsentation temporelle du signal passe-bande, des parties relle (Re[Low Pass Signal]) et imaginaire
(Im[Low Pass Signal]) du signal passe-bas, et des parties relle (Re[Symbols]) et imaginaire (Im[Symbols]) des
symboles BPSK
Consignes pour toute la suite des travaux pratiques : toujours donner une lgende explicite aux figures.
Diffrencier les courbes par un trac diffrent (en prvision dune ventuelle impression sur une imprimante
noir et blanc).
Remarque importante : les tracs MATLAB prsents dans cet nonc ont t retravaills (zoom sur les axes,
travail sur le type de trait, les markers , etc). Cela permet de prsenter au lecteur les parties intressantes
des figures. Pour le compte-rendu, lenjeu nest pas de restituer exactement les figures de cet nonc mais de
trouver une manire personnelle de prsenter les informations contenues dans une figure.
Reprsentation en frquence
Sur une mme figure, tracer les trois densits spectrales de puissance correspondant aux trois signaux
suivants :
4
2
Les densits spectrales de puissance seront traces entre Fs/2 et Fs/2 en utilisant une chelle logarithme . Un
exemple est fourni Figure 4. Expliquer la forme des densits spectrales de puissance.
2
10
0
10
-2
10
-4
10
Symbols
Low-Pass Signal
-6 Band-Pass Signal
10
-0.5 0 0.5
Normalized Frequency
Figure 4 : reprsentation en frquence des symboles (Symbols), du signal passe-bas (Low-Pass Signal) et du signal
passe-bande (Band-Pass Signal) relatifs une modulation BPSK
Indication : la densit spectrale de puissance dun signal donn est obtenue par le module au carr de sa
3
transforme de Fourier (FFT pour Fast Fourier Transform) , divis par le nombre de points de la transforme de
Fourier, not nfft.
Rappel : le calcul de la FFT est optimis pour un nombre de points de la FFT multiple dune puissance de deux. Il
est conseill de choisir, comme nombre de points de la FFT, la puissance de deux immdiatement suprieure
au nombre dchantillons du signal. Pour cela, utiliser la fonction nextpow2().
2
Le recentrage des densits spectrales de puissance se fait par la fonction fftshift().
3
Sous MATLAB, lchelle des frquences par dfaut va de 0 1 o 1 reprsente la frquence dchantillonnage.
5
Remarque : lenjeu consiste ici expliquer les diffrences qui peuvent exister entre deux courbes. Sil est
possible dexpliquer la valeur absolue de chaque reprsentation, cest un plus mais ce nest pas obligatoire.
Compte-rendu : explications de linfluence des diffrents paramtres sur les reprsentations en temps et en
4
frquence des signaux . Cette partie nest pas obligatoire mais optionnelle.
2 2 2
10 10 10
0 0 0
10 10 10
-2 -2 -2
10 10 10
-4 -4 -4
10 10 10
2 2 2
10 10 10
0 0 0
10 10 10
-2 -2 -2
10 10 10
-4 -4 -4
10 10 10
Figure 5 : densits spectrales de puissance simules (trait continu) et thoriques (ronds) des symboles BPSK (figures
de gauche), du signal passe-bas (figures centrales) et du signal passe-bande (figures de droite) pour un nombre de
statistiques Nstat gal 1 (figures du haut) et 100 (figures du bas)
4
Il est possible dajouter des figures pour illustrer une explication.
6
Gnralisation aux cas des modulations M-PSK (M-ary Phase Shift Keying)
(programme MATLAB intitul modulateur_v2.m)
Adapter le programme modulateurv0.m pour quil puisse traiter toutes les modulations de phase (cas
M>2). Une tape de codage supplmentaire sera rajoute entre la gnration de bits et la gnration de
symboles dans la mesure o les fonctions pskmod() nacceptent, en entre, que des entiers entre 0 et M-1.
Indication : initialiser la valeur de M en dbut de programme et faire ensuite un programme qui fonctionne
pour tout M.
Compte-rendu : programme MATLAB, intitul modulateur_v2.m, lensemble des figures pour le cas M=4 et
explications.
Indication : initialiser le choix de la modulation en dbut de programme et faire ensuite un programme qui
fonctionne quelle que soit la modulation.
Compte-rendu : programme MATLAB, intitul modulateur.m, et lensemble des figures pour une
modulation 16-QAM et 4-PAM. Expliquer les changements observs au niveau des reprsentations en temps et
en frquence des signaux par rapport la modulation BPSK.
Rcepteurs numriques
Lobjectif de cette partie consiste tudier la conception des rcepteurs numriques (voir Figure 6).
Larchitecture optimale sera implante et simule. Des versions sous-optimales seront galement implantes
afin dapprhender les lments caractristiques de larchitecture optimale.
7
SIGNAL PASSE-BAS
SIGNAL
PASSE-BANDE
BITS
+
BRUIT
SYMBOLES
Le caractre optimal fait rfrence ici au fait que le rcepteur propos permet datteindre les meilleures
performances en termes de BER (Bit Error Rate) dans le cas dun canal AWGN (Additive White Gaussian Noise).
Ds lors que le canal nest plus AWGN, le rcepteur propos nest plus optimal.
Cette partie commence par ltude de la modulation de phase BPSK. Puis, ltude est gnralise aux autres
modulations de phase, aux modulations damplitude et aux modulations damplitude en quadrature, et enfin
aux modulations de frquence.
Filtrage adapt
La suppression des lobes secondaires en -2fc et +2fc, et la maximisation du rapport signal bruit aux instants de
prise de dcision, se font grce au filtrage adapt. Filtrer le signal (aprs le retour en bande de base) par un
filtre adapt au filtre de mise en forme.
NB : si le filtre de mise en forme est de la forme g(t), alors la rponse impulsionnelle du filtre adapt est de la
forme g*(-t) o est un retard servant compenser les retards introduits par les diffrentes tapes de filtrage.
La permutation de la rponse impulsionnelle du filtre de mise en forme se fait par la fonction fliplr().
NB : tenir compte des retards introduits par le filtre de mise en forme et le filtre adapt (temps de propagation
de groupe) dans le sous-chantillonnage. En particulier, le premier chantillon, tant donn la forme des filtres
(porte de dure T), est prlev linstant T.
Le peigne de Dirac reprsentant les symboles mis (seules les amplitudes non-nulles seront traces),
Le signal passe-bas la sortie du filtre de mise en forme,
Le signal passe-bas la sortie du filtre adapt,
8
Le peigne de Dirac correspondant aux chantillons prlevs toutes les priodes T (seules les
amplitudes non-nulles seront traces).
Laxe temporel sera exprim en multiples de la priode symbole T. Une figure sera pour la partie relle et une
autre pour la partie imaginaire. Un exemple est fourni Figure 7. Commenter.
10
Re[Low Pass Signal]
Re[Symbols]
Re[Matched Filter Output]
0
Re[Output Samples]
-10
0 2 4 6 8 10
[t/T]
10
Im[Low Pass Signal]
Im[Symbols]
Im[Matched Filter Output]
0
Im[Output Samples]
-10
0 2 4 6 8 10
[t/T]
Figure 7 : reprsentation temporelle de la partie relle et imaginaire du signal passe-bas (Low Pass Signal), des
symboles mis (Symbols), de la sortie du filtre adapt (Matched Filter) et des chantillons prlevs la sortie du filtre
adapt (Output Samples)
Reprsentation en frquence
Superposer, sur la mme figure, les densits spectrales de puissance, entre Fs/2 et Fs/2, du signal passe-bande
mis et du signal la sortie du filtre adapt (voir Figure 8). Commenter.
2
10
0
10
-2
10
-4
10
Figure 8 : reprsentation en frquence du signal passe-bande mis (Band-Pass Signal) et du signal la sortie du filtre
adapt (Matched Filter Output)
9
Impact de diffrents paramtres sur la rception des signaux numriques
Impact de la synchronisation
Changer la frquence porteuse du rcepteur. Pour cela, rajouter 1% la valeur de fc. Sur une mme figure,
superposer deux signaux : lun correspondant la sortie du filtre adapt pour une frquence porteuse la
rception gale fc, et lautre avec la nouvelle valeur de la frquence porteuse. Commenter et conclure.
Changer la phase de la porteuse du rcepteur. Faire varier la phase de 0 /2 par pas de /8 (la valeur de la
frquence porteuse doit reprendre sa valeur initiale). Sur une mme figure, superposer les signaux la sortie
du filtre adapt, correspondants aux quatre valeurs de la phase. Commenter et conclure.
Compte-rendu : programme MATLAB demodulateur_v1.m, les figures pour le cas M=4 et explications.
10
Compte-rendu : programme MATLAB demodulateur.m et les explications relatives au facteur de
normalisation.
Lhypothse de blancheur du bruit signifie que les chantillons de bruit ne sont pas corrls les uns aux autres.
Cette hypothse est justifie par le fait que la valeur dun chantillon dune source de bruit particulire nest
pas conditionne par les valeurs des chantillons des autres sources de bruit. Par exemple, si les dfauts de
connectique de lmetteur sont modliss par une source de bruit, les valeurs alatoires de cette source de
bruit ne sont pas conditionnes par les valeurs de la source de bruit modlisant les perturbations lies au
milieu de transmission, et inversement.
Test : tracer la fonction dauto-corrlation dun bruit blanc gaussien de moyenne nulle, de variance unit, de
1000 chantillons. Un exemple est fourni dans la figure suivante. Commenter.
0.8
0.6
0.4
0.2
-0.2
-1000 -500 0 500 1000
Lhypothse relative au caractre gaussien de la source de bruit vient du thorme de la limite centrale : toute
combinaison de N variables alatoires indpendantes et identiquement distribues tend vers une loi
gaussienne lorsque N tend vers linfini. Il faut galement tenir compte du fait que toute combinaison linaire de
variables alatoires gaussiennes suit une loi gaussienne. Ainsi, lorsque toutes les sources de bruit sont
11
rfrences et classes selon leur densit de probabilit, il est possible dobtenir, pour chaque classe, une
variable alatoire globale suivant une loi gaussienne qui, une fois combine aux variables alatoires globales
associes aux autres classes, donne galement une loi gaussienne.
Test : tracer lhistogramme dune variable alatoire uniformment rpartie sur [0,1[ comprenant 1000
ralisations. Puis, tracer lhistogramme de la somme de 2, 3 et 10 variables alatoires. Commenter. Tracer
lhistogramme dun bruit AWGN. Un exemple est fourni dans la figure suivante.
Histogram of an AWGN
60
50
40
30
20
10
0
-4 -3 -2 -1 0 1 2 3 4
Figure 10 : histogramme d'un bruit AWGN de variance unit (1000 chantillons gnrs)
Compte-rendu : les courbes dauto-corrlation, lhistogramme dun bruit blanc et les commentaires.
Es =
Eg 2
2
a avec
+
Eg = g(t ) dt et
2
[ ]
a2 = E ak
2
o Eg est lnergie du filtre de mise en forme et a est la variance des symboles, condition que ces derniers
soient centrs et o les ak sont les symboles de la constellation. Lnergie moyenne Eb dun bit est de la forme
Es/log2(M). Le bruit blanc gaussien doit tre gnr avec une variance N0/2.
Reprsentation en temps
Sur une premire figure, superposer le signal reu sans bruit et le signal reu avec bruit. Sur une autre figure,
superposer le signal la sortie du filtre adapt dans le cas sans bruit et dans le cas avec bruit. Commenter.
Reprsentation en frquence
Sur une mme figure, reprsenter les densits spectrales de puissance des trois signaux suivants : le signal
mis, le signal reu et le bruit. Commenter. La visualisation dpend de la ralisation des processus observs.
Pour visualiser le comportement moyen, moyenner 100 ralisations et reprsenter les trois densits spectrales
sur la mme figure. Commenter. Un exemple est fourni dans la figure suivante.
12
1
10
0
10
-1
10
Received Signal
Band-Pass Signal
-2 AWGN
10
-0.5 0 0.5
Normalized Frequency
Figure 11 : reprsentation en frquence du signal mis (Band-Pass Signal), du signal reu (Received Signal) et du
bruit AWGN (AWGN)
Recommencer les reprsentations en temps et en frquence pour un rapport Eb/N0 de 10 dB. Commenter.
1 E
Pb = erfc b
2 N0
Commenter. Modifier le code MATLAB pour obtenir des courbes compltes. Superposer aux courbes
prcdentes les courbes dcart-type sur lestimation de la probabilit derreur par le BER, cest--dire
BER BER . Commenter.
Indication : le BER est un estimateur de la probabilit derreur Pb avec les caractristiques suivantes
mBER = E [BER] = Pb
Pb (1 Pb )
2
BER [ ]
= E BER 2 mBER
2
=
Nb
La prcision de la mesure de la probabilit derreur est quantifie par lerreur relative quadratique moyenne,
, dfinie par
2
BER
= 2 =
2 (1 Pb ) avec 2 =
1
si Pb << 1
mBER NbPb NbPb
13
Indication : la quantit NbPb reprsente le nombre derreurs commises au cours de la simulation.
Le calcul du BER consistait jusqu prsent envoyer un nombre de bits arbitrairement grand pour toutes les
valeurs du rapport Eb/N0 et observer un certain nombre derreurs. Ainsi, par rapport lexpression suivante
du BER
Ne
BER =
Nb
o Ne reprsente le nombre derreurs observes, la procdure prcdente consistait fixer Nb et compter Ne.
La procdure va maintenant tre modifie. Au lieu de fixer Nb, cest Ne qui va tre fix. Expliquer limpact de
cette procdure sur la prcision de la mesure.
Etablir une nouvelle courbe de BER avec la mthode suivante. Pour chaque valeur du rapport Eb/N0, envoyer
des blocs de Nb bits, avec Nb assez faible (exemple : Nb=100), et accumuler ces erreurs dans une variable
intermdiaire Nerror. Envoyer autant de blocs de Nb bits jusqu ce que Nerror soit suprieur un seuil Ne fix ici
100 (Ne=100). Le BER estim sera alors gale au nombre derreurs observes, divis par Nb, multipli par le
nombre de blocs envoys.
Commenter les rsultats et expliquer les avantages de cette mthode par rapport la prcdente en se
fondant sur la notion de prcision de la mesure. Tracer galement les courbes BER+/- BER o BER est
maintenant calcul en fonction de . Expliquer et commenter. Un exemple est fourni dans la figure suivante.
-1
10
-2
10
Figure 12 : courbe de BER thorique (Theoretical BER), courbe de BER simul (Simulated BER), courbes de BER
avec les carts-type de lestimateur suivant que le nombre de bits mis est fixe (Fixed Nb) ou que le nombre derreurs
observes est fixe (Fixed Ne)
Compte-rendu : le code MATLAB, nomm perf_BPSK.m, les courbes de BER, les commentaires et
explications.
14
Scatter plot
15
10
Quadrature 0
-5
-10
-15
-15 -10 -5 0 5 10 15
In-Phase
Figure 13 : dispersion des chantillons d'une modulation QPSK avant la prise de dcision
15
en termes de BER et en prsence dISI seront tablies. Puis, la chane de base sera adapte ce nouveau
scnario. Des filtres en racine de cosinus surlevs seront utiliss. Les performances seront de nouveau tablies
avec ces nouveaux filtres.
Diagramme de lil
Tracer le diagramme de lil la sortie du filtre adapt pour une modulation QPSK grce la fonction
eyediagram(). La visualisation sera faite sans bruit. Un exemple est donn dans la figure suivante.
Commenter.
10
Amplitude
-10
-20
-20 -10 0 10 20
Time
Eye Diagram for Quadrature Signal
20
10
Amplitude
-10
-20
-20 -10 0 10 20
Time
Figure 14 : diagramme de l'il d'une modulation QPSK la sortie du filtre adapt pour un canal AWGN
NB : faire attention modifier correctement la variance du bruit AWGN de faon ce que le rapport signal
bruit soit identique au cas prcdent.
16
Eye Diagram for In-Phase Signal
20
10
Amplitude
0
-10
-20
-20 -10 0 10 20
Time
Eye Diagram for Quadrature Signal
20
10
Amplitude
-10
-20
-20 -10 0 10 20
Time
Figure 15: diagramme de l'il d'une modulation QPSK la sortie du filtre adapt pour un canal bande limite
Compte-rendu : le code MATLAB intitul chaine_w_ISI.m ainsi que les courbes de BER et les digrammes de
lil et les commentaires.
Compte-rendu : le code MATLAB intitul chaine_wo_ISI.m ainsi que les courbes de BER et les digrammes
de lil et les commentaires.
17