Amine GUISOUS Mémoire Fin D'étude

REPUBLIQUE ALGERIENNE DEMOCRATIQUE ET POPULAIRE
MINISTERE DE LENSEIGNEMENT SUPERIEUR ET DE LA RECHERCHE

SCIENTIFIQUE

Universit de Mohamed El-Bachir El-Ibrahimi - Bordj Bou Arrridj

Institut des Sciences et de la technologie
Dpartement dElectronique

Mmoire
Prsent pour obtenir

LE DIPLOME DE MASTER

FILIERE : ELECTRONIQUE

Spcialit : Rseaux et technologie de tlcommunication

Par

Mr. Amine GUISSOUS & Melle. Sarah BOUTEHRA

Intitul

Soutenu le : 26/06/2013

Anne Universitaire 2012/2013

Etude comparative des diffrents codes dEtalement.
Applications ltude de linter-corrlation et du
temps dacquisition dans les Systmes GPS/Galileo.
Encadrs par :
Dr. K. ROUABAH
Mr. S. ATIA
Mr. M. FLISSI
Devant les jurys :
- Dr. S. AIDEL Prsident
- Pr. D.CHIKOUCHE Membre
-

REPUBLIQUE ALGERIENNE DEMOCRATIQUE ET POPULAIRE
MINISTERE DE LENSEIGNEMENT SUPERIEUR ET DE LA RECHERCHE
SCIENTIFIQUE

Universit de Mohamed El-Bachir El-Ibrahimi - Bordj Bou Arrridj

Institut des Sciences et de la technologie
Dpartement dElectronique

Mmoire
Prsent pour obtenir

LE DIPLOME DE MASTER

FILIERE : ELECTRONIQUE

Spcialit : Rseaux et technologie de tlcommunication

Par

Mr. Amine GUISSOUS & Melle. Sarah BOUTEHRA

Intitul

Soutenu le : 25/06/2013

Anne Universitaire 2012/2013
Applications ltude de linter-corrlation et du
temps dacquisition dans les Systmes GPS/Galileo.
Encadrs par :
Dr. K. ROUABAH
Mr. S. ATIA
Mr. M. FLISSI
Devant les jurys :
- Dr. S. AIDEL Prsident
- Pr. D.CHIKOUCHE Membre
-

Chapitre 0

Applications { ltude de linter-corrlation et du temps
dacquisition dans les Systmes GPS/Galileo.
Chapitre 0
Introduction gnrale

es codes sont les lments fondamentaux dans n'importe quel systme d'accs
multiple par rpartition des codes tel que les systmes de navigation par
satellites GPS et Galileo. Ces codes reprsentent l'outil qui permet un rcepteur
GNSS de distinguer et didentifier un satellite parmi d'autres dans la constellation.
Malgr leur grande importance, aucune grande innovation n'a t faite dans le monde
de la navigation par satellites concernant la gnration de ce type de codes. Le
systme GPS ancienne gnration emploi des codes de type Gold depuis plus de 30
ans.
1- Problmatique
L'obtention de meilleures performances en termes de rsistance vis--vis du
bruit et des Multitrajets ncessite l'utilisation des codes trs longs. Cela implique que
l'espace de recherche sera beaucoup plus grand en comparaison avec les structures
classiques. Par consquent, l'utilisation des processus d'acquisition et de poursuite
traditionnels pour l'accrochage de ce type de code peut tre trs pnalisante en
termes dau moins un facteur qui est la vitesse d'acquisition. Les caractristiques de
ce type de code obligent les concepteurs de rcepteurs de dvelopper de nouvelles
techniques d'acquisition rapide ou d'optimiser ce type de code. Le dveloppement de
nouvelles techniques d'acquisition et de poursuite ncessite un changement matriel
ce qui peut compliquer l'architecture du rcepteur. Par consquent l'optimisation des
codes longs, qui seront utiliss dans les applications Galileo et GPS nouvelle
gnration, s'avre une solution meilleure. C'est justement dans ce cadre que s'inscrit
ce travail de recherche. En effet, depuis quelques annes, un travail important a t
ralis afin de dvelopper de nouveaux algorithmes permettant la gnration de
codes optimiss. De ce fait, le systme Galileo, annoncera pour la premire fois cette
opration dans un espace fortement multidimensionnel, pour amliorer la prcision
de localisation. De ce fait, Le systme Galileo emploie des codes alatoires et le GPS
nouvelle gnration emploie son tour de nouvelles structures de codes bass sur
des squences de type Legender.
Le but de ce travail, serait donc l'tude et la comparaison des diffrents codes
d'talement pour les applications GPS/Galileo. Le compromis entre les inter-
corrlations des diffrents codes utiliss dans la mme application et la dure
d'intgration optimale sera considr pour effectuer cette tude.
L

Notre mmoire est donc organis de la manire suivante :
Dans le premire chapitre, nous allons dcrire le principe du systme de navigation
GNSS ainsi que les diffrentes techniques d'acquisition des codes. Lobjectif du 2me
chapitre de ce mmoire est de prsenter les diffrentes techniques d'accs multiple.
Nous allons baser ainsi notre tude sur la technique d'accs multiple par rpartition
de code. Dans les chapitres 3 et 5, nous allons tudier les diffrents types de codes
pseudo-alatoires utiliss dans le CDMA d'une faon gnrale et dans les applications
GNSS comme cas exceptionnel. Dans le dernier chapitre, des tudes comparatives
peuvent tre envisages pour en arriver finalement choisir le code le plus
appropri. Nous allons clturer ce mmoire par des rsultats de simulation suivies
dune conclusion dordre gnrale et des perspectives.

I. Chapitre I
Chapitre I
Les systmes de navigation
GNSS

Plan du chapitre
1. Introduction
2. Prsentation du systme GNSS
3. Principe de la localisation
4. Actualits
5. Le rcepteur GNSS
5.1. Acquisition des signaux GNSS
5.1.1. La recherche srie
5.1.2. La recherche parallle
5.1.3. La recherche hybride
5.2. Les modes de dtection
6. Les sources derreurs de la localisation
7. Conclusion

Rsum

Ce chapitre est une prsentation gnrale des systmes de navigation par
satellite GNSS pour "Global Navigation System Satellite". Il fait l'objet aussi
d'une description dtaille du principe de localisation par satellites. Un aperu
sur lactualit de ces systmes de navigation ainsi que les nouvelles
gnrations seront galement abords. Dans une partie finale, les lments
principaux des rcepteurs savoir les diffrents types des structures
dacquisition et de recherche tel que, la recherche srie, parallle et hybride
seront prsents avec une tude de leurs sensibilits vis--vis les sources
d'erreurs.

CHAPITRE I
Les systmes de navigation GNSS
1
I.1 Introduction
es systmes de navigation par satellite GNSS (Global Navigation Satellite System)
[1], tel que GPS (Global Positioning System), et Galileo(Le nouveau systme de
navigation par satellites) prsentent de nombreuses applications dans les deux
domaine civil et militaire[2]. Les signaux mis par les satellites de la constellation
permettent { tout utilisateur (mobile ou immobile) quip dun rcepteur adquat
capable de recevoir quatre signaux avec une puissance suffisante pour se positionner,
pour connaitre son temps et mme sa vitesse { nimporte quel endroit du globe (sur
la terre, la mer, dans laire et mme dans lespace proche de la terre) par rapport {
une rfrence de temps et un repre de coordonns dans lespace [3]. Les avances en
recherches scientifique favorisent lutilisation des GNSS dans des environnements de
plus en plus contraints, comme le centre-ville [4]. Toutefois, elles ne permettent pas
dassurer la prcision de la localisation quel que soit lenvironnement.
Les systmes de navigation permettent un positionnement autonome du porteur, trs
prcis court terme ; mais cette prcision va se dgrader de plus en plus au cours du
temps : la vitesse de cette drive dpend de plusieurs paramtres (environnement,
type de rcepteur et la qualit des services fournis). On peut citer deux exemples
dapplications de ces services : les applications militaires avec une erreur de lordre
de quelques millimtre, par contre les applications civiles permettent une localisation
avec une erreur pouvant atteindre plus de 10 m.
Plusieurs solutions, particulirement dans le domaine de traitement des signaux
GNSS permettent damliorer tout dabord la qualit de lmission et la rception par
lintgration des nouveaux concepts pour la gnration, la modulation et lacquisition
des signaux, en suite la conception des rcepteurs optimiss.
Vu du cot utilisateur, le systme GNSS prend la forme dun rcepteur, autonome ou
intgr, fournissant une information de positionnement et de datation. Derrire

1 (Diggelen, 2009)
2 (Bradford, 1996)
3 (BONIN.G)
4 (Rouabah, September 2009)
L

CHAPITRE I
2
chacun de ces rcepteurs se cache une architecture et un principe de positionnement.
C'est ce que nous allons voir dans ce chapitre.
I.2 PRSENTATION DU SYSTME GNSS
Les systmes Galileo et GPS reposent sur le mme principe de fonctionnement :
Une constellation de satellites en orbite autour de la Terre :
Chaque satellite de la constellation diffuse en permanence un signal vers
lensemble des zones visibles de la Terre.
Chaque satellite inclut dans son signal les informations prcisant sa position dans
lespace.
Un nombre illimit de rcepteurs utilisateurs, avec certaines proprits :
rception des signaux, mesure de distances et calcule de la position.

Figure I-1 - La constellation des satellites autour de la terre.
En plus de lutilisation de la technique dtalement de spectre, Galileo et GPS utilisent
de nouvelles structures de signaux et de modulations pour chaque bande de
frquence afin dassurer une grande couverture ainsi qu'une grande prcision de
localisation [5]. La Figure I-2 prsente le mcanisme de la rception des signaux
GNSS.

5 (Braasch M.S, Jan.1990)

CHAPITRE I
3

Figure I-2 - un mcanisme simple de la rception des signaux GNSS.
Comme l'illustre cette figure, on trouve une antenne spcifique permettant de capter
les signaux en provenance des satellites. On trouve aussi le bloque de dmodulation
pour liminer la porteuse et transformer le signal en bande de base. Le rcepteur
GNSS effectue une opration de corrlation pour faire la synchronisation
(Acquisition+ poursuite) dans le temps du code reu depuis les satellites visibles et
celui gnr localement [6].
I.3 Principe de la localisation
En gnral, tous les systmes GNSS reposent sur le mme principe de
fonctionnement .Trois segments (spatial, contrle et utilisateurs) sont ncessaires
pour dterminer la position et mme la vitesse dun objet quelconque (quil soit
mobile ou immobile), en se basant sur les phnomnes de la propagation des ondes
EM et les lois de la gomtrie.
Ces systmes sont bass sur des mesures de distances entre le rcepteur et un certain
nombre de satellites dont les positions sont connues tout instant. Ces distances sont
dduites du temps de propagation (Time-delay) du signal radiolectrique sur un
trajet simple, entre le point de rfrence du satellite metteur et le rcepteur. Afin de
calculer ce temps de propagation, le rcepteur utilise les mesures de pseudo-distance
de code et ventuellement aussi de phase de la porteuse.
La Figure I-3 illustre le principe de mesure de la phase entre le code reu et le code
gnrer par le rcepteur.

6 (Jean-Marc, 2006)

CHAPITRE I
4

Figure I-3 - Le principe pour dterminer la distance.
Les satellites envoient des signaux qui se propagent jusquaux rcepteurs { la vitesse
de la lumire. Lors de la rception, le dcalage en temps du code reu, permet de
dterminer la distance entre le satellite et le rcepteur.
Pour dterminer les coordonnes exactes du porteur (longitude latitude et altitude)
On rsout un system d'quations trois inconnus.
{
( )
( )
( )
( )
( )
( )
( )
( )
( )

( )
Pour se faire, on utilise 3 satellites pour dterminer les distances D1, D2 et D3 plus un
4
me
satellite pour synchroniser les 3 satellites.

Figure I-4 - Le principe de la triangulation.
La rsolution de ce systme dquations est effectue { chaque instant par le
rcepteur. Cest ce quon appelle le calcul des coordonnes (x, y, z), qui donne la
position, la vitesse et mme la direction du porteur lorsque il sagit dun mobile. Le

CHAPITRE I
5
calcul donne galement accs au paramtre , qui reprsente lcart du temps du
rcepteur par rapport au temps de rfrence; pour les services autoriss un
paramtre qui permet de corriger les erreurs du positionnement pour avoir une
meilleure prcision { laide du coefficient .
I.4 Actualits
Les gnrations modernes du GPS et Galileo qui sont devenues plus performantes,
intgrent des applications avances et plus sres par rapport aux systmes de
navigation par satellites anciens. Pour pouvoir satisfaire les besoins futurs et
grandissant en matire de navigation, un systme combinant les diffrents systmes
de navigations devrait introduit. La combinaison de ces systmes de navigation va
permettre une grande couverture et une grande prcision de localisation. Cest une
technologie de pointe qui rsulte de l'mission, { partir dune constellation de plus de
60 satellites. Ceci permet chaque personne ayant un rcepteur adquat de recevoir
les signaux de localisation des diffrents systmes de navigations, par exemple les
travaux [7](le nombre des satellites = 31 GPS+ 4 SBAS + 3 QZSS + 22 GLONASS= 60
satellite).
I.5 Le rcepteur GNSS
Ltude du rcepteur nest pas le but principal de notre travail mais la comprhension
de son fonctionnement est une tape ncessaire et importante. En effet, le rcepteur
est compos dune chane de traitement des signaux satellitaires reus. Son rle est de
calculer la position de lantenne de rception grce aux mesures du temps de
propagation des signaux satellitaires [8]. Pour cela, un rcepteur spcifique
dtermine les distances entre son antenne et chacun des satellites reus. Ces
distances sont par la suite appeles pseudo-distances pour signifier le caractre
erron de ces mesures. Le rcepteur a pour principales fonctions de recevoir les
signaux radios mis par les satellites et de diffrencier les signaux issus des diffrents
satellites en visibilit directe. Il doit galement dcoder le message de navigation

7 (FRANK VAN DIGGELEN, 2011)
8 (El-Rabbany, 2002)

CHAPITRE I
6
pour dterminer la position du satellite, sa vitesse et les paramtres de lhorloge.
Durant ce processus, le rcepteur effectue les oprations suivantes :
1) Une recherche et une identification des satellites visibles.
2) Acquisition du signal.
3) Dmodulation.
4) Poursuite du code.
5) Poursuite de phase de la porteuse et extraction des donnes.
6) Lecture du message de navigation.
Lors de la poursuite du code par le rcepteur, celui-ci effectue la corrlation
[9] du code reu du satellite avec sa rplique qui est gnr localement au niveau du
rcepteur. Lorsque les deux codes sont en phase (align parfaitement),
lautocorrlation est maximale. Le dcalage temporel appliqu { la rplique
reprsente la mesure du temps de propagation de la porteuse. En plus de la poursuite
de code, la plupart des rcepteurs peuvent galement poursuivre la phase de la
porteuse. Lobservation de la phase, permet dobtenir une information de distance
plus prcise. La boucle verrouillage de phase, mesure la diffrence de phase entre le
signal reu du satellite et un signal gnr par le rcepteur la mme frquence. La
valeur recueillie est en fait une mesure de phase instantane qui prend en compte le
nombre entier de cycles de la phase de battement perus depuis le dbut des
mesures. Cependant, la valeur obtenue est ambigu puisque elle est tronque dun
nombre entier de cycles, correspondant au nombre entier (inconnu) de longueurs
donde sparant le satellite du rcepteur { linitialisation de ce dernier. Ce nombre
entier est appel ambigit de la phase. Si laccrochage sur le signal est maintenu,
lambigit est constante au cours du temps et peut tre estime grce des
mthodes diffrentielles.
I.5.1 Acquisition des signaux GNSS
Ltape de lacquisition est une tape primordiale pour la navigation, elle
permet de chercher et de dtecter les satellites visibles, de synchroniser le rcepteur
avec les satellites (synchronisation grossire) et de dterminer quelques paramtres

9 (D. V. Sarwate, 1980)

CHAPITRE I
7
essentiels de la navigation, tel que le dcalage en temps et en phase du code reu
ainsi que la variation de la frquence centrale qui est cause par leffet de Doppler.
Plusieurs approches sont possibles, les rcepteurs peuvent raliser une acquisition
dun seul satellite ou bien de plusieurs satellites en parallle. Lune des
caractristiques principales de lacquisition est la rapidit de la recherche.
Plusieurs solution sont proposes pour diminuer le temps de la recherche, la plus
importante cest la recherche intelligente qui est base sur des processus trs
complexes pour identifier [10] des satellites qui ont une position connue (une
constellation contrle). C'est--dire le rcepteur peut diminuer le nombre des
satellites cherchs et aussi minimiser lnergie par les deux modes de fonctionnement
aveugle et almanach. En gnral, lacquisition dpend de la qualit du rcepteur. Afin
que ce rcepteur dtecte les satellites, il fait une recherche en deux axes, le premier
cest laxe des frquences et lautre c'est laxe de temps. La figure suivante illustre le
principe de la recherche bidimensionnelle.

Figure I-5 - Le principe de la recherche bidimensionnelle.
En gnral, il existe trois schmas pour la recherche :

10 (Barkat, 2005)

CHAPITRE I
8
I.5.1.1 La recherche srie
Le premier schma est la recherche srie (Serial search). Dans cette mthode
le circuit dacquisition essaie de changer progressivement la phase de la squence du
code local et de tester toutes les phases possibles une par une jusqu{ la dtection
d'un alignement de phase. Le dcalage a un pas fixe de longueur , Polydoros et
Weber [11] ont tudi les performances de tel systmes pour un canal AWGN et
Chang et Lee lont tudi pour un canal { trajet multiples [12].
La figure suivante reprsente le mcanisme de la recherche srie pour lacquisition.

Figure I-6 - le mcanisme de la recherche srie.
Lavantage de cette mthode est la simplicit. Par contre le temps dacquisition est
trs lent.
I.5.1.2 La recherche parallle
La mthode de la recherche parallle permet de tester simultanment toutes les
phases possibles. la Figure I-7 reprsente cette mthode.
Ce schma dacquisition reste une mthode destimation qui consiste choisir la plus
grande valeur de la variable estimer.
La recherche parallle sappelle aussi lacquisition rapide base sur le calcul de la FFT.
Elle garantit un gain trs important du temps de la recherche. En contrepartie, le
nombre de ressources matriels devient trs excessif pour des codes longs.

11 (Weber, May 1984.)
12 (Moon, pp. 926-928)

CHAPITRE I
9

Figure I-7 - le mcanisme de la recherche parallle (FFT).
I.5.1.3 La recherche hybride
La recherche hybride est une combinaison de la recherche srie et de la recherche
parallle, elle permet de faire un compromis entre la vitesse dacquisition et la
simplicit du system.
I.6 Les modes de la dtection
Le dtecteur joue un rle fondamental dans le processus de lacquisition. Il permet de
faire la dtection, avec un degr de prcision lev. Trois approches de base sont
possibles : la dtection cohrente, la dtection non cohrente et la dtection par la loi
carre. Ces types de dtecteurs sont reprsents sur la Figure I-8 :
Le premier mode (a) consiste multiplier le signal reu avec une porteuse
() (
). Lintgrateur permet de corrler le code dmodul ( ) avec

un code local ( ) dont la valeur de dcalage est contrle. La sortie est attaque
par un dtecteur(en gnral un comparateur). Avec le mme principe, les deux modes
(b) et (c) fonctionnent. La seule dfrence est le type de lestimateur de la corrlation.
Pratiquement la sortie des corrlateurs est une variable alatoire avec une
distribution prsente dans les figures Figure I-9 et Figure I-10.

CHAPITRE I
10

Figure I-8 - Les diffrents modes de dtection.
Dans les trois cas, la variable de dcision y est compare un seuil T [13].Si y dpasse
la valeur du seuil T, le dtecteur dclare que les codes sont ventuellement en phase
bien sr avec certaine probabilit derreur. Les deux fausses alarmes [12] (erreurs)
possibles sont :
Dtecter un satellite qui nest pas visible P (1|0) ou ne pas dtecter un satellite visible
P(0|1). En gnral, une fausse alarme gnre une augmentation dans le temps
dacquisition. En effet, lopration de la poursuite [14] du code sera alors active mais
le systme va vite sapercevoir quil sagit dune fausse acquisition. Dans ce cas, il rend
la main au processus de lacquisition pour reprendre la recherche aprs un temps de
pnalit.

13 (Landry, 1998)
14 (Mattos, April 2006)

CHAPITRE I
11

Figure I-9 - Schma synoptique de la dtection des signaux binaires.
Ce modle permet destimer les symboles
partir des chantillons () chaque

instant.
Alors les chantillons (), sont une variable alatoire avec une densit de
probabilit :
() (
(|
) ( )

Figure I-10 - la distribution des symboles et la prise de dcision pour le cas binaire.
La densit de probabilit conditionnelle de () scrit comme :

(|
) ) ( ()) ( ()) ( )
- La probabilit derreur par symbole est dfinie par :
) .

CHAPITRE I
12
- La probabilit derreur par bit est dfinie par :
)
Il existe une relation entre
et
tell que :

( )
La valeur minimale de
est atteinte quand un bit erron conduit au symbole

erron.
I.7 Les sources derreurs de la localisation
Des modles mathmatiques des erreurs de la localisation (ionosphre, troposphre)
qui sont implants dans les rcepteurs, les informations transmises dans les
phmrides du message de navigation (La dgradation de lhorloge du satellite)
ainsi que les mthodes de traitement du signal permettent tous de diminuer les
sources derreurs. La prcision des mesures de pseudo-distance de code aprs ces
corrections est caractrise par lUERE (User Equivalent Range Error) : cette
grandeur est estime en calculant la somme quadratique des composantes derreur
rsiduelle affectant la mesure. Le Tableau I-1 illustre le bilan derreur des mesures de
pseudo-distance avec le code C/A.
Segment Source de lerreur Erreur (m) Erreur en %
Espace Stabilit de lhorloge de satellite 3 14,56311
Perturbation du satellite 1 4,854369
Autre (radiation thermique, etc.) 0,5 2,427184
Control Erreur de prdiction des phmrides 4,2 20,38835
Autre (performance du propulseur, etc.) 0,9 4,368932
U
t
i
l
i
s
a
t
e
u
r

Retard ionosphrique 5 24,27184
Retard troposphrique 1,5 7,281553
Bruit de rcepteur et de rsolution 1,5 7,281553
Multitrajet 2,5 12,13592
Autre (biais inter-canal, etc.) 0,5 2,427184
UERE Total 20,6 100
Tableau I-1- Budget derreur affectant les mesures de pseudo-distance de code pour
le code C/A.

CHAPITRE I
13
Conclusion
Dans ce chapitre nous avons donn une brve description du principe de la
navigation par satellites. En effet, nous avons introduit les notions de base de la
navigation et de la localisation par rcepteurs GNSS. Les fonctions principales de la
rception, de la recherche et de la dtection des satellites, par le concept de
lacquisition des signaux capts, ont t aussi prsentes.
Comme il a t dmontr, les diffrent types dacquisition, srie, parallle est
hybride sont caractriss par un temps de recherche et une architecture soit simple
pour les modle de recherche en srie ou complexe pour dautres. Toutes ces
structures sont bases sur la corrlation entre le signal reu et un code gnr
localement. Nous avons prsents dans une dernire phase plusieurs approches pour
corrler les signaux savoir, la corrlation cohrente qui ne prend pas en compte le
dcalage en phase, et la corrlation non cohrente.

CHAPITRE I
14
I.8 Bibliographie
[1] Diggelen, F. v. (2009). A-GPS: Assisted GPS, GNSS, and SBAS. BOSTON LONDON:
ARTECH HOUSE.
[2] Bradford, W. P. (1996). Global Positioning System: Theory end applications.
Spikler. American Institute of Aeronautics and Astronautics .
[3] BONIN.G. (2010). Systme GPS de positionnement par satellites. Techniques de
lingnieur.
[4] Rouabah, K. (September 2009). GPS/Galileo Multipath Detection and Mitigation
Using Closed-Form Solutions. Mathematical Problems in Engineering Vol. 2009.
[5] Braasch M.S. (Jan.1990). GPS receiver architectures and measurements A.J.
Proceedings of the IEEE Volume 87, Issue 1, 48 64.
[6] Jean-Marc, P. (2006). GPS et Galileo Systmes de navigation par satellites.
EYROLLES. Paris .
[7] FRANK VAN DIGGELEN, C. A. (2011, Mars/Avril). GPS,GLONASS,QZSS and SBAS in
a Singlechip. GNSS Inside Mobile Phones.
[8] El-Rabbany, A. (2002). Introduction to GPS The Global Positioning System.
Boston .London: Artech House Mobile Communications Series.
[9] D. V. Sarwate, M. B. (1980). Crosscorrelation Properties of Pseudorandom and
Related Sequences. IEEE Trans. Commun. Technol., vol. COM-68, pp. 594-596.
[10] Barkat, M. (2005). Signal Detection and Estimation. Second Edition. Boston
london: ARTECH HOUSE.
[11] Weber, A. P. (May 1984.). A unified approach to serial search spread spectrum
code acquisitionPart II: matchedfilter receiver. EE Trans. on Communications vol.
32, pp. 550-560.
[12] Moon, B. L. (pp. 926-928). Detection and false alarm probability of PN code
acquisition in DS-CDMA system. Electronics Letters, vol. 33, November 1997.
[13] Landry, R. (1998). Techniques dabaissement des seuils dacquisition et de
poursuite pour les rcepteurs GPS. Paris, France: Editions HERMS.
[14] Mattos, P. G. (April 2006). Galileo L1C Acquisition complexity: cross correlation
benefits, sensitivity discussions on the choice of pure pilot, secondary code, or
something different. IEEE/ION Position Location, and Navigation Symposium, pp.
845852.

II. Chapitre II
Chapitre II
Ltalement de spectre
par squence directe

Plan du chapitre
1. Introduction
2. Mthodes daccs multiples
2.1. Accs multiples par la rpartition de frquences (FDMA)
2.2. Accs multiples par la rpartition dans le temps (TDMA)
2.3. Accs multiple par la rpartition dans les codes (CDMA)
3. Ltalement de spectre
3.1. Introduction
3.2. Historique
3.3. Le secret de ltalement de spectre
3.4. Les principes de ltalement de spectre
4. Modlisation du signal au niveau de lmetteur et de rcepteur
5. Conclusion

Rsum

Les systmes de navigation GNSS reposent sur la CDMA (Code Division
Multiple Access). Ce chapitre dveloppe une description gnrale des
techniques daccs multiples utilises dans les transmissions numriques et
spcialement la technique CDMA base sur ltalement de spectre. Dans ce
chapitre nous dcrivons lhistorique de cette technique ainsi que sa dfinition
et ces caractristiques. Ensuite, nous donnons la modlisation du signal au
niveau de lmetteur dans un systme CDMA.

CHAPITRE II
Ltalement de spectre par squence directe
16
II.1 Introduction (un calcul pour une introduction)
Supposons que le signal transmettre soit un signal analogique pralable quantifi en
amplitude sur q niveaux reprsents par des nombres binaires de n bits, avec
naturellement
. Il conviendra donc de transmettre n bits pour chantillon (si la

valence m=2), soit
impulsion par seconde, si
est la frquence
dchantillonnage. Le dbit de moments est donc :

() ( )
La bande passante ncessaire { la transmission PCM dun signal analogique dont le
spectre est born par valeur suprieure
est donc, daprs la relation de

Nyquist [15]:

() ( )
Avec
frquence maximum prsente dans le spectre du signal analogique (au sens

de Shannon), correspondant usuellement la frquence de coupure du filtre anti
repliement du convertisseur analogique numrique.
On constate ainsi que la bande passante ncessaire la transmission du signal
analogique numris est n fois plus importante que celle ncessaire sa transmission
analogique en bande de base. Pour contourner cette pnalit intrinsque, deux voies
sont traditionnellement exploites :
La rduction du nombre de bits transmettre, avec les techniques de compression
de donnes.
Laugmentation de la rapidit pour transmettre plus de deux bits pendant la dure
dun symbole. Dans la suite de ce chapitre, nous allons en particulier dvelopper
et analyser ce principe en commenant par les mthodes daccs.
II.2 METHODES DACCES MULTIPLES
Les communications dans les systmes radio-mobiles utilisent une bande de
frquence alloue au systme, par des organismes de rgulation, dont la largeur est
limite. Cette bande de frquence doit tre utilise d'une faon plus judicieuse afin
dcouler le maximum de communications. Elle est partage en canaux qui sont

15 (Marvin K. Simon, 2002)

CHAPITRE II
17
allous la demande aux mobiles pour permettre lchange dinformations dun
terminal mobile avec le rseau ou avec dautres mobiles. La dfinition des canaux de
communication dpend de la mthode daccs multiples retenue.
Les trois principales techniques daccs multiple sont les suivantes :
Accs multiples par une rpartition de frquences (Frequency Division
Multiple Access, FDMA).
Accs multiples par une rpartition dans le temps (Time Division Multiple
Access, TDMA).
Accs multiples par une rpartition des codes (Code Division Multiple Access,
CDMA).
II.2.1 Accs multiples par la rpartition de frquences (FDMA)
Cest la mthode daccs multiples la plus ancienne. Elle est utilise principalement
dans les systmes analogiques et elle est combine la mthode TDMA dans la
majorit des systmes numriques. Elle consiste subdiviser la bande alloue (canal)
en petites bandes continues de frquences (porteuses). Chaque porteuse est utilise
pour vhiculer un appel unique et dans un seul sens la fois (sens montant ou sens
descendant, Uplink ou Downlink). En fonction de la capacit du systme et ses
besoins en signalisation, un ou plusieurs canaux de contrle sont utiliss comme le
reprsente la Figure II-1.

Figure II-1 : Partage des canaux en FDMA

CHAPITRE II
18
Les principales caractristiques de la mthode daccs multiples FDMA sont les
suivantes :
Un circuit unique par porteuse : chaque canal FDMA est dfini pour vhiculer une
seule communication.
Transmission continue : quand les canaux de communication dans les deux sens
ont t attribus, les deux extrmits mettent en continu et de faon simultane.
Faible largeur de bande : les canaux FDMA sont relativement troits, en gnral
30 kHz ou moins, puisquils ne vhiculent quun circuit par porteuse.
Faible complexit du terminal mobile : en effet, la transmission en mode FDMA ne
ncessite pas dgalisation ou de tramage complexe et de synchronisation associe
la transmission de bursts (rafales ou salves), comme dans le cas des systmes TDMA,
puisque les informations sont mises et reues sans interruption de faon synchrone.
Ncessit dutiliser un duplexeur : comme lmetteur et le rcepteur doivent
fonctionner de manire simultane, le mobile doit utiliser un duplexeur permettant
dviter les brouillages entre lmetteur et le rcepteur du mobile. Cet quipement
entrane galement des cots supplmentaires.
II.2.2 Accs multiples par la rpartition dans le temps (TDMA)

La technique TDMA est la premire alternative la technique FDMA. Elle est mise en
uvre dans les systmes numriques comme alternance principale { la technique
CDMA. Elle permet de transmettre des dbits dinformations plus importants quun
systme FDMA. La porteuse (frquence radio) est partage en N intervalles de temps
(slots) et peut tre utilise par N terminaux, chacun utilisant un slot particulier. Le
nombre de slots par canal est choisi en fonction de plusieurs facteurs tels que la
technique de modulation, la bande de frquence disponible, le dbit recherch, etc.
La Figure II-2 reprsente la rpartition de la plage de temps :

CHAPITRE II
19

Figure II-2 : Partage des canaux en TDMA

Les diffrentes caractristiques dun systme TDMA sont :

Plusieurs circuits par porteuse : tous les systmes TDMA multiplexent au moins
deux circuits par porteuse. Le systme GSM par exemple multiplexe 8 circuits (ou 8
canaux) par porteuse.
Transmission par bursts : la transmission venant dun mobile nest pas continue.
A chaque instant, seule une fraction des mobiles connects au systme et en mission,
do limpact sur le niveau dinterfrences sur le canal qui vari de faon importante
dun slot { lautre.
Bande large ou troite : les largeurs de bandes ncessaires pour les systmes
TDMA sont de lordre de quelques dizaines { quelques centaines de . La largeur
de bande est dtermine par la technique de modulation choisie avec des dbits de
quelques .
Dbits binaires et dures symboles : les dbits du canal les plus levs crent des
interfrences inter-symboles bien plus importantes que dans un systme FDMA. Des
retards observs en milieu urbain. Linterfrence inter-symbole a donc un impact
important dans un systme TDMA.
Non ncessit dun duplexeur : comme les missions et les rceptions ont lieu sur
des slots diffrents, le recours { un duplexeur nest pas ncessaire.

CHAPITRE II
20
II.2.3 Accs multiple par la rpartition des codes (CDMA)
Larchitecture CDMA repose sur la technique de modulation talement de spectre
(spread spectrum, SS). Ltalement de spectre est un moyen de transmission selon
lequel les donnes occupent une largeur de bande plus large que celle ncessaire au
transfert des donnes dune communication. La technique dtalement de spectre
comme son nom lindique exige que le signal { transmettre occupe une largeur de
bande beaucoup plus importante que le cas des systmes FDMA et TDMA, comme le
montre la Figure II-3.

Figure II-3 Partage des canaux en CDMA.
Donc, la technique CDMA est une mthode daccs o chaque usager est autoris
dutiliser toute la bande (le cas du TDMA) durant toute la dure dappel (le cas du
FDMA).

Les caractristiques principales de ce type de systme sont :
Nombre de circuit par porteuse trs lev : les systmes CDMA utilisent un canal
unique ou trs peu de canaux frquentiels. Ils peuvent thoriquement transporter des
dizaines dappels sur chaque porteuse.
Largeur de bande : les systmes CDMA ncessitent des largeurs de bande trs
importantes.

CHAPITRE II
21
Dbits binaires et dures symboles : cause des dbits binaires levs, la dure
de symbole est trs courte. Avec un dbit de , chaque symbole a une dure
environ (dans le cas dune modulation BPSK par exemple). Cette proprit
permet damliorer la rsolution frquentielle qui est alors proportionnelle 1/NT au
lieu de 1/T. Ceci est intressant pour la mesure de la distance entre le mobile et la
station de base, par exemple pour la rcupration des signaux issus des trajets
multiples.
Complexit au niveau du mobile : le traitement des informations reues et mises
est beaucoup plus important que dans les autres types de systmes puisquil faut
rajouter un niveau de codage supplmentaire.
Ncessit du contrle de puissance : lun des principaux inconvnients de la
mthode daccs CDMA est la ncessit de disposer un mcanisme de contrle de
puissance trs performant. Le contrle de puissance rapide est le seul moyen
permettant de maximiser le nombre dutilisateurs communicants en mme temps
dans le rseau.

CHAPITRE II
22
II.3 LETALEMENT de spectre
II.3.1 Introduction
Toutes les mthodes de modulation se rfrent au mme principe : une certaine
bande passante (ou canal) la plus troite possible est affecte une voie de
communication. Les interfrences entre canaux sont empches en vitant les
superpositions de bandes passantes. Une transmission par talement de spectre
abandonne ce schma : Les bandes passantes se superposent. En dautres termes,
cest par dizaines ou centaines que des voies de communication partagent les mmes
gammes de frquence.
II.3.2 Historique
Le principe dtalement de spectre est connu depuis la fin des annes quarante, mais
le secret a longtemps t gard et son usage rserv { des fins militaires jusquau
dbut des annes quatre-vingts. Dans ce contexte, les buts recherchs ont t la lutte
contre le brouillage intentionnel de canaux de transmission dune part, la
confidentialit des communications dautre part. Les deux techniques corres-
pondantes mises en uvre sont ltalement de spectre par vasion de frquence et
ltalement de spectre par squence directe [16].
Les applications civiles, la radio navigation par exemple, tirent profit aujourdhui de
ces mthodes. Le vocabulaire a bien sur chang. Lutilisation de la ressource spectrale
sen trouve par ailleurs optimise : il devient inutile de rserver des intervalles vides
entre deux canaux. Cest dailleurs un avantage, puisque lefficacit spectrale se
trouve rduite par lopration dtalement de spectre.
En particulier, Ltalement de spectre est une technique dveloppe historiquement
pour le cryptage : Un signal tal par une technique adapte ne peut tre repr par
balayage des frquences, il ne peut tre brouill par mission dun message qui
interfrerait, et de plus se confond avec le bruit naturel dune transmission (si lon
utilise une technique dtalement fonde sur des squences pseudo alatoires).
On pourra citer le premier brevet dpos en 1941 sous le nom de Secret

16 (Scholtz)

CHAPITRE II
23
Communication System sur la technique dtalement par lactrice Hedy Lamarr avec
laide du compositeur Georges Antheil qui a dcrit en 1941 une liaison hertzienne
bloque pour commander des torpdos (brevet 2.292.387 dpos le 10 juin 1941 en
tats-Unis)Figure II-4. Bien que l'ide n't pas prise au srieux au dbut, elle est
mme oublie. Hedy et Georges ont tait trs intresss la ralisation d'un systme
permettant la transmission de donnes de manire sr. Celui-ci a t utilis par les
militaires ds 1960.

Figure II-4 - Le schma propos sous le nom Secret Communications System.
La communaut scientifique la redcouvert en 1957 Sylvania Electronic Systems
Division. Quelques annes plus tard larme la utilis sur le terrain { cause de ses
capacits de scurisation des communications et sa rsistance aux brouilleurs (anti-
jam).
En 1948 Shannon et Hartley publient une quation tablissant que la capacit
maximale dun canal de communication (c'est--dire ; le dbit maximal dun canal)
pour transmettre une information sans erreur est proportionnelle la bande
passante de ce canal et au logarithme du rapport signal bruit exprim en terme de
puissance, sous rserve dun procd de codage adapt (voir quation II. 3).
Shannon publia en 1950 les principes de ltalement de spectre qui seront mis en
application, { titre dexemple, dans les annes 70 par larme amricaine dans la
technique GPS (Global Positioning System) [17].
Il fallait attendre les annes 90, avec lavnement des tlcommunications mobiles,

17 (Sass)

CHAPITRE II
24
pour voir lutilisation de ltalement de spectre pour le partage de ressources (laccs
multiples) : limplmentation du CDMA pour les communications cellulaires est
tudie, puis standardise sous la norme IS-95 en 1993. Ce protocole, associ la
tlphonie mobile de seconde gnration. De nombreux systmes utilisent les
proprits du CDMA pour lapplication aux communications cellulaires, en particulier
le W-CDMA (pour Wideband Code Division Multiple Access) qui sera reconnu par
lITU (International Telecommunication Union) comme standard de la tlphonie de
troisime gnration de mme quun standard concurrent le CDMA2000. Le premier
rseau commercial utilisant le W-CDMA est mis en place en 2001 au Japon par la
socit NTT DoCoMo. Paralllement, lEurope dveloppe son projet de tlphonie de
troisime gnration sous lappellation UMTS (Universal Mobile Telecommunications
System), qui utilise la technique daccs W-CDMA, tout en tant { lorigine
incompatible avec le systme japonais [18]
,
[19].
II.3.3 Le secret de ltalement de spectre
Considrons le thorme de Shannon et Hartley concernant la capacit dun canal de
communication :

) ( )
Dans cette quation, C reprsente la capacit maximale dun canal en bits par seconde
(bit/s ou bps), cest le taux de transfert maximum pour un taux derreur binaire (Bit
Error Rate, BER) nul, { condition quun procd de codage adquat puisse tre
trouv. B tant la bande passante du canal en Hertz et
, le rapport de puissance
signal/bruit [20].
On peut donc augmenter la capacit maximale en agissant sur la largeur de bande de
faon linaire ou en agissant sur le rapport signal bruit de faon logarithmique.
Les erreurs peuvent tre soit tout simplement ignores soit corriges par lutilisation
de protocoles de transmission de niveau suprieur. Au niveau de la formule, en

18 (McGillem)
19 (Jabbari)
20 (Proakis, 1999)

CHAPITRE II
25
fonction du type de bruit et du procd de codage/dcodage, on peut intgrer le BER
sous la forme de laddition dune constante au rapport signal sur bruit.
Dans le cas du CDMA, le bruit est constitu principalement par les autres utilisateurs
dont on cherchera { augmenter le nombre. Il en rsulte quen rgle gnrale un
systme CDMA opre sur des rapports signal bruit faibles, voire trs faibles. Par
changement de base des logarithmes (base 2 vers base e), lquation (II. 3) devient :
()
(
) (
) ( )

Si la puissance du signal est infrieure la puissance du bruit, on peut simplifier et
linariser lexpression (II.3), en appliquant le dveloppement en srie de Mac Laurin
de ( ) :
] ( )
Puisque ltalement du spectre permet un rapport
trs faible et que la puissance du

signal utile pouvant tre infrieure au niveau du bruit. Pour un rapport
,
lquation (II.5) devient alors :
) ( )
Et par approximation on obtient :
( )
La dpendance capacit/rapport signal bruit est approximativement linaire. La
bande tale permet donc la transmission de signaux perturbs par dautres signaux
considrs alors comme du bruit, cest { dire la transmission de signaux sur le mme
support. Le nombre de canaux utiliss un instant donn pourra varier de faon
souple puisque laugmentation du nombre dutilisateurs se traduira simplement par
une augmentation, pour tous, du taux derreur. Ceci permet en tlphonie de
maintenir une qualit de service sensiblement gale pour tous, (plutt quune
dprciation totale pour un utilisateur) ajustable et relativement facile.

CHAPITRE II
26
II.3.4 Les principes de ltalement de spectre
Le principe de ltalement du spectre, quelle que soit la mthode utilise, repose sur
le codage de linformation { transmettre avec une squence pseudo alatoire
(Pseudo-Noise, PN), connue uniquement par lmetteur et le rcepteur.
La Figure II-5 et Figure II-6 rsume le principe de base qui consiste multiplier
chaque symbole informatif (data) de priode symbole
par une squence pseudo-

alatoire, dont les lments ou chips ont une priode

(priode chip) beaucoup
plus faible que
. La longueur
de la squence PN et la priode de chip

sont
deux caractristiques importantes des transmissions DS-SS.
Dans la Figure II-5, la dure de la squence est gale la priode symbole. Cependant,
il existe des systmes de transmission o la dure totale de la squence PN est plus
grande que

[21]. Toutefois, nous nous limitons aux cas o la dure de la squence
PN est gale la priode symbole.

Figure II-5 - Le principe de la CDMA.

21 (R. L. Pickoltz, May 1982)

CHAPITRE II
27

Figure II-6 - la chaine de communication par codage CDMA.
Grce aux proprits des squences pseudo alatoires utilises, l'nergie utilise
transmettant le signal de l'information de bande de base tant tal sur une bande
plus large [22] comme illustr sur les figures (a, b et c) de Figure II-7.

Figure II-7 - L'opration de l'talement et le ds-talement de spectre, le signal tal
est devenue sous le niveau du bruit avec une bande de frquence plus large.
Le signal transmis se comporte alors comme un bruit vis--vis des autres utilisateurs
qui travaillent en bande troite ou de ceux qui ne possdent pas le mme code.

Figure II-8 - Le codage et le dcodage des signaux informatif.
La Figure II-8, illustre lopration du codage et du dcodage de linformation
transmet, linformation rcupr ce nest quune estimation pour linformation
origine.

22 (Ziemer, 1985)

CHAPITRE II
28
II.4 Modlisation du signal au niveau de Lmetteur et de
rcepteur
Le principe de ltalement de spectre semble simple et vident, mais son excution est
plus complexe [23]. Afin d'accomplir ces objectifs, diffrents systmes de ltalement
de spectre sont disponibles [24], mais ils ont tous le mme principe : ils
perfectionnent lopration de ltalement et de destalement au moyen d'un code
pseudo alatoire [25] fix un utilisateur de la voie de transmission. La faon
d'insrer ce code dans la chane de transmission dfinit les techniques particulires
de ltalement de spectre, la Figure II-9 montre en dtail un exemple de ces
techniques.
Selon la faon d'insrer le code dans la chane, on distingue deux mthodes :
Etalement de spectre par squence direct - Direct Sequence Spread Spectrum
(DS-SS)- quand le code est insr au niveau de linformation. En pratique, la squence
pseudo-alatoire est mlange ou multiplie avec le signal de l'information. C'est le
cas pour le system GPS et Galile.
Etalement de spectre par saut de frquence -Frequency Hopping Spread
Spectrum (FH-SS)- quand le code agit au niveau de la porteuse. Les codes de FHSS
force la porteuse changer ou sauter selon cette squence pseudo-alatoire. C'est le
cas pour la technologie Bluetooth.
La Figure II-9 illustre la structure de base de lmetteur, du canal et du rcepteur pour
une transmission DS-SS :

Figure II-9 - Structure gnrale dune chaine de transmission DS-SS.

23 (Garg, 2007)
24 (Turin, March 1980)
25 (Popovic, June 1999)

CHAPITRE II
29
Pour une meilleure comprhension, on suppose quun seul utilisateur, dans notre
systme DS-CDMA, nous utilisons la modulation de phase (BPSK) avec des squences
de codes PN bipolaires prenant des valeurs de lensemble { } pour ltalement
spectral. Le signal transmis est alors donn par [26]:
()
()() (
) ( )
O

est la puissance moyenne, () le signal de donnes, ()le signal de la squence
dtalement PRN,
la porteuse et la phase initiale du signal. La relation entre la

dure du bit de la donne et la dure dun chip est :
nous pouvons alors

exprimer ()et ()de la manire suivante:
()
()
) ( )

()
()
) ( )

O
()

()
prennent des valeurs de lensemble { 1, +1}
() et
()sont des
impulsions rectangulaires damplitude 1 et de dure

et
respectivement.
Le rcepteur est modlis par lexpression ( ).
Le signal { lentre du rcepteur est donn par la relation
()
( )( ) (
) () ( )
O () est un bruit blanc gaussien(AWGN) de moyenne nulle et une densit
spectrale de puissance
, cest le retard relatif entre lmetteur et le rcepteur

due la transmission, est le dphasage qui est modlis par une variable alatoire
indpendante uniformment distribues sur , -
Le rcepteur utilise un filtre adapt, alors la sortie du rcepteur est donne par :
() ()
() (
) ( )

26 (Viterbi, 1979.)

CHAPITRE II
30
II.5 Conclusion
Dans ce chapitre on a essay de prsenter, le plus explicitement possible, une
technique puissante dans les transmissions numriques et spcialement dans le
domaine de la navigation par satellites. Il sagit en fait de la technique CDMA qui est
base sur ltalement de spectre et qui possde des proprits trs importantes pour
lutter contre le bruit et les interfrences grce aux caractristiques des codes allous
cette technique.
Dans les chapitres suivants, on va tudier les codes afin de dterminer les
critres permettant de choisir les codes les plus performants pour cette technologie
et en particulier pour les applications de la navigation GNSS.

CHAPITRE II
31
II.6 Bibliographie
[15] Marvin K. Simon, J. K. (2002). SPREAD SPECTRUM COMMUNICATIONS
HANDBOOK. New York : McGraw-Hill, Inc.Electronic Edition.
[16] Scholtz, R. A. (May 1982). The origins of spread spectrum communications. IEEE
Trans., vol. 30, pp. 90-92.
[17] Sass, P. F. (July 1983). Why is the army interested in spread spectrum? IEEE
Communications Magazine, vol. 21, pp. 23-25.
[18] McGillem, G. R. (1986). Modern communications and spread spectrum. Inc., USA:
McGrawHill.
[19] Jabbari, E. H. (September 1998). Spreading codes for direct sequence CDMA and
wideband CDMA cellular networks. IEEE Communications Magazine.
[20] Proakis, J. G. (1999). Digital communications, Third edition. New York, USA: Mc
Graw-Hill.
[21] R. L. Pickoltz, D. L. (May 1982). Theory of spread-spectrum. IEEE Trans. on
Communications.
[22] Ziemer, R. E. (1985). Digital Communications and Spread Spectum Systems. New
York, NY.: MacMillan.
[23] Garg, V. K. (2007). WIRELESS COMMUNICATIONS AND NETWORKING.
Amsterdam Boston Heidelberg: Morgan Kaufmann ELSEVIER.
[24] Turin, G. L. (March 1980). Introduction to spread-spectrum antimultipath
techniques and their application to urban digital radio. Proceeding of the IEEE, vol.
68, pp. 328-354.
[25] Popovic, B. M. (June 1999). Spreading sequences for multicarrier cdma systems.
IEEE Trans. on Communications, vol. 47 pp. 918-926.
[26] Viterbi, A. (1979.). Principles of Digital Communication and Coding. McGrawHill

III. Chapitre III
Applications ltude de linter-corrlation et du temps
Chapitre III
Les codes de ltalement
et leurs performances

Plan du chapitre

1 Introduction
2. Les squences binaires pseudo-alatoires
2.1. Gnration dune squence pseudo-alatoire
2.2. Proprits statistiques
3. Les Squences Pseudo Alatoires PN (Pseudo Noise)
4. Les codes de Gold
5. Les codes de Kasami
6. Les codes de Walsh-Hadamard
7. Les codes OVSF
8. Les codes de Golay
9. Les codes de Barker
10. Codes de Weil et les Squences de Legendre
11. Les codes mmoires
12 Conclusion

Rsum

Les codes sont les acteurs principaux dans les systmes CDMA pour
ltalement de spectre. Cependant, les codes dtalement dans les systmes de
navigation ont un autre objectif essentiel qui est la synchronisation entre les
diffrents satellites et les rcepteurs afin doffrir une meilleure prcision dans
la localisation. Cette prcision dpend directement du type de code utilis et
de ces performances en termes de corrlation. Dans ce chapitre on va tudier
les diffrents types de codes dtalement en se basant sur les mthodes
arithmtiques de gnration que leurs niveaux de corrlation.

CHAPITRE III
Les codes de ltalement et leurs performances
33
III.1 Introduction
Dans les applications des tlcommunications modernes et spcialement
pour les applications de la navigation GNSS, le nombre infini dapplications exige une
recherche de codes plus longs et plus efficaces. Actuellement, les codes jouent un rle
exceptionnel dans les systmes de navigation GNSS comme GPS et Galileo.
Ces codes sont des squences binaires qui doivent tre aussi longs et aussi alatoires
que possibles. C'est dire ils doivent apparatre comme un bruit pour certains
rcepteurs et dterministes pour dautre. Cette caractristique est la plus importante
pour ltalement et le ds-talement qui dpend directement de la qualit du code.
Pratiquement, les codes alatoires sont les meilleurs codes en termes de
performances et doptimisation. Mais ce type de codes a un temps dexcution
relativement long. Le plus important cest que les codes doivent tres
reproductibles ; autrement dit, le rcepteur devrait tre capable de rgnrer les
codes dsirs pour dtecter les satellites et synchroniser ces codes avec les codes
reus mais en temps rel. C'est la raison pour laquelle les squences pseudo-alatoire
sont le choix le plus favorable.
Alors comment peut-on choisir ces codes et comment juger si un code est meilleur ou
non ? Pour rpondre { cette question, on va dabord tudier les diffrents types de
codes et voir leurs caractristiques.
Comme nous avons dj prcis, pour la plupart des codes, ce sont uniquement les
codes alatoires qui ncessitent une mmoire pour le stockage de toutes les
squences. Nous avons besoin d'un gnrateur de squences pseudo-alatoires et
dun algorithme qui produit des squences qui ne sont pas vraiment alatoires, mais
priodiques avec des priodes trs longues. Sur la notion pseudo alatoire des codes,
John Von Neumann (1951) a mmorablement nonc, Anyone who considers
arithmetical methods of producing random digits is, of course, in a state of sin.
Donc, nimporte quelle squence gnre par une arithmtique nest pas une
squence alatoire.
Dans la section suivante on va citer les arithmtiques de gnration de quelques
types de codes dtalement ainsi que leurs caractristiques en termes de corrlation.
Les squences binaires pseudo-alatoires

CHAPITRE III
34
Une squence pseudo-alatoire est une suite dvnements cre de faon
dterministe mais qui apparait alatoire tout observateur non averti (deux
vnements successifs ne sont pratiquement pas corrls). Les squences pseudo
alatoires prsentent lavantage dtre centres [27] et de possder un spectre riche.
Les plus utilises sont les squences binaires pseudo alatoires (SBPA) appeles aussi
PN (Pseudo Noise) ou m-squences pour (maximum length sequences MLS-), mais
on peut concevoir dautres types comme les squences ternaires pseudo-alatoires et
les squences pseudo alatoire base de registre dcalage non linaire.
III.1.1 Gnration dune squence pseudo-alatoire
Pour gnrer une SBPA, on utilise les registres dcalage linaire (Figure - III-1)
boucls par un OU-exclusif ().
Pour un ensemble fix de valeurs des paramtres *
+,
*+ la squence
construite est priodique de longueur maximale
car

implique une sortie constamment binaire. Pour tudier les proprits de ces
squences, on dfinit leur polynme caractristique. Daprs le fonctionnement des
registres dcalage, on a :
*+ ( )
O R reprsente loprateur de dcalage. Ce qui peut scrire, en tenant compte de :

( )
En dfinissant le polynme caractristique par :
()
( )

27 (Willett, April 1976)

CHAPITRE III
35

Figure - III-1 : Gnration d'une SBPA { base dun registre { dcalage linaire LFSR.
La relation ()
. La langueur de la squence est fixe par le choix des
, et
elle est maximale si le polynme () dordre N est irrductible [28], [29]. Des
exemples de tels polynmes pour N donn sont indiqus dans le Tableau III-1:
N P(R)
3 1+R+R
3
4 1+R+R
4
5 1+R
2
+R
5
6 1+R+R
6
7 1+R
3
+R
7
8 1+R
2
+R3+R
4
+R
8
9 1+R
4
+R
9
10 1+R
3
+R
10
12 1+R+R
4
+R
6
+R
12

14 1+R+R
6
+R
10
+R
14

16 1+R+R
3
+R
12
+R
16

18 1+R
2
+R
18

20 1+R
3
+R
20

22 1+R+R
22

24 1+R+R
2
+R
7
+R
24

26 1+R+R
2
+R
6
+R
26

28 1+R
3
+R
28

30 1+R+R
2
+R
23
+R
30

32 1+R+R
2
+R
22
+R
32

Tableau III-1 : Les polynmes irrductibles correspondant aux diffrentes tailles du
registre gnrateur.

28 (McEliece, 1972)
29 (N. HADJ-SAID, 2009)

CHAPITRE III
36
La SBPA est obtenue partir des valeurs binaires successives du bit
, sur lesquelles
on applique une modulation en bande de base NRZ par la transformation :
{

( )
Sur une priode . On peut ainsi obtenir une SBPA continue (analogique) en la
considrant constante sur , ( )- (bloqueur) ou une SBPA discrte
(chantillonne) en ne considrant que les instant (Figure - III-2).

Figure - III-2 : Signal binaire pseudo-alatoire. (a) reprsentation continue, (b)
reprsentation discrte.
III.1.2 Proprits statistiques
Une squence de priode contient
et
. A partir de cette
remarque, on peut calculer, sur une priode, des estimations [30] des fonctions de
corrlations de ces signaux :
a. Cas discret :
- La moyenne :
*+
()
( )
- La fonction dautocorrlation :
()
() ( )
( )

30 (N. E. Bekir, 1978)

CHAPITRE III
37

()
( )
de on a
couples (() ( )) tel que les couples

soient de mme signe, on obtient donc :
()
( )
b. Cas continu
- La moyenne :
*+
()
( )
- La fonction dautocorrlation :
()
() ( ) ( )
() ( ) ( )

Si || :
Supposons, dans un premier temps . Soient deux intervalles de temps
conscutifs
et
. Sil ny a pas de changement de signe de entre
et
alors :
() ( )
( )

Et sil y a un changement de signe, alors :
() ( )
) ( )

Pour une squence de longueur , avec un changement supplmentaire de signe, on
obtient :
()
( )

) ( )
On applique ce raisonnement pour , on obtient finalement :
()

||
) ( )
Si || :

CHAPITRE III
38
De on a
intervalles , ( )-, il y a un nombre pair de

changement de signe entre () ( ) ( , ( )- , on a alors :
() ( )
( )
()

On obtient ainsi pour toute valeur de || :
()
( )
Les fonctions de corrlations obtenues sont reprsentes par la (Figure - III-3).
Dans le cas o est petit et suffisamment grand, ces fonctions de corrlations
peuvent tre assimiles des impulsions de Dirac :

Figure - III-3 : Les fonctions de corrlations (a) continue, (b) discrte.
Continu :
()
( )
( ) ( )

Discret :
()
( )

Ce qui permet dassimiler les squences binaires pseudo-alatoires un bruit blanc
presque centr.

(a) (b)

CHAPITRE III
39
III.2 Les squences pseudo alatoires PN (Pseudo Noise) :
Les codes PN sont une classe de squence binaire pseudo alatoires avec de bonnes
proprits de la corrlation priodiques. Ces codes sont appels aussi les codes
longueur maximale (m-squences). Les proprits des m-squences sont rcapitules
dans les articles [31], [32].
Des registres dcalage linaire binaire (LFSR) [33] de n tages sont dcrits par des
polynmes dordre n, produisent une squence binaire priodique de priode N, qui
correspond aux polynmes primitifs caractristique, o N = 2
n
- 1. La reprsentation
binaire pour le registre qui correspond au polynme :
()
( )
Le vecteur binaire (
) avec une taille de , o

qui reprsente les deux prises de connexion qui correspondent au rebouclage du n
me

tage et 1
er
tage. Pour sil y a une prise de connexion vers ladditionneur
modulo 2 alors
,si non
. La squence de la sortie satisfait la relation de

rcurrence pour chaque :
( )
O cest laddition modulo 2. Le chargement initial de registre est spcifi par le
mot binaire
.
La figure suivante prsente un schma dun gnrateur de code PN { laide dun
registre LFSR de 7 tages.

Figure III-4 - Schma d'un gnrateur de code PN (7 3 0).
La Figure III-4 montre un gnrateur de SBPA avec trois prises de slection de
rtroaction. N'importe quel gnrateur de ce type est dcrit par un polynme primitif

31 (D. V. Sarwate, 1980)
32 (Fredricsson, 1975)
33 (Sloane, December 1976)

CHAPITRE III
40
qui dpend de ces raccordements. Le nombre de m-squence qui existe pour une
longueur de registre donne se dveloppe trs rapidement, mais cest extrmement
difficile de trouver les squences optimales. Le Tableau III-2 prsente le nombre des
m-squences et leurs tailles en fonction du nombre dtages du registre (de 13 17).
Taille du registre LFSR
N
Taille de la squence
N=2
n
-1
Nombre des squences
13 8191 630
14 16383 756
15 32767 1800
17 131071 7710
Tableau III-2 - La taille et le nombre des squences en fonction de la taille du registre
gnrateur.
Puisque les squences gnres par ce type de gnrateurs sont priodiques et font
un cycle pour chaque valeur binaire possible, le registre peut tre initialis
n'importe quel tats, sauf le vecteur nul, qui s'appelle le vecteur d'initialisation
bloqueur.
Les proprits des squences longueur maximale.
Les proprits des m-squences ont t formules par S. Golomb [34] dans son article
Shift Register Sequences . Celles-ci incluent :
Proprit d'quilibre. Le nombre des bits 1 du code est gal au nombre des
bits 0 ou plus grand dun seul bit.
Proprit de la gnration. Le vecteur de linitialisation permet de gnrer
plusieurs squences en utilisant le mme registre dcalage avec la mme
configuration.
Proprit de corrlation. L'autocorrlation et linter-corrlation du m-
squence est priodique.
Les squences de longueur maximale peuvent tre modlises par la
transformation de Hadamard.
La proprit la plus importante des codes PN cest la fonction dautocorrlation
(FAC).
Une FAC dune squence de code PN de taille
est :
() {

34 (Scholtz, 1978)

CHAPITRE III
41
Pour le cas de la fonction dinter corrlation entre deux squences x et y, elle prend 3
valeurs :
(
(())
)

(())

Les figures Figure III-5, Figure III-6 reprsentent respectivement la FAC et la FIC dun
code PN avec une taille de
.
Les deux squences utilises ont la configuration suivante :
Paramtres Polynme caractristique Etat initial du registre
Squence 01 [7 6 0] [0 0 0 0 0 0 1]
Squence 02 [7 3 0] [0 0 0 0 0 0 1]
Tableau III-3 - La configuration de la gnration des deux squences PN (N=127).

Figure III-5 - La fonction d'auto corrlation du code PN (N=127).

CHAPITRE III
42

Figure III-6 - La fonction d'inter corrlation entre deux codes PN (N=127).
Les valeurs de la FAC et de la FIC sont indiques dans le Tableau III-4 :
Niveaux FAC, FIC
PN (N=127)
Max Min
FAC 127 -1
FIC 15 -17
Tableau III-4 - Les maximums de FAC et FIC pour les codes de PN (N=127).
III.3 Les codes de Gold
Actuellement, les codes de Gold sont la famille de codes la plus populaire dans le
monde de la navigation par satellites. Les codes de Gold sont nomms du nom de leur
fondateur Robert Gold, qui les a prsents dans deux articles [35], [36]. Une famille
de code de Gold est base sur le produit de XOR de deux squences de PN de
longueurs identiques avec une fonction dinter corrlation qui prend trois valeurs
possibles :
() {
()
()

Avec ()
(())
, o n reprsente le nombre dtages des registres {
dcalages. Les codes de Gold de priode
existent pour tous les nombres

35 (Gold, October 1967)
36 (Gold, Jan 1968 )

CHAPITRE III
43
entiers positifs n. Chaque classe des codes de Gold de priode
se
compose de toutes les squences qui peuvent tre produites par un registre
dcalage linaire particulier qui a tages rpartis sur deux registres. Chaque
registre correspond une classe particulire de polynme primitif de degr n. Le
polynme qui produit les squences de Gold peut tre crit comme produit d'une
paire prfre [37] de polynmes primitifs de mme ordre n, o la multiplication et
l'addition de coefficient sont effectues par l'arithmtique modulo 2.

Figure III-7 - Un gnrateur de code de Gold 2n=20 tages.
La figure III-7, est un exemple de gnrateur de codes de Gold constitu de 2 registres
10 tages modliss par une paire prfr de polynmes :

()
()

Deux configurations sont utilises pour linitialisation des deux registres et une phase
de (1,1).
Avec cette arithmtique d'une paire prfre de polynmes, un premier code de Gold
est produit, le jeu de phases de XOR permet de gnrer plusieurs squences avec le
mme gnrateur et la mme configuration. Pratiquement la taille globale de la
famille de code est limite mais elle possde une bonne fonction dinter corrlation
dans le cas o aucun dcalage en frquence de Doppler n'est considr.
Les deux figures (Figure III-8 etFigure III-9) prsentent la fonction dautocorrlation
et de linter corrlation dun code de Gold avec n=10 ;
La fonction analytique de lauto corrlation est donne par la relation suivante.

37 (Chu, 1972)

CHAPITRE III
44
() {
(
(())
)

(())

Figure III-8 - La fonction d'autocorrlation d'un code de Gold avec n=10.

Figure III-9 - La fonction d'inter corrlation de deux squences de Gold avec n=10.
Avec ce modle de gnrateur de codes de Gold on peut gnrer des squences
binaires priodiques de taille et de priode
. Ce type de codes est

caractris par de bonnes fonctions dauto et dinter corrlation. Le tableau suivant
indique les maximums et les minimums de linter corrlation pour cette
configuration.

CHAPITRE III
45
Niveaux de FIC
Gold (n=10)
Max (dB) Min (dB)
Mode cohrent -11.9696 -30.0988
Mode non cohrent -11.9696 -12.1054
Tableau III-5 - Les maximums de l'inter corrlation pour les codes de Gold n=10.
III.4 Les codes de Kasami
Si
() est un polynme primitif de degr n, et
() un polynme primitif bien

choisi ([28], [31]) de degr , alors le produit () ()() permet de
gnrer un ensemble de
squences de priode avec une corrlation

priodique
()

Nous rfrons un tel ensemble de squences par la classe Small Kasami [38],
puisqu'il y a galement plusieurs classes des squences qui ont t dcouvertes par
Kasami (voient [39] pour plus de rsultat de Kasami).
Par exemple, pour , on peut multiplier le polynme
()
du et le polynme
()
pour obtenir le polynme ()
. Un registre dcalage linaire de 12

tages correspondant ce polynme produit un ensemble de 16 squences de
priode 255 avec une limite de corrlation priodique de 17.
Si on pose un nombre pair et une squence PN de priode
gnre
par (), alors on peut considrer la squence ,()- ,
-. Daprs
les relations de Kasami, toutes les squences ont une forme de (
) et donc chaque squence y de priode

gnre par h(x) appartient
() dfini par :

() {
} ( )
O reprsente l'addition modulo 2 du vecteur et une version dcale de
Cette formule indique que le nombre des squences du Small-Kasami est
certainement trs limit d la rptition du vecteur et .
Le maximum de linter-corrlation de ce type de code est compris entre les trois
valeurs [40] de lquation III.22, cette proprit se comporte comme une
caractristique pour des codes optimiss.

38 (Kasami, 1969)
39 (Fredricsson, 1975)

CHAPITRE III
46
() {
( )
Limplmentation des codes de Small set Kasami est plus complique puisque le
processus de la gnration utilise la dcimation dune squence par lui-mme ce qui
exige des horloges plus rapides. La solution pour rduire cette complexit est de
produire une nouvelle classe des codes de Kasami qui sappelle Kasami Larg set, cette
dernire est gnre la base des codes de Gold :

( )
Avec et tant une paire prfre des m-ordres et tant la version dcime d'un
des deux m-squences. La fonction dinter-corrlation de ce type de codes comporte
cinq valeurs :
() {
()

( )
O est le nombre de prises de rtroaction du registre et
indiquent les deux

codes qui produisent la squence de Larg Kasami.
La figure suivante prsente un exemple de la FAC dun code de type Kasami Small set
avec n=8 (
);
Le nombre des squences gnres par cette configuration est de 15 squences, les
trois niveaux de linter-corrlation pour ces codes sont :
() {

( )
Comme montre la Figure III-10, la FAC dun code de type Kasami est un pic centr
avec une dure de 2 bits et avec un niveau maximal de
. Le reste des
niveaux sont des niveaux qui varient entre (15,-1 et -17) en fonction de la valeur de ,
cest le mme cas pour la FIC sauf que pour ce cas, il ny a que les trois niveaux (15,-1
et -17) cause de labsence du pic de corrlation puisque les deux squences de
linter corrlation sont diffrentes comme montr dans la Figure III-11.

40 (Godfrey, 1966)

CHAPITRE III
47

Figure III-10 - La fonction d'autocorrlation du code Kasami pour n=8.

Figure III-11 - La fonction d'inter corrlation de deux codes de Kasami pour n=8.
Le tableau suivant prsente les maximums et les minimums de linter corrlation
pour deux types de rception cohrente et non cohrente.
Niveau de FIC
Kasami (8 4 3 2 1)
Max (dB) Min (dB)
Mode cohrent -11.7609 -24.0654
Mode non cohrent -11.7609 -12.3045
Tableau III-6 - Les maximums de l'inter corrlation pour les codes de Kasami n=8.

CHAPITRE III
48
III.5 Les codes de Walsh-Hadamard
Un autre type de codes dtalement avec une meilleure proprit dinter corrlation
sont les codes de Walsh Hadamard. Dans les prochaines sections, on va tudier les
proprits de corrlation des codes de Walsh-Hadamard et galement leurs
performances.
Pour ce type de codes les
codes de Walsh-Hadamard acquis par la matrice de

Walsh-Hadamard [41] (
), seulement codes sont orthogonaux pour toutes

les phases possibles.
III.5.1 Les proprits des codes de Walsh-Hadamard
Les codes de Walsh-Hadamard sont des codes orthogonaux [42] utiliss comme des
codes dtalement dans les systmes CDMA, pour obtenir cette proprit
dorthogonalit, il faut bien choisir ces codes { partir des ranges ou des colonnes
dune matrice de Hadamard.
Une matrice de Hadamard est une matrice carrs de (
) dans les lments

prennent les valeurs * +.
Pour gnrer une grande matrice de Hadamard, on utilise le processus de Walsh-
Hadamard, on fait une rptition de la matrice de dpart H(0) jusqu' atteindre la
taille de (
) :
() [
( ) ( )
( ) ( )
] ( )
Voici quelques proprits de la matrice de Hadamard :
1. Toutes les vectrices lignes ou colonnes de la matrice de Hadamard sont
orthogonaux, si :
(
) (
) ( )
Donc, X et Y sont deux colonnes ou lignes spares de la matrice de Hadamard. Alors,
la corrlation entre ces deux vecteurs est nulle :
()
( )

41 (Beer, May 1981)
42 (Beauchamp, 1984)

CHAPITRE III
49
Lquation prcdente est un cas particulier de la fonction de la corrlation
discrte qui est dfinie par:
()
( )

Alors, on peut dduire que linter corrlation de toutes les squences possibles de la
matrice de Hadamard et pour une phase est toujours gale 0.
2. Pour le cas o , on a :
( )

Vu la dfinition de la FAC en ( 31), la fonction d'autocorrlation phase nulle de
chaque range ou colonne d'une matrice de Hadamard est :
()
( )
Considrant cette proprit, tous les codes de Walsh-Hadamard ont une mme
nergie [43].
3. La matrice transpose de la matrice de Hadamard

O :
( ) {

Par consquent, les lments de la diagonale de la matrice Z sont la fonction
d'autocorrlation des diffrentes ranges et colonnes de la matrice H pour une phase
0, les autres lments de la matrice Z sont la FIC des ranges ou colonnes spares de
la matrice H pour une phase 0.
4. Proprits de la corrlation des codes de Walsh-Hadamard
Comme on a vu, lune des proprits des codes de Walsh-Hadamard est la valeur nulle
des FICs. La Figure III-12 montre cette proprit pour une matrice de
Hadamard (
). D'ailleurs, la Figure III-12 montre les FICs de tous les codes de

Walsh avec une longueur de 64 bits. On voit que dans toutes les situations, la valeur
de la FIC est nulle pour le cas :

43 (Abouei, 2006)

CHAPITRE III
50

( )
(()

.

Figure III-12 La FIC entre les 64 lignes de la matrice H (64x64).
Le problme de ce type de code est les faibles performances en termes de la FAC, les
niveaux maximaux du lobe latral sont trs levs. La Figure III-13 reprsente la
fonction dautocorrlation pour index=64 de la matrice H(64,64).

Figure III-13 La FAC linaire pour index 64 du code Walsh Hadamard H(64x64).

CHAPITRE III
51
Les niveaux maximums et minimums du lobe latral de la FAC sont rsums dans le
tableau :
Niveau de FAC
H(64,64)
Corrlation linaire Corrlation circulaire
Max (dB) Min (dB) Max (dB) Min (dB)
Mode cohrent -3.0103 -3.0103 -3.0103 -3.0103
Mode non cohrent -2.0412 -18.0618 -3.0103 -12.0412
Tableau III-7 - Les maximums de l'autocorrlation pour les codes de Walsh-
Hadamard H(64,64).
III.6 Les codes OVSF :
La transmission multiple spectre tal comporte deux oprations. La premire est
une opration de synchronisation (canalisation), qui transforme chaque symbole de
donnes en un certain nombre de morceaux (chip). Le nombre des chips par symbole
de donnes est appel le facteur de ltalement. Les codes Orthogonal Variables
Spreading Factor sont utiliss comme codes de synchronisation pour assurer
l'orthogonalit entre les diffrents canaux de liaison montante et descendante. La
deuxime est une opration de brouillage, o un code de brouillage est appliqu un
signal dtalement. Les codes de Gold ou les squences pseudo alatoires avec une
inter-corrlation limite sont utiliss en tant que codes de brouillage pour prserver
l'aspect alatoire entre deux utilisateurs diffrents.
Dans les systmes CDMA, la gnration des codes OVSF [44] ce fait par des
assembleurs qui contiennent des circuits numriques base de bascules et de portes
logiques est un codeur de plusieurs coefficients qui sont produits par la matrice
suivante :
()
(
()
(
()
(
)
]
( )
: Reprsente ltiquette du code et qui est cod par la table de Gray (Tableau 1).
: Lordre du code dans les diffrents niveaux de larbre gnrateur.

44 (K. Okawa, 1998)

CHAPITRE III
52
L'arbre gnrateur des codes entiers d'OVSF est produit en stockant les
o n est
lindex le plus grand de i pour chaque niveau qui est dfini par son facteur
dtalement FE tell que
est la nouvelle taille de la squence fille. L'arbre des

codes d'OVSF qui prservent l'orthogonalit [45] entre les diffrentes squences filles
est montre dans la Figure III-14.

Figure III-14 - l'arbre gnrateur de code OVSF pour un facteur FE= 8.
Dans la Figure III-14, les codes OVSF sont indexs par
() o FE est le facteur de
ltalement et i le numro du code dans chaque niveau avec K possibilit, avec la
condition . Chaque niveau de larbre dfini une squence de code de
taille FE. Par consquent, il y a
squences de codes dtalement avec la longueur de
au i
me
niveau, On a
une matrice de la dimension
avec N ensembles
de codes binaires de (0,1) avec une longueur N pour chaque niveau.

() . Alors, la matrice gnratrice
peut tre produite partir

de la matrice prcdente
selon la relation ci-dessus :
()
()
()
()
( )
( )]
()
()
()
()
()
()
()
()
( )
( )
( )
( )
( )
tant donn que
()
est le complment binaire, par consquent, deux codes

quelconques de la mme couche dans l'arbre gnrateur sont orthogonaux. En outre,

45 (T. Minn, 2000)

CHAPITRE III
53
deux codes quelconques de diffrentes couches sont orthogonaux { lexception du
cas o lun deux est un code mre de l'autre.
En gnral, il existe plusieurs classes pour la gnration des codes OVSF[46]. La plus
utilise est celle qui produit des codes en combinant deux gnrateurs OVSF { laide
dune opration bien dtermine telle que la somme modulo 2, avec un dcalage
circulaire dont le pas est fix entre les deux squences des diffrents gnrateurs,
comme prsent dans la figure suivante :

Figure III-15 - La gnration des nouveaux codes OVSF.
Tous les codes produits de lalgorithme OVSF ne peuvent pas tre assigns d'autres
utilisateurs de la mme bande de frquence. Cette restriction permet de maintenir la
proprit de l'orthogonalit. Voici un exemple de la gnration des codes OVSF avec
une taille de 64 bit, les codes gnrs par lalgorithme de OVSF sont une matrice
carre (64 x 64) bit, Linter-corrlation entre toutes les squences de cette matrice est
toujours nulle. La figure suivante prsente la fonction dinter-corrlation entre une
squence A (index=5) et la somme des squences B (index 5).

Figure III-16 - La fonction d'inter-corrlation entre les codes OVSF (64 bit).

46 (Chen, 2005)
0 10 20 30 40 50 60 70
-1
-0.8
-0.6
-0.4
-0.2
0
0.2
0.4
0.6
0.8
1
N
i
v
e
a
u

d
'
i
n
t
e
r
-
c
o
r
r
l
a
t
i
o
n
Tho
L'inter-corrlation entre deux codes OVSF de tailles 64 bit

CHAPITRE III
54
Cette caractristique de linter-corrlation permet de classer les codes OVSF avec les
codes orthogonaux pour les applications de la transmission synchronise [47].
Le Tableau III-8 est une table de conversion du code binaire en code de Gray afin de
classer les diffrentes squences gnres par lalgorithme dOVSF.
Dcimale
Code en Binaire Code en Gray
X1 X2 X3 X4 Y1 Y2 Y3 Y4
0 0 0 0 0 0 0 0 0
1 0 0 0 1 0 0 0 1
2 0 0 1 0 0 0 1 1
3 0 0 1 1 0 0 1 0
4 0 1 0 0 0 1 1 0
5 0 1 0 1 0 1 1 1
6 0 1 1 0 0 1 0 1
7 0 1 1 1 0 1 0 0
8 1 0 0 0 1 1 0 0
9 1 0 0 1 1 1 0 1
10 1 0 1 0 1 1 1 1
11 1 0 1 1 1 1 1 0
12 1 1 0 0 1 0 1 0
13 1 1 0 1 1 0 1 1
14 1 1 1 0 1 0 0 1
15 1 1 1 1 1 0 0 0
Tableau III-8 - Table de conversion du code Binaire en code de Gray.
III.7 Les code de Golay
Les codes de Golay sont des codes correcteurs derreurs et aussi utiliss pour
ltalement de spectre. En 1948, Hamming a propos un code correcteur d'erreurs
trs important [48], Il a propos une matrice de contrle de parit pour dtecter et
corriger les erreurs de la transmission, cependant, cette dtection et cette possibilits
de correction est limite, elle ne peut dtecter que les erreurs de 2 bit. En gnral, la
matrice Golay (23.12) peut corriger les erreurs de chaque 3 bits. Le code de Golay
sappelle aussi le code parfait [49].
Le code de Golay est un code linaire binaire (23,12) dont la distance minimale est
gale 7. Ce code est le seul code parfait binaire avec les codes de Hamming. En effet,
pour N = 23, K = 12 et e = 3 on a lgalit :

47 (R.G. Cheng, 2000)
48 (Lin, 1983)
49 (Golay, 1949)

CHAPITRE III
55

( )
Les codes de Golay sont aussi des codes cycliques. La matrice gnratrice
systmatique du code de Golay (23,12) est la suivante :
()
(
( )
Cette matrice peut tre modlise par le polynme caractristique suivant :
()
( )
Cette matrice permet de coder chaque 3 bits de la squence informative pour
obtenir un mot de 23 bits.
Voici un exemple du ce type de codage :
On a un ensemble de 6 bits de linformation x=101001 qui peut tre reprsent par
les deux valeurs entires 5,1 alors la sortie du codeur sera :
Y=
A partir du code de Golay (23,12), il est possible de construire le code de Golay
tendu (24,12) en ajoutant un bit de parit. La distance minimale de ce code est gale
8. La matrice gnratrice non systmatique du code de Golay tendu (24,12) est la
suivante :
()
(
( )

CHAPITRE III
56
Le polynme gnrateur de cette matrice est le suivant :
()

Parmi les matrice gnres par Golay, la matrice de control M(20,8) prsente une
matrice de codage trs prcise et qui permet de corrger le maximum de bits de
linformation. En termes dtalement, cette matrice permet de coder chaque 8 bits de
linformation sur 20 bits o les 12 bits ajouts sont des bits de parit selon la matrice
M(20,8)qui est donne par :
()
(
( )
Le polynome gnrateur de cette matrice est :
()

Pour gnrer une squence code en Golay de taille L, on a besoin dutiliser un
ensemble de bits de linformation avec un nombre de bits gal ()
Le tableau suivant reprsente la configuration utilise pour gnrer une squence de
Golay avec une taille de 1040 bits.
Paramtres Valeurs
Taille de linformation 156 bit
Taille de la squence 1040 bit
Tableau III-9 - la configuration pour la gnration du code Golay 1040.
La FAC et la FIC de cette squence sont reprsentes dans Figure III-17 et la Figure
III-18.

CHAPITRE III
57

Figure III-17 - La fonction d'auto corrlation du code Golay 1040.

Figure III-18 - La fonction d'inter corrlation de deux codes Golay 1040.
Les valeurs maximales et minimales de la FIC sont indiques dans le tableau suivant :
Les niveaux Min (dB) Max (dB)
Lobe latral FAC -8.2391 -4.2375
FIC -7.8150 -4.4185
Tableau III-10 - Les valeurs max et min de la corrlation pour le code de Golay 1040.
III.8 Les codes de Barker
Les codes de Barker sont des squences pour le codage de la phase [50], utiliss dans
les systmes radar [51] pour dtecter et identifier la phase du signal. Les codes de
Barker sont un ensemble de 9 squences reprsentes dans la table suivante :

50 (Golomb N. C., 1994)

CHAPITRE III
58
Index La squence
2 [1 1]
2 [-1 1]
3 [1 1 -1]
4 [1 1 -1 1]
4 [1 1 1 -1]
5 [1 1 1 -1 1]
7 [1 1 1 -1 -1 1 -1]
11 [1 1 1 -1 -1 -1 1 -1 -1 1 -1]
13 [1 1 1 1 1 -1 -1 1 1 -1 1 -1 1]
Tableau III-11 - Liste des 9 codes principaux de Barker
Les codes de Barker ont des FACs avec des lobes latraux identique : ce sont les seuls
codes pour lesquels les lobes latraux de la FAC 0 Doppler ont le niveau de N=
(longueur de code). Malheureusement, ils sont trop courts pour quelques
applications - la plus grande rduction du lobe latral, pour N = 13, est de -22.3dB.
La figure suivante reprsente la FAC du code de Barker pour index =13.

Figure III-19 - La fonction d'auto corrlation pour le code de Barker 13.
Niveaux max et min
de FAC pour le
Code Barker 13
Lobe latral de la FAC
Max (dB) Min (dB)
-8.4510 -161.3548
Tableau III-12 - Les niveaux max et min de la FAC pour le code 13 de Barker.

51 (Poor, 1998)

CHAPITRE III
59
Cependant il est possible de combiner des codes de Barker pour crer de plus longs
codes avec de bonnes proprits [52].
Lors de la combinaison des codes de Barker lune des proprits importantes avec
deux combinaisons possibles (B[m,n], B[n,m]) qui permet de gnrer deux
squences
et
diffrentes.
Par exemple, pour gnrer une squence avec une taille de 15, on peut combiner les
deux squences 3 et 5 de Barker suivantes: B3=[1 1 -1] et B5=[1 1 1 -1 1]. La
squence gnre peut prendre les deux cas possibles :
, - ( )
, - ( )
Cette combinaison des codes sappelle le produit de Kronecker de deux vecteurs
.
On utilise ce principe pour gnrer des squences dune taille de 1183 en combinant
les deux squences B13 et B7. Pour ce choix on a 3 cas possibles:
et
. Pour comparer ces trois rsultats, on fait

une comparaison des FACs des diffrentes squences.
Les figures Figure III-20, Figure III-21 et Figure III-22, reprsentent les fonctions
dautocorrlation pour les trois configurations :

Figure III-20 - La FAC pour Barker (13 13 7).

52 (Friese, 1996)

CHAPITRE III
60


Le tableau suivant prsent les valeurs max et min du lobe latral de la FAC pour
chaque cas :
Le Code Max (dB) Min (dB)
Barker (13 13 7) -3.6798 -8.4510
Barker (13 7 13) -2.5147 -5.6110
Barker (7 13 13) -16.9020 -8.4510
Tableau III-13 - Les max et les min de la FAC pour les squences de Barker 1183.

CHAPITRE III
61
Rellement, cette arithmtique de la gnration des codes de Barker contient une
certaine priodicit dans les squences gnres ; pour voir cette priodicit, on
calcule la FAC bidimensionnelle des squences B (3 5 7), B(3 7 5), B(5 3 7), B(5 7 3),
B(7 3 5) et B(7 5 3) avec une taille de 105 bits qui est base sur la loi suivante :
( ) ,() ( ) ( )--
Les rsultats pour le code Barker (3 5 7) sont reprsents dans la Figure III-23 , les
max et les min des lobes latraux de la FAC pour les diffrentes squences de chaque
combinaison sont gaux.
Le Code Max (dB) Min (dB)
Barker (105) -4.7712 -6.9897
Tableau III-14 - Le maximum et le minimum de la FAC pour les squences Barker 105.
Les FAC des diffrent cas possible des squences de Barker (105) sont indiques dans
la figure III-23.

CHAPITRE III
62

Figure III-23 - Les FAC des diffrentes squences pour Barker (L=105).

CHAPITRE III
63

Figure III-24 - FAC bidimensionnelle du code Barker (7 5 3).
La Figure III-24 est une reprsentation qui permet de dterminer la priodicit dans
la squence gnre par Barker (7 5 3), on peut constater clairement quil y a deux
types de priodicit dans cette squence, la premire priodicit est celle qui est
indique par les deux zones entoures par des cercles en rouge, la deuxime est celle
situe { lintrieur de chaque cercle.
III.9 CODES DE WEIL ET LES SEQUENCES LEGENDEr
Les codes de Weil ont t engendrs par Andr Weil. Ils sont bass sur les squences
de Legendre qui sont obtenues partir du XOR de deux squences de Gold [35] [36],
de longueur
, et sur des registres dcalage linaires - Linear Feedback Shift

Registers - (LFSR), les squences de Legendre doivent tre de longueur primaires.
Cependant la densit des nombres primaires dans l'ensemble est plus leve que
cela de
qui est le problme des codes de Gold.

Les squences de type Legender sont bases sur les rsidus quadratiques d'un
nombre primaire (voir [53] [54] [55] [56] [57]). Pour un nombre primaire et un

53 (Kitabayashi S., 1992)
54 (Jensen, January 1988)
55 (Paterson, July 1997)

CHAPITRE III
64
nombre entier positif alatoire ( ), on dit que est le rsidu quadratique de
si l'quation
a une solution . On peut dterminer si est un rsidu

quadratique ou pas en calculant la valeur de l'expression ()
()

En effet, cette expression peut prendre seulement les deux valeurs +1 et -1 et la
squence alatoire de Legender est forme de la faon suivante :
,() () ( )-
La squence obtenue a des proprits trs importantes, elle est priodique de priode
, quilibre et caractrise par une relation symtrique au centre de la squence.
() ()
( ) ( )
Avec

( )
( )
Une autre manire pour gnrer des squences de Legendre, vient du fait que les
squences de Legendre existent pour toutes les longueurs qui sont primaires. Donc,
on peut construire les squences de Legender, en utilisant les symboles de la
squence (i/L) qui est dfinie comme :
(
) {

( )
Une squence de Legender est forme des symboles
tel que
) ( )
Selon la taille de la squence on distingue deux classes de squences de Legender :
Classe 1 : avec une FAC priodique :
() {

( )
Classe 2 : avec une FAC priodique :
() {

()
La fonction d'autocorrlation dune squence de code Legender de la classe 1 est
reprsente dans la Figure III-25.

56 (Ding C., May 1998)
57 (Green D.H., Marc 2001)

CHAPITRE III
65

Figure III-25 - La fonction d'autocorrlation d'une squence de code Legender de la
classe 1
Laddition en modulo 2(XOR) d'une squence de Legender avec une version dcale
de la mme squence permet de crer une squence W dune nouvelle famille de
codes nomme codes Weil.
() {
/
} ()
Les FACs et les FIC des squences de code Weil sont reprsentes dans la Figure
III-26 et la Figure III-27, Les valeurs max et min de ces corrlations sont inclues dans
le Tableau III-15:
0 1000 2000 3000 4000 5000 6000 7000 8000 9000
0
0.2
0.4
0.6
0.8
1
1.2
Temps d'autocorrlation
A
m
p
l
i
t
u
d
e

n
o
r
m
a
l
i
s
La fonction d'autocorrlation d'une squence de Legendre de Classe 1

CHAPITRE III
66

Figure III-26 - La fonction d'autocorrlation d'une squence de code Weil (L=1013)

Figure III-27 - La fonction d'inter corrlation de deux squences de code Weil (1013).
Niveau de FC
Weil 1013
FAC (Lob latrale) FIC
Max (dB) Min (dB) Max (dB) Min (dB)
Mode cohrent -12.0627 -12.2028 -12.0627 -12.2028
Mode non cohrent -12.0621 -30.0561 -12.0627 -30.0561
Tableau III-15 - Les max et min de FAC et de FIC pour les codes de Weil (1013).

CHAPITRE III
67
On note que les codes de Weil sont des codes dont les performances sont quivalentes
aux codes de Gold pour des tailles maximales, mais lavantage des code de Weil
comme on a dj{ cit est que larithmtique de Weil permet de gnrer des codes
avec des tailles relativement non limit par rapport aux tailles des codes de Gold, PN
(
), Kasami (
), etc. La figure suivante reprsente une

comparaison entre les tailles possibles gnres par Gold et Weil entre 1000 et
30 000.

Figure III-28 - Les diffrentes tailles gnrer par Weil et Gold entre [1000 et 30 000].
Larithmtique de Weil est la plus importante en termes de tailles des codes gnrs.
Le nombre des tailles (entre 1000 et 30 000) que Weil peut gnrer est de 3077
tailles. Cependant, Gold peut gnrer uniquement 5 tailles.
III.10 Les codes mmoires
Les codes alatoires ou dite les codes mmoire puisque ce sont des codes gnrs par
des mthodes statistiques et qui doivent tre enregistrs dans les mmoires internes
des rcepteurs.
Le problme des squences pseudo alatoires gnres par les arithmtiques
prcdentes est la limite du choix de la taille des codes, en plus, mme lide de
tronquer les squences pseudo alatoires de quelques chips dtruit compltement
leurs proprits de lautocorrlation et de linter corrlation ; en voici un exemple :

CHAPITRE III
68
On prend une squence de code de Weil (1013) avec 2 scnarios, le premier est de
tronquer cette squence en liminant 7 chips et le deuxime scnario est dajouter {
la squence originale 7 chips. Ensuite on calcule la FAC et la FIC pour voir la
dgradation cause par cette dformation.

Figure III-29 - l'effet de tronquer une squence de Weil 1013 de 7chips.

Figure III-30 - l'effet d'ajouter 7chips une squence de Weil 1013.
FC pour code de Weil
1013 tronqu
MAX (dB) Min (dB)
FAC
Cohrent -9.7729 -2.2898 -10.4836 -1.7192
Non Cohrent -9.7729 -2.2898 -27.0157 -3.0404
FIC
Cohrent -10.1137 -1.9490 -10.4836 -1.7192
Non Cohrent -10.1137 -1.9490 -27.0157 -3.0404
Tableau III-16 - les valeurs de la dgradation des niveaux max et min de la FAC et de
la FIC.
L'ide dutiliser des codes alatoires est expose par J. Winkel pour lobjectif de
produire une famille de codes qui accomplit les proprits dun bruit alatoire aussi

CHAPITRE III
69
possible pour nimporte quelle longueur de code. Les codes peuvent tre conduits
pour accomplir les proprits spciales telles que l'quilibre des squences o la
probabilit de 0 et 1 ne doit pas tre identique mais dans une marge bien dfinie et o
la proprit du lobe latral de la FAC est nulle (ASZ Autocorrelation Side lobe Zero).
La proprit de lASZ garantit que les valeurs d'autocorrlation de chaque code se
corrlant avec une version du mme code dcal dun chip est nulle.
Parmi les mthodes utilises pour gnrer des codes alatoires, linsertion des bits
1 ou 0 dune faon bien dfinie pour amliorer la FAC et la FIC des squences
pseudo alatoires. Le nombre de choix pour insrer le bit 0 ou 1 pour une famille
entire de code est
(n : est la longueur du code et M le nombre des codes) qui est

un nombre inimaginable. Par consquent, des algorithmes spciaux doivent tre
appliqus pour produire des codes alatoires plus efficaces et moins coteux.
La figure suivante reprsente les valeurs du lobe latral de la FAC pour un code
alatoire optimis et un autre non optimis avec une taille de 1023.

Figure III-31 - Les niveaux de la FAC d'une squence de code alatoire optimis et non
optimis.

100 200 300 400 500 600 700 800 900 1000
-0.2
-0.15
-0.1
-0.05
0
0.05
0.1
0.15
0.2
0.25
0.3
La FAC pour un code alatoire (1023) optimis et non optimis

Code Alatoire non optimis
Code Alatoire optimis

CHAPITRE III
70
III.11 Conclusion
Les codes sont des lments fondamentaux dans nimporte quel systme
CDMA, rellement un grand nombre de codes sont appliqus pour des diffrentes
applications bases sur cette technique. Dans ce chapitre, on a essay dtudier les
arithmtiques de la gnration dun ensemble de codes qui sont les codes
dtalement. Notre choix a t fait selon des caractristiques qui rpondent aux
besoins des systmes de navigation GNSS.
On a bas sur la complexit de la gnration, les tailles des codes gnrs et les
niveaux maximums et minimums de la FAC et la FIC.
Daprs les rsultats, certains codes tels que PN, Gold et Kasami ont de meilleures
performances en termes de la FAC, c'est--dire une meilleure dtection et aussi une
simple architecture en hardware et software dans la conception et limplmentation.
Mais dun autre ct, ils ne peuvent pas couvrir toutes les applications des GNSS
puisquils sont limits dans le choix de la taille des squences.
Lun des objectifs principaux des systmes GNSS est dobtenir des codes qui ont de
bonnes FICs car les interfrences dans les systmes de la navigation restent toujours
le problme principal qui gne la dtection des satellites et qui est devenu la bte
noire de ces systmes. Les codes de Walsh, Hadamard et OVSF, malgr quils sont des
codes parfaits avec une meilleure caractristique de linter corrlation (puisque ce
sont des codes orthogonaux qui possdent un niveau dinter corrlation nul) mais
avec une mauvaise FAC.
Les codes Golay et Barker sont des codes qui peuvent amliorer la dtection mme
dans les mauvais cas, puisque ces deux codes sont des codes correcteurs, o Golay
permet de corriger les erreurs de la transmission de linformation et Barker permet
damliorer la dtection en prsence du leffet Doppler. Malheureusement la taille de
ces codes est trs limite pour les utiliser seuls dans les systmes GNSS.
Donc, le choix dun meilleur code pour les applications GNSS nest pas une tche
facile, plusieurs critres et conditions sont imposs pour tester les performances des
diffrents codes. Dans la suite de ce travail, on va dterminer les critres quon va
utiliser pour choisir les meilleurs codes.

CHAPITRE III
71
III.12 Bibliographie
[27] Willett, M. (April 1976). Characteristic m-sequences. Mathematics of
Computation vol. 30, pp. 306-311 .
[28] McEliece, L. D. (1972). Weights of irreducible cyclic codes. IEEE Inform. Cont vol.
20, pp. 158-175.
[29]N. HADJ-SAID, A. A.-P. (2009). Mthodes de Dterminations des Polynmes
Primitifs. IEEE ,SETIT th 5 International Conference: Sciences of Electronic, 1-4.
[30] N. E. Bekir, R. A. (1978). Partial-period correlation properties of PN sequences.
National Telecommunications Conf Rec vol. 3, pp. 35.1.1-4.
[31]D. V. Sarwate, M. B. (1980). Crosscorrelation Properties of Pseudorandom and
Related Sequences. IEEE Trans. Commun. Technol., vol. COM-68, pp. 594-596.
[33] Sloane, M. W. (December 1976). Pseudo-random sequences and arrays.
Proceedings of the IEEE, vol. 64, pp. 1715-1729.
[34] Scholtz, W. G. (1978). Group characters: Sequences with good correlation
properties. IEEETrans. Inform. Theory, vol. VOI. IT-24 , pp. 537-545.
[35]Gold, R. (October 1967). Optimal binary sequences for spread spectrum
multiplexing. IEEE Transactions on Information Theory, vol. IT-13, pp. 619-621.
[36]Gold, R. (Jan 1968 ). Maximal recursive sequences with 3 valued recursive
crosscorrelation function(corresp). IEEE Transactions on Information Theory
Volume14 Issue 1 , pp 154-156.
[37]Chu, D. C. (1972). Poly-phase codes with good periodic correlation properties.
IEEE Trans. Inform. Theory, vol. IT-18, pp. 531-432.
[38]Kasami, T. (1969). Weight distribution of Bose Chaudhuri-Hocquenghem codes.
Combingtorial Mathematics and Its Applications. University of North Carolina Press:
Chapel Hill, NlC.
[39]Fredricsson, S. (1975). Pseudo-randomness properties of binary shift register
sequences. IEEE Truns. Inform. Theory, vol. IT-2 1, pp. 115-120.
[40]Godfrey, K. R. (1966). Three-level m-sequences. IEEE Electron. Lett. vol. 2, pp.
241-243.
[41]Beer, T. (May 1981). Walsh transforms. American journal of physics, Vol. 49,
Issue 5.

CHAPITRE III
72
[42]Beauchamp, K. (1984). Applications of Walsh and related functions - with an
introduction to sequency theory . Academic Press.
[43]Abouei, J. (2006). A Set of Cyclic Orthogonal Codes Acquired from Walsh-
Hadamard Matrix. Yazd University,Iran Electrical Engineering Department.
[44]K. Okawa, F. A. (1998). Orthogonal forward link using orthogonal multi spreading
factor codes for coherent DS-CDMA mobile radio. IEICE Trans. Commun. E81-B (4),
pp. 778779.
[45]T. Minn, K. S. (2000). Dynamic assignment of orthogonal variable-spreading-
factor codes in W-CDMA. IEEE J. Sel. Areas Commun. 18 (8), pp. 14291440.
[46]Chen, W.-S. E. (2005). Optimum OVSF code reassignment in Wideband CDMA
forward radio link. ELSEVIER Information Sciences 174 , pp 81101.
[47] R.G. Cheng, P. L. (2000). OVSF code channel assignment for IMT-2000.
Proceedings IEEE Vehicular Technology Conference (VTC00), vol. 3, pp. 21882192.
[48]Lin, S. a. (1983). Error Control Coding: Fundamentals and Applications. Prentice
Hall, Inc.
[49]Golay, M. (1949). Notes on Digital Coding. IRE.
[50]Golomb, N. C. (1994). On n-phase Barker sequences. IEEE Trans on Information
Theory, Vol 40/4, pp 1251-1253.
[51] Poor, X. W. (1998). Blind Multiuser Detection: a Subspace Approach. IEEE Trans
on Information Theory, Vol 44 , pp 677-691.
[52]Friese, M. (1996). Polyphase Barker sequences up to length 36. IEEE Trans on
Information Theory, 42/4, pp 1248-1250.
[53]Kitabayashi S., O. T. (1992). Prop-erty of the Legendre Subsequence. Singapore
ICCS ll SLTA .
[54]Jensen, H. E. (January 1988). Determination of the Merit Factor of Legendre
Sequences. IEEE Transactions on Information Theory, Vol. 34,No. 1.
[55]Paterson, K. G. (July 1997). New Binary Sequence Sets with Favourable
Correlations. ISLT 1997, Ulm, Germany.
[56]Ding C., H. T. (May 1998). On the Linear Complexity of Legendre Sequences. IEEE
Transactions on Information Theory, Vol. 44, No.3.
[57]Green D.H., G. P. (Marc 2001). Polyphase related-prime sequences. IEE Proc.-
Comput-Digit.Tech, Vol. 148, No. 2.

IV. Chapitre VI

Applications ltude de linter-corrlation et du temps
Chapitre IV
Les codes choisis par les
systmes GNSS

Plan du chapitre
1 Introduction
2 Les codes des systmes GNSS
3 La gnration des codes de GNSS
3.1 Les codes du systme GPS
3.1.1 Les codes CA de GPS
3.1.2 Les code P et P(Y) de GPS
3.2 Les codes de Galileo
3.3 Les Tiered-codes
3.4 Les nouveaux codes E1 OS de Galileo et L1C de GPS
4 Les critres des codes
4.1 La fonction dautocorrlation (FAC)
4.2 La fonction dinter-corrlation (FIC)
5 Conclusion
6 Bibliographie

Rsum

Les systmes de navigation GNSS sont bass sur des codes dtalement de
types spciaux. Ce chapitre est une description dtaille des diffrents codes
implments dans les systmes GNSS, en particulier dans les systmes GPS et
Galileo ancienne et nouvelle gnrations. Les critres de choix des nouveaux
codes pour les applications GNSS modernes ont t bien dcrits.

74
CHAPITRE IV
Les codes choisis par les systmes GNSS
IV.1 Introduction
Les systmes GNSS ont t incroyablement russis, avec des applications dpassant
ceux que les concepteurs imagins en 1973. Certaines de ces nouvelles applications
prolongent les frontires auxquelles on s'attend ce que ces systmes fonctionnent.
certains d'autres trouvent de nouvelles utilisations pour les signaux GNSS. Quelques
applications augmentent le systme pour amliorer la disponibilit ou la prcision
dans les environnements difficiles et d'autres visent simplement amliorer
l'exactitude dans des conditions standard. Nanmoins, en tout cas il est ncessaire de
dvelopper les signaux qui sont prsents, y compris le signal de code C/A, qui est le
signal dominant actuellement en service.
Malheureusement, le code C/A a t dvelopp un moment o les
techniques de traitement des signaux nont pas t ralisables avec la technologie
actuelle. Pour la conception du signal GNSS un compromis entre les conditions
contradictoires d'une acquisition rapide qui exige une longueur de code courte et des
bonnes proprits de corrlation qui exigent une longueur de code plus longue.
Actuellement, les codes C/A ne sont pas appropris aux applications modernises du
GNSS.
Le dveloppement de nouvelles techniques avec de nouveaux matriaux,
associ la conception assiste par les processeurs de plus en plus utiliss par les
rcepteurs GNSS, rend ncessaire de dvelopper des codes optimiss plus rapides qui
permettent une synchronisation rapide et prcise, durant ltape de lacquisition.
cette dernire fournit des informations sur la visibilit et la disponibilit des
satellites, sur la frquence de Doppler et sur la phase des codes des satellites
permettant ainsi de rduire le temps de recherche pendant le processus de poursuite
des signaux reus.
Dans ce chapitre nous allons faire une prsentation des codes proposs pour
les applications GPS et Galileo et qui permettent une meilleure rsistance vis--vis du
bruit, des Multitrajets et des interfrences.

75
CHAPITRE IV
IV.2 Les codes des systmes GNSS
Dans cette section, on va dcrire les diffrents codes utiliss par les systmes GNSS.
Pour le cas GPS, les satellites transmettent leurs signaux sur deux frquences
porteuses, L1 et L2. La porteuse L2 est tale par un code de prcision appel le code
P(Y) transmis 10,23 M chip/s pour une dure de 7jours. La porteuse L1 est modul
par deux codes, le code P(Y) et un code (C/A : Coarse /Acquisition) qui est transmis
un taux de 1,023 M chip/s et qui a une dure de 1ms. Le code P (Y) est disponible
pour les utilisateurs PPS, alors que le code de C/A est disponible pour les utilisateurs
SPS.
Un programme de modernisation du GPS permettra dtendre ses possibilits { de
nouveaux services par la mise en place de signaux supplmentaires.
SPS (Standard Positioning Service)
PPS (Precise Positioning Service)
Le systme GPS est donc dot progressivement de nouveaux codes, par exemple le
code M, sur les signaux L1 et L2, qui est rserv aux usagers militaires et les codes
L1C sur les nouveaux signaux L5 pour un usage civil.
Galileo est conu pour offrir six services aux usagers dots de rcepteurs. Ces services
sont fournis par lintermdiaire des signaux diffuss par lensemble des satellites de
la constellation. Chaque service Galileo diffuse des signaux (E5, E6 et L1) qui sont
spcifis par une modulation et un ensemble de codes spcifiques.
IV.3 La gnration des codes de GNSS
Actuellement les diffrents codes des systmes GNSS sont des codes produits par des
gnrateurs de Gold bass sur des registres LFSR.
IV.3.1 Les codes du systme GPS
Le grand succs des applications militaires et civiles du GPS a conduit les concepteurs
et les chercheurs dvelopper un programme de modernisation du systme GPS
permettant la fourniture de signaux et services supplmentaires tout en assurant le
maintien des services SPS et PPS existants.
Le code P(Y) du service PPS introduit certaines limitations par rapport aux besoins
des utilisateurs militaires du GPS, en particulier concernant la rsistance aux

76
CHAPITRE IV
brouillages, la vitesse dacquisition du code et du service et la capacit distinguer les
diffrents utilisateurs.
IV.3.1.1 Les codes CA de GPS
Les codes C/A sont des codes pseudo alatoire priodique de priode 1ms avec un
dbit de 1.023 M chip/S, gnrer par un gnrateur de Gold deux registres LFSR de
10 tage selon la configuration indiqu dans le Tableau IV-1.
Les paramtres La configuration
Polynme gnrateur 1

Etats initial de gnrateur 1 [1 1 1 1 1 1 1 1 1 1]
Polynme gnrateur 2

Etats initial de gnrateur 2 [1 1 1 1 1 1 1 1 1 1]
Les slecteurs de phases (S1, S2) diffrents pour chaque satellite

Tableau IV-2 - La configuration d'un gnrateur de code CA de GPS.
Pour les slecteurs de phases, chaque satellite a ca propre configuration. Ces
slecteurs sont un pair de valeurs entieres qui varient de 1 10, exemple : (2, 6) pour
Sat 1 et (3, 8) pour Sat 31.
Le principe de la gnration du code CA est montr dans la Figure IV-1. Ce mme
schma est utilis pour gnrer les 32 autres squences qui correspondent chaque
satellite de la constellation GPS. Les diffrents satellites sont identifis par la phase du
code
.
Figure IV-1 - Schma d'un gnrateur de code CA du system GPS.

77
CHAPITRE IV
Lavantage du code CA est la simplicit de son architecture, mais ce type de code nest
pas scuris, c'est--dire, nimporte le qui peut dterminer la squence gnrer par ce
modle, en ce basons sur les proprits des registres LFSR et un nombre ditrations
pour dterminer ltat initial des 10 tages.
Les figures suivantes prsentent, respectivement la squence module en NRZ, la
fonction dautocorrlation et la densit spectrale de puissance pour le satellite 11.

Figure IV-2 - la squence binaire gnre par le satellite 11.

Figure IV-3 - La FAC correspond la squence gnre par le satellite 11.

78
CHAPITRE IV

Figure IV-4 - La densit spectrale de puissance du code CA (11).
Daprs les figures Figure IV-3 et Figure IV-4 on peut dduire que le code CA ou les
codes de Gold, en gnral ont une meilleure performance en termes de fonction
dautocorrlation, ce qui permet une dtection parfaite pour un seuil minimal. Les
caractristiques de ces squences sont indiques dans le tableau suivant.
Code CA11 Squence FAC PSD
Paramtre
Taille
(chip)
Dure
(ms)
Nombre
bits 0
Nombre
bits 1
Max
(dB)
Min
(dB)
Fc
(MHz)
Bande
(MHz)
Valeurs 1023 1 512 511 -24.21 -23.93 5.11 0.8870
Tableau IV-3 - les caractristiques du code GPS CA (11).
IV.3.1.2 Les code P et P(Y) de GPS
Le gnrateur du code P et trs complexe, il est bas sur des gnrateurs pseudos
alatoires constitus de 4 registres LFSR configurs comme suit :
Les paramtres La configuration
Polynme gnrateur 1

Etats initial de gnrateur 1 [0 0 1 0 0 1 0 0 1 0 0 0]
Polynme gnrateur 2

Polynme gnrateur 2

Polynme gnrateur 2


Tableau IV-4 - la configuration du gnrateur de code P du GPS.

79
CHAPITRE IV

Figure IV-5 - Schma d'un gnrateur 4 registres de code P.
Les squences engendres par ce type de gnrateurs sont des squences binaires
pseudo alatoire non priodiques dont la dure est de 23 017 555.5 (s). Les
premires 32 squences sont utiliss par les satellites de la constellation GPS, les
autres sont rserves aux applications internes du system GPS. Lexploitation de ces
codes ncessite laccs { des informations particulires et du matriel (hardware)
adapt.
En plus des 4 registres, le gnrateur de code P utilise des processeurs numrique
pour contrler et coder la sortie de chaque registre, chaque registre est contrl par
un contrleur dhorloge CH. Pour chaque cycle dhorloge les phases des quatre
gnrateurs se change pour gnrer des squences plus longues.
Malheureusement, on na pas pu gnrer la taille maximale du code P, mais on a
gnr une squence avec une dure de 8 ms, les caractristiques de ce code sont
reprsentes dans le Tableau IV-5.

80
CHAPITRE IV
Les caractristiques du code P de GPS gnr
Code P Squence FAC
Paramtre
Taille
(chip)
Dure
(ms)
Nombre
bits 0
Nombre
bits 1
Max
(dB)
Min
(dB)
Bande
(MHz)
Valeurs 81900 8 41152 40748 -37.06 -36.88 1.57
Tableau IV-5 - les caractristiques de base du code P de GPS.

Figure IV-6 - les densits spectral de puissance des codes CA et P de GPS.

Figure IV-7 - l'enveloppe de la densit spectrale de puissance pour les deux code CA
et P de GPS.

81
CHAPITRE IV
Comme illustre la figure IV -7, la bande de frquence du code P et plus large que la
bande du code CA, ce dernier est caractris par une densit spectrale de puissance
troite.
Le Code P(Y) est une version crypte du code P, cod par un codage convolutionel
[58].
IV.3.2 Les codes de Galileo
Les codes utiliss dans le systme europen Galileo sont utiliss sur trois bandes E1,
E6 et E5 chaque bande a un ensemble de codes pour les diffrents signaux data et
pilotes.
Le tableau suivant contient les diffrentes caractristiques de ces codes.
Galileo
Frquence Signal Type Modulation
F_chip
(MHz)
Taille du
code
Dure
(ms)
E1
1575.42
A PRS CBOC (15,2.5) 1.023 25575*1 10
B DATA BOC(1,1) 1.023 4092*1 4
C Pilot BOC(1,1) 1.023 4092*25 100
E6
1278.750
A PRS CBOC (15, 5) 5.115 5115*10 10
B DATA BPSK(5) 5.115 5115*1 1
C Pilot BPSK(5) 5.115 10230*100 100
E5
1191.795
1176.450
1207.140
a_I DATA AltBOC(15,10) 10.23 10230*20 20
a_Q Pilot AltBOC(15,10) 10.23 10230*100 100
b_I DATA AltBOC(15,10) 10.23 10230*4 4
b_Q Pilot AltBOC(15,10) 10.23 10230*100 100
Tableau IV-6 - - les caractristiques des diffrents signaux Galileo.
La configuration des gnrateurs de quelques codes pseudo alatoires du systme
Galileo est donne dans les tableaux suivants :
Galileo E1 B BOC(1,1) (4092x1) chip
dure 4 (ms)
Polynme 1

Initialisation [1 1 1 1 1 1 1 1 1 1 1 1 1]
Polynme 2

Initialisation [1 1 1 1 1 1 1 1 1 1 1 0 0]
Tableau IV-7 - la configuration du gnrateur E1C Galileo.

58 (Wolf, Feb 1986)

82
CHAPITRE IV
Galileo E1 C BOC(1,1) (4092x4) chip
dure 100 (ms)
Polynme 1

Initialisation [0 1 0 0 0 1 0 1 1 1 1 1 1]
Polynme 2

Initialisation [1 1 1 1 1 1 1 1 1 1 1 1 1]
Tableau IV-8 - la configuration du gnrateur E1 B Galileo
Galileo E5 a_I AltBOC(15,10) (10230x20) chip
dure 20 (ms)
polynme 1

Initialisation [1 1 1 1 1 1 1 1 1 1 1 1 1 1]
polynme 2

Initialisation [1 1 1 0 1 0 1 0 1 1 1 1 1 1]
Tableau IV-9 - la configuration du gnrateur E5a I de Galileo
Galileo E5 a_Q AltBOC(15,10) (10230x100) chip
dure 100 (ms)
polynme 1

Initialisation [1 1 1 1 1 1 1 1 1 1 1 1 1 1]
polynme 2

Initialisation [1 0 0 0 1 1 1 0 1 0 1 1 0 0]
Tableau IV-10 - la configuration du gnrateur E5a Q de Galileo
Galileo E5 b_I AltBOC(15,10) (10230x4) chip
dure 4 (ms)
polynme 1

Initialisation [1 1 1 1 1 1 1 1 1 1 1 1 1 1]
polynme 2

Initialisation [0 1 0 1 0 0 0 0 0 1 0 1 1 1]
Tableau IV-11 - la configuration du gnrateur E5b I de Galileo
Galileo E5 b_Q AltBOC(15,10) (10230x4) chip
dure 4 (ms)
polynme 1

Initialisation [1 1 1 1 1 1 1 1 1 1 1 1 1 1]
polynme 2

Initialisation [1 1 1 0 0 1 0 0 1 0 1 1 1 0]
Tableau IV-12 - la configuration du gnrateur E5b Q de Galileo

83
CHAPITRE IV
IV.3.3 Les Tiered-codes
En principe le systme Galileo adopte une nouvelle technique pour la gnration des
codes. En effet, il utilise ce quon appelle les Tiered-codes, dont le principe est
dutiliser un code secondaire priodique bien dtermin pour crypter (moduler) le
code primaire comme l'illustre la figure suivante.
Un code primaire est constitu dun nombre N de priodes des codes pseudo
alatoires, le code secondaire est un code dans la frquence de la gnration
. Il permet de multiplier chaque priode (poche) du code primaire par un bit de 1

ou 0 (par addition modulo 2).

Figure IV-8 - Principe de la gnration des Tiered-Codes.
Le tableau suivant indique les tailles des codes secondaires correspond chaque
signal de Galileo.
Signal
La priode du
Tiered-codes (ms)
La taille du code (chip)
Code primaire Code secondaire
E5a-I 20 10230 20
E5a-Q 100 10230 100
E5b-I 4 10230 4
E5b-Q 100 10230 100
E1-B 4 4092 /
E1-C 100 4092 25
Figure IV-9 - les tailles des codes secondaire dans la configuration Galileo.

84
CHAPITRE IV
Les codes secondaires sont des codes dont la taille est constitue de 100 chips, leur
rle est lidentification des signaux Galileo permettant ainsi daugmenter la vitesse de
lacquisition par le concept (la recherche indexe ou la recherche raccourci).
IV.3.4 Les nouveaux codes E1 OS de Galileo et L1C de GPS
Les codes E1 OS et L1C emploient la mme modulation dans les bandes E1/L1. Mais,
comme le montre le Tableau IV-13 les systmes Galileo et GPS ont diffrentes ,
configurations, ayant de ce fait, comme consquence, diffrentes longueurs des codes
primaires.
Les parameters de code Galileo E1 OS GPS L1 C
Primary code length 4092 10 230
Secondary code length 25 1800
Number of codes 50 2 = 100 210 2 = 420
Symbol rate 250 symb/s 100 symb/s
Tableau IV-13 - Les codes E1 OS de Galileo et L1C de GPS.
La famille des codes de L1C comprend les 420 codes suivants :
63 codes pour les signaux (data) de L1C
63 codes pour les signaux pilotes de L1C
56 codes pour les signaux futurs (data) de GBAS
56 codes pour les signaux futurs (pilote) de GBAS
39 codes pour SBAS (data)
39 codes pour SBAS (pilote)
52 codes pour les signaux (data) pour d'autres applications de GNSS
52 codes pilotes pour d'autres applications de GNSS
Tandis que la famille des codes E1OS de Galileo se compose de :
50 codes pour les donnes de E1 OS
50 codes pour le pilote de E1 OS
Un signal est dit pilot s'il ne contient pas de donnes. Il est data s'il contient des
donnes.
Daprs le Tableau IV-13, nous pouvons voir que la famille des codes L1C pour le
systme GPS utilise des codes plus longs, ce qui nous permet de dduire directement
que les codes L1C ont de meilleures proprits de corrlation.

85
CHAPITRE IV
IV.4 Les critres des codes
Le choix dun meilleur est bas sur la slection de celui qui prsente de meilleures
performances par rapport lapplication et aux exigences demandes. Ce n'est pas
une tche facile, par exemple, le dveloppement des codes pour Galileo dcrit dans
[59] montre que la qualit d'un code devrait indpendante de la ralisation et de
limplmentation des rcepteurs. En effet, concevoir des codes optimiss pour toutes
les applications potentielles de GNSS est pratiquement impossible. La loi de Welch a
connu une grande importance ces dernires annes en tant que solution approprie
pour valuer les codes pseudo alatoires.
IV.4.1 La bande de Welch
Le thorme de Welch pour une bande minimale est dfini [60] par :

( )
Cette valeur caractrise le minimum thorique de la valeur maximale de linter
corrlation qui peut tre obtenue pour un code de longueur dans un ensemble de
codes de . Cette expression peut tre encore simplifie quand le nombre de
squences est relativement haut. Dans ce cas, la bande de Welch peut tre
approxime par
, ou bien aprs normalisation de la fonction

dautocorrlation par rapport la taille du code qui est lamplitude du pic principal
par:
. Des expressions plus prcises sur cette limite ont t obtenues

dans [61]. Il claire que la bande de Welch dpend directement de la longueur du code.
Pour l'optimisation des codes de Galileo, diffrents algorithmes ont t proposs pour
recouvrir les diffrentes utilisations de Galileo l'avenir [62].

59 (Soualle F, July 2005)
60 (L.Welch, May 1974)
61 (Tang X. H., March 2000)
62 (G.W. Hein, April 2006)

86
CHAPITRE IV
IV.4.2 La fonction dautocorrlation (FAC)
()

( )
Avec :
{
[
( )
]

La taille de la squence du code.
Le k
me
chip du code.
La frquence centrale de Doppler.
La frquence dchantillonnage.
Pour un nombre M de codes :
((
* +
La FAC pour les deux modes paire et impaire est un moyen efficace pour estimer la
corrlation entre le signal reu et le signal local. Pratiquement pour les signaux pilot,
une FAC paire offre une bonne estimation, puisquil ny a pas de changement de signe
tout le long du signal reu (N priodes). Pour les signaux data moduls par un signal
pilot conduisant ainsi un changement de signe entre deux priodes successives, ce
changement de signe nest corrig que par le mode impair.
IV.4.3 La fonction dinter-corrlation (FIC)
()

( )
La fonction d'autocorrlation et dinter-corrlation sont deux reprsentations qui
contiennent beaucoup d'informations sur les proprits du code et qui permettent de
juger leurs performances. Ceci est la base pour l'optimisation des codes L1C GPS et E1
OS Galileo.
En outre, plusieurs recherches scientifiques ont galement proposs, de nouveaux
facteurs pour identifier les meilleurs codes. Tel que le facteur dantibrouillage tudi

87
CHAPITRE IV
dans [63],ainsi que des travaux sur des squences des codes avec de grandes zones de
corrlation nulles- (ZCZ) qui reste toujours un sujet qui peut tre exploits par les
recherches dans le domaine GNSS.
En se basant sur ces critres les concepteurs ont propos d'utiliser les codes de Weil
et les codes m+-sequences pour les applications GPS. En outre, ils ont propos
d'utiliser les codes alatoires et les Teired codes pour les applications Galileo.
On rsume dans le tableau suivant les types des codes allous aux systmes GNSS
nouvelles gnrations.

Signaux Galileo GPS
E1/L1 Les codes alatoires Les codes de Weil
E5/L5 Les Teired codes m- sequences
Tableau IV-14 - les nouveaux codes choisis pour les nouveaux services de GNSS.

IV.5 Conclusion
Dans ce chapitre nous avons prsent les diffrents codes utiliss pour couvrir les
besoins des systmes GNSS. Les codes de Gold sont les codes les plus utiliss dans les
deux systmes GPS et Galileo. Les proprits souhaitables des codes de Gold
permettent une meilleure performance pour la navigation. Les besoins croissants
des tailles aussi longs pour les applications modernes a conduit les concepteurs de
trouver de nouvelles familles de codes.
Les solutions proposes pour les nouveaux signaux E1 OS Galileo et L1C GPS est
dutiliser les codes alatoires et les codes de type Legender (Weil).

63 (MacMullan S. J.)

88
CHAPITRE IV
IV.6 Bibliographie
[58] Wolf, H. M. (Feb 1986). On Tail Biting convolutional codes. IEEE Transactions on
Communications, Vol. COM-34, No. 2, pp. 104-11.
[59] Soualle F, S. ,.-R. (July 2005). Spreading Code Selection Criteria for the Future
GNSS Galileo. Proceedings of GNSS 2005,Munich 19-22 .
[60] L. Welch. (May 1974). Lower bounds on the maximum cross correlation of
signals. IEEE Transactions on Information Theory, Vol. IT-20, No. 3, pages 397-399.
[61] Tang X. H., F. P. (March 2000). Lower bounds on correlation of spreading
sequence set with low or zero correlation zone. Electronic Letters, Vol. 36, No. 6, pp
210-218.
[62] G.W. Hein, J. A.-R. (April 2006). MBOC: The New Optimized Spreading Modula-
tion Recommended for GALILEO L1 OS and GPS L1C. Proceedings of IEEE/ION PLANS
2006 4-27.
[63] MacMullan S. J., C. O. (s.d.). A Comparison of Known Codes, Random Codes, and
the Best Codes. IEEE Transactions on Information Theory.

V. Chapitre V
CHAPITRE
Chapitre V
Simulation et rsultats

Plan du chapitre
1. Introduction
2. L'organigramme des tapes suivies pour la simulation
3. Les conditions sur les codes dtalement
4. leffet de la taille sur les niveaux du lobe latrale pour chaque
type de code
5. Leffet de la modulation sur les performances des codes
6. Leffet du bruit
7. Lestimation moyenne de la FC (Accumulation)
8. Les statistiques de la dtection
9. Comparaison entre les codes dtalement
10. Leffet du mode de la dtection
11. Amlioration de performances des codes :
12. La fonction dinter-corrlation des codes CA (GPS L1)
13. Ltude de lintra systmes
14. Conclusion

Rsum

Dans ce chapitre on va faire une tude dtaille sur les codes en termes de
gnration et en caractristiques de performances dans les applications GNSS.
Les rsultats de simulations des performances seront introduits et discuts, on
value par simulation, les performances des codes en termes de FAC, FIC et
dtection pour les applications GNSS ainsi ltude des intra-systmes L1C de
GPS et E1 OS de Galileo.

CHAPITRE V
Rsultats et simulations
91
V.1 Introduction
Dans ce chapitre on va faire une tude dtaille sur les codes en termes de leurs
gnrations, leurs caractristiques et de leurs performances dans les applications
GNSS. L'tudes des performances, base sur les rsultats de simulations, sera
introduite et discute, en termes de FAC, de FIC et de dtection.
V.2 L'organigramme des tapes suivies pour la simulation
En gnral, les tapes suivies sont prsentes dans lorganigramme suivant :

Figure V-1 - L'organigramme qui rsume les tapes suivie pour la simulation.
Dans ce qui suit nous allons dtailler les diffrentes dmarches pour l'tude ralise.
Pour effectuer les comparaisons entres les diffrents codes, nous allons
d'abord les deviser en plusieurs familles. Par la suite, nous allons employer les
fonctions d'autocorrlation et d'inter-corrlation pour caractriser la dfrence entre
ces familles de codes. En effet, la comparaison directe des diffrents codes n'est pas
possible cause de la diffrence en terme de taille.
Pour vrifier les performances des codes, une tude statistique des codes,
dans des conditions diffrentes, est ralise. Cette tude est base sur les quantits
suivantes:

- Variance,
- La STD,

CHAPITRE V
92
- Le RMSE,
- La probabilit.
D'autres critres de base sont utiliss pour tudier linfluence des autres systmes de
golocalisation (Intra-systmes) sur le systme considr.
V.3 Les conditions sur les codes dtalement
Pour que les codes rpondent aux besoins des applications GNSS, ils doivent avoir les
proprits suivantes :
- Ils doivent tre dterministes.
- Le niveau de la corrlation entre deux codes diffrents (Inter corrlation)
doit tre trs faible.
- La corrlation dun code avec sa version retarde (autocorrlation) doit tre
faible partout sauf au point de lalignement parfait.
- Puisque les codes utiliss pour les systmes GNSS sont priodiques, les
longueurs de ces derniers devront avoir une grande importance.

CHAPITRE V
93
V.4 Leffet de la taille sur les niveaux du lobe latrale pour
chaque type de code
Lors de la gnration des codes, la taille du code est un paramtre essentiel. Avons de
faire la comparaison des niveaux dinter corrlation des diffrents codes, on propose
ltude de leffet des longueurs sur les performances des codes.
Pour raliser cette tape, on doit gnrer plusieurs codes pour diffrentes valeurs
des tailles. Puisque chaque code est caractris par sa propre taille, alors l'ide est de
gnrer plusieurs codes avec des tailles trs proches l'une par rapport l'autre. Les
tableaux suivants illustrent respectivement les diffrentes tailles choisis et leur
configuration de gnration.
Index 1 2 3 4 5 6 7
PN 15 63 255 1023 4095 8191 16383
Gold 15 63 255 1023 4095 8191 16383
Alatoire 16 64 256 1024 4096 8192 16384
Weil 17 61 257 1021 4093 8191 16381
Tableau V-1 - l'ensemble des tailles de rfrences.
Index 1 2 3 4 5 6 7
Polynme
s
, - , - , - , - , - , - , -
Taille 15 63 255 1023 4095 8191 16383
Tableau V-2 - La configuration de la gnration de diffrentes tailles du code PN.
Pour la gnration des codes de Gold et PN, le problme principal est de trouver les
paires prfres des polynmes caractristiques, ce qui limite le choix de la taille.
La proprit intressante des codes alatoires, cest quon peut les gnrer avec
nimporte quelle taille. En revanche, les codes de Weil ont des tailles primaires.
Pour voir leffet de la taille, on procde de la manire suivante:
On calcule les niveaux maximum et minimum des lobes secondaires pour chaque
taille. Les rsultats sont illustrs dans les figures suivantes.

CHAPITRE V
94

Figure V-2 - les niveaux maximum cohrents du lobe latral des FAC pour diffrents
type de codes.

Figure V-3 - les niveaux minimum cohrents du lobe latral des FAC pour diffrents
type de codes.

CHAPITRE V
95

Figure V-4 - les niveaux maximum non cohrents du lobe latral des FAC pour
diffrents type de codes.

Figure V-5 - Figure 4 - les niveaux minimum non cohrents du lobe latral des FAC
pour diffrents type de codes.
Comme l'illustrent ces figures, plus la taille des codes est longue, plus les niveaux
maximum et minimum des lobes latraux sont rduits. Par consquent, on peut

CHAPITRE V
96
dduire que le choix des codes avec des tailles longues s'avre important pour avoir
une meilleure prcision.
V.5 Leffet de la modulation sur les performances des codes
En gnral, les systmes GNSS transmettent leurs signaux sur des mme bandes de
frquence, ce qui provoque une inter corrlation entre les diffrents signaux. Pour
sparer les signaux de chaque systme et ceux des diffrentes applications, les
systmes GNSS utilisent des techniques de modulation diffrentes comme la
modulation BPSK et la modulation BOC.
Pour tudier l'effet de la modulation, on a choisi un code pseudo alatoire avec une
taille de 4095 chip. Ce code est modul diffrents types de modulation savoir BPSK,
BOC(1,1), MBOC (6,1,1/11) et TMBOC(6,1,4/33).
Les niveaux max et min des lobes latraux de la fonction dautocorrlation, pour des
configurations cohrente et non cohrente, sont prsents dans le tableau suivant :
Modulation
Max coh
(dB)
min coh
(dB)
Max non coh
(dB)
min non coh
(dB)
BPSK -26.27 -25.68 -25.6775 -inf
BOC(1,1) -26.27 -25.68 -25.6775 -inf
MBOC (6,1,1/11) -27.84 -27.25 -27.25 -345.98
TMBOC (6,1,4/33) -28.35 -27.76 -27.76 -141.86
Tableau V-3 - les niveaux max et min du lob latrale FAC pour diffrent type de
modulations.
Comme l'illustre ce tableau, on remarque que les codes moduls avec les modulations
BPSK et BOC(1,1) prsentent les mmes nivaux des lobes latraux. En revanche, les
modulations MBOC (6,1,1/11) et TMBOC(6,1,4/33) prsentent des rsultats qui sont
attirants. En effet, les niveaux max et min des lobes latraux sont rduits ce qui offre
une meilleure performance par apport aux modulations BPSK et BOC(1,1).
La figure suivante indique linfluence du type de modulation sur le pic de la FAC.

CHAPITRE V
97

Figure V-6 - l'effet de la modulation sur le pic principale de la FAC.
V.6 Leffet du bruit
La transmission par les ondes radios est toujours affect par la prsence du bruit
dans la canal de transmission. Dans ce scnario on va faire l'tude des performances
en prsence d'un bruit blanc gaussien aditif AWGN.
Les paramtres de la simulation sont indiqus dans le tableau suivant. les rsultats de
cette simulation sont illustrs sur les figures V-7 et V-8.
Paramtres Valeurs
Frquence dchantillonnage Hz
Frquence du chip Hz
Type de code Code pseudo alatoire PN
Taille du code
chip
Nombre de priodes 1 priode
SNR -30 dB
Tableau V-4- - Les paramtre de configuration pour tudier l'effet du bruit sur la
dtection

Figure V-8 - la FAC du code PN en absence du bruit.

CHAPITRE V
98

Figure V-9 - L'effet du bruit sur la qualit de la dtection d'un code pseudo alatoire.
Il est vident que la dtection du pic principal dans ces condition (SNR=-30 dB) est
impossible, les niveaux max et min des fonctions de corrlation sont rsums dans le
tableau V-5.
Mesures FAC FIC Erreur
paramtres Max (dB) Min (dB) Max (dB) Min (dB) (m)
Sans bruit -30.0988 dB - -9.7245 -10.3215 0
Avec bruit 0.76663 0.03572 0.47379 0.31288 3577
Tableau V-6 - Les niveaux de la corrlation d'un code pseudo alatoire en prsence et
en absence u bruit.
V.7 Lestimation moyenne de la FC (Accumulation)
Parmi les solutions intressantes pour amliorer la qualit de la transmission en
prsence du bruit est dutiliser un filtre moyenneur. Ce scnario a pour objectif la
mise en vidence de leffet de laccumulation des priodes du code pour l'estimation
de la fonction dauto corrlation. Les paramtres de cette simulation sont rsums
dans le tableau suivant :
Paramtres Valeurs
Frquence du chip Hz
Taille du code
chip
Nombre de priodes 10 priodes
SNR -30 dB
Tableau V-7 - la configuration pour tudier l'effet de l'accumulation.
Les rsultats de la simulation sont indiqus dans la Figure V-10.

CHAPITRE V
99

Figure V-10 - l'effet de l'accumulation sur la fonction d'auto corrlation.
Comme l'illustrent ces deux figures, pour les mmes conditions prcdentes,
laccumulation des priodes nous a permis damliorer la qualit de la fonction de
corrlation et par consquent la position du pic principal. Cette amlioration est
rsume dans le tableau suivant.
mesures FAC (side lobe) FIC
paramtres Max (dB) Min (dB) Max (dB) Min (dB)
1 priode 0.76663 0.03572 0.47379 0.31288
10 priodes -4.2299 -4.331 -4.4355 -4.3556
Tableau V-8 - les rsultats de simulation pour un nombre de priode = 10
V.8 Performances de dtection du pic principal
Pour tudier les performances de dtection du pic principal, on propose le scnario
suivant. Pour se faire, on suppose une transmission dun signal pilot (code PN seul)
dans un canal de transmission AWGN avec un rapport SNR qui varie de -35 dB -10
dB. Les rsultats de simulation, en termes Variance, RMSE, STD et probabilit
derreur, sont reprsents dans les figures suivantes.
-4 -3 -2 -1 0 1 2 3 4
x 10
-4
-0.5
0
0.5
1
1.5
La fonction d'auto corelation signal sans bruit | Code PN 10 Accimulations
temps s
A
m
p
l
i
t
u
d
e
-4 -3 -2 -1 0 1 2 3 4
x 10
-4
-0.5
0
0.5
1
1.5
La fonction d'auto corelation signal bruit ( -30 dB ) | Code PN 10 Accimulation
temps s
A
m
p
l
i
t
u
d
e

CHAPITRE V
100
Les paramtres de cette simulation sont indiqus dans le tableau suivant :
Paramtres Valeurs
Frquence du chip Hz
Taille du code
chip
Nombre ditration 900
SNR [-30 :2 :-10] dB
Tableau V-9 - les paramtres de simulation pour l'tude de performances de la
dtection.
Dans cette simulation le nombre ditrations est choisis gal { 900.

Figure V-11 la variance de la dtection pour un code PN 1023 chip.
Comme l'illustre la Figure V-11, le systme bas sur des codes pseudo alatoires
atteint sa meilleure prcision partir de la valeur -25 dB de SNR.

CHAPITRE V
101

Figure V-12 la STD de la dtection pour un code PN 1023 chip.
Le mme rsultat peut s'observer pour la cas de la STD.
Les deux mesures (VAR et STD) sont utilises dans le cas de la prsence d'un seul
phnomne perturbateur (1 seul processus alatoire). En prsence de plusieurs
processus alatoires, on mesure le RMSE comme illustr dans la figure suivante.

Figure V-13 - L'RMSE de la dtection pour un code PN 1023 chip.
La probabilit derreur pour ce type de code est prsente dans la Figure V-14

CHAPITRE V
102

Figure V-14 - la probabilit de la dtection pour un code PN 1023 chip.
Pour faire une comparaison entre les diffrents types de codes, on va mesurer la VAR,
STD, RMSE et la probabilit derreur pour chaque code.

V.9 Comparaison entre les codes dtalement
La mme configuration du scnario prcdent est prise en compte. En premier lieu
l'exprience est ralise pour comparer les codes PN avec le code CA GPS.

Figure V-15 - la variances de la dtection pour les codes PN et CA (1023 chip).

CHAPITRE V
103

Figure V-16 La probabilit derreur de la dtection pour les codes PN et CA (1023
chip).
Nous remarquons, d'aprs cette figure, que la probabilit derreur pour les codes PN
approche celle du code CA-GPS. La mme observation est valable pour le cas de la
variance (Figure V-17).
Les figures "Figure V-18", "Figure V-19", "Figure V-20" et "Figure V-21" illustrent les
rsultats de comparaison entre les codes Gold, Kasami, Hadamard, Walsh, OVSF,
Barker et Weil. Les rsultats sont effectues en terme de variance, STD, RMSE et
probabilits derreur. Les courbes sont reprsentes en fonction du SNR qui vari
entre -35dB et -10dB.
-35 -33 -31 -29 -27 -25 -23 -21 -19 -17 -15 -13 -11
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1

Code CA
Code PN

CHAPITRE V
104

Figure V-18 - la variance des diffrents types de codes.

Figure V-19 - La STD des diffrents types de codes.

CHAPITRE V
105

Figure V-20 - Les statistiques de la dtection pour les diffrents types de codes.
La probabilit derreur des diffrents type de codes est prsente par :

Figure V-21 - la probabilit d'erreur pour les diffrents types de codes (1023 chip).
Comme l'illustrent toutes ces figures, on remarque que les codes de Walsh et
Hadamard sont meilleurs. Le code de Weil prsente une faible rsistance vis--vis du
bruit.

CHAPITRE V
106
Comparaison entre les deux codes CA et P de GPS
Les figures "Figure V-22", "Figure V-23"et "Figure V-24" illustrent les rsultats de
comparaison entre les codes C/A et P pour le systme GPS. Les rsultats sont
effectues en terme de variance, STD, RMSE et probabilits derreur. Les courbes sont
reprsentes en fonction du SNR qui varie entre -35 dB et -16 dB.

Figure V-22 - La STD de la dtection pour les deux codes CA et P de GPS.

CHAPITRE V
107

Figure V-23 Les RMSE de la dtection pour les deux codes CA et P de GPS.

Figure V-24 - Les probabilit derreur de la dtection pour les deux codes CA et P de
GPS.
Il est vident, d'aprs ces courbes, que le code P-GPS est meilleur que celui du code
CA-GPS.
Le mme scnario est utilis pour faire la comparaison entre les codes de nouvelles
gnrations GPS et Galileo savoir le code de Weil, le code mmoire et le code bas

CHAPITRE V
108
sur la FFT. Les rsultats des probabilits de dtections sont illustres dans les figures
Figure V-25 et Figure V-26.

Figure V-25 - La probabilit d'erreur pour les codes de Weil et les codes alatoires.
En comparons les rsultats, on peut dire que les codes alatoires sont les codes les
plus performants.

Figure V-26 - La probabilit d'erreur pour les codes de Weil, les codes alatoires et les
nouveaux codes gnrs.

CHAPITRE V
109
De plus, les codes bass sur lalgorithme de la FFT prsentent des performances trs
proches des codes alatoires.
V.10 Leffet du mode de la dtection
Dans cette section nous allons comparer les diffrents modes de dtection. En effet,
les deux modes simuls sont les modes "dtecteur seuil fixe" et "dtecteur seuil
adaptatif".
La dtection seuil fixe, consiste a compar les niveaux de corrlation avec un niveau
de rfrence fixe.
La dtection adaptative est une dtection dynamique qui dpend de la fonction de la
corrlation. Le seuil est calcul par la formule suivante :

(())
(())

La figure suivante illustre ce principe.

Figure V-27 - le principe de la dtection adaptatif.
La configuration de ce scnario est :
Paramtres Valeurs
Frquence du chip Hz
Taille du code
chip
Nombre de priodes 1 priode
Seuils ([0.1 0.2 0.3 0.4 0.5], adaptatif)
Tableau V-10 - Les paramtres de simulation de l'effet du mode de la dtection.

CHAPITRE V
110

Figure V-28 - les probabilits d'erreur pour un code PN pour diffrent seuils de
dtection.
Nous remarquons, d'aprs cette figure, que la dtection adaptative prsente le
meilleur choix.
V.11 Amlioration de performances des codes :
Pratiquement, les performances des codes gnrer par des oprations arithmtiques
ne sont pas des codes optimiss. Il existent plusieurs algorithmes pour amliorer ces
codes afin de les utiliss dans les applications GNSS.
La figure suivante prsente lorganigramme des deux algorithmes utiliss pour
amliorer les performances des codes en terme de fonction de corrlation.

Figure V-29 - l'algorithme 1 pour amliorer les performances des codes pseudo
alatoires.

CHAPITRE V
111

Figure V-30 - l'algorithme 2 pour amliorer les performances des codes pseudo
alatoires.
Les rsultats pour un exemple des codes mmoires est prsent dans la figure
suivante.

Figure V-31 - l'amlioration des performances dun code de mmoire (1023 chip).
La figure V-29 montre les niveaux max du lobe latral de la fonction dauto
corrlation en fonction du nombre doprations. Comme le montre cette figure,
l'algorithme utilis permet de minimiser le niveau max jusqu{ un niveau de -27.3 dB.
0 500 1000 1500 2000 2500
-27.4
-27.2
-27
-26.8
-26.6
-26.4
-26.2
Nombre d'oprations
M
A
X
(
L
o
b

l
a
t
r
a
l

F
A
C

)

e
n

d
B
Amlioration des squences pseudo alatoires

CHAPITRE V
112
Dans le deuxime exemple nous prsentons le rsultat de l'amlioration des codes de
Barker.

Figure V-32 - l'amlioration des performances du code de Barker (1023 chip).
Comme l'illustre cette figure, lamlioration apporte aux codes de Barker est trs
intressante. En effet, l'algorithme en question permet damliorer jusqu{ 73% de la
fonction dauto corrlation.
V.12 La fonction dinter-corrlation des codes CA (GPS L1)
En gnral un rcepteur GPS peut recevoir des signaux de diffrents satellites.
Lobjectif de ce scnario est l'tude de leffet des interfrences de plusieurs satellites
sur un satellite de la constellation. On calcule ainsi la fonction dinter-corrlation
entre le code local et le code reu qui est le rsultat de la superposition de plusieurs
signaux.
En premier lieu, nous prsentons le cas d'une rception en prsence de quatre
signaux de diffrents satellites de la constellation GPS. Le rsultat est illustr dans les
quatre plans de la Figure V-33.
0 500 1000 1500 2000 2500
-24
-22
-20
-18
-16
-14
-12
-10
-8
Nombre d'opration
M
A
X
(
l
o
b
e

l
a
t
r
a
l

F
A
C
)

e
n

d
B
Amlioration des codes de Barker (1023 chip)

CHAPITRE V
113

Figure V-33 - les fonctions d'inter corrlation des 4 satellites GPS [1 2 3 4].
On peut clairement constater, d'aprs cette figure, que le niveau de linter-corrlation
pour les diffrents cas ne dpasse pas le seuil de 35%. Cela implique qu'il ny a aucune
dpendance entre les codes choisis. Le tableau suivant prsente les niveaux max et
min de linter-corrlation.
Le cas Max (dB) Min (dB) ()
SAT 1 || 2 3 4 -16.8875 -16.5750 -51.3874
SAT 5 ||6 7 8 -16.4729 -15.9778 -46.2791
SAT 9 || 10 11 12 -17.5435 -17.4311 -62.3735
Tableau V-11 - Les niveaux d'inter-corrlation pour 4 stellites de GPS.
Dans le deuxime cas on suppose la rception de 8 signaux de 8 satellites diffrents
de la bande L1 GPS. Les rsultats de simulation sont prsents dans la Figure V-34 et
le Tableau V-12.

CHAPITRE V
114

Figure V-34 - les fonctions d'inter corrlation des 8 satellites GPS.
Le cas Max (dB) Min (dB) ()
1 -11.8429 -13.5834 -30.1847
2 -14.1304 -13.5834 -42.9144
3 -12.9214 -12.0365 -36.7240
Tableau V-12 - Les niveaux d'inter-corrlation pour 8 stellites de GPS.
Daprs les rsultats indiqus dans le Tableau V-12, on observe que les performances
sont proportionnelles aux nombres de satellites visibles. Plus le nombre de signaux
dinterfrence augmente, plus la dgradation est importante.
V.13 Ltude de lintra systmes
Pour le calcul de la fonction dinter corrlation intra-systme nous prsentons quatre
scnarios diffrents,
Le GPS interfra Galileo
Galileo interfra GPS
Le GPS interfra GPS
Galileo interfra Galileo
Pour chaque scnario on calcule les niveaux max et min de linter corrlation ainsi
que le niveau du lobe latral de la fonction dautocorrlation.
-4 -2 0 2 4
x 10
-4
0
0.5
1
Fonction d'inter-correlation satelitte : 9 || 10 11 12 13 14 15 16
-4 -2 0 2 4
x 10
-4
0
0.5
1
-4 -2 0 2 4
x 10
-4
0
0.5
1
-4 -2 0 2 4
x 10
-4
0
0.5
1
-4 -2 0 2 4
x 10
-4
0
0.5
1
-4 -2 0 2 4
x 10
-4
0
0.5
1
-4 -2 0 2 4
x 10
-4
0
0.5
1
-4 -2 0 2 4
x 10
-4
0
0.5
1

CHAPITRE V
115
Le calcul du maximum de linter-corrlation inter systmes prsent dans ce travail
est juste une approche. Des investigations plus prcises doivent tre effectues
renforcer les rsultats. Nanmoins, les rsultats de notre travail sont tout fait
reprsentatifs car cette valuation bas sur lauto corrlation.
Lapproche suivie pour raliser les scnarios prcdents est la suivante:
Pour calculer les FAC des diffrents scnarios, on utilise le principe suivant, les deux
tailles des codes GPSL1 et Galileo E1 sont dffirents, donc on utilise le principe de
lencapsulation.
,
Scnario 1 : Galileo interfre GPS

Figure V-35 - Intra - systmes Galileo E1 OS -> GPS L1C.
Figure V-36 - Les niveaux de la corrlation (a) paire et (b) impaire du scnario 1

(a)

(b)
0 2000 4000 6000 8000 10000 12000
-80
-70
-60
-50
-40
-30
-20
M
a
x
i
m
u
m

e
n

d
B

Les niveaux maximums de lobe latral FAC pair Scnario 1 0 Hz

FC pair scenario 1
FAC pair Weil 10230
0 2000 4000 6000 8000 10000 12000
-90
-80
-70
-60
-50
-40
-30
-20
M
a
x
i
m
u
m

e
n

d
B

Les niveaux maximums de lobe latral FAC impair Scnario 1 0 Hz

FC impair scenario 10
FAC impair Weil 10230

CHAPITRE V
116
Scnario
Doppler 0 Hz Doppler 1000 Hz
MAX FAC MAX FAC
Pair Pair Pair Impair
Scnario 1 -28,6 -29,1 -24,5 -25,2
Weil 10230 -31 -30,2 -28,2 -27,5
( ) -2,4 -1,1 -3,7 -2,3
Tableau V-13 Les rsultats du scnario 1 Galileo interfre GPS.
Scnario 2 : GPS interfre Galileo

Figure V-37 - Intra - systmes GPS L1C -> Galileo E1 OS.

(a)

(b)
Figure V-38 - Les niveaux de la corrlation (a) paire et (b) impaire du scnario 2
Scnario
MAX FAC (dB) MAX FAC (dB)
Scnario 2 -25,5 -26,1 -23,2 -22,5
L1 4092 -27,9 -29,4 -26,2 -25
( ) -2,4 -3,3 -3 -3,1
Tableau V-14 Les rsultats du scnario 2 GPS interfre Galileo.
0 500 1000 1500 2000 2500 3000 3500 4000 4500
-65
-60
-55
-50
-45
-40
-35
-30
-25
-20
M
a
x
i
m
u
m

e
n

d
B


FC pair scenario 2
FAC pair L1 4092
0 500 1000 1500 2000 2500 3000 3500 4000 4500
-80
-70
-60
-50
-40
-30
-20
M
a
x
i
m
u
m

e
n

d
B


FAC impair L1 4092

CHAPITRE V
117
Scnario 3 GPS interfre GPS :

Figure V-39 - Intra - systmes GPS L1C -> GPS L1C.

(a)

(b)
Figure V-40 - Les niveaux de la corrlation (a) paire et (b) impaire du scnario 3.
Scnario
Scnario 3 -29 -29 -25,2 -24,4
Weil 10230 -31 -30,2 -28,2 -27,5
() -2 -1,2 -3 -3,1
Tableau V-15 Les rsultats du scnario 3 GPS interfre GPS.

0 2000 4000 6000 8000 10000 12000
-75
-70
-65
-60
-55
-50
-45
-40
-35
-30
-25
M
a
x
i
m
u
m

e
n

d
B


FC pair scenario 3
FAC pair Weil 10230
0 2000 4000 6000 8000 10000 12000
-90
-80
-70
-60
-50
-40
-30
-20
M
a
x
i
m
u
m

e
n

d
B


FAC impair Weil 10230

CHAPITRE V
118
Scnario 4 Galileo interfre Galileo :

Figure V-41 - Intra - systmes Galileo E1 OS -> Galileo E1 OS.

(a)

(b)
Figure V-42 - Les niveaux de la corrlation (a) paire et (b) impaire du scnario 4.
Scnario
Scnario 4 -25,9 -27,3 -22,7 -21,2
L1 4092 -27,9 -29,4 -26,2 -25
( ) -2 -2,1 -3,5 -3,8
Tableau V-16 Les rsultats du scnario 4 Galileo interfre Glileo.
Discussion :
Pour les quatre scnarios, une dgradation considrable dans les performances caus
par linterfrence entres les codes des diffrents systmes (intra systmes).
Pour le cas o F doppler gale 0 Hz, le maximum de la dgradation ne dpasse pas la
valeur de 2 dB.
Les scnarios 2,3 et 4 ont une dgradation dans les performances de la FAC pair (F
Doppler = 1000 Hz) avec une marge de 3,8 dB.
Le seul cas o la dgradation des performances de la FAC impair est grande, cest le
cas du premier scnario.
0 500 1000 1500 2000 2500 3000 3500 4000 4500
-65
-60
-55
-50
-45
-40
-35
-30
-25
-20
M
a
x
i
m
u
m

e
n

d
B


FC pair scenario 4
FAC pair L1 4092
0 500 1000 1500 2000 2500 3000 3500 4000 4500
-80
-70
-60
-50
-40
-30
-20
M
a
x
i
m
u
m

e
n

d
B


FAC impair L1 4092

CHAPITRE V
119
V.14 Conclusion

Nous nous sommes intresss dans ce chapitre { ltude de performances de
certaines familles de codes des simulations pour tudier leurs caractristiques et
dterminer les critres de comparaison pour choisir les meilleurs codes. Les objectifs
de ce chapitre sont :
Gnrer le maximum de codes dtalement de spectre ainsi que les nouvelles
squences qui sont dveloppes.
Etudier et interprter les caractristiques de ces codes (les niveaux de la FAC, de
lFIC et aussi loccupation spectrale).
Etudier leffet de la taille des codes sur les performances des systmes CDMA en
prsence du bruit.
Afin damliorer les performances de ces codes, une tude sur les statistiques des
squences en terme de la variance, lSTD, RMSE et la probabilit derreur.
Comme dernire tape, une implmentation (par simulation) des codes alatoire
(4092 chip) et des codes de Weil (10230) dans les systmes de navigation GNSS
(GPS/Galileo) tout en essayant dvaluer les performances de ces codes pour les
diffrents scnarios de lintra-systmes.
D'une part, les codes de Weil sont capables de construire une grande famille de codes
avec une bonne proprit de corrlation, mais ils ne sont pas aussi flexibles que les
codes alatoires. En ce qui concerne limplmentation, les codes Weil peuvent tre
mis en application facilement tous comme les codes bass sur les LFSR.
Pour linter corrlation intra-systme, quatre scnarios ont t prsentes, les effets
inter-systmes sont compatibles aux effets alatoires de code.
En gnral, L'utilisation des nouveaux codes dtalement prsent de nombreux
avantages, notamment en termes de prcision, de rapidit et de rsistance aux
interfrences.

Table des matires Liste des figures Liste des tableaux
i
Table des matires i
Table des matires ii
Liste des figures iii
Liste des tableaux iv
Liste des abrviations v
Chapitre 0 INTRODUCTION GENERALE .............................................................. 4
I. Chapitre I Les systmes de navigation GNSS ....................................................... 1
I.1 INTRODUCTION ........................................................................................................................................................ 1
I.2 PRSENTATION DU SYSTME GNSS ......................................................................................................... 1
I.3 PRINCIPE DE LA LOCALISATION ............................................................................................................................ 3
I.4 ACTUALITES .............................................................................................................................................................. 5
I.5 LE RECEPTEUR GNSS ............................................................................................................................................. 5
I.5.1 ACQUISITION DES SIGNAUX GNSS 6
I.6 LES MODES DE LA DETECTION .............................................................................................................................. 9
I.7 LES SOURCES DERREURS DE LA LOCALISATION............................................................................................. 12
CONCLUSION ................................................................................................................................................................ 13
I.8 BIBLIOGRAPHIE ..................................................................................................................................................... 14
II. Chapitre II Ltalement de spectre par squence directe ..................................... 15
II.1 INTRODUCTION (UN CALCUL POUR UNE INTRODUCTION) ......................................................................... 16
II.2 METHODES DACCES MULTIPLES .......................................................................................................... 16
II.2.1 ACCES MULTIPLES PAR LA REPARTITION DE FREQUENCES (FDMA) 16
II.2.2 ACCES MULTIPLES PAR LA REPARTITION DANS LE TEMPS (TDMA) 18
II.2.3 ACCES MULTIPLE PAR LA REPARTITION DES CODES (CDMA) 19
II.3 LETALEMENT DE SPECTRE .......................................................................................................................... 22
II.3.1 INTRODUCTION 22
II.3.2 HISTORIQUE 22
II.3.3 LE SECRET DE LETALEMENT DE SPECTRE 24
II.3.4 LES PRINCIPES DE LETALEMENT DE SPECTRE 25
II.4 MODELISATION DU SIGNAL AU NIVEAU DE LEMETTEUR ET DE RECEPTEUR ........................................ 28
II.5 CONCLUSION ........................................................................................................................................................ 30
II.6 BIBLIOGRAPHIE ................................................................................................................................................... 31
III. Chapitre III Les codes de ltalement et leurs performances ................................ 32
III.1 INTRODUCTION .................................................................................................................................................. 33
III.2 LES SEQUENCES BINAIRES PSEUDO-ALEATOIRES ....................................................................................... 33
III.2.1 GENERATION DUNE SEQUENCE PSEUDO-ALEATOIRE 34
III.2.2 PROPRIETES STATISTIQUES 36
III.3 LES SEQUENCES PSEUDO ALEATOIRES PN (PSEUDO NOISE) : ................................................................ 39
III.4 LES CODES DE GOLD.......................................................................................................................................... 42
III.5 LES CODES DE KASAMI ..................................................................................................................................... 45

ii
III.6 LES CODES DE WALSH-HADAMARD ............................................................................................................. 48
III.6.1 LES PROPRIETES DES CODES DE WALSH-HADAMARD 48
III.7 LES CODES OVSF : ............................................................................................................................................ 51
III.8 LES CODE DE GOLAY ......................................................................................................................................... 54
III.9 LES CODES DE BARKER..................................................................................................................................... 57
III.10 CODES DE WEIL ET LES SEQUENCES LEGENDER ...................................................................... 63
III.11 LES CODES MEMOIRES ................................................................................................................................... 67
III.12 CONCLUSION .................................................................................................................................................... 70
III.13 BIBLIOGRAPHIE ............................................................................................................................................... 71
IV. Chapitre VI Les codes choisis par les systmes GNSS ...................................... 73
IV.1 INTRODUCTION .................................................................................................................................................. 74
IV.2 LES CODES DES SYSTEMES GNSS ................................................................................................................... 74
IV.3 LA GENERATION DES CODES DE GNSS .......................................................................................................... 75
IV.3.1 LES CODES DU SYSTEME GPS 75
IV.3.2 LES CODES DE GALILEO 81
IV.3.3 LES TIERED-CODES 83
IV.3.4 LES NOUVEAUX CODES E1 OS DE GALILEO ET L1C DE GPS 84
IV.4 LES CRITERES DES CODES................................................................................................................................. 84
IV.4.1 LA BANDE DE WELCH 85
IV.4.2 LA FONCTION DAUTOCORRELATION (FAC) 86
IV.4.3 LA FONCTION DINTER-CORRELATION (FIC) 86
IV.5 CONCLUSION ...................................................................................................................................................... 87
IV.6 BIBLIOGRAPHIE ................................................................................................................................................. 88
V. Chapitre V Simulation et rsultats .......................................................... 89
V.1 INTRODUCTION .................................................................................................................................................... 91
V.2 L'ORGANIGRAMME DES ETAPES SUIVIES POUR LA SIMULATION ............................................................... 91
V.3 LES CONDITIONS SUR LES CODES DETALEMENT ......................................................................................... 91
V.4 LEFFET DE LA TAILLE SUR LES NIVEAUX DU LOBE LATERALE POUR CHAQUE TYPE DE CODE ............ 93
V.5 LEFFET DE LA MODULATION SUR LES PERFORMANCES DES CODES ........................................................ 95
V.6 LEFFET DU BRUIT ............................................................................................................................................... 97
V.7 LESTIMATION MOYENNE DE LA FC (ACCUMULATION) ............................................................................. 98
V.8 LES STATISTIQUES DE LA DETECTION ............................................................................................................. 99
V.9 COMPARAISON ENTRE LES CODES DETALEMENT ..................................................................................... 102
V.10 LEFFET DU MODE DE LA DETECTION ........................................................................................................ 109
V.11 AMELIORATION DE PERFORMANCES DES CODES : .................................................................................. 110
V.12 LA FONCTION DINTER-CORRELATION DES CODES CA (GPS L1) ...................................................... 112
V.13 LETUDE DE LINTRA SYSTEMES ................................................................................................................. 114
V.14 CONCLUSION ................................................................................................................................................... 119
Chapitre VI. Conclusion gnrale 118
Rsum 120

iii
Liste des figures
FIGURE I-1 - LA CONSTELLATION DES SATELLITES AUTOUR DE LA TERRE. ........................................................................ 2
FIGURE I-2 - UN MECANISME SIMPLE DE LA RECEPTION DES SIGNAUX GNSS. .................................................................. 3
FIGURE I-3 - LE PRINCIPE POUR DETERMINER LA DISTANCE. ................................................................................................... 4
FIGURE I-4 - LE PRINCIPE DE LA TRIANGULATION. ...................................................................................................................... 4
FIGURE I-5 - LE PRINCIPE DE LA RECHERCHE BIDIMENSIONNELLE. ....................................................................................... 7
FIGURE I-6 - LE MECANISME DE LA RECHERCHE SERIE. .............................................................................................................. 8
FIGURE I-7 - LE MECANISME DE LA RECHERCHE PARALLELE (FFT). ..................................................................................... 9
FIGURE I-8 - LES DIFFERENTS MODES DE DETECTION. ............................................................................................................. 10
FIGURE I-9 - SCHEMA SYNOPTIQUE DE LA DETECTION DES SIGNAUX BINAIRES. ............................................................ 11
FIGURE I-10 - LA DISTRIBUTION DES SYMBOLES ET LA PRISE DE DECISION POUR LE CAS BINAIRE. ......................... 11
FIGURE II-1 : PARTAGE DES CANAUX EN FDMA ....................................................................................................................... 17
FIGURE II-2 : PARTAGE DES CANAUX EN TDMA ...................................................................................................................... 19
FIGURE II-3 PARTAGE DES CANAUX EN CDMA. ....................................................................................................................... 20
FIGURE II-4 - LE SCHEMA PROPOSE SOUS LE NOM SECRET COMMUNICATIONS SYSTEM. ..................................... 23
FIGURE II-5 - LE PRINCIPE DE LA CDMA..................................................................................................................................... 26
FIGURE II-6 - LA CHAINE DE COMMUNICATION PAR CODAGE CDMA. .............................................................................. 27
FIGURE II-7 - L'OPERATION DE L'ETALEMENT ET LE DES-ETALEMENT DE SPECTRE, LE SIGNAL ETALE EST
DEVENUE SOUS LE NIVEAU DU BRUIT AVEC UNE BANDE DE FREQUENCE PLUS LARGE. .................................... 27
FIGURE II-8 - LE CODAGE ET LE DECODAGE DES SIGNAUX INFORMATIF........................................................................... 27
FIGURE II-9 - STRUCTURE GENERALE DUNE CHAINE DE TRANSMISSION DS-SS. ......................................................... 28
FIGURE - III-1 : GENERATION D'UNE SBPA A BASE DUN REGISTRE A DECALAGE LINEAIRE LFSR. ...................... 35
FIGURE - III-2 : SIGNAL BINAIRE PSEUDO-ALEATOIRE. (A) REPRESENTATION CONTINUE, (B) REPRESENTATION
DISCRETE...................................................................................................................................................................................... 36
FIGURE - III-3 : LES FONCTIONS DE CORRELATIONS (A) CONTINUE, (B) DISCRETE. ..................................................... 38
FIGURE III-4 - SCHEMA D'UN GENERATEUR DE CODE PN (7 3 0). ........................................................................................ 39
FIGURE III-5 - LA FONCTION D'AUTO CORRELATION DU CODE PN (N=127). .................................................................. 41
FIGURE III-6 - LA FONCTION D'INTER CORRELATION ENTRE DEUX CODES PN (N=127). ............................................ 42
FIGURE III-7 - UN GENERATEUR DE CODE DE GOLD A 2N=20 ETAGES. ............................................................................. 43
FIGURE III-8 - LA FONCTION D'AUTOCORRELATION D'UN CODE DE GOLD AVEC N=10. .............................................. 44
FIGURE III-9 - LA FONCTION D'INTER CORRELATION DE DEUX SEQUENCES DE GOLD AVEC N=10.......................... 44
FIGURE III-10 - LA FONCTION D'AUTOCORRELATION DU CODE KASAMI POUR N=8. .................................................... 47
FIGURE III-11 - LA FONCTION D'INTER CORRELATION DE DEUX CODES DE KASAMI POUR N=8. .............................. 47
FIGURE III-12 LA FIC ENTRE LES 64 LIGNES DE LA MATRICE H (64X64). .................................................................... 50
FIGURE III-13 LA FAC LINEAIRE POUR INDEX 64 DU CODE WALSH HADAMARD H(64X64). ............................... 50
FIGURE III-14 - L'ARBRE GENERATEUR DE CODE OVSF POUR UN FACTEUR FE= 8. ..................................................... 52
FIGURE III-15 - LA GENERATION DES NOUVEAUX CODES OVSF......................................................................................... 53
FIGURE III-16 - LA FONCTION D'INTER-CORRELATION ENTRE LES CODES OVSF (64 BIT). ........................................ 53
FIGURE III-17 - LA FONCTION D'AUTO CORRELATION DU CODE GOLAY 1040. ............................................................... 57

iv
FIGURE III-18 - LA FONCTION D'INTER CORRELATION DE DEUX CODES GOLAY 1040. ................................................ 57
FIGURE III-19 - LA FONCTION D'AUTO CORRELATION POUR LE CODE DE BARKER 13. ................................................ 58
FIGURE III-20 - LA FAC POUR BARKER (13 13 7). ................................................................................................................... 59
FIGURE III-23 - LES FAC DES DIFFERENTES SEQUENCES POUR BARKER (L=105). ...................................................... 62
FIGURE III-24 - FAC BIDIMENSIONNELLE DU CODE BARKER (7 5 3). ................................................................................ 63
FIGURE III-25 - LA FONCTION D'AUTOCORRELATION D'UNE SEQUENCE DE CODE LEGENDER DE LA CLASSE 1 .. 65
FIGURE III-26 - LA FONCTION D'AUTOCORRELATION D'UNE SEQUENCE DE CODE WEIL (L=1013) ......................... 66
FIGURE III-27 - LA FONCTION D'INTER CORRELATION DE DEUX SEQUENCES DE CODE WEIL (1013). ................... 66
FIGURE III-28 - LES DIFFERENTES TAILLES GENERER PAR WEIL ET GOLD ENTRE [1000 ET 30 000]. ..................... 67
FIGURE III-29 - L'EFFET DE TRONQUER UNE SEQUENCE DE WEIL 1013 DE 7CHIPS........................................................ 68
FIGURE III-30 - L'EFFET D'AJOUTER 7CHIPS A UNE SEQUENCE DE WEIL 1013................................................................. 68
FIGURE III-31 - LES NIVEAUX DE LA FAC D'UNE SEQUENCE DE CODE ALEATOIRE OPTIMISE ET NON OPTIMISE. 69
FIGURE IV-1 - SCHEMA D'UN GENERATEUR DE CODE CA DU SYSTEM GPS. .................................................................. 76
FIGURE IV-2 - LA SEQUENCE BINAIRE GENEREE PAR LE SATELLITE 11. ............................................................................ 77
FIGURE IV-3 - LA FAC CORRESPOND A LA SEQUENCE GENEREE PAR LE SATELLITE 11. ............................................ 77
FIGURE IV-4 - LA DENSITE SPECTRALE DE PUISSANCE DU CODE CA (11). ...................................................................... 78
FIGURE IV-5 - SCHEMA D'UN GENERATEUR A 4 REGISTRES DE CODE P. ........................................................................... 79
FIGURE IV-6 - LES DENSITES SPECTRAL DE PUISSANCE DES CODES CA ET P DE GPS. ................................................. 80
FIGURE IV-7 - L'ENVELOPPE DE LA DENSITE SPECTRALE DE PUISSANCE POUR LES DEUX CODE CA ET P DE GPS.
......................................................................................................................................................................................................... 80
FIGURE IV-8 - PRINCIPE DE LA GENERATION DES TIERED-CODES. ..................................................................................... 83
FIGURE IV-9 - LES TAILLES DES CODES SECONDAIRE DANS LA CONFIGURATION GALILEO. ...................................... 83
FIGURE V-1 - L'ORGANIGRAMME QUI RESUME LES ETAPES SUIVIE POUR LA SIMULATION. ........................................ 91
FIGURE V-2 - LES NIVEAUX MAXIMUM COHERENTS DU LOBE LATERAL DES FAC POUR DIFFERENTS TYPE DE
CODES............................................................................................................................................................................................ 94
FIGURE V-3 - LES NIVEAUX MINIMUM COHERENTS DU LOBE LATERAL DES FAC POUR DIFFERENTS TYPE DE
CODES............................................................................................................................................................................................ 94
FIGURE V-4 - LES NIVEAUX MAXIMUM NON COHERENTS DU LOBE LATERAL DES FAC POUR DIFFERENTS TYPE
DE CODES. .................................................................................................................................................................................... 95
FIGURE V-5 - FIGURE 4 - LES NIVEAUX MINIMUM NON COHERENTS DU LOBE LATERAL DES FAC POUR
DIFFERENTS TYPE DE CODES. ................................................................................................................................................ 95
FIGURE V-6 - L'EFFET DE LA MODULATION SUR LE PIC PRINCIPALE DE LA FAC. ........................................................... 97
LES PARAMETRES DE LA SIMULATION SONT INDIQUES DANS LE TABLEAU SUIVANT. LES RESULTATS DE CETTE
SIMULATION SONT ILLUSTRES SUR LES FIGURES V-7 ET V-8. ................................................................................... 97
FIGURE V-8 - LA FAC DU CODE PN EN ABSENCE DU BRUIT. ................................................................................................. 97
FIGURE V-9 - L'EFFET DU BRUIT SUR LA QUALITE DE LA DETECTION D'UN CODE PSEUDO ALEATOIRE. ................. 98
FIGURE V-10 - L'EFFET DE L'ACCUMULATION SUR LA FONCTION D'AUTO CORRELATION. .......................................... 99

v
FIGURE V-11 LA VARIANCE DE LA DETECTION POUR UN CODE PN 1023 CHIP............................................................. 100
FIGURE V-12 LA STD DE LA DETECTION POUR UN CODE PN 1023 CHIP. ....................................................................... 101
FIGURE V-13 - L'RMSE DE LA DETECTION POUR UN CODE PN 1023 CHIP. .................................................................. 101
FIGURE V-14 - LA PROBABILITE DE LA DETECTION POUR UN CODE PN 1023 CHIP. ................................................... 102
FIGURE V-15 - LA VARIANCES DE LA DETECTION POUR LES CODES PN ET CA (1023 CHIP). .................................. 102
FIGURE V-16 LA PROBABILITE DERREUR DE LA DETECTION POUR LES CODES PN ET CA (1023 CHIP). ........ 103
NOUS REMARQUONS, D'APRES CETTE FIGURE, QUE LA PROBABILITE DERREUR POUR LES CODES PN APPROCHE
CELLE DU CODE CA-GPS. LA MEME OBSERVATION EST VALABLE POUR LE CAS DE LA VARIANCE (FIGURE
V-17). ......................................................................................................................................................................................... 103
FIGURE V-18 - LA VARIANCE DES DIFFERENTS TYPES DE CODES. ..................................................................................... 104
FIGURE V-19 - LA STD DES DIFFERENTS TYPES DE CODES. ............................................................................................... 104
FIGURE V-20 - LES STATISTIQUES DE LA DETECTION POUR LES DIFFERENTS TYPES DE CODES. ........................... 105
FIGURE V-21 - LA PROBABILITE D'ERREUR POUR LES DIFFERENTS TYPES DE CODES (1023 CHIP). ...................... 105
FIGURE V-22 - LA STD DE LA DETECTION POUR LES DEUX CODES CA ET P DE GPS. .............................................. 106
FIGURE V-23 LES RMSE DE LA DETECTION POUR LES DEUX CODES CA ET P DE GPS. ....................................... 107
FIGURE V-24 - LES PROBABILITE DERREUR DE LA DETECTION POUR LES DEUX CODES CA ET P DE GPS. ...... 107
FIGURE V-25 - LA PROBABILITE D'ERREUR POUR LES CODES DE WEIL ET LES CODES ALEATOIRES. ................... 108
FIGURE V-26 - LA PROBABILITE D'ERREUR POUR LES CODES DE WEIL, LES CODES ALEATOIRES ET LES
NOUVEAUX CODES GENERES. ............................................................................................................................................. 108
FIGURE V-27 - LE PRINCIPE DE LA DETECTION ADAPTATIF. ................................................................................................ 109
FIGURE V-28 - LES PROBABILITES D'ERREUR POUR UN CODE PN POUR DIFFERENT SEUILS DE DETECTION. ...... 110
FIGURE V-29 - L'ALGORITHME 1 POUR AMELIORER LES PERFORMANCES DES CODES PSEUDO ALEATOIRES. .... 110
FIGURE V-30 - L'ALGORITHME 2 POUR AMELIORER LES PERFORMANCES DES CODES PSEUDO ALEATOIRES. .... 111
FIGURE V-31 - L'AMELIORATION DES PERFORMANCES DUN CODE DE MEMOIRE (1023 CHIP)............................... 111
FIGURE V-32 - L'AMELIORATION DES PERFORMANCES DU CODE DE BARKER (1023 CHIP). .................................... 112
FIGURE V-33 - LES FONCTIONS D'INTER CORRELATION DES 4 SATELLITES GPS [1 2 3 4]........................................ 113
FIGURE V-34 - LES FONCTIONS D'INTER CORRELATION DES 8 SATELLITES GPS. ........................................................ 114
FIGURE V-35 - INTRA - SYSTEMES GALILEO E1 OS -> GPS L1C.................................................................................... 115
FIGURE V-36 - LES NIVEAUX DE LA CORRELATION (A) PAIRE ET (B) IMPAIRE DU SCENARIO 1 .............................. 115
FIGURE V-37 - INTRA - SYSTEMES GPS L1C -> GALILEO E1 OS.................................................................................... 116
FIGURE V-38 - LES NIVEAUX DE LA CORRELATION (A) PAIRE ET (B) IMPAIRE DU SCENARIO 2 .............................. 116
FIGURE V-39 - INTRA - SYSTEMES GPS L1C -> GPS L1C................................................................................................. 117
FIGURE V-40 - LES NIVEAUX DE LA CORRELATION (A) PAIRE ET (B) IMPAIRE DU SCENARIO 3.............................. 117
FIGURE V-41 - INTRA - SYSTEMES GALILEO E1 OS -> GALILEO E1 OS........................................................................ 118
FIGURE V-42 - LES NIVEAUX DE LA CORRELATION (A) PAIRE ET (B) IMPAIRE DU SCENARIO 4.............................. 118

vi

Liste des tableaux

TABLEAU I-1- BUDGET DERREUR AFFECTANT LES MESURES DE PSEUDO-DISTANCE DE CODE POUR LE CODE
C/A. ............................................................................................................................................................................................... 12
TABLEAU III-1 : LES POLYNOMES IRREDUCTIBLES CORRESPONDANT AUX DIFFERENTES TAILLES DU REGISTRE
GENERATEUR. ............................................................................................................................................................................. 35
TABLEAU III-2 - LA TAILLE ET LE NOMBRE DES SEQUENCES EN FONCTION DE LA TAILLE DU REGISTRE
GENERATEUR. ............................................................................................................................................................................. 40
TABLEAU III-3 - LA CONFIGURATION DE LA GENERATION DES DEUX SEQUENCES PN (N=127). ............................ 41
TABLEAU III-4 - LES MAXIMUMS DE FAC ET FIC POUR LES CODES DE PN (N=127). ................................................. 42
TABLEAU III-5 - LES MAXIMUMS DE L'INTER CORRELATION POUR LES CODES DE GOLD N=10................................ 45
TABLEAU III-6 - LES MAXIMUMS DE L'INTER CORRELATION POUR LES CODES DE KASAMI N=8. ............................ 47
TABLEAU III-7 - LES MAXIMUMS DE L'AUTOCORRELATION POUR LES CODES DE WALSH-HADAMARD H(64,64).
......................................................................................................................................................................................................... 51
TABLEAU III-8 - TABLE DE CONVERSION DU CODE BINAIRE EN CODE DE GRAY. .......................................................... 54
TABLEAU III-9 - LA CONFIGURATION POUR LA GENERATION DU CODE GOLAY 1040. ................................................. 56
TABLEAU III-10 - LES VALEURS MAX ET MIN DE LA CORRELATION POUR LE CODE DE GOLAY 1040.................... 57
TABLEAU III-11 - LISTE DES 9 CODES PRINCIPAUX DE BARKER .......................................................................................... 58
TABLEAU III-12 - LES NIVEAUX MAX ET MIN DE LA FAC POUR LE CODE 13 DE BARKER. ......................................... 58
TABLEAU III-13 - LES MAX ET LES MIN DE LA FAC POUR LES SEQUENCES DE BARKER 1183. ................................. 60
TABLEAU III-14 - LE MAXIMUM ET LE MINIMUM DE LA FAC POUR LES SEQUENCES BARKER 105. ....................... 61
TABLEAU III-15 - LES MAX ET MIN DE FAC ET DE FIC POUR LES CODES DE WEIL (1013). ....................................... 66
TABLEAU III-16 - LES VALEURS DE LA DEGRADATION DES NIVEAUX MAX ET MIN DE LA FAC ET DE LA FIC. .. 68
LES CODES C/A SONT DES CODES PSEUDO ALEATOIRE PERIODIQUE DE PERIODE 1MS AVEC UN DEBIT DE 1.023 M
CHIP/S, GENERER PAR UN GENERATEUR DE GOLD A DEUX REGISTRES LFSR DE 10 ETAGE SELON LA
CONFIGURATION INDIQUE DANS LE TABLEAU IV-1. ..................................................................................................... 76
TABLEAU IV-2 - LA CONFIGURATION D'UN GENERATEUR DE CODE CA DE GPS........................................................... 76
TABLEAU IV-3 - LES CARACTERISTIQUES DU CODE GPS CA (11). ..................................................................................... 78
TABLEAU IV-4 - LA CONFIGURATION DU GENERATEUR DE CODE P DU GPS. .................................................................. 78
TABLEAU IV-5 - LES CARACTERISTIQUES DE BASE DU CODE P DE GPS............................................................................ 80
TABLEAU IV-6 - - LES CARACTERISTIQUES DES DIFFERENTS SIGNAUX GALILEO. ......................................................... 81
TABLEAU IV-7 - LA CONFIGURATION DU GENERATEUR E1C GALILEO. ............................................................................ 81
TABLEAU IV-8 - LA CONFIGURATION DU GENERATEUR E1 B GALILEO ............................................................................ 82
TABLEAU IV-9 - LA CONFIGURATION DU GENERATEUR E5A I DE GALILEO .................................................................... 82
TABLEAU IV-10 - LA CONFIGURATION DU GENERATEUR E5A Q DE GALILEO ................................................................ 82
TABLEAU IV-11 - LA CONFIGURATION DU GENERATEUR E5B I DE GALILEO .................................................................. 82
TABLEAU IV-12 - LA CONFIGURATION DU GENERATEUR E5B Q DE GALILEO ................................................................ 82

vii
TABLEAU IV-13 - LES CODES E1 OS DE GALILEO ET L1C DE GPS. .................................................................................. 84
TABLEAU IV-14 - LES NOUVEAUX CODES CHOISIS POUR LES NOUVEAUX SERVICES DE GNSS. ............................... 87
TABLEAU V-1 - L'ENSEMBLE DES TAILLES DE REFERENCES. ................................................................................................. 93
TABLEAU V-2 - LA CONFIGURATION DE LA GENERATION DE DIFFERENTES TAILLES DU CODE PN. ........................ 93
TABLEAU V-3 - LES NIVEAUX MAX ET MIN DU LOB LATERALE FAC POUR DIFFERENT TYPE DE MODULATIONS.
......................................................................................................................................................................................................... 96
TABLEAU V-4- - LES PARAMETRE DE CONFIGURATION POUR ETUDIER L'EFFET DU BRUIT SUR LA DETECTION .. 97
IL EST EVIDENT QUE LA DETECTION DU PIC PRINCIPAL DANS CES CONDITION (SNR=-30 DB) EST IMPOSSIBLE,
LES NIVEAUX MAX ET MIN DES FONCTIONS DE CORRELATION SONT RESUMES DANS LE TABLEAU V-5. .... 98
TABLEAU V-6 - LES NIVEAUX DE LA CORRELATION D'UN CODE PSEUDO ALEATOIRE EN PRESENCE ET EN
ABSENCE U BRUIT...................................................................................................................................................................... 98
TABLEAU V-7 - LA CONFIGURATION POUR ETUDIER L'EFFET DE L'ACCUMULATION. .................................................... 98
TABLEAU V-8 - LES RESULTATS DE SIMULATION POUR UN NOMBRE DE PERIODE = 10 ................................................ 99
TABLEAU V-9 - LES PARAMETRES DE SIMULATION POUR L'ETUDE DE PERFORMANCES DE LA DETECTION. ...... 100
TABLEAU V-10 - LES PARAMETRES DE SIMULATION DE L'EFFET DU MODE DE LA DETECTION. ............................ 109
TABLEAU V-11 - LES NIVEAUX D'INTER-CORRELATION POUR 4 STELLITES DE GPS. ................................................ 113
TABLEAU V-12 - LES NIVEAUX D'INTER-CORRELATION POUR 8 STELLITES DE GPS. ................................................ 114
TABLEAU V-13 LES RESULTATS DU SCENARIO 1 GALILEO INTERFERE GPS. ............................................................... 116
TABLEAU V-14 LES RESULTATS DU SCENARIO 2 GPS INTERFERE GALILEO. ............................................................... 116
TABLEAU V-15 LES RESULTATS DU SCENARIO 3 GPS INTERFERE GPS......................................................................... 117
TABLEAU V-16 LES RESULTATS DU SCENARIO 4 GALILEO INTERFERE GLILEO. ......................................................... 118

Conclusion gnrale
Les systmes GNSS sont les systmes de localisation les plus utiliss qui sont capables
datteindre une meilleure prcision surtout lorsquon lassocie avec des codes de
meilleurs performances.
Pour plusieurs annes, il y avait des efforts de recherche significatifs pour ltude et
lamlioration des performances des codes dtalement. Nous avons constat que
plusieurs paramtres peuvent affects dune manire directe ou indirecte, parmi ces
paramtres on peut citer : le type de code, la taille du code et aussi les mthodes
dacquisitions. Le but de ce travail est ltude comparative des codes de ltalement et
lanalyse des performances de ces codes dans les systmes de navigation GNSS.
Dans le cadre de ce mmoire, nous prsentons dans un premier chapitre quelques bases
concernant les transmissions numriques et la chaine de communication pour ce type de
communication et on a introduit les rcepteurs optimaux. Ensuite, on prsente les
notions de la probabilit.
Dans le deuxime chapitre, on a prsent d'un point de vue gnral, les systmes de
navigation par satellites en termes de signaux et principes de navigation pour cela, on a
pris l'exemple d'un rcepteur GPS.
Le troisime chapitre prsente une tude thorique de ltalement de spectre de la
technique CDMA employ dans les systmes de navigation. Ce chapitre prsente aussi
l'historique sur lorigine de ltalement de spectre, ainsi l'importance de ce dernier pour
le codage du canal avec des codes bien dtermins.
On a termin par un dernier chapitre, o on a prsent les rsultats de simulation, les
conditions sur les caractristiques des codes et les critres de la comparaison qui seront
suivis.
Lobjectif principal de notre travail est d'tudier les performances de codes, en termes
de dtection, temps dacquisition et de la synchronisation des codes qui est une tape
trs cruciale dans les systmes de navigation par satellites.
Pour atteindre cet objectif, on a suivi un modle simple et logique. Ltude thorique des
algorithmes de gnration des codes nous a permis de crer une base de donnes dun
grand nombre de codes, pour classer et tester les performances des codes selon les
critres des fonctions dautocorrlation et de linter-corrlation ainsi que les statistiques
des squences binaires.

En particulier, on a tudi leffet de la taille des codes ce qui nous a permis dobtenir
des meilleurs performances, en terme de dtection et dimmunit contre le bruit, trs
significatives pour des tailles de codes.
Nous avons calcul et trac les fonctions dautocorrlation et de linter-corrlation
des signaux GPS et Galileo afin dillustrer linfluence de la taille des codes sur les
niveaux de linter corrlation de code et du lobe latral de la fonction
dautocorrlation. Tous les rsultats sont exploits afin dvoluer les performances
des systmes GNSS.
Les rsultats de simulation montrent que les codes alatoires optimiss sont les
meilleurs codes pour les applications GNSS, mais les codes de Weil sont les plus
efficaces, ils peuvent couvrir le grand nombre dapplications GNSS.
Par consquence, des tudes et des recherches pour trouver des moyens pour
optimiser les codes de Weil.
Les proprits exceptionnelles des diffrents codes dtalement, ont suscit un grand
intrt mais aussi des questions, la principale est ce que on peut combiner diffrents
types de codes pour amliorer leurs performances.

PERSPECTIVES
Comme continuit de ce travail, il sera intress dimplmenter ces codes ainsi que
les algorithmes doptimisation des codes dans des systmes de traitement
numriques en temps rel tel que les systmes de DSP ou bien FPGA. Qui peuvent
tre par la suite intgr dans un systme dacquisition et de poursuite des codes dans
les systmes de navigation.
Rsum

es squences pseudo alatoires occupent une place importante dans les
systmes de tlcommunication modernes.
Le travail prsent dans ce mmoire porte sur les caractristiques communes des
codes pseudo alatoires pour les applications GNSS nouvelles gnrations (de GPS
moderniss ainsi que Galileo). La poursuite de tels codes est base sur la dtection de
dcalages relatifs entre le rcepteur et les satellites. Les caractristiques des
nouveaux signaux ont oblig les concepteurs de rcepteurs de dvelopper de
nouvelles techniques dacquisition rapide ou de modifier le type de code { utiliser. Le
dveloppement de nouvelles techniques dacquisition et de poursuite ncessitent un
changement matriel ce qui peut compliquer larchitecture du rcepteur. La
modification des codes utiliss savre une solution meilleure. Lobjectif de ce travail
est damliorer les performances dacquisition des signaux pour le but de faire une
tude comparative entre les diffrents codes dtalements afin de dduire les
meilleurs codes pour les applications GPS/Galileo selon des critres fixs par les
besoins des applications GNSS (au terme de fonction dautocorrlation, fonction
dinter-corrlation , bande passant et rapidit dacquisition). Pour cela, on a suivi cinq
tapes indispensables, premirement une tude sur les systmes de navigation GNSS,
les principes de la navigation, les techniques dacquisitions des codes par rcepteurs
GNSS et les sources derreurs de la localisation. Dans un deuxime chapitre, les
techniques daccs et spcialement la CDMA ainsi le principe de ltalement de
spectre. Le troisime chapitre est une tude dtaille des arithmtiques de la
gnration des diffrents types de codes et leurs caractristiques en termes de la
corrlation. Le quatrime chapitre dcrit les caractristiques des codes choisis par les
systmes GNSS et les critres de base pour choisis ces codes.
Et enfin, nous terminons ce manuscrit par un chapitre cinq, consacr aux rsultats de
simulation, un classement de ces codes selon leurs performances concernant les
fonctions dautocorrlation et dinter-corrlation par la suite une tude de leurs
caractristiques (probabilit et statistique) tel que la variance, STD, RMSE et la
probabilit derreur pour chaque type de code ensuite une tape trs importante qui
est lamlioration des codes, o on a tudi lquilibre des squences. De nouveaux
algorithmes pour amliorer les performances des codes. Ltape finale est une tude
comparative entre les deux systmes de navigation GPS et Galileo moderniss en
termes dacquisition des nouveaux codes alatoire et de Weil.
L

Amine GUISOUS Mémoire Fin D'étude

Transféré par

Droits d'auteur :

Formats disponibles

Amine GUISOUS Mémoire Fin D'étude

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Amine GUISOUS Mémoire Fin D'étude

Transféré par

Droits d'auteur :

Formats disponibles

REPUBLIQUE ALGERIENNE DEMOCRATIQUE ET POPULAIRE

MINISTERE DE LENSEIGNEMENT SUPERIEUR ET DE LA RECHERCHE

). Lintgrateur permet de corrler le code dmodul ( ) avec

partir des chantillons () chaque

est atteinte quand un bit erron conduit au symbole

. Il conviendra donc de transmettre n bits pour chantillon (si la

impulsion par seconde, si

est donc, daprs la relation de

frquence maximum prsente dans le spectre du signal analogique (au sens

trs faible et que la puissance du

par une squence pseudo-

de la squence PN et la priode de chip

la porteuse et la phase initiale du signal. La relation entre la

nous pouvons alors

, cest le retard relatif entre lmetteur et le rcepteur

. La langueur de la squence est fixe par le choix des

couples (() ( )) tel que les couples

. Sil ny a pas de changement de signe de entre

intervalles , ( )-, il y a un nombre pair de

) avec une taille de , o

. La squence de la sortie satisfait la relation de

existent pour tous les nombres

. Ce type de codes est

() est un polynme primitif de degr n, et

() un polynme primitif bien

squences de priode avec une corrlation

pour obtenir le polynme ()

. Un registre dcalage linaire de 12

) et donc chaque squence y de priode

indiquent les deux

codes de Walsh-Hadamard acquis par la matrice de

), seulement codes sont orthogonaux pour toutes

) dans les lments

). D'ailleurs, la Figure III-12 montre les FICs de tous les codes de

est la nouvelle taille de la squence fille. L'arbre des

squences de codes dtalement avec la longueur de

une matrice de la dimension

() . Alors, la matrice gnratrice

peut tre produite partir

selon la relation ci-dessus :

est le complment binaire, par consquent, deux codes

. Pour comparer ces trois rsultats, on fait

, et sur des registres dcalage linaires - Linear Feedback Shift

qui est le problme des codes de Gold.

a une solution . On peut dterminer si est un rsidu

La fonction d'autocorrlation d'une squence de Legendre de Classe 1

), etc. La figure suivante reprsente une

(n : est la longueur du code et M le nombre des codes) qui est

. Il permet de multiplier chaque priode (poche) du code primaire par un bit de 1

, ou bien aprs normalisation de la fonction

. Des expressions plus prcises sur cette limite ont t obtenues

La frquence centrale de Doppler.

Vous aimerez peut-être aussi