Si, au cours de l'épreuve, un candidat repère ce qui lui semble être une erreur d'énoncé, il le signale sur sa copie et poursuit sa composition en expliquant les raisons des initiatives qu'il est amené à prendre.
Etude du fonctionnement du GPS : des satellites au récepteur
Figure 1 : illustration de la constellation des satellites GPS.
Étude du fonctionnement du GPS : des satellites au récepteur
1 Présentation de l'étude
1.1 Introduction au GPS
Le GPS (Global Positioning System) est un système de localisation par satellites conçu par les Américains dans les années 1980 pour des applications militaires. Les satellites GPS tournent en orbite autour de la Terre, comme illustré sur la figure 1, et émettent en permanence des signaux vers celle-ci, signaux qui sont reçus et décodés par le récepteur que porte l'utilisateur. Ce récepteur GPS fournit alors à ce dernier une solution dite «PVT»:
P, pour Position : position de l'antenne GPS du récepteur (trois coordonnées);
V, pour Vitesse : vitesse de l'antenne GPS du récepteur (trois coordonnées);
T, pour Temps : référence de temps du récepteur, synchronisée avec les horloges atomiques extrêmement précises des satellites GPS.
Figure 1 - Illustration de la constellation des satellites GPS
Depuis sa mise en service, les domaines d'application du GPS n'ont cessé de se développer et de s'étendre dans les applications militaires et surtout dans les applications civiles : transports terrestres, aériens et maritimes, opérations de secours, agriculture, géologie, prospection pétrolière, communications, et bien évidemment en tant que système intégré au téléphone mobile dans la vie de tous les jours.
L'Europe n'est pas en reste et développe son propre système de localisation par satellite : Galileo. Ce système permettra tout d'abord à l'Europe d'assurer son indépendance dans ce domaine stratégique (le GPS étant un système américain) et également de bénéficier de performances et de fonctionnalités supérieures à celles offertes actuellement par le GPS.
Cinq caractéristiques essentielles du GPS (et bientôt également de Galileo) ont permis son immense succès dans de multiples domaines :
la solution PVT fournie par le récepteur est disponible n'importe où dans le monde;
le récepteur GPS fournit une solution PVT suffisamment précise pour de nombreuses applications; en particulier l'erreur est de l'ordre de la dizaine de mètres pour la position et de quelques centimètres par seconde pour la vitesse;
le récepteur GPS fournit une position et une vitesse en trois dimensions;
le récepteur GPS fournit une solution PVT en continu et avec un taux de rafraîchissement suffisamment rapide, de l'ordre de quelques Hertz;
le nombre d'utilisateurs est illimité.
Q1. Donner 2 exemples concrets d'utilisation des informations fournies par un récepteur GPS.
1.2 Objectif de l'étude
L'objectif de l'étude est de comprendre comment le GPS fonctionne, c'est-à-dire de comprendre comment une position, une vitesse et un temps peuvent être estimés par un récepteur à partir de signaux émis par des satellites en orbite autour de la Terre.
Pour ce faire, le fonctionnement du GPS est modélisé comme une succession de 4 étapes, telle qu'illustrée sur la figure 2:
Étape 1 : génération et émission des signaux GPS par les satellites en orbite;
Étape 2 : propagation de ces signaux des satellites vers l'antenne du récepteur GPS;
Étape 3 : réception et décodage de ces signaux par le récepteur;
Étape 4 : estimation de la solution PVT par le récepteur.
Figure 2 - Vue d'ensemble du fonctionnement du GPS
Outre la partie 1 de présentation de l'étude, le sujet se décompose en 5 parties:
Partie 2 : détermination de la solution PVT à travers une approche simplifiée;
Partie 3 : modélisation de la trajectoire des satellites et analyse de la génération des signaux;
Partie 4 : étude de la propagation des signaux GPS du satellite vers l'antenne du récepteur;
Partie 5 : étude du décodage des signaux par le récepteur et mise en relation avec l'étude menée dans la partie 2;
Partie 6 : synthèse de l'étude.
Même si certaines parties peuvent être traitées séparément, il est préférable de suivre l'ordre de questionnement proposé.
2 Détermination de la solution PVT
On considère satellites , en orbite autour de la Terre qui émettent en permanence un signal se propageant à la vitesse de la lumière. Ces satellites possèdent des horloges atomiques extrêmement précises et synchronisées entre elles, permettant de définir ainsi une horloge universelle dénommée «temps GPS »>. Tous les signaux sont émis à un même instant connu. Chaque satellite inclut dans son signal les informations donnant son identité (numéro de satellite) ainsi que les paramètres de son orbite nécessaires au récepteur pour que ce dernier puisse estimer précisément et en temps réel la position ( ) du satellite dans le référentiel géocentrique .
Ce référentiel, défini sur la figure 3, a pour origine le centre de la Terre et les axes de la base ( ) sont définis à partir d'étoiles lointaines. La position d'un satellite de la constellation est définie par ses coordonnées sphériques ou par ses coordonnées cartésiennes .
On considère également un récepteur GPS à une position de coordonnées inconnues ( ) que l'on cherche à déterminer. Ce récepteur capte les signaux des satellites . D'après le paragraphe précédent, grâce au décodage des signaux, le récepteur est capable d'identifier l'identité du satellite générateur de chaque signal et d'en déduire la position de chaque satellite . On définit la distance entre le récepteur et le satellite , comme illustré sur la figure 4.
On prend tout d'abord comme hypothèse que l'horloge du récepteur est synchronisée avec le temps GPS des satellites.
Q2. Démontrer que le récepteur est capable de déterminer la distance . Donner l'expression de en fonction des coordonnées de et de .
Q3. Donner les deux conditions nécessaires sur les satellites afin que le récepteur puisse déterminer les coordonnées du point .
En réalité, l'horloge du récepteur est beaucoup moins précise que celles des satellites. L'horloge du récepteur n'est donc pas synchronisée avec le temps GPS. Ce biais entre le référentiel de temps du récepteur et celui du GPS est inconnu ( ). De ce fait, le récepteur ne mesure plus directement la distance mais une «pseudo-distance» . On suppose que le temps du récepteur est en retard sur le temps GPS.
Q4. Donner l'expression de en fonction de et .
Q5. Donner les deux conditions nécessaires sur les satellites afin que le récepteur puisse déterminer les coordonnées du point .
Pour la suite de cette partie, on considère que les satellites visibles par le récepteur remplissent les deux conditions précédemment déterminées. Les questions suivantes détaillent l'algorithme permettant au récepteur de déterminer sa position à partir des mesures de pseudo-distances . Cet algorithme doit donc permettre de résoudre un système de équations à 4 inconnues : . L'algorithme considéré est celui des moindres carrés, dans lequel la minimisation de la somme des erreurs au carré est réalisée par itérations successives.
On désigne par le vecteur représentant les pseudo-distances le vecteur représentant les inconnues à déterminer, le vecteur représentant les positions du satellite calculées par le récepteur grâce au décodage des informations transmises par le signal, et la transposée d'une matrice quelconque.
Le système considéré possède donc équations à 4 inconnues. Afin de le résoudre, la fonction suivante, représentant la somme des erreurs au carré, est minimisée:
La fonction représente le calcul de la pseudo-distance à partir du vecteur d'inconnues et des données . Comme cette fonction est non linéaire, la minimisation de conduit au système suivant, linéarisé autour des variations du vecteur :
Dans cette équation, est la matrice jacobienne donnée par et est l'évolution des valeurs du vecteur entre la et la itération. Le problème revient donc à résoudre, par la méthode de Newton-Raphson, le schéma numérique suivant :
Afin de résoudre le schéma numérique, la méthode du pivot de Gauss avec recherche partielle du pivot est utilisée pour la résolution linéaire du problème local.
Q6. Le problème local est équivalent au système linéaire . Expliciter la matrice , le vecteur et le vecteur du problème en fonction de et de .
Q7. Donner la relation que doit vérifier la matrice pour que le système admette une solution unique.
La première étape de la résolution consiste à trouver la forme échelonnée réduite par lignes de la matrice par la méthode du pivot de Gauss avec recherche partielle du pivot. À la fin de cette étape, un système linéaire équivalent au système linéaire est obtenu.
Nota :
Une matrice est dite échelonnée par lignes si elle vérifie les deux propriétés suivantes:
si une ligne est nulle, toutes les lignes suivantes le sont aussi ;
à partir de la deuxième ligne, dans chaque ligne non nulle, le premier coefficient non nul à partir de la gauche est situé à droite du premier coefficient non nul de la ligne précédente;
le «pivot» est le premier coefficient non nul de chaque ligne;
une matrice échelonnée par lignes est dite échelonnée réduite par lignes si tous ses pivots sont égaux à 1 et sont les seuls éléments non nuls de leur colonne.
Q8. Expliquer comment obtenir la forme échelonnée réduite par lignes en précisant les opérations à effectuer sur les coefficients de la matrice ou sur ses lignes et en énonçant le principe de la recherche partielle du pivot. Déterminer également les relations afin de garder l'équivalence entre le système linéaire et le système linéaire .
Q9. Proposer une fonction écrite dans le langage Python permettant de résoudre le système linéaire où la matrice est une matrice échelonnée réduite en ligne. Cette fonction devra prendre comme arguments d'entrée la matrice et le vecteur associé . Elle devra alors retourner la valeur solution du système . Les coefficients de la matrice sont stockés dans une liste de liste telle que est le terme de la matrice .
On suppose dans la suite du sujet que le schéma de Newton-Raphson converge. Nous nous intéressons au coût de la résolution de ce schéma numérique. Pour cela, les hypothèses suivantes sont prises :
l'affectation d'une variable vaut une unité de temps, c'est-à-dire 1 UT ;
les opérations sur les nombres et les booléens valent chacune 1 UT;
l'accès à la valeur d'une variable pour les nombres et les booléens vaut chacun UT;
les comparaisons, les tests d'égalité ou de différence valent chacun 1 UT;
l'accès à l'élément d'une liste, l'accès à la longueur d'une liste, l'insertion en fin de liste valent chacun 5 UT ;
la suppression d'un élément d'une liste de taille , l'insertion en début de liste de taille , ainsi que le parcours de la liste de taille sont en temps linéaire et valent chacun 5.n UT ;
la concaténation de deux listes de taille et est en temps linéaire et coûte UT.
Q10. À partir des coûts d'opérations proposées, déterminer la complexité temporelle au pire des cas de la fonction proposée à la question précédente en fonction de , où est la dimension de la matrice .
L'unité de temps est d'environ 10 ns . Le coût temporel du pivot de Gauss est d'environ .
Q11. Conclure quant au coût de la résolution du système par la méthode numérique proposée.
Après la résolution du pivot de Gauss, il est possible d'écrire la suite du schéma numérique de façon plus explicite :
Q12. Proposer un critère d'arrêt pour le schéma numérique de la suite .
Il a été démontré dans cette partie que la détermination du vecteur , c'est-à-dire des solutions P et T , est possible grâce à la réception et au décodage des signaux émis par les satellites.
L'estimation de la vitesse du récepteur se fait de façon similaire, grâce à la mesure de «pseudovitesses » (car le biais de temps subit en réalité quelques variations temporelles) et la connaissance de la vitesse des satellites via les informations transmises dans les signaux. Ainsi, à partir des signaux GPS émis par les satellites, le récepteur est capable de fournir la solution PVT à l'utilisateur.
3 Modélisation du segment spatial
3.1 Etude de la constellation des satellites GPS
Afin d'assurer une bonne couverture du réseau satellite GPS, la constellation des satelittes GPS doit garantir un nombre minimum de 21 unités opérationnelles dans du temps. Ces satellites sont répartis dans 6 plans orbitaux décalés de entre eux avec un décalage de par rapport au plan de l'équateur. Actuellement, la constellation comporte 31 unités opérationnelles ce qui permet pour l'utilisateur d'avoir un nombre moyen de 8 satellites visibles (ce nombre de satellites visibles est compris entre 4 et 12).
Afin d'obtenir une certaine stabilité de la couverture du réseau, un satellite de la constellation doit retrouver une position identique, dans le référentiel géocentrique (voir figure 3), deux fois par jour. Dans ce cas, son orbite terrestre correspond à une période deux fois plus petite qu'un jour sidéral ( ). Un point quelconque de l'espace est représenté par le vecteur position en coordonnées sphériques:
Le satellite de la constellation est représenté par le point situé à une distance de l'origine :
3.1.1 Caractéristiques de l'orbite idéale d'un satellite de la constellation GPS
Dans cette première approche, nous allons considérer uniquement la Terre et un satellite de la constellation GPS: l'influence des autres astres est négligée. Ce satellite, vu comme une masse ponctuelle, n'est donc soumis qu'au champ de gravitation de la Terre. Le mouvement du satellite est décrit dans le référentiel géocentrique. À l'instant initial, le satellite est à une altitude du sol terrestre et possède une vitesse initiale . Notre planète est considérée comme une boule solide de masse et de rayon . Son champ de gravitation est donné en tout point de l'espace par l'équation suivante :
avec le potentiel de gravitation au point .
Q13. Montrer que le mouvement du satellite par rapport au référentiel géocentrique est plan.
Les coordonnées polaires ( ) sont utilisées pour définir la position du point dans un plan de normale , où l'axe ( ) est fixe par rapport au référentiel géocentrique . L'angle est défini par rapport à un axe ( ) de référence tel que :
Figure 5 - Coordonnées polaires
Q14. Déterminer la norme du moment cinétique du satellite en fonction de et . Montrer également que l'on a la relation suivante:
Q15. À partir du principe fondamental de la dynamique, montrer que le mouvement du satellite par rapport au référentiel géocentrique décrit la trajectoire d'une courbe dont l'équation dans le plan polaire vérifie l'équation différentielle:
où est une constante du problème à déterminer en fonction de et .
Q16. Préciser pour quelle(s) condition(s) le satellite ne se libère pas de l'attraction de la Terre.
Q17. En supposant que le satellite tourne en orbite circulaire autour de la Terre, avec une période 2 fois plus petite qu'un jour sidéral, déterminer l'altitude de l'orbite du satellite en fonction de et . Donner également l'expression de sa vitesse de satellisation en fonction de et .
3.1.2 Évaluation des perturbations de positionnement du satellite dues au couple Soleil-Lune
Afin d'obtenir une modélisation plus précise de la trajectoire des satellites, l'influence de la Lune, de masse , et du Soleil, de masse , sur leur mouvement est désormais prise en compte. Pour cela, nous considérons le système formé par les trois astres: Terre, Soleil et Lune, dont les centres d'inertie sont respectivement et . La Terre tourne autour de l'axe ( ) où le point est le centre d'inertie du système . La Terre a également un mouvement de rotation sur elle-même autour d'un axe ( ) incliné d'un angle de par rapport ( ). Nous considérons également le mouvement de rotation de la Lune autour de l'axe ( ) incliné d'un angle de par rapport à l'axe ( ).
Figure 6 - Système Soleil-Terre-Lune
Les hypothèses simplificatrices suivantes sont prises:
seules les attractions du Soleil, de la Terre et de la Lune sont prises en compte;
toutes les orbites sont circulaires et on ne considère que les distances moyennes entre les éléments célestes ;
le référentiel galiléen est le référentiel lié au centre d'inertie des trois astres considérés et dont les directions restent identiques à celles du référentiel géocentrique .
Les notations suivantes sont utilisées:
Q18. Déterminer les caractéristiques du mouvement du référentiel par rapport au référentiel galiléen . Déterminer alors la vitesse du point par rapport au référentiel en fonction de la vitesse du point par rapport au référentiel et de la vitesse du point par rapport au référentiel .
Q19. Déterminer l'accélération du centre de la terre par rapport au référentiel galiléen .
Q20. Déterminer l'accélération du satellite par rapport au référentiel galiléen en fonction de l'accélération du satellite par rapport à et des constantes et .
Q21. Faire l'inventaire des actions mécaniques appliquées au satellite. Expliciter chacune de ces actions mécaniques en fonction des données du problème.
L'application au satellite du principe fondamental de la dynamique fait apparaître de nouveaux termes par rapport au modèle idéal étudié dans le §3.1.1. Ces termes peuvent être considérés comme des forces perturbant l'orbite du satellite.
Q22. À partir du principe fondamental de la dynamique appliqué au satellite, déterminer et les perturbations dues à la présence respective du Soleil et de la Lune.
Q23. Évaluer la norme maximale de et de . Les résultats seront exprimés en fonction de , et de .
3.1.3 Évaluation des perturbations de positionnement du satellite dues à la non-sphéricité de la Terre
Dans l'analyse qui suit, la Terre n'est plus considérée comme étant de forme sphérique mais comme une ellipsoïde avec une forme plus aplatie au niveau des pôles.
Figure 7 - Géoïde de référence
Les hypothèses simplificatrices suivantes sont prises:
le potentiel de gravitation de la Terre s'exprime par , le développement zonal d'ordre 2 du potentiel précédent :
le satellite est suffisamment haut pour considérer que est toujours porté par .
Pour rappel, le gradient en coordonnées sphériques s'écrit comme suit :
Q24. Déterminer la nouvelle force de gravitation de la Terre issue du potentiel . Déterminer également , l'écart de force par rapport au modèle simplifié du §3.1.1.
Q25. Évaluer alors la norme maximale de . Le résultat sera exprimé en fonction de et .
Les applications numériques des résultats obtenus au §3.1.2 et au §3.1.3 donnent
une perturbation due au couple luni-solaire de l'ordre de par rapport à l'attraction terrestre;
une perturbation due à la non-sphéricité du champ de gravitation terrestre de l'ordre de par rapport au potentiel de gravitation d'une Terre hypothétiquement ronde.
Ces perturbations mécaniques, bien que numériquement faibles en valeurs relatives, ont une influence non négligeable sur l'orbite des satellites. En effet, par intégration temporelle, si elles ne sont pas prises en compte dans les calculs de trajectographie, elles générent des erreurs d'estimation de position des satellites trop fortes pour obtenir un système GPS performant. Par conséquent, les signaux envoyés par les satellites contiennent tous les paramètres et correctifs correspondants nécessaires pour que le récepteur GPS puisse précisément évaluer en temps réel la position des satellites.
3.2 Génération des signaux GPS
3.2.1 Expression et représentation des signaux GPS
On considère satellites , en orbite autour de la Terre. Les horloges embarquées à bord de ces satellites sont des horloges atomiques extrêmement précises et synchronisées entre elles, définissant ainsi le «temps GPS». Chaque satellite de la constellation émet en permanence un signal vers la Terre. Ces
signaux sont tous émis au même instant, mais leur contenu est propre à chaque satellite puisque s'y trouvent les informations sur sa trajectoire, informations essentielles pour que le récepteur puisse ensuite déterminer la position du satellite. Ces informations sont appelées « message de navigation».
Dans la suite du sujet sont utilisées les deux formulations suivantes:
signal «binaire»: succession de 0 et de 1 ;
signal «binaire NRZ» (Non Retour à Zéro) : succession de -1 et de 1 ; ce signal étant créé à partir d'un signal binaire, les valeurs -1 et 1 correspondent respectivement aux 0 et 1 du signal binaire.
Le signal émis par le satellite , noté , est la combinaison de plusieurs signaux élémentaires :
le message de navigation , signal binaire NRZ propre à chaque satellite, de fréquence ;
le code pseudo-aléatoire dit «code » (Coarse/Acquistion ou acquisition grossière en français) , code binaire NRZ propre à chaque satellite, de fréquence ;
la porteuse dite «porteuse L1» , signal sinusoïdal identique pour tous les satellites, de fréquence et d'amplitude .
La génération de est représentée sur la figure 8 . Du fait de l'utilisation de la représentation binaire NRZ, l'opération notée sur la figure 8 est une multiplication classique de signaux : si les 2 entrées sont et , la sortie de cet opérateur est avec .
Figure 8 - Génération du signal
Nota : compte tenu de la représentation binaire NRZ, la porteuse subit des sauts de phase de en fonction de la valeur du signal utile ( -1 ou 1). La modulation représentée sur la figure 8 est donc de type BPSK (Binary Phase-Shift Keying).
Q26. Donner l'expression du signal de sortie en fonction de .
Un exemple de signaux et est donné sur la figure A du document réponse, ainsi que le début du signal correspondant.
Nota : pour permettre une bonne visualisation graphique, les rapports de fréquences des différents signaux n'ont pas été conservés sur cette figure.
Q27. Compléter la forme des signaux et .
3.2.2 Etude de la génération du code
Tous les signaux émis par les satellites ont la même porteuse de fréquence . Afin que le récepteur puisse distinguer les différents signaux qu'il reçoit, un code binaire NRZ est utilisé, comme défini précédemment. Ce code possède les propriétés suivantes:
P1 : chaque code est propre au satellite ;
P2 : tous les codes ont la même fréquence et ont le même nombre de bits NRZ avant de se répéter;
P3 : la séquence de bits NRZ des codes n'est pas confidentielle: la structure de ces codes est connue du grand public;
P4 : les codes sont dits «orthogonaux» (codes de Gold), c'est-à-dire que leur intercorrélation est quasi-nulle, où l'intercorrélation entre deux codes et est définie par l'équation suivante: pour tout entiers naturels non nuls et tout entier relatif ;
P5 : l'autocorrélation d'un code est quasi-nulle pour un retard d'au moins 1 bit NRZ, où l'autocorrélation d'un code est définie par l'équation suivante:
Nota : comme il le sera démontré dans la partie 5, la propriété P4 permet au récepteur de séparer les différents signaux tandis que la propriété P5 permet au récepteur de déterminer la pseudo-distance séparant l'antenne du récepteur et le satellite émettant le signal décodé.
La puissance d'un code , notée , est définie par . Numériquement, on considère dans cette étude que les propriétés P4 et P5 s'écrivent : pour tout entiers naturels non nuls et tout entier relatif; pour tout entier naturel non nul et tout entier relatif.
Q28. Déterminer la valeur numérique de . Cette valeur dépend-elle du code considéré?
Chaque satellite possède à son bord un générateur de code qui permet de créer un code possédant les caractéristiques précédentes. La figure 9 illustre la génération du code utilisé par le satellite , les autres codes étant générés de façon similaire. Pour ce faire, deux registres et sont utilisés, chaque registre étant un emplacement de mémoire interne à un processeur. Chacun de ces registres est constitué de 10 bits, ceux-ci étant notés et , pour entier naturel compris entre 1 et 10 .
Figure 9 - Générateur du code
À l'initialisation du processus de génération du code , les registres et sont composés uniquement de bits valant 1 . Ensuite, au pas de calcul, les opérations suivantes sont réalisées successivement afin de déterminer les bits des registres et , opérations à l'issue desquelles le nouveau bit du code binaire est calculé :
registre
registre
sortie
Enfin, le code binaire est transformé en code binaire NRZ par simple passage du bit 0 à -1 .
L'opération notée dans la figure 9 est un sommateur bit à bit modulo 2 , c'est-à-dire un «OU exclusif » dont les règles sont énoncées dans le tableau 1 .
Table 1- Règle pour l'opérateur
Entrée E1
Entrée E2
Sortie S=E1 E2
0
0
0
0
1
1
1
0
1
1
1
0
Nota: les autres codes binaires pour les autres satellites sont générés en sélectionnant des bits du registre différents des bits et choisis ici.
Sur la figure B du document réponse sont représentés 20 pas de temps de et parmi les 1023 au total.
Q29. Compléter sur la figure du document réponse les séquences des registres et .
Sur la figure C du document réponse sont représentés les 20 bits NRZ du code correspondant aux 20 pas de temps des registres et de la figure B .
Q30. Compléter sur la figure du document réponse la séquence du code .
Afin de pouvoir mener les calculs de façon simple, on considère désormais que les codes possèdent au total bits NRZ au lieu des 1023 en réalité. La valeur numérique de la puissance , déterminée à la question Q28, doit être recalculée en conséquence.
Q31. Calculer les autocorrélations pour . En exploitant le fait que le code possède par hypothèse la propriété , vérifier que le code satisfait bien la propriété .
Sur la figure 10 sont représentés 20 pas de temps correspondants à 3 codes binaires NRZ. L'un d'entre eux est celui du code .
Code a)
1
-1
-1
-1
-1
1
-1
1
1
-1
1
1
-1
1
1
1
1
-1
1
-1
Code b)
1
-1
1
1
-1
1
-1
1
1
-1
1
1
1
1
-1
1
1
-1
1
-1
Code c)
1
-1
-1
1
1
-1
-1
1
-1
1
1
-1
1
1
1
1
1
1
-1
-1
Figure 10 - Codes potentiels pour
Q32. En vous appuyant sur la propriété , déterminer lequel des 3 signaux correspond à un code valide.
4 Etude de la propagation du signal entre le satellite et le récepteur
Un signal radio-électrique se propageant dans le vide se déplace à vitesse constante (célérité de la lumière) et en ligne droite, si on fait abstraction des effets relativistes. Vis-à-vis de ces hypothèses, l'évolution du signal à travers les différentes couches de l'atmosphère va provoquer un allongement du temps de propagation du signal entre le satellite et le récepteur du fait de la diminution de la vitesse de propagation et de la déviation du signal. Le récepteur GPS estime la distance le séparant des satellites via la durée de propagation du signal : cet allongement du temps de propagation doit par conséquent être modélisé finement pour obtenir une estimation précise de cette distance.
Dans notre étude, nous considérons que les deux seules couches atmosphériques pouvant entraîner un allongement du temps de propagation sont la troposphère et l'ionosphère:
La troposphère est la couche la plus basse de l'atmosphère. Elle s'élève jusqu'à 15 km d'altitude environ. Le retard troposphérique dépend essentiellement des conditions climatiques (température, taux d'humidité, pression de l'air). Le retard troposphérique engendre une erreur sur la pseudo-distance pouvant aller jusqu'à 30 mètres.
Figure 11 - Couches atmosphériques prises en compte
L'ionosphère est la couche atmosphérique située de 60 km à 700 km d'altitude environ. C'est un milieu ionisé par les radiations solaires. Cette couche atmosphérique peut être considérée comme un plasma neutre. Elle est constituée de cations de charge positive d'une densité volumique de charge d'environ ainsi que d'électrons de charge de même densité volumique . Afin de déterminer le comportement du signal électrique dans l'ionosphère, on élabore un modèle dont les hypothèses sont les suivantes:
les cations sont immobiles du fait de leur grande masse devant celle d'un électron;
l'agitation thermique, les collisions entre les particules ainsi que l'interaction entre particules chargées sont négligées ;
la vitesse des porteurs de charges reste faible devant la vitesse de la lumière.
On se propose d'étudier dans ce milieu la propagation d'ondes du type
où la phase s'exprime en fonction de la pulsation et du vecteur d'onde par , et . Dans le plasma, les électrons ne sont soumis qu'à la force de Lorentz due aux champs électromagnétiques de l'onde.
La présente étude ne concerne que le retard ionosphérique. Le retard dû à la troposhpère est quant à lui compensé par le récepteur via des paramètres correctifs présents dans le message envoyé par le satellite.
Q33. Déterminer l'expression simplifiée de l'accélération de l'électron pour les vitesses faibles.
Q34. En supposant que l'amplitude de déplacement des porteurs de charges reste faible devant la longueur d'onde des charges, déterminer l'expression de défini par la relation : , où est le vecteur densité de courant. Cette expression doit être fonction de et .
Q35. Écrire les équations de Maxwell à l'intérieur de l'ionosphère. Déterminer l'expression de dans le plasma. Mettre en évidence une fréquence caractéristique dite fréquence plasma : donner son expression.
Numériquement, la fréquence plasma est évaluée à .
Q36. Les fréquences de l'ensemble des signaux GPS étant comprises entre 1000 MHz et 1600 MHz , justifier ce choix de gamme de fréquences.
Q37. Déterminer l'expression de l'indice du plasma défini par : où est la vitesse de phase, et en conclure si le milieu est dispersif ou non. Donner également l'expression de la vitesse de groupe à laquelle se déplace le signal électrique émis par le satellite dans l'ionosphère en fonction de et de .
L'ionosphère a une épaisseur et le nombre d'électrons reste constant dans cette épaisseur. On note l'angle d'élévation de l'onde électromagnétique entre la normale à la surface et la direction de propagation (voir la figure 12). En traversant l'ionosphère, la différence de vitesse retarde la réception du signal d'une durée , ce qui engendre une erreur de mesure de la distance entre le satellite et le récepteur selon la relation :
Q38. Ecrire le retard sous la forme , où est la fréquence de la porteuse du signal, en exprimant en fonction de la fréquence plasma , de la célérité , de l'angle d'élévation et de l'épaisseur de l'ionosphère traversée .
Figure 12 - Influence de l'ionosphère sur le trajet du signal
Dans la pratique, le satellite émet le même signal sur deux porteuses de fréquences et .
Q39. En exploitant la forme de de la question précédente, démontrer comment l'utilisation de ces deux porteuses permet de corriger les erreurs dues à la traversée de l'ionosphère.
5 Etude des traitements effectués par le récepteur
Soit un signal émis par un satellite . On suppose que l'horloge du récepteur et celles des satellites sont synchrones : elles constituent même référence temporelle. D'après les parties précédentes, le signal reçu au niveau du récepteur peut s'écrire : où
est le temps dans le référentiel temporel du récepteur (identique au temps GPS par hypothèse);
est le temps de propagation du signal dans l'espace;
est l'amplitude du signal reçu;
est le message de navigation contenant les informations sur la trajectoire du satellite , de fréquence et prenant comme valeurs (signal binaire NRZ);
est le code de Gold propre au satellite , de fréquence et comprenant une succession de bits NRZ, i.e. de valeur ;
est la fréquence de la porteuse sinsuoïdale, .
5.1 Analyse de la puissance des signaux reçus
Au niveau du récepteur, plusieurs éléments génèrent du bruit qui s'ajoute au signal reçu , entre autres les câbles de liaison et les différents circuits électroniques, tels que les amplificateurs linéaires intégrés. Au final, la puissance totale du bruit générée s'écrit , où est la densité spectrale de bruit totale et la largeur de bande de fréquence utile du signal. La fréquence du code étant d'environ 1 MHz , nous obtenons numériquement . De plus, .
Nota : si est une puissance exprimée en Watts, elle s'exprime en dBW selon la relation avec .
Q40. En prenant comme approximation , calculer la valeur numérique de la puissance du bruit , en .
Q41. Sachant que la puissance des signaux reçus est , calculer le rapport signal sur bruit . Commenter.
5.2 Analyse des traitements effectués par le récepteur
Les différents traitements effectués par le récepteur sur le signal reçu sont représentés sur la figure 13 .
Le signal reçu est tout d'abord transposé en fréquence afin de supprimer la porteuse et donc de récupérer un signal utile qui sera exploité dans la suite des traitements. Cette transposition en fréquence s'effectue grâce aux utilisations successives
d'un multiplieur, qui multiplie avec un signal sinusoïdal d'expression ;
d'un filtre passe-bas de fréquence de coupure .
Figure 13 - Traitements effectués par le récepteur
Q42. Préciser quelle relation doivent vérifier et . Les valeurs numériques de et permettent-elles de vérifier cette relation? En supposant le filtre passe-bas idéal, écrire alors l'expression du signal de sortie du filtre passe-bas sous la forme , où sera exprimé en fonction de et .
Un second multiplieur est ensuite utilisé, prenant comme deuxième entrée un générateur de codes de Gold identique à celui embarqué à bord des satellites GPS. L'horloge du récepteur étant supposée synchronisée aux satellites, le code généré par le récepteur est synchrone de celui généré par les satellites. Ne connaissant a priori ni l'identité du satellite émettant le signal reçu ni le temps de propagation du signal, le code généré par le récepteur est décalé de : les valeurs de et de sont choisies arbitrairement par le récepteur. La sortie du multiplieur est ensuite sommée sur les 1023 bits NRZ du signal résultant de cette multiplication pour obtenir une valeur .
Q43. Compte tenu de la différence de fréquence entre les signaux et , donner l'expression de en fonction de et .
Q44. Indiquer pour quelles valeurs de et la valeur maximale de la corrélation est atteinte. Donner alors l'expression de dans ce cas-là en fonction de , et .
Q45. En déduire les étapes pour que le récepteur soit capable d'identifier l'identité du satellite émetteur du signal et d'estimer le temps de propagation de celui-ci.
On retrouve ainsi les hypothèses de l'étude réalisée dans la partie 2, qui a mis en évidence la possibilité de déterminer une solution PVT.
6 Synthèse de l'étude
Q46. En exploitant la démarche de modélisation menée dans le sujet ainsi que les résultats de l'étude réalisée, synthétiser les éléments clés du fonctionnement du GPS.
Fin de l'énoncé.
Document-réponse à rendre avec votre copie
Figure A - Chronogramme des signaux générés par le satellite
1
0
1
0
1
0
0
1
1
0
0
0
0
1
0
1
0
1
1
1
1
1
0
1
0
1
0
0
1
1
0
0
0
0
1
0
1
0
1
1
1
1
1
0
1
0
1
0
0
1
0
0
0
0
0
1
0
1
0
1
0
1
1
1
0
1
0
1
0
0
1
0
0
0
0
0
1
0
1
0
1
0
1
1
1
0
1
0
1
0
1
1
0
0
0
0
0
1
0
1
1
1
0
1
1
1
0
1
0
1
1
1
1
0
0
0
0
0
1
0
1
1
1
0
1
1
1
0
1
0
1
1
1
1
0
0
0
0
0
1
1
1
1
1
0
1
1
1
0
1
1
1
1
1
1
0
0
0
0
0
0
1
1
1
1
0
1
1
1
0
0
1
1
1
1
1
0
0
0
0
1
0
1
1
1
1
0
1
1
1
1
0
1
1
1
1
1
0
0
0
0
1
0
1
1
1
1
0
1
1
0
1
0
1
1
1
1
1
0
0
1
0
1
0
1
1
1
1
0
1
1
0
1
0
1
1
1
1
1
0
0
1
0
1
0
1
1
1
1
0
0
1
0
1
0
1
1
1
1
1
0
0
1
0
1
0
1
1
1
1
1
0
1
0
1
0
1
1
1
1
Figure B - Générateur de code pas de temps pour (gauche) et (droite)
1
1
1
-1
-1
1
-1
-1
1
1
1
1
1
-1
Figure C - Générateur de code pas de temps du code
X ENS Modélisation PSI 2016 - Version Web LaTeX | WikiPrépa | WikiPrépa
