UNIVERSITE
ABDELMALEK ESSAADI
ECOLE NATIONALE
DES SCIENCES APPLIQUEES
TANGER
TAVEAUX PRATIQUES
OFDM ET CDMA
ANNEE
2008/2009
Prof : mohamed moussaoui
ENSA Tanger, Route Ziaten, BP 1818
Tanger principale
[email protected] Prof : mohamed moussaoui
ENSA Tanger, Route Ziaten, BP 1818
Tanger principale
[email protected] ECOLE NATIONALE DES SCIENCES APPLIQUEES DE TANGER
PARTIE A :
Modulations Multi-porteuses OFDM
Lcart frquentiel entre deux porteuses orthogonales
II
Transmission dun signal OFDM
III Etude du cas de la Transmission OFDM dans le Wifi IEEE 802.11a
IV Etude de la transmission et la rception dun signal OFDM
Prof: Mohamed Moussaoui
TRAVAUX PRATIQUES OFDM ET CDMA
ECOLE NATIONALE DES SCIENCES APPLIQUEES DE TANGER
I. Lcart frquentiel entre deux porteuses orthogonales
Notre objectif est de calculer lcart minimal entre deux sinusodes orthogonales de
frquences f 1 et f 2 , de dure T .
On
considre que la diffrence de phase entre les deux sinusodes,
avec 0, 2 .
Pour vrifi lorthogonalit :
cos(2f 1 t + ) cos(2f 2 t )dt = 0
On peut simplifie lquation:
sin(2 ( f 1 + f 2 )T ) sin(2 ( f 1 f 2 )T )
+
cos( )
2 ( f 1 f 2 )
2 ( f 1 + f 2 )
cos(2 ( f 1 + f 2 )T ) 1 cos(2 ( f 1 f 2 )T ) 1
+ sin( )
+
=0
2 ( f 1 + f 2 )
2 ( f 1 f 2 )
(1)
On note:
sin(n ) = 0 et cos(2n ) = 1 , avec n est un entier.
On assume que ( f 1 + f 2 )T est un entier
Donc
sin(2 ( f 1 + f 2 )T ) = 0 et cos(2 ( f 1 + f 2 )T ) = 1 .
Finalement on peut crire lquation (1) sous la forme:
cos( )
sin(2 ( f 1 f 2 )T )
cos(2 ( f 1 f 2 )T ) 1
+ sin( )
= 0.
2 ( f 1 f 2 )
2 ( f 1 f 2 )
(2)
Questions:
1. Donner la condition de lorthogonalit dans le cas ou 0, 2
2. Donner la condition de lorthogonalit dans le ou = 0
3. crire un programme matlab qui explique les deux cas
On note
La dure de deux sinusodes T = 1.
La frquence dchantillonnage fs = 100
La fonction matlab rand donne une distribution uniforme
Prof: Mohamed Moussaoui
TRAVAUX PRATIQUES OFDM ET CDMA
ECOLE NATIONALE DES SCIENCES APPLIQUEES DE TANGER
II. Transmission dun signal OFDM
Lcart frquentiel entre deux sinusodes de phases arbitraires est
1
. Avec T est la
T
priode du symbole
En OFDM, on utilise des sinusodes dcart frquentiel
1
.
T
On dfini les sinusodes dun signal OFDM sous la forme :
g k (t ) =
2kt
T
exp j
rect T t
2
T
T
Avec
T
rect T t : un signal rectangulaire entre [0, T )
2
k = 0, 1, ..., K 1 correspond la frquence de la sinusode
Chaque sinusoide g k est multiple par sa propre information a k , et la somme des
sinusodes modules forme un signal OFDM s(t).
Questions
1. Donner la forme gnrale dun signal OFDM s(t).
2. Trouver la relation entre un signal OFDM et la transformation de Fourier discrte
a0
a1
g 0 (t )
a K 1
g K 1 (t )
g1 (t )
s (t )
Figure 1. Ralisation possible dun modulateur OFDM
Prof: Mohamed Moussaoui
TRAVAUX PRATIQUES OFDM ET CDMA
ECOLE NATIONALE DES SCIENCES APPLIQUEES DE TANGER
III Etude du cas de la Transmission OFDM dans le Wifi IEEE 802.11a
Les spcifications techniques de IEEE 802.11
Paramtre
valeur
FFT size. nFFT
Number of used subcarriers. nDSC
FFT Sampling frequency
Subcarrier spacing
Used Subcarrier index
Cyclic prefix duration, Tcp
Data symbol duration, Td
Total Symbol duration, Ts
64
52
20MHz
312.5kHz
-26 to -1, +1 to 26
0.8us
3.2us
4us
La dure symbole T = 3.2 s , signifi quon peut utiliser les sous-porteuses :
k
= 312.5kHz , 625kHz .......
T
La bande passante de 20Mhz est divise entre les 64 sous-porteuses dindices :
k = [ 32,.....31] .
Le nombre des sous-porteuses utilises est 52 dindices k = 26,.. 1, +1,..., +26 .
Elles sont utilises pour transmettre les squences dinformation a1 a 52 .
Chaque bit du symbole est assign son entre IFFT.
Lopration IFFT nous permet de prsenter le signal dans le domaine temporel
Intervalle de garde (cyclic prefix) : Il est possible dajouter un espace entre les
symboles OFDM, dune dure suprieure ltalement des retards du canal. Ainsi les
derniers chos du symbole OFDM auront lieu durant cet intervalle dit "de garde", et
le symbole OFDM suivant ne sera plus perturb par le prcdent.
En pratique on recopie le segment de la fin vers le dbut du symbole.
la Simulation
On note :
la taille de la IFFT est 64
52 sous porteuses utilises dindices k = 26,.. 1, +1,..., +26
Le nombre de bits par symbole OFDM est 52
Le nombre de bits totale 2500
1. gnrer une squence binaire alatoire de longueur nBit
2. convertir les bits en symboles BPSK ( 0 1, 1 +1 )
3. prsenter le signal BPSK sous forme de matrice :
Nombre de lignes= nombre de symboles OFDM
Nombre de colonnes= nombre de bit par symbole OFDM)
Pour ce but utiliser la fonction matlab reshape (voir help)
4. Pour chaque symbole OFDM, On assigne les bits a1 a 52 aux sous frquences
dindices k = 26,.. 1, +1,..., +26
5. gnrer le signal OFDM s (t ) avec lintervalle de garde
6. Tracer la densit spectrale de puissance DSP du signal OFDM
Prof: Mohamed Moussaoui
TRAVAUX PRATIQUES OFDM ET CDMA
ECOLE NATIONALE DES SCIENCES APPLIQUEES DE TANGER
IV Etude de la transmission et la rception dun signal OFDM
Dans cette section, on sintresse
tudier un simple metteur rcepteur OFDM
trouver la relation entre E b / N 0 (Bit to Noise ration) et E s / N 0 (Signal to Noise
ratio)
Calculer la probabilit derreurs dans le cas de la modulation BPSK
Cyclic prefix
Dans le cas du canal AWGN, on sait que lintervalle de garde ne porte pas dune
information extra.
Lnergie du signal est tale sur Td + Tcp , alors que lnergie du bit est tale sur Td .
c--d :
E s Td + Tcp = E bTd E s =
Td
Eb .
Td + Tcp
Ltalement frquentiel (Frequency spread)
En transmission OFDM, les sous-porteuses disponibles ne sont pas toutes utilises pour
transmettre les donnes.
Par exemple le IEEE802.11a :
la bande passante est 10 MHz , 10 MHz , mais on utilise seulement la bande de
8.1250MHz (-26/64*20MHz) +8.1250MHz (+26/64*20MHz).
Lnergie du signal est tale sur une bande passante de 16.250MHz, alors que le
bruit est tal sur la bande passante de 20MHz (-10MHz to +10MHz).
c-d :
20 Mhz . E s = 16.25. E b .
Pour simplifier
Es =
nDSC
Eb .
nFFT
La relation entre E b / N 0 et E s / N 0
En combinant les deux aspects, la relation entre lnergie du symbole et lnergie du
bit est la suivante:
Es
E
= b
N0 N0
nDSC Td
nFFT Td + Tcp
Td
Es
E
nDSC
= b + 10 log 10
+ 10 log 10
N 0 dB N 0 db
nFFT
Td + Tcp
Prof: Mohamed Moussaoui
TRAVAUX PRATIQUES OFDM ET CDMA
ECOLE NATIONALE DES SCIENCES APPLIQUEES DE TANGER
La simulation
On note:
nFFT = 64;
nDSC = 52
nBitPerSym = 52
nSym = 10000
EbN0dB = [0:10]
Pour chaque EbN0dB
1. gnrer une squence binaire alatoire longueur nBit
2. gnrer le signal BPSK
3. Pour chaque symbole OFDM, assigner les bits a1 a 52 aux sous frquences
4.
5.
6.
7.
8.
dindices k = 26,.. 1, +1,..., +26 , et gnrer le signal OFDM transmis avec
lintervalle de garde
ajouter le bruit AWGN
diviser le vecteur reu en multiples symboles, et garder les sous porteurs utiles
dmoduler et convertir aux bits
calculer le nombre de bits errons
tracer la porbabilit drreur en fonction EbN0dB
On donne la probabilit derreur thorique de la modulation BPSK :
Pb , BPSK =
Eb
1
erfc
N
2
0
Prof: Mohamed Moussaoui
TRAVAUX PRATIQUES OFDM ET CDMA
ECOLE NATIONALE DES SCIENCES APPLIQUEES DE TANGER
PARTIE B :
LAccs Multiple de Division de Code
(CDMA)
I. Le modle de CDMA asynchrone
II. La dtection conventionnelle
III. Simulation de CDMA dans un canal AWGN
IV.Simulation de CDMA asynchrone dans un scnario
Multi-Utilisateurs
Prof: Mohamed Moussaoui
TRAVAUX PRATIQUES OFDM ET CDMA
ECOLE NATIONALE DES SCIENCES APPLIQUEES DE TANGER
I. Le modle de CDMA asynchrone
Dans cette section nous prsentons un modle mathmatique de systme DS-CDMA
asynchrone pour la liaison montante (Uplink) avec une modulation BPSK.
Un systme gnral de CDMA multi-utilisateurs est illustr dans la figure 2.
b1 ( t )
P1 c1 ( t )
b2 ( t )
r (t )
P2 c 2 ( t )
n(t )
bK (t )
PK c K (t )
Fig. 2 - Modle de lmetteur pour le systme CDMA asynchrone
Dans ce modle, les K utilisateurs partagent les mmes mdias de communication
et les signaux transmis par les utilisateurs passent dans des canaux spars et
indpendants.
Les sorties des canaux sont ajoutes un bruit commun aditif blanc gaussien.
Chaque utilisateur transmet un flux de donnes binaires modul comme suit
bk ( t ) =
N b 1
b
n=0
(t nT ) ,
k ,n
avec T la priode de symbole,
bk ,i le bit i de lutilisateur k, bk , i { 1,1}et (t ) la
forme d'onde de symbole.
Le flux des donnes de chaque utilisateur est tal par son propre code de taux
lev, sa forme donde est exprime comme :
ck (t ) =
N c 1
m =0
( t mTc )
k ,m
Avec c k ,m est le chip m de l'utilisateur k, Tc est la priode de chip, N c =
T
est le
Tc
gain de traitement (ou facteur dtalement : Spreading Factor), et (t ) la forme
d'onde de chip. On suppose que c k (t ) est relle (l'analyse peut tre simplement
gnralise dans le cas complexe) et normalise, donc,
Prof: Mohamed Moussaoui
TRAVAUX PRATIQUES OFDM ET CDMA
ECOLE NATIONALE DES SCIENCES APPLIQUEES DE TANGER
c k ( t ) = 0 si t [0, T ) et
( t ) dt = 1 , les chips sont binaires, c--d c k ,m { 1,1} . La
forme donde de la signature est priodique de priode T.
Le signal mis par chaque utilisateur aprs ltalement est
s k (t ) =
Pk bk ( t )c k ( t ) ,
avec Pk la puissance en bande de base de l'utilisateur k.
Il est suppos aussi que le systme CDMA tudi ici est asynchrone et que les
retards sont uniformment distribus dans lintervalle k [0, T ) k .
Le systme CDMA est dit synchrone si les retards sont gaux (et, ainsi, normalisables
zro), c.--d., 1 = 2 = ... = K = 0 et
quasi-synchrone
si les retards sont petits
durant l'intervalle de symbole.
Le signal reu est :
K
r ( t ) = s k ( t k ) + n( t ) =
k =1
N b 1 K
n = 0 k =1
Pk bk , n c k ( t nT k ) + n( t ) ,
Le bruit n(t ) est de moyenne nulle ergodique et stationnaire AWGN complexe avec
un module de spectre de puissance N 0 .
II. La dtection conventionnelle
Le dtecteur conventionnel est une batterie de K corrlateurs, comme montr dans
la figure 3.
r (t )
1
T
1
T
1
T
t = iT
Dcision
b1
t = iT
Dcision
b2
c1 ( t )
T
c2 (t )
T
t = iT
Dcision
b K
c K (t )
FIG. 3 - Le dtecteur conventionnel de DS-CDMA : une batterie de corrlateurs
Prof: Mohamed Moussaoui
TRAVAUX PRATIQUES OFDM ET CDMA
10
ECOLE NATIONALE DES SCIENCES APPLIQUEES DE TANGER
Ici, chaque code est rgnr et corrl avec le signal reu dans une branche
spare de dtecteur. Les sorties des corrlateurs sont chantillonnes au temps de
bit, pour donner des estimations "souples" (soft) ( Z m ) des donnes transmises.
En final, des dcisions "dures" (hard) ( bm ) des donnes sont donnes selon le signe
des estimations souples.
La sortie du corrlateur de l'utilisateur
Z m ,i =
1
T
( i +1)
iT
avec Am = 1
{r ( t )c1 (t )}dt = Am +
( i +1)
k =1
km
k ,m
Pm bm ( t )c m ( t )c m ( t ) dt =
T iT
1444442444443
m pendant lintervalle de bit i est
+ ,
Pm bm,i
L' utilisateu r _ Dsir
Im =
k =1
km
k ,m
1 ( i + 1 )T K
=
Pk bk ( t k )c k ( t k ) c m ( t )dt
iT
T
kk =1m
443
4444424444444
14444
MAI
( i + 1 )T
1
{n( t )c m ( t )}dt .
T iT
14442444
3
Bruit
Avec Am , i k , m et les contributions de l'utilisateur dsir, linterfrence de l'utilisateur
k, et le bruit AWGN, respectivement
Dans lhypothse SGA (approximation gaussienne standard), la probabilit d'erreur
moyenne approximative de symbole dans un canal asynchrone est donne sous la
forme :
E (Z m , i )
PSGA = Q
var (Z )
m ,i
Prof: Mohamed Moussaoui
2
Pk
k =1
1
1 k m
= Q
+
.
Pm
2 E b1 3 N
N
0
TRAVAUX PRATIQUES OFDM ET CDMA
11
ECOLE NATIONALE DES SCIENCES APPLIQUEES DE TANGER
III. Simulation de CDMA dans un canal AWGN
On note:
le nombre de bits simuler : x_nm = 10000;
La longueur du code : N = 64;
La dure de chip : Tc=1 ;
La densit spectrale du bruit est fixe 2
Lintervalle du rapport signal/bruit : EbNodB = -20:2:10;
Pour chaque EbN0dB
1. Convertir le rapport signal/bruit en linaire
2. Calculer la puissance dmission
3. Gnrer les symboles BPSK
4. Gnrer le code Hadamard
5. Gnrer le signal tal plus le bruit AWGN
6. Ds-taler le signal et rcuprer les symboles reus
7. Tracer la probailit derreur en fonction du rapport signal/bruit
Prof: Mohamed Moussaoui
TRAVAUX PRATIQUES OFDM ET CDMA
12
ECOLE NATIONALE DES SCIENCES APPLIQUEES DE TANGER
IV.Simulation de CDMA asynchrone dans un scnario
Multi-Utilisateurs
On note:
le nombre de bits simuler : x_nm = 300000;
le nombre dutilisateurs: K = 4
La longueur du code : N = 32
La dure de chip : Tc=1
La densit spectrale du bruit est fixe 2
Lintervalle du rapport signal/bruit : EbNodB = -20:3:50,
1- gnrer les retards en chips entre les utilisateurs
2- pour chaque EbNodB
Convertir le rapport signal/bruit en linaire
Calculer la puissance dmission
Gnrer les symboles BPSK
Gnrer le code Hadamard
prendre le code 1 de la matrice comme le code de lutilisateur 1 et
gnrer son signal tal
ajouter les signaux des autres utilisateurs au signal de lutilisateur 1
ajouter le bruit AWGN
Ds-taler le signal de lutilisateur 1 et rcuprer les symboles reus
2-a) Tracer la probabilit derreur en fonction du rapport signal/bruit
2-b) tracer la probabilit derreur en fonction du nombre dutilisateurs (La capacit)
Prof: Mohamed Moussaoui
TRAVAUX PRATIQUES OFDM ET CDMA
13
