FR3066027A1 - Procede de positionnement par gnss - Google Patents
Procede de positionnement par gnss Download PDFInfo
- Publication number
- FR3066027A1 FR3066027A1 FR1753907A FR1753907A FR3066027A1 FR 3066027 A1 FR3066027 A1 FR 3066027A1 FR 1753907 A FR1753907 A FR 1753907A FR 1753907 A FR1753907 A FR 1753907A FR 3066027 A1 FR3066027 A1 FR 3066027A1
- Authority
- FR
- France
- Prior art keywords
- gnss
- biases
- receiver
- distances
- pseudo
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 239000000654 additive Substances 0.000 claims abstract description 5
- 230000000996 additive effect Effects 0.000 claims abstract description 5
- 238000005259 measurement Methods 0.000 claims description 27
- 238000012937 correction Methods 0.000 claims description 3
- 238000004590 computer program Methods 0.000 claims description 2
- 230000001419 dependent effect Effects 0.000 abstract 1
- 239000000243 solution Substances 0.000 description 17
- 239000011159 matrix material Substances 0.000 description 13
- 238000001514 detection method Methods 0.000 description 9
- 230000014509 gene expression Effects 0.000 description 9
- 230000008569 process Effects 0.000 description 7
- 230000000694 effects Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000012360 testing method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000000873 masking effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 230000007480 spreading Effects 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000005054 agglomeration Methods 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000002790 cross-validation Methods 0.000 description 1
- 230000001627 detrimental effect Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- LIWAQLJGPBVORC-UHFFFAOYSA-N ethylmethylamine Chemical compound CCNC LIWAQLJGPBVORC-UHFFFAOYSA-N 0.000 description 1
- 230000007717 exclusion Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000003116 impacting effect Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 235000015145 nougat Nutrition 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/22—Multipath-related issues
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Navigation (AREA)
Abstract
Le procédé de positionnement comprend l'estimation d'un ensemble de biais multitrajets sur la base d'un modèle théorique des écarts entre les pseudo-distances mesurées et les pseudo-distances estimées. Le modèle comprend un terme de géométrie dépendant de la position du récepteur GNSS et des émetteurs GNSS et un terme additif représentant l'ensemble desdits biais. Les inconnues du modèle comprennent la position courante du récepteur GNSS et lesdits biais. L'estimation des biais repose sur une minimisation des résidus entre les écarts et le modèle théorique, en fonction des inconnues, sous la contrainte que l'ensemble des biais doit être parcimonieux. Les pseudo-distances mesurées sont corrigées sur base des biais estimés et une position courante du récepteur GNSS est prédite sur la base des pseudo-distances mesurées ainsi corrigées. L'invention permet de réduire l'impact des multi-trajets sur la solution de positionnement et donc d'en améliorer l'exactitude.
Description
Domaine Technique [0001] De manière générale, l’invention concerne un procédé de positionnement par GNSS (acronyme de « Global Navigation Satellite System », système de positionnement par satellites), plus particulièrement un procédé adapté aux environnements caractérisés par la présence de multi-trajets. Un autre aspect de l’invention concerne un navigateur GNSS configuré pour la mise en œuvre du procédé.
Arrière-plan technologique [0002] L’invention s’intéresse à la problématique des multi-trajets. Le problème se pose surtout en milieu urbain, mais pas exclusivement. Le terme « multi-trajet » vise les effets qui résultent de la réception de signaux GNSS qui, sur leur parcours entre l’émetteur et le récepteur, ont été réfléchis par une surface. Ces signaux réfléchis peuvent se superposer aux signaux ayant suivi la ligne de mire (en anglais : « line of sight » ou « LOS ») entre leur émetteur et le récepteur ou, en cas de masquage, se substituer à ces signaux. La topographie des villes et des agglomérations favorise l’occurrence de multi-trajets et de masquage de satellites, notamment à cause des façades verticales en bordure des espaces fréquentés par les utilisateurs des GNSS.
[0003] Le problème des multi-trajets est d’autant plus important le plus grand nombre d’utilisateurs de GNSS (comme p.ex. GPS, GLONASS, COMPASS et GALILEO) se trouvent en milieu urbain. En outre, la densité géographique des chemins empruntables par les utilisateurs fait qu’une erreur de positionnement peut facilement se répercuter au niveau du navigateur et être perceptible par l’utilisateur.
[0004] Différentes stratégies existent pour assurer la détection des multi-trajets. Le brevet US 5 917 445 propose un test de détection de multi-trajets qui combine des mesures de la phase porteuse et des pseudo-distances réalisées avec plusieurs antennes GNSS espacées entre elles. L’ensemble des mesures disponibles est hautement redondant en cas de fonctionnement normal. Si un multi-trajet est détecté pour un signal particulier, les mesures en rapport avec ce signal ne sont pas prises en considération pour le positionnement. Cette solution nécessitant plusieurs antennes et des traitements indépendants des différents signaux reçus, elle n’est pas applicable aux récepteurs GNSS destinés au grand public.
[0005] L’article de M. Irsig 1er et G. W. Hein, Development of a Real-Time Multipath Monitor Based on Multi-Correlator Observations, Proceedings of the 18th International Technical Meeting of the Satellite Division of the Institute of Navigation, ION GNSS 2005, propose une méthode de détection des multi-trajets basée sur l’observation de la forme du pic de corrélation. Plusieurs points de corrélation sont mesurés et à partir de ceux-ci divers métriques sont établies. Pour chacune de ces métriques un test de détection est mis en place pour vérifier si cette dernière s’écarte trop de sa valeur théorique obtenue en absence de multi-trajets.
[0006] L’article d’O. M. Mubarak et A. G. Dempster, Exclusion of Multipath-Affected Satellites Using Early Late Phase, Journal of Global Positioning System, 2010, Vol. 9, No.2 :145-155, propose un test de détection de multi-trajets utilisant la différence de phase des points de corrélation « Early » (en avance) et « Late » (en retard), nommée ELP (« Early Late Phase »). Le test utilise les voies I et Q des corrélateurs et suppose que la mesure ELP a une amplitude non-nulle en présence de multi-trajets.
[0007] D’autres procédés de détection de multi-trajets exploitent le fait que la réflexion change la polarisation des signaux GNSS.
[0008] Les signaux de radionavigation émis par des satellites (ou pseudolites) d'un GNSS se présentent sous la forme d'une porteuse modulée par une forme d'onde d'étalement contenant un code binaire pseudo-aléatoire. La modulation de la porteuse provoquant l'étalement du spectre autour de la fréquence de la porteuse, les signaux de radionavigation sont souvent appelés « en spectre étalé ». Les codes pseudoaléatoires représentent un identificateur du signal et donc du satellite émetteur. Connus des récepteurs, ils permettent à ceux-ci un Accès Multiple à Répartition de Code (AMRC). A l’exception des signaux « pilotes », les signaux GNSS transportent en outre des données utiles (p. ex. le message de navigation) sous forme d'une séquence binaire (à rythme nettement moins élevé que le code pseudo-aléatoire) modulée en plus sur la porteuse.
[0009] La réception d'un signal de radionavigation comprend sa démodulation à l'aide d'une réplique de la porteuse générée dans le récepteur par un oscillateur piloté par une boucle de poursuite de porteuse (ou boucle de poursuite de phase) et d'une réplique de la forme d'onde d'étalement (ou réplique de code) produite par un générateur de forme d'onde piloté par une boucle de poursuite de code. Le signal d'écart de phase entre la porteuse du signal entrant et la réplique de porteuse fournit à chaque pas de temps une première observable (l'observable de phase). Le signal de décalage temporel entre la forme d'onde d'étalement du signal entrant et la réplique de code produit à chaque pas de temps représente une deuxième observable (l'observable de code). Les observables de phase et de code sont ensuite utilisées pour estimer la pseudo-distance et la delta-distance (en anglais « delta range »). Une delta-distance équivaut à une estimation indépendante de la dérivée de la pseudodistance correspondante. Les étapes sont réalisées en parallèle sur les différents canaux du récepteur pour les différents signaux. L’ensemble des pseudo-distances et des delta-distances sont ensuite amenées dans le navigateur qui les combine et en dérive une estimation de la position, de la vitesse, du temps, etc.
[0010] Les techniques de détection de multi-trajets identifiées ci-dessus sont réalisées au niveau du récepteur GNSS et nécessitent l’accès aux boucles de poursuite ou au moins la connaissance de certains paramètres internes du récepteur. L’implémentation des techniques connues demande donc des modifications de l’architecture du récepteur, ce qui implique un effort considérable en termes de développement. En outre des adaptations sont nécessaires pour chaque type de récepteur.
[0011] La présente invention a pour objectif de proposer une solution qui peut être adaptée plus facilement à plusieurs types de récepteurs GNSS.
Description générale de l’invention [0012] Un premier aspect de l’invention concerne un procédé de positionnement par GNSS, comprenant: o l’obtention de mesures GNSS réalisées par un récepteur GNSS, les mesures comprenant un ensemble de pseudo-distances entre le récepteur GNSS et des émetteurs GNSS sur des satellites ; o l’estimation d’un ensemble de pseudo-distances prédites sur la base d’une position précédente du récepteur GNSS ; o l’estimation d’un ensemble de biais susceptibles d’être causés par des multi-trajets (ci-après aussi appelés « biais multi-trajets ») sur la base d’un modèle théorique des écarts entre les pseudo-distances mesurées et les pseudodistances estimées, le modèle comprenant un terme de géométrie dépendant de la position du récepteur GNSS et des émetteurs GNSS et un terme additif représentant l’ensemble desdits biais, les inconnues du modèle comprenant la position courante du récepteur GNSS et lesdits biais ; l’estimation des biais reposant sur une minimisation des résidus entre les écarts et le modèle théorique des écarts, en fonction des inconnues du modèle, sous la contrainte que l’ensemble des biais doit être parcimonieux (« contrainte de parcimonie ») ; o la correction des pseudo-distances mesurées sur base des biais estimés ; et o la prédiction d’une position courante du récepteur GNSS ayant réalisé les mesures GNSS sur la base des pseudo-distances mesurées ainsi corrigées.
[0013] Le procédé repose sur le fait que la présence de multi-trajet sur un signal GNSS se répercute dans la pseudo-distance mesurée associée comme un biais additif. Pour réduire l’effet des multi-trajets, ces biais peuvent être estimés et soustraits des pseudo-distances et delta-distances concernées. L’article A. Giremus, J.-Y. Tourneret et V. Calmettes, « A particle filtering approach for joint detection/estimation of multipath effects on GNSS measurements, » IEEE Trans. Signal Processing, vol. 55, no. 4, pp. 1275-1285, avril 2007, décrit des méthodes Monte-Carlo séquentielles (ou filtre particulaire) pour la réalisation de cette estimation. Toutefois, ces méthodes étant très coûteuses en termes d’efforts de calcul, leur implémentation en temps réel est très compliquée sinon impossible dans des applications pratiques. Le procédé selon le premier aspect de l’invention impose une contrainte de parcimonie sur les composantes du vecteur représentant les biais additifs. En d’autres mots, le procédé admet que la majorité des biais soient nuis. Cette hypothèse revient à supposer que la plupart des signaux GNSS reçus ne sont pas affectés par des multi-trajets ou que leurs effets peuvent être négligés pour la plupart des signaux GNSS.
[0014] La contrainte de parcimonie peut prévoir que plus que le quart, préférablement plus que la moitié, encore plus préférablement plus de deux tiers, toujours plus préférablement plus que trois quarts et le plus préférablement plus que quatre cinquièmes, des biais susceptibles d’être causés par des multi-trajets soient nuis (dans la mesure de la tolérance).
[0015] Un avantage intéressant du procédé selon le premier aspect de l’invention réside dans le fait qu’un traitement des multi-trajets au niveau du récepteur GNSS (c.-à-d. au niveau des boucles de poursuite) n’est pas nécessaire et éventuellement même contre-indiquée. Le présent procédé intervient en aval du récepteur GNSS, au niveau de l’étage normalement appelé « le navigateur. » Une modification de l’architecture (matérielle) du récepteur GNSS ou l’accès à des paramètres internes du récepteur ne sont donc pas nécessaires. Un autre avantage est qu’il existe un standard pour le format des données numériques produites en sortie du récepteur. La plupart des récepteurs respectent en effet le format RINEX (du terme anglais « Receiver Independent Exchange Format », format d’échange indépendant du récepteur). On note que le traitement des données produites par le récepteur se fait généralement au niveau logiciel (à l’aide d’un processeur à usage général ou d’un processeur spécialisé). L’implémentation du procédé de cette invention est donc relativement facile, p.ex. par une mise à jour du microprogramme (« firmware ») des jeux de puces (« chipsets ») GNSS existants.
[0016] La minimisation des résidus repose de préférence sur le critère des moindres carrés (la somme des carrés des résidus, c.-à-d. la norme f2, est minimisée), la contrainte de parcimonie étant implémentée comme une pénalisation qui augmente avec le nombre de biais multi-trajets non nuis.
[0017] Selon un mode de réalisation du procédé, l’estimation de l’ensemble des biais utilise un opérateur LASSO (acronyme du terme anglais « least absolute shrinkage and sélection operator », opérateur de sélection et de rétrécissement par moindre valeur absolue). L’opérateur LASSO permet, dans le cas de ce procédé, la régression avec le double objectif de minimiser les résidus et d’effectuer une sélection des biais multi-trajets directement lors de l’estimation de la position.
[0018] Selon un mode de réalisation préféré du procédé, l’estimation de l’ensemble des biais multi-trajets comprend une pénalisation de signaux en provenance d’émetteurs GNSS sur des satellites à basse élévation par rapport à des signaux en provenance d’émetteurs GNSS sur des satellites à haute élévation. La pénalisation d’un signal est choisie de sorte à augmenter la probabilité qu’un biais multi-trajet soit détecté sur ce signal si le signal provient d’un satellite ayant une élévation basse. En d’autres mots, on applique une pondération des signaux qui favorise la détection de multi-trajets sur les signaux qui proviennent de satellites éloignés du zénith.
[0019] De manière alternative ou additionnellement, l’estimation de l’ensemble des biais peut comprendre une pénalisation de signaux présentant un rapport C/No faible par rapport à des signaux présentant un rapport C/No élevé. Dans ce cas, la pénalisation d’un signal est donc choisie de sorte à augmenter la probabilité qu’un biais multi-trajet soit détecté sur ce signal si le signal présentant un rapport C/No élevé. Il est possible de combiner la pénalisation basée sur l’élévation et celle qui est basée sur le rapport C/No.
[0020] Il est également possible que l’estimation de l’ensemble des biais comprenne une pénalisation de signaux préalablement affectés par des multi-trajets. Il est possible de combiner cette pénalisation avec les pénalisations basées sur l’élévation et le rapport C/No.
[0021] Les choix de pénalisation sont justifiés par les constats que des signaux en provenance de satellites proches de l’horizon ont généralement une probabilité plus élevée d’être réfléchis sur le chemin parcouru entre le satellite et le récepteur. On observe également qu’un multi-trajet a tendance à aller de pair avec un rapport C/No faible. Finalement, si un multi-trajet est détecté sur un signal GNSS donné à un instant, il y a de fortes chances que le signal reste affecté de multi-trajet l’instant d’après.
[0022] Selon un mode de réalisation du procédé, l’estimation des biais utilise un paramètre de régularisation pour mettre en œuvre la contrainte de parcimonie, le paramètre de régularisation étant ajustable pour relâcher ou donner plus d’importance à ladite contrainte (c.-à-d. pour tolérer que plus ou moins de biais multi-trajet soient non nuis). Le paramètre de régularisation pourrait être ajusté de manière dynamique, p.ex. en fonction de la position du récepteur GNSS. Dans des zones cartographiées, on pourrait, p.ex. utiliser un paramètre reflétant l’angle solide que le ciel ouvert occupe à une position donnée, un paramètre reflétant la hauteur moyenne des bâtiments, un degré d’urbanisation (p.ex. l’indicateur DEGURBA). Le paramètre de régularisation pourrait également être ajusté en fonction des bruits C/No. On notera toutefois que la contrainte de parcimonie ne saurait être relâchée au point d’être dénudée de sens : le nombre de biais multi-trajet non nuis doit rester bien en-deçà du nombre des mesures.
[0023] De préférence, les mesures GNSS réalisées par le récepteur GNSS sont obtenues au format RINEX.
[0024] De préférence, l’estimation des pseudo-distances prédites repose non seulement sur une position précédente mais également sur une vitesse précédente et un biais récepteur du récepteur GNSS. (On désigne par « biais récepteur » l’écart entre l’horloge interne du récepteur et une horloge de référence, p.ex. le temps GNSS.) Dans ce cas, le terme de géométrie du modèle utilisé pour l’estimation des biais multitrajet dépend également de la vitesse et du biais récepteur du récepteur GNSS ainsi que de la position et de la vitesse des émetteurs GNSS (c.-à-d. des satellites), et les inconnues du modèle comprennent en outre la vitesse et le biais récepteur actuels du récepteur GNSS.
[0025] Selon un mode de réalisation avantageux de l’invention, la prédiction d’une position courante du récepteur GNSS ayant réalisé les mesures GNSS est réalisée par un filtre de Kalman (étendu) utilisant les pseudo-distances mesurées et corrigées pour mettre à jour l’état prédit.
[0026] Selon un mode de réalisation préféré, le procédé fonctionne en temps réel. En d’autres termes, le procédé est réalisé de telle façon à respecter strictement des contraintes de temps prédéfinies, en particulier le retard maximal (p. ex. 1 seconde) que les résultats peuvent accumuler avant de parvenir à l’utilisateur (p. ex. un automobiliste).
[0027] Un deuxième aspect de l’invention concerne un programme d’ordinateur comprenant des instructions de code de programme pour l’exécution des étapes du procédé tel que décrit, lorsque ledit programme est exécuté sur un ordinateur.
[0028] Un troisième aspect de l’invention concerne un navigateur GNSS, comprenant une mémoire et un processeur, ainsi que des instructions de code de programme enregistrés dans la mémoire et accessibles par ledit processeur, le processeur effectuant le procédé décrit lorsqu’il exécute les instructions de code de programme.
[0029] Le navigateur GNSS peut comprendre ou être relié à un récepteur GNSS configuré pour obtenir les mesures GNSS et les mettre à disposition du processeur lorsque les instructions de code de programme sont exécutées.
Brève description des dessins [0030] D'autres particularités et caractéristiques de l'invention ressortiront de la description détaillée de certains modes de réalisation avantageux présentés ci-dessous, à titre d'illustration, avec référence aux dessins annexés qui montrent :
Fig. 1: le schéma d’un appareil de positionnement selon un premier mode de réalisation ;
Fig. 2: le schéma d’un appareil de positionnement selon un deuxième mode de réalisation ;
Fig. 3: le schéma d’un procédé de positionnement selon un aspect de l’invention.
Description détaillée d’un mode de réalisation de l’invention [0031] La tâche principale d’un navigateur GNSS consiste à trouver une solution au problème de navigation. Une solution comprend la position, le temps et, optionnellement, le vecteur vitesse du récepteur. Si la vitesse est estimée, on parle aussi de solution « PVT » (position, vitesse et temps). Afin de trouver la solution au problème de navigation, le navigateur GNSS se base sur des mesures GNSS, comprenant des pseudo-distances, fournies par le récepteur GNSS.
[0032] Un récepteur GNSS fournit, pour chaque signal GNSS qu’il poursuit, une pseudo-distance définie par : (Eq. 1) où :
l’exposant i identifie le satellite qui a envoyé le signal GNSS (i = 1, ..., N, N étant le nombre de satellites « visibles » par le récepteur), o p1 représente la pseudo-distance mesurée entre le récepteur (RX) et le satellite i par une mesure dite « de code » ; o x1 est la position (connue) du satellite i à l’émission du signal GNSS reçu par le récepteur : x1 = (xYyYzf7 ; o xRX est la position du récepteur (qui doit être estimée) à la réception du signal : xRX = (x>y>z)T ; o |H|2 est la distance géométrique (norme o b est le biais récepteur (inconnu a priori), exprimé en unités de longueur ; et ο ει est un terme d’erreur associé au signal en provenance du satellite i : ce terme inclut notamment le délai troposphérique, le délai ionosphérique, le biais émetteur i (c.-à-d. l’écart entre l’horloge interne du satellite i et une horloge de référence, p.ex. le temps GNSS), l’incertitude de la position du satellite i, l’effet du multi-trajet et tout bruit non modélisé.
[0033] On note que p' diffère de la distance géométrique par les termes b et ε\ d’où la désignation de « pseudo-distance ».
[0034] On peut estimer une position courante du récepteur en linéarisant l’Eq. 1 autour d’une position précédemment estimée par le navigateur GNSS (p.ex. la dernière position estimée par le navigateur). On trouve une expression de la forme :
(Eq. 2) où : o yp (e IRW) est la différence entre la pseudo-distance mesurée et la pseudodistance précédemment estimée ; o G (g IRWx4) est la matrice jacobienne associée au système linéarisé ; o x (e IR4) est la différence entre l’état (position et biais récepteur) estimé à l’instant précédent et l’état actuel ; o mp (g IRw) est un vecteur de biais multi-trajet ; et o np (e IRW) est un vecteur de bruit gaussien.
[0035] La matrice G peut être exprimée comme :
(Eq. 3) où:
(Eq. 4) x0 étant la position de référence autour de laquelle l’expression Eq. 1 a été linéarisée. [0036] Par différentiation de l’expression Eq. 1 par rapport au temps, on obtient :
(Eq. 5) [0037] Il convient de noter que la grandeur // est mesurée indépendamment de la pseudo-distance p' ; on l’appelle « observable de phase », « observable de porteuse », « delta-distance » ou encore « observable Doppler ». On a, en effet, p1 = -^porteuse où Δ/' est le décalage Doppler mesuré pour le signal i et kporteuse est
la longueur d’onde nominale de la porteuse du signal. Par linéarisation de l’Eq. 5, on obtient une expression similaire à l’Eq. 2 :
(Eq. 6) où: o yr (e IRW) est la différence entre l’observable de phase mesurée et l’observable de phase précédemment estimée ; o x (e IR4) est la différence entre l’état (vecteur vitesse et dérive de l’horloge du récepteur) estimé à l’instant précédent et l’état actuel ; o mr (g IRw) est un vecteur de biais multi-trajet ; et o nr (e IRW) est un vecteur de bruit gaussien.
[0038] Selon l’invention la résolution du problème de navigation se fait sous la contrainte que le vecteur des biais multi-trajet
soit un vecteur parcimonieux.
[0039] Pour simplifier les expressions, une autre notation sera appliquée par la suite. Un vecteur d’état x sera exprimé de la façon suivante : x = (x,x,y,ÿ,z,z,b, è) , (Eq. 7) où (x,y,z)T est le vecteur position du récepteur, (x,ÿ,z)T le vecteur vitesse du récepteur, b le biais récepteur et b la dérive de l’horloge récepteur.
[0040] Si on note Ôy = [^p] et n = [™p], les équations 2 et 6 deviennent :
(Eq. 8) où H est la matrice jacobienne (ou matrice des cosinus directeurs) et δχ est la différence entre le vecteur d’état à l’instant actuel et le vecteur d’état à l’instant précédent.
[0041] Pour estimer les biais multi-trajet impactant les mesures de code et de phase, on minimise l’expression suivante :
(Eq. 9) où ôyth(ôx,m) = Ηδχ - m est un modèle théorique de l’écart entre l’état mesuré et l’état estimé, W est une matrice diagonale de pénalisation ayant la dimension 2N x
2N, Il Il ! est la norme f1 et λ est un paramètre de régularisation. Le dernier terme de l’Eq. 9 traduit la contrainte que l’ensemble des biais doit être parcimonieux.
[0042] La matrice diagonale IV peut pénaliser les signaux en provenance d’émetteurs GNSS sur des satellites à basse élévation par rapport à des signaux en provenance d’émetteurs GNSS sur des satellites à haute élévation. Les éléments diagonaux de IV pourraient alors prendre la forme suivante :
(Eq. 10) où Elevj est l’élévation de l’émetteur GNSS lors de l’émission du signal GNSS.
[0043] Une autre possibilité serait de définir les éléments diagonaux de IV par :
(Eq. 10’) où Elevseuii est une élévation-seuil ; p.ex., Elevseuii = 5° ou 10°. Cette expression a été proposée dans l’article E. Realini et. al, Meas. Sci. Technol. 24, 115010 (2013).
[0044] La matrice diagonale IV peut alternativement ou additionnellement pénaliser les signaux présentant un rapport C/No faible par rapport à des signaux présentant un rapport C/No élevé. Les éléments diagonaux de IV pourraient alors prendre la forme suivante :
(Eq. 11) [0045] Une autre possibilité serait de définir les éléments diagonaux de IV par (voir l’article E. Realini et. al) :
(Eq. 11’) pour i = 1 ,.,.,Ν et où £ = C/Noi, T est un seuil, a, A et F sont des paramètres de réglages. P. ex. :(T, a, F, A) = (45, 80, 20,1/30).
[0046] La matrice diagonale IV peut aussi pénaliser de signaux préalablement affectés par du multi-trajet précédemment. Les éléments diagonaux de IV pourraient alors prendre la forme suivante :
(Eq. 12)
où hj est le biais estimé précédemment affectant le signal provenant de l’émetteur GNSS i.
[0047] Selon un mode de réalisation préféré, les pénalisations décrites ci-dessus sont combinées. Les éléments diagonaux de IV pourraient alors prendre la forme suivante :
(Eq. 13) qui correspond à la multiplication des éléments diagonaux des Eqs. 10, 11 et 12. Alternativement, les éléments diagonaux donnés aux Eqs. 10’, 11’et 12, Eqs. 10, 11’ et 12 ou Eqs. 10’, 11 et 12 pourraient être multipliés.
[0048] Il sera apprécié qu’une pénalisation différente pourrait être affectée à l’observable de code et à l’observable de phase.
[0049] Le paramètre de régularisation λ peut être fixe ou ajusté de manière dynamique en fonction de la position du récepteur GNSS, p. ex. sur base de l’indicateur DEGURBA, et/ou par validation croisée. Le paramètre de régularisation permet de relâcher ou donner plus d’importance à ladite contrainte, c.-à-d. pour tolérer que plus ou moins de biais multi-trajets soient non nuis.
[0050] La solution de minimisation de l’Eq. 9 selon Sx admet une expression explicite : Sx = (H' H^H'f Sy - ni). Cette expression correspond à la solution, sans contrainte, au sens des moindres carrés pour des mesures (Sy - m). Après injection de cette solution pour Sx dans l’Eq. 9, la minimisation peut être reformulée de la façon suivante : (Eq. 14) où I2N est la matrice identité ayant une dimension 2N x 2N et P = Η^Η)-1!!7. De plus, en définissant : =(I2„-P)<Sy, (Eq. 15) H = (I2JV - PM1 , (Eq. 16) Q = Wm, (Eq. 17) il est finalement possible de réécrire le problème de minimisation défini à l’Eq. 9 de la façon suivante :
(Eq. 18) [0051] Dans l’hypothèse où m est un vecteur parcimonieux, θ est bien un vecteur parcimonieux lui aussi car la matrice W est diagonale. Le problème de minimisation de l’Eq. 18 a une forme permettant d’utiliser l’algorithme LASSO pour trouver une solution. Après avoir trouvé une solution Θ, il suffit d’inverser l’Eq. 17 pour obtenir le vecteur m des biais multi-trajets estimés : m = M7_1Ô.
[0052] La Fig. 1 illustre système de positionnement 10 selon un mode de réalisation de l’invention. Le système 10 comprend une antenne radiofréquence 12 configurée pour recevoir des signaux radiofréquence dans l’espace (p. ex. les signaux dans la bande de fréquences L1 et/ou L5). L’antenne 12 est connectée à un processeur 14 de signaux GNSS (p.ex. un jeu de puces GNSS). Le processeur 14 traite les signaux GNSS de l’antenne 12 afin d’obtenir une solution PVT du problème de navigation. Le traitement des signaux par le processeur 14 est réalisé en plusieurs étapes. Un étage de tête radiofréquence (« RF front end » en anglais) est connecté à l’antenne 12 et convertit les signaux GNSS reçus par l’antenne 12 en signaux numérisés en bande de base (signaux dits « IQ » pour « in-phase » et « quadrature » en anglais). Un processeur bande de base 18 (« baseband processor » en anglais) ou « récepteur GNSS » acquiert et poursuit les signaux GNSS en provenance des différents satellites et fournit les mesures GNSS correspondantes. Les mesures GNSS comprennent les observables GNSS (notamment les observables de code, respectivement de phase), le rapport C/No. Le processeur bande de base 18 peut également extraire les données des messages de navigation (éphémérides, données almanach, corrections d’horloge émetteur, etc.) véhiculés par les signaux GNSS. (De manière alternative, ces données peuvent être obtenues par d’autres canaux, p.ex. via un service de radiodiffusion terrestre, une liaison Internet ou autre.) Les mesures GNSS et les données GNSS peuvent être présentées au format RINEX.
[0053] Les mesures GNSS et les données GNSS sont finalement traitées par un navigateur GNSS 20. Le logiciel du navigateur GNSS 20 estime dans un premier temps un ensemble de biais multi-trajets et corrige les observables GNSS associées. L’estimation de l’ensemble de biais causés par des multi-trajets est réalisée selon les lignes développées ci-dessus. Dans un second temps, le logiciel résout les équations de navigation, p. ex. grâce à un filtre de Kalman étendu, sur base des observables
GNSS corrigées. Les données de sortie du système de positionnement 10 comprennent la solution PVT des équations de navigation, p.ex. au format CSV ou NMEA.
[0054] Le mode de réalisation illustré à la Fig. 1 est particulièrement intéressant lorsqu’il est implémenté par le constructeur du système de positionnement 10. En effet, le constructeur a la possibilité de mettre à jour le logiciel du navigateur 20.
[0055] La Fig. 2 illustre un système de positionnement 10 selon un second mode de réalisation de l’invention. Il est particulièrement adapté lorsqu’une personne tierce souhaite tenir compte des multi-trajets lors de la résolution des équations de navigation mais n’a pas la possibilité de modifier le logiciel présent dans le navigateur GNSS 20.
[0056] Dans ce mode de réalisation, le processeur bande de base 18 met les données de mesure GNSS à disposition de l’utilisateur. Les téléphones portables utilisant, p. ex., le système d’exploitation « Android Nougat », version 7 ou plus élevée, fournissent à l’utilisateur les données de mesures GNSS brutes. Il est alors possible d’implémenter un navigateur GNSS personnalisé 22 qui réalise le procédé de positionnement par GNSS tel que décrit dans ce document.
[0057] Il est à noter que cette configuration peut être intéressante pour comparer la solution PVT fournie par un navigateur conventionnel (ne corrigeant pas les multi-trajets) avec la solution PVT fournie par le procédé de la présente invention.
[0058] La Fig. 3 illustre un procédé réalisé dans un navigateur GNSS, selon un mode Γηί
de réalisation de l’invention. A chaque pas de temps, le vecteur d’observables y = U est tout d’abord obtenu d’un récepteur GNSS 18. Ensuite, le navigateur GNSS détermine le vecteur Sy = y- ypréd (étape 24), où ypréd est le vecteur d’observables prédit sur la base de la solution PVT précédente. Sy représente le vecteur des écarts entre les pseudo-distances mesurées y et les pseudo-distances estimées yprèd. Si on marque les grandeurs à un instant k par l’indice correspondant, on a, pour yprèd à l’instant k :
(Eq. 19) où et Xfc_! sont, respectivement, la matrice jacobienne et le vecteur d’état à l’instant précédent.
[0059] Le navigateur GNSS estime ensuite les biais multi-trajets (c.-à-d. le vecteur m, voir Eq. 18) sont ensuite déterminés 26 et corrige ôy, le vecteur des écarts entre les pseudo-distances mesurées y et les pseudo-distances estimées ypréd :
(Eq. 20) [0060] Le vecteur ôycorr est injecté dans un filtre de Kalman étendu 28 à chaque pas de temps en tant qu’observation. Le filtre de Kalman 28 utilise un modèle tenant compte de certaines données GNSS 30, comme, p.ex., de la position et de la vitesse des satellites, de la position et de la vitesse précédemment estimées du récepteur, des erreurs relativistes, des dérives et décalages des horloges des satellites et du récepteur, d’erreur sur l’estimation du temps d’arrivée des signaux GNSS, du retard de groupe, des biais ionosphériques (p.ex. sur base du modèle Klobuchar), des biais troposphériques (p. ex. sur base de la fonction de Niell), etc.
[0061] Le filtre de Kalman 28 fournit en sortie, à l’instant k, le vecteur variation d’état δχ , qui est la variation entre l’état précédemment estimé (k-1 ) et l’état actuel (k), ainsi que le vecteur d’état x = xk. Le nouvel état x = xk peut être trouvé par addition du vecteur variation d’état δχ à l’état précédent x^. D’autre part, le filtre de de Kalman 28 met également à jour la matrice H, ainsi que p.ex les décalages et dérives des horloges, les erreurs relatives au temps d’arrivée des signaux GNSS, la matrice de propagation (de transition), etc, pour la prochaine itération (k <- k + 1).Il a été dit plus haut que le paramètre de régularisation λ (voir Eq. 9) peut être fixe ou ajusté de manière dynamique. Un choix possible est λ = 1. Une valeur plus élevée de λ (tous les autres paramètres restant inchangés) donne plus d’importance à la contrainte de parcimonie, c.-à-d. fait ressortir un vecteur des biais multi-trajets plus parcimonieux (davantage de composantes de la solution m sont égales à 0). Une valeur de λ moins élevée relâche la contrainte de parcimonie. En cas d’ajustement dynamique de A, il est possible de définir une plage de valeurs dans lesquelles λ pourra être choisi en fonction de certains paramètres qui sont surveillés. De manière alternative, des valeurs discrètes peuvent être définies parmi lesquelles λ pourra être choisi.
[0062] Il convient de noter que le procédé fonctionne également si le récepteur GNSS ne fournit pas les delta-distances p1 (i = Ι,.,.,Ν). Dans ce cas, le vecteur des biais multi-trajets m à estimer se réduit à mp
[0063] Alors que des modes de réalisation particuliers viennent d’être décrits en détail, l’homme du métier appréciera que diverses modifications et alternatives à ceux-là puissent être développées à la lumière de l’enseignement global apporté par la présente divulgation de l’invention. Par conséquent, les agencements et/ou procédés spécifiques décrits ci-dedans sont censés être donnés uniquement à titre d’illustration, sans intention de limiter la portée de l’invention.
Claims (16)
- Revendications1. Procédé de positionnement par GNSS, caractérisé par: l’obtention de mesures GNSS (y) réalisées par un récepteur GNSS, les mesures (y) comprenant un ensemble de pseudo-distances entre le récepteur GNSS et des émetteurs GNSS sur des satellites ; l’estimation d’un ensemble de pseudo-distances prédites sur la base d’une position précédente du récepteur GNSS ; l’estimation (26) d’un ensemble de biais susceptibles d’être causés par des multi-trajets sur la base d’un modèle théorique des écarts entre les pseudo-distances mesurées et les pseudo-distances estimées, le modèle comprenant un terme de géométrie (H) dépendant de la position du récepteur GNSS et des émetteurs GNSS et un terme additif représentant l’ensemble desdits biais, les inconnues du modèle comprenant la position courante du récepteur GNSS et lesdits biais ; l’estimation (26) des biais reposant sur une minimisation des résidus entre les écarts et le modèle théorique des écarts, en fonction des inconnues du modèle, sous une contrainte de parcimonie de l’ensemble des biais ; la correction (26) des pseudo-distances mesurées sur base des biais estimés ; et la prédiction d’une position courante (x) du récepteur GNSS ayant réalisé les mesures GNSS (y) sur la base des pseudo-distances mesurées ainsi corrigées.
- 2. Procédé selon la revendication 1, dans lequel la contrainte de parcimonie prévoit que plus que le quart, préférablement plus que la moitié, encore plus préférablement plus de deux tiers, toujours plus préférablement plus que trois quarts et le plus préférablement plus que quatre cinquièmes, des biais susceptibles d’être causés par des multi-trajets soient nuis.
- 3. Procédé selon la revendication 1 ou 2, dans lequel l’estimation (26) de l’ensemble des biais utilise un algorithme LASSO.
- 4. Procédé selon l’une quelconque des revendications 1 à 3, dans lequel l’estimation (26) de l’ensemble des biais comprend une pénalisation de signaux en provenance d’émetteurs GNSS sur des satellites à basse élévation par rapport à des signaux en provenance d’émetteurs GNSS sur des satellites à haute élévation.
- 5. Procédé selon l’une quelconque des revendications 1 à 4, dans lequel l’estimation (26) de l’ensemble des biais comprend une pénalisation de signaux présentant un rapport C/No faible par rapport à des signaux présentant un rapport C/No élevé.
- 6. Procédé selon l’une quelconque des revendications 1 à 5, dans lequel l’estimation (26) de l’ensemble des biais comprend une pénalisation de signaux préalablement affectés par des multi-trajets.
- 7. Procédé selon l’une quelconque des revendications 1 à 6, dans lequel l’estimation (26) des biais utilise un paramètre de régularisation pour mettre en œuvre la contrainte de parcimonie, le paramètre de régularisation étant ajustable pour relâcher ou donner plus d’importance à ladite contrainte.
- 8. Procédé selon la revendication 7, dans lequel le paramètre de régularisation est ajusté de manière dynamique.
- 9. Procédé selon la revendication 8, dans lequel le paramètre de régularisation est ajusté en fonction de la position (x) du récepteur GNSS.
- 10. Procédé selon l’une quelconque des revendications 1 à 9, dans lequel les mesures GNSS (y) réalisées par un récepteur GNSS sont obtenues au format RINEX.
- 11. Procédé selon l’une quelconque des revendications 1 à 10, dans lequel l’estimation des pseudo-distances prédites repose également sur une vitesse précédente et un biais récepteur du récepteur GNSS et dans lequel le terme de géométrie (H) du modèle utilisé pour l’estimation des biais susceptibles d’être causés par des multi-trajets dépend également de la vitesse et du biais récepteur du récepteur GNSS ainsi que de la position et de la vitesse des émetteurs GNSS, les inconnues du modèle comprenant en outre la vitesse et le biais récepteur actuels du récepteur GNSS.
- 12. Procédé selon l’une quelconque des revendications 1 à 11, dans lequel la prédiction d’une position courante (x) du récepteur GNSS ayant réalisé les mesures GNSS (y) est réalisée par un filtre de Kalman étendu (28) utilisant les pseudo-distances mesurées et corrigées pour mettre à jour l’état prédit.
- 13. Procédé selon l’une quelconque des revendications 1 à 11, fonctionnant en temps réel.
- 14. Programme d’ordinateur comprenant des instructions de code de programme pour l’exécution des étapes du procédé selon l’une quelconque des revendications 1 à 13, lorsque ledit programme est exécuté sur un ordinateur.
- 15. Navigateur GNSS (22), comprenant une mémoire et un processeur, ainsi que des instructions de code de programme enregistrés dans la mémoire et accessibles par ledit processeur, le processeur effectuant le procédé selon l’une quelconque des revendications 1 à 13 lorsque les instructions de code de programme sont exécutées par le processeur.
- 16. Navigateur GNSS (22) selon la revendication 15, comprenant un récepteur GNSS configuré pour obtenir les mesures GNSS (y) et les mettre à disposition du processeur lorsque les instructions de code de programme sont exécutées.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1753907A FR3066027B1 (fr) | 2017-05-03 | 2017-05-03 | Procede de positionnement par gnss |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1753907 | 2017-05-03 | ||
FR1753907A FR3066027B1 (fr) | 2017-05-03 | 2017-05-03 | Procede de positionnement par gnss |
Publications (2)
Publication Number | Publication Date |
---|---|
FR3066027A1 true FR3066027A1 (fr) | 2018-11-09 |
FR3066027B1 FR3066027B1 (fr) | 2020-07-31 |
Family
ID=59811423
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR1753907A Active FR3066027B1 (fr) | 2017-05-03 | 2017-05-03 | Procede de positionnement par gnss |
Country Status (1)
Country | Link |
---|---|
FR (1) | FR3066027B1 (fr) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112415541A (zh) * | 2019-08-22 | 2021-02-26 | 安波福技术有限公司 | 用于gnns反射表面映射与位置定位估计的系统与方法 |
FR3136285A1 (fr) | 2022-06-02 | 2023-12-08 | Commissariat à l'Energie Atomique et aux Energies Alternatives | Procédé de localisation d'une unité de navigation |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120127029A1 (en) * | 2010-11-22 | 2012-05-24 | The Charles Stark Draper Laboratory, Inc. | Coded filter |
-
2017
- 2017-05-03 FR FR1753907A patent/FR3066027B1/fr active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120127029A1 (en) * | 2010-11-22 | 2012-05-24 | The Charles Stark Draper Laboratory, Inc. | Coded filter |
Non-Patent Citations (1)
Title |
---|
FEI XIANG ET AL: "A multipath mitigation discriminator for GPS receiver", AEU - INTERNATIONAL JOURNAL OF ELECTRONICS AND COMMUNICATIONS, ELSEVIER, AMSTERDAM, NL, vol. 67, no. 10, 1 October 2013 (2013-10-01), pages 839 - 847, XP001585661, ISSN: 1434-8411, DOI: HTTP://DX.DOI.ORG/10.1016/J.AEUE.2013.04.007 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112415541A (zh) * | 2019-08-22 | 2021-02-26 | 安波福技术有限公司 | 用于gnns反射表面映射与位置定位估计的系统与方法 |
FR3136285A1 (fr) | 2022-06-02 | 2023-12-08 | Commissariat à l'Energie Atomique et aux Energies Alternatives | Procédé de localisation d'une unité de navigation |
EP4290267A1 (fr) | 2022-06-02 | 2023-12-13 | Commissariat à l'énergie atomique et aux énergies alternatives | Procédé de localisation d'une unité de navigation |
Also Published As
Publication number | Publication date |
---|---|
FR3066027B1 (fr) | 2020-07-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
FR3076354B1 (fr) | Procede de controle de l'integralite de l'estimation de la position d'un porteur mobile dans un systeme de mesure de positionnement par satellite | |
CA2682278C (fr) | Procede de traitement de signaux de radionavigation | |
CN111025346A (zh) | 一种快速估计gnss精密卫星钟差的方法及存储介质 | |
EP2335085B1 (fr) | Traitement de signaux de radionavigation utilisant une combinaison widelane | |
EP2245479A1 (fr) | Systeme de navigation a hybridation par les mesures de phase | |
EP1654559A2 (fr) | Procede de generation de corrections d'horloges pour systeme gps differentiel a grande distance ou mondial | |
EP3223038B1 (fr) | Procédé de geopositionnement par satellites et terminal associé | |
EP2755050B1 (fr) | Dispositif de poursuite de signaux de radionavigation | |
EP2511734A2 (fr) | Récepteur bi-fréquences de positionnement par satellites et procédé de réception associé | |
FR3066027A1 (fr) | Procede de positionnement par gnss | |
EP2221628B1 (fr) | Procédé de lever d'ambiguïté, procédé de localisation d'un récepteur par radionavigation comprenant une étape de lever d'ambiguïté et récepteur de localisation | |
CA2754178C (fr) | Procede de geopositionnement utilisant des donnees d'assistance | |
EP2442137B1 (fr) | Système d'augmentation de la disponibilité et des performances d'un système de géolocalisation par satellite | |
FR2918765A1 (fr) | Procede de determination d'une erreur d'asservissement dans une boucle de poursuite d'un code pseudo-aleatoire. | |
FR3020687A1 (fr) | Procede de determination de la position d'un recepteur de systeme de navigation par satellites, et systeme associe | |
EP1493041B1 (fr) | Procede et dispositif de determination de la position relative de deux points | |
FR2764708A1 (fr) | Perfectionnements aux procedes et systemes de localisation radiosatellitaire en temps reel, notamment de type gps | |
EP4107549B1 (fr) | Procédé de navigation d'un porteur à l'aide d'un filtre de kalman estimant un état de navigation du porteur | |
FR3137972A1 (fr) | Méthode MAC de monitoring, avec compensation de biais commun, de l’intégrité d’un procédé de positionnement ponctuel par balises virtuelles | |
WO2022090157A1 (fr) | Procédé de détection d'un masquage d'un ou plusieurs satellites, dispositif électronique de détection et produit programme d'ordinateur associé | |
FR2914999A1 (fr) | Procede de mitigation du phenomene de multitrajet dans un recepteur de positionnement par satellites |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PLFP | Fee payment |
Year of fee payment: 2 |
|
PLSC | Publication of the preliminary search report |
Effective date: 20181109 |
|
PLFP | Fee payment |
Year of fee payment: 3 |
|
PLFP | Fee payment |
Year of fee payment: 4 |
|
PLFP | Fee payment |
Year of fee payment: 5 |
|
PLFP | Fee payment |
Year of fee payment: 6 |
|
PLFP | Fee payment |
Year of fee payment: 7 |
|
PLFP | Fee payment |
Year of fee payment: 8 |