Le domaine de l'invention est celui de la réception d'un signal de radionavigation provenant d'un système de positionnement par satellites tel que le système GPS (acronyme de l'expression anglo-saxonne « Global Positioning System »).

Dans la suite, on désigne par le terme générique récepteur GNSS tout type de récepteur GPS, GLONASS, GALILEO, etc.

Pour bien fonctionner, les récepteurs actuels nécessitent généralement une réception à vue directe des satellites. Le positionnement se détériore rapidement notamment en termes de précision et de temps d'acquisition, lorsque la réception est perturbée comme c'est le cas à l'intérieur d'un bâtiment ou plus généralement en milieu dégradé (obstacles, multi-trajets, interférence, intérieur d'un bâtiment, etc.).

Lors de conditions de réception dégradées la récupération du signal de navigation à faible rapport signal sur bruit est possible mais nécessite un temps d'intégration plus long pour filtrer le bruit.

La difficulté réside dans le fait que pour intégrer le signal reçu sur un temps long et récupérer de l'information il faut que la dynamique apparente au niveau du récepteur soit très faible pendant ce temps. Lorsque le récepteur (ou son antenne) est fixe il n'y a plus de dynamique due au mouvement. Par contre il reste la dynamique apparente due à la dérive de l'horloge du récepteur.

On va préciser en quoi consiste cette dérive en rappelant le principe de fonctionnement d'un récepteur GNSS. On a représenté figure 1a un code émis par un satellite à t _ém ission et figure 1 b, ce code arrivant au récepteur à tréception, soit Δt _pr op = tréception - témission secondes plus tard, Δt _pr op étant le temps de propagation du signal entre le satellite et le récepteur. Ce temps de propagation multiplié par la vitesse de la lumière dans le vide donne une mesure de distance ; il est mesuré de la manière suivante. Une réplique du code émis par le satellite et représentée figure 1c peut être générée localement par le récepteur, en phase. Le décalage entre le code reçu (figure 1 b) et le code local généré (c'est-à-dire la réplique représentée figure 1c) correspond à l'écart de temps Δtprop recherché. Ce décalage représenté figure 1d est mesuré en mettant en

phase le code reçu et le code local ; le critère de mise en phase correspond à la maximisation de la fonction de corrélation des deux codes. Si les horloges du satellite et du récepteur sont parfaitement synchronisées, le décalage qui maximise la fonction de corrélation des deux codes fournit le temps Δtprop-

Dans la pratique ce n'est pas le cas. On a représenté figure 2 le temps réellement mesuré ; témission est le temps vrai auquel le signal est émis par le satellite, tréception le temps vrai auquel le signal atteint le récepteur, t _h ce temps indiqué par l'horloge et Δt _h le biais de l'horloge du récepteur. Finalement le temps de propagation mesuré est : Δt _pr op mesuré = Δt _pr op + Δt _h . Dans le cas idéal d'une horloge parfaite, Δt _h = 0.

De plus lorsque le temps d'intégration de la corrélation entre le signal local et le signal reçu est long, le biais Δt _h de l'horloge du récepteur peut varier. En fait dans le cas des applications grand public les horloges à faible coût utilisées dans les récepteurs ont des dérives d _h (t) très importantes et sont très instables. On rappelle que la dérive de l'horloge est par définition la dérivée par rapport au temps du biais d'horloge : dh(t) = d (Δt h (t)) /dt.

Pour la réception du signal à l'intérieur d'un bâtiment, cette dérive d'horloge peut empêcher la récupération du signal utile et le positionnement. Une solution consiste à réduire le temps d'intégration cohérente de la corrélation mais cela augmente les pertes dues au bruit et donc aussi le temps nécessaire pour obtenir le positionnement; une autre solution éventuellement combinée à la précédente consiste à tester plusieurs hypothèses de dérives d'horloge en parallèle, mais cela alourdit la charge de calcul. De plus lorsque la dérive est très instable la tâche devient très difficile, voire impossible.

Un but important de l'invention est donc de compenser la dynamique apparente due à la dérive de l'horloge du récepteur, de manière à permettre d'utiliser un temps d'intégration de la corrélation suffisamment long pour pouvoir distinguer le signal du bruit, malgré des conditions de réception dégradées.

Pour atteindre ce but, l'invention propose un dispositif de positionnement, qui comporte une horloge locale présentant une dérive, et

un premier récepteur d'un signal de radionavigation, principalement caractérisé en ce qu'il comprend en outre un deuxième récepteur d'un signal radioélectrique provenant d'une balise fixe qui présente une horloge stable, l'horloge locale étant commune au premier et au deuxième récepteurs, et en ce qu'il comprend relié au premier récepteur, un élément de calcul de la dérive de l'horloge locale, à partir du signal radioélectrique provenant de la balise fixe.

Ce dispositif permet ainsi de compenser la dérive d'horloge d'un récepteur GNSS dans le but d'intégrer le signal reçu sur un temps long. On peut alors déterminer la position du récepteur dans un environnement difficile.

Le deuxième récepteur est par exemple un dispositif mobile de télécommunication.

Le signal de radionavigation provient d'un système de positionnement par satellites et/ou pseudolites.

L'invention concerne également un procédé de positionnement à partir d'un signal de radionavigation et d'un signal d'horloge local, caractérisé en ce qu'il comprend une étape de réception d'un signal provenant d'une balise fixe qui présente une horloge stable, et de calcul de la dérive du signal d'horloge local à partir du signal de la balise fixe.

D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de la description détaillée qui suit, faite à titre d'exemple non limitatif et en référence aux dessins annexés dans lesquels : la figure 1 déjà décrite, illustre schématiquement le principe de fonctionnement d'un récepteur GNSS, la figure 2 déjà décrite, illustre schématiquement le temps réellement mesuré, la figure 3 représente schématiquement des exemples de courbes du temps t vrai et du temps t _h (t) indiqué par l'horloge à l'instant t vrai, la figure 4 illustre schématiquement des exemples de résultats de la corrélation du signal reçu avec le signal local en fonction de plusieurs hypothèses de temps, τ étant le retard du code local représenté figure 1d par rapport au code reçu représenté figure 1 b,

la figure 5 représente schématiquement des exemples de variation sur le temps d'intégration, d'une part du temps t _reçu (t) auquel est émis par le satellite le code C(t _re çu(t)) reçu par le récepteur à l'instant t vrai et des différentes hypothèses de temps t _re çu estimé(t) reçu estimé dans le cas où l'on ne compense pas la dérive d'horloge, la figure 6 représente schématiquement un récepteur GNSS selon l'art antérieur, la figure 7 représente schématiquement des exemples de variation sur le temps d'intégration, d'une part du temps t _re çu(t) reçu vrai et des différentes hypothèses de temps t _re çu estimé(t) reçu estimé dans le cas où l'on compense la dérive d'horloge, la figure 8 représente schématiquement un dispositif de positionnement selon l'invention, la figure 9 représente schématiquement un récepteur non GNSS permettant de déterminer la dérive d'horloge locale, la figure 10 représente schématiquement les principaux éléments d'un dispositif de positionnement selon l'invention.

On va exprimer plus en détail les temps indiqués figures 1 et 2. On a : t _re çu(t) = t - Δt _prop (t), t _re çu(t) et t étant des temps vrais.

Le signal GNSS émis par le satellite est de la forme :

Sémis(t) = Code( témis(t) ) . sin( φ _é mise(t) ) = Code( t ) . sin( ω. t ) avec témis(t) = t φémise(t) = ω. t

Le signal GNSS reçu par le récepteur est de la forme: Sreçu(t) = Code( treçu(t) ) . sin( φ _reÇ ue(t) ) avec treçu(t) = t - Δt _prop (t) φreçue(t) = ω ( t - Δt _prop (t) )

Par ailleurs dans la mesure où l'horloge présente un biais Δt _h (t), on a : th(t) = t + Δth(t), th étant le temps indiqué par l'horloge du récepteur à l'instant t.

On a typiquement des dérives Δth(t) de 10 ^"5 seconde par seconde pour des horloges bas de gamme, soit une variation du biais d'horloge de 10 microsecondes sur un intervalle de 1 seconde, c'est-à-dire 10 fois la durée d'un chip de code, ce qui est rédhibitoire. On a représenté figure 3 les courbes illustrant t et th(t).

Mais on ne connaît pas t ni par conséquent Δt _h (t). On a en fait :

testimé(t) = t _h (t) - Δt _h estimé(t) ( 1 )

De même, on ne connaît pas t _re çu(t) ni Δt _pr op(t). On a en fait : treçu estimé(t) ⁼ t _e stimé(t) - Δtprop estimé(t) (2)

On a d'après (1 ) et (2), treçu estimé(t) = th(t) - Δth estimé(t) - Δtprop estimé(t) (3)

Lors de l'acquisition du signal par le récepteur, et notamment lors du calcul de corrélation du signal reçu avec le signal généré localement, le récepteur ne peut faire qu'une estimation grossière du temps : cette estimation t _re çu estimé(t) n'est pas assez précise pour faire suffisamment coïncider le signal local avec le signal reçu.

Pour cela le récepteur tente plusieurs corrélations avec différentes hypothèses de biais d'horloge estimés Δt _h estimé(t) et donc aussi de temps t _reÇ u estimé(t) décalées entre elle d'une quantité T telle qu'au moins une des hypothèses soit suffisamment proche du temps t _re çu(t) vrai pour que les signaux soient corrélés, comme illustré figure 4. Le résultat de la corrélation dépend de τ avec : t ⁼ treçu(t) - treçu estimé(t)

T = ( t - testimé(t) ) " ( Δtprop(t) - Δtprop estimé(t) )

L'erreur due à l'estimation du temps de propagation est considérée comme nulle car on suppose que l'on connaît la position du satellite avec une précision suffisante. On peut donc écrire : τ ≡ Δth(t) - Δt _h estimé(t)

Sur la figure 5, on a représenté les variations sur le temps d'intégration, d'une part du temps t _re çu(t) vrai et des différentes hypothèses

de treçu estimé(t) dans le cas où l'on ne compense pas la dérive d'horloge, c'est- à-dire où l'on suppose que Δt _h (t) et donc Δt _{h e} stimé(t) sont constants.

Les courbes de la figure 5 qui représentent t _re çu estimé(t) diffèrent les unes des autres en prenant comme différentes hypothèses pour Δt _h estimé(t), nT, (n+1 ) T, (n+2) T, etc.

En raison de la dérive de l'horloge, aucune hypothèse de temps reçu estimée ne coïncide avec le temps t _reçu (t) vrai sur la durée d'intégration. De ce fait le récepteur est incapable de déterminer quelle est l'hypothèse qui coïncide avec le signal reçu, même en augmentant le temps d'intégration. On a représenté figure 6 un exemple d'un tel récepteur 10. Il comprend de manière classique, un générateur 11 d'une porteuse locale de la forme e' ^ωt , un démodulateur 18 de porteuse, un corrélateur apte à multiplier le signal reçu démodulé par le code local issu d'un générateur 12 de code local et à l'intégrer au moyen d'un intégrateur par intervalles 13, un dispositif d'asservissement 14 pour l'acquisition et la poursuite des phases des signaux reçus élaborant des corrections en vitesse de la boucle de code et de la boucle de porteuse, une horloge locale 15 générant un signal puisé périodique de fréquence F _h , un NCO (acronyme de l'expression anglo- saxonne « Numerically Controlled Oscillator ») code 16 et un NCO porteuse 17 intégrant sur un intervalle T _in tégration ces corrections de vitesse et des corrections permettant de prendre en compte l'effet Doppler, pour fournir respectivement des estimées de phase (ou de temps) de code et de porteuse. On a représenté sur cette figure les éléments pour un canal de réception correspondant à un satellite. Il y a généralement dans un récepteur plusieurs canaux de réception en parallèle.

On rappelle qu'un NCO code réalise l'opération : treçu estimé(t) = th(t) - Δth estimé(t) + I[th(to), th(t)] V _CO de (x) dτ avec Vcode = - Vitesse radiale du satellite / c + Correction de code en vitesse (ou en fréquence) c : vitesse de la lumière.

Le terme « - Vitesse radiale du satellite / c » correspond à la prise en compte de l'effet Doppler. On a donc d'après (3) :

J[th(to), th(t)] Vcode (τ) dτ = - Δtprop estimé(t)

Selon une variante, le NCO réalise l'opération suivante : treçu estimé(t) = th(t) - Δth estimé(t) + J [th(to), th(t)] V _CO de(x) dτ

On a : Δth estimé(t) = Δth estimé(to) + I[th(to), m)] dh estimée(τ) dτ par définition

Et : th(t) = th(to) + I[th(to), th(t)] 1 dτ

On obtient finalement : treçu estimé(t) = th(to) - Δth estimé(tθ) + I[th(to), th(t)] (1 " dh estimée + V _CO de)(x)dτ

La mise en œuvre de cette variante est représentée figure 10. Un NCO porteuse réalise la même opération en multipliant en outre le résultat par la pulsation ω.

Selon l'invention le récepteur 10 de type GNSS est associé à un récepteur 20 d'un signal de communication non GNSS (désigné récepteur non GNSS) qui présente un rapport signal sur bruit élevé et qui est utilisé pour identifier la dérive de l'horloge locale 15 commune aux deux récepteurs et la compenser, comme représenté figure 8.

Le signal non GNSS est émis par une balise fixe 30 ayant une horloge 31 ou plus généralement une référence de temps stable, typiquement avec une dérive inférieure en valeur absolue à 10 ^"10 seconde par seconde.

On a représenté figure 9, un exemple de récepteur 20 non GNSS permettant de déterminer cette dérive d'horloge. Il comprend notamment une boucle à verrouillage de phase (« PLL » ou Phase Lock Loop en anglais) qui comporte de manière classique un démodulateur 22 de porteuse, un dispositif d'asservissement 24 incluant un filtre de boucle, un NCO intégrateur 27 relié à une horloge 15 et un générateur 21 de porteuse locale de la forme e' ^ωt . La sortie du NCO intégrateur de phase fournit la phase locale φ'iocaie-

La récupération du signal au niveau du récepteur non GNSS permet d'estimer la dérive de l'horloge locale de la manière suivante.

Le signal non GNSS émis démodulé est de la forme : Sémis(t) = sin( φ'émis(t) ) avec φ'émis(t) = ω'. t

Le signal non GNSS reçu : Sreçu(t) = sin( φ' _re çu(t) ) avec φ' _re çu(t) = ω'. ( t - Δt' _P ro _P (t) ) (4)

Δt'prop(t) étant le temps de propagation entre la balise fixe et le récepteur non GNSS, constant car le récepteur et la balise sont censés être immobiles.

Le signal non GNSS local :

Siocal I (t) = Sin( φ'iocale(t) ) Siocal û(t) = COS( φ'bcale(t) ) avec φ'iocale(t) = φ'reçue estimée(t) = ω'. ( t' _re çu estimé(t) ) (5)

Le récepteur non GNSS peut par exemple mettre en œuvre une boucle à verrouillage de phase (« PLL ») pour récupérer le signal de communication non GNSS.

Quand la PLL a convergé les phases reçue et locale coïncident, et on obtient : φ'iocaie = φ'reçue + n x 2π + ε ε « π/2 (erreur de phase) soit en vertu de (4) et (5) : t'reçu estimé(t) = t - Δt'prop + n X 2π/ ⁽ θ' + ε(t) /ω' (6)

Par ailleurs on a :

Δt _h (t) = t _h (t) - t soit en vertu de (6) :

Δth(t) = t _h (t) - t'reçu estimé(t) + Δt'prop - n X 2π/(θ' - ε(t)/(θ'

Si on néglige l'erreur de phase ε(t) et lorsque δt _pr op(t) est constant car le récepteur non GNSS ne bouge pas, on obtient finalement :

Δth(t) - Δth(to) ≡ ( t _h (t) - th(to) ) - (t' _re çu estimé(t) - t'reçu estimé(tθ))

Δth(t) - Δth(to) ≡ ( t _h (t) - th(to) ) - (φ'iocaie(t) /ω' - φ'iocaie(to) /ω').

On a :

Δth(t) - Δth(to) = I[to, t] d _h (τ) dτ par définition,

th(t) - th(to) = I _[th (to), th (t)] 1 dτ par définition, φ'ι∞aie(t) /ω' - φ'iocaie(to) /ω' = Wo), m)] (1 + Vp _L ι_(τ) ) dτ en sortie du NCO du récepteur non GNSS.

On peut donc écrire :

I[to, t] dh(τ) dτ ≡ J[th (to), th <t)] 1 dτ - Imto), th (t)] (1 + Vp _L ι_(τ) ) dτ

D'où : d _h (t) ≡ - VPLLW

Ainsi on est capable de déterminer la dérive d'horloge instantanée d _h (t), à partir de la correction de porteuse en vitesse de la PLL qui est envoyée au NCO de porteuse de la PLL.

Selon l'invention, cette horloge locale est commune aux deux récepteurs c'est-à-dire au récepteur GNSS et au récepteur non GNSS. La dérive d'horloge estimée en seconde/seconde est ajoutée aux corrections en vitesse (issues de l'asservissement et de la compensation de l'effet Doppler) des signaux locaux GNSS pour la corrélation, à savoir le code local et la porteuse locale. Sur la figure 7, on a représenté des résultats ainsi obtenus c'est-à- dire les variations sur le temps d'intégration d'une part du temps t _re çu(t) vrai et des différentes hypothèses dans le cas où l'on compense la dérive d'horloge. Dans ce cas lorsque les hypothèses de temps reçu sont suffisamment proches, il y en a une qui coïncide avec le temps reçu vrai, ce qui permet une détection d'énergie en sortie de la corrélation et donc l'estimation d'un temps reçu.

On a représenté figure 10 un exemple de dispositif de positionnement selon l'invention. Il comprend une horloge locale 15, un récepteur GNSS 10, un récepteur 20 d'un signal radioélectrique non GNSS provenant d'une balise fixe présentant une référence de temps stable, l'horloge locale 15 étant commune aux deux récepteurs.

Il comprend également, relié au récepteur GNSS, un élément de calcul de la dérive d'horloge à partir du signal radioélectrique non GNSS provenant de la balise fixe, le calcul étant par exemple tel que décrit précédemment. Cet élément est relié aux NCO (les NCO 16 et 17 sur la figure 10) des signaux locaux GNSS, à savoir le code local et la porteuse

locale, de manière à ajouter cette dérive d'horloge instantanée aux corrections existantes. Cet élément de calcul est par exemple un microprocesseur ou un circuit comprenant notamment un NCO pour effectuer l'intégration. Dans l'exemple de la figure 10, la dérive égale à -Vp _L ι_(t) est calculée par le dispositif d'asservissement 24 de la PLL du récepteur non GNSS.

Le récepteur GNSS 10 comprend comme dans le cas de la figure 6, une horloge locale 15, un générateur 1 1 d'une porteuse locale de la forme e' ^ωt , un démodulateur 18 de porteuse, un corrélateur apte à multiplier le signal reçu par le code local issu d'un générateur 12 de code local et l'intégrer au moyen d'un intégrateur par intervalle 13, un dispositif d'asservissement 14 pour l'acquisition et la poursuite des phases des signaux reçus élaborant des corrections en vitesse de la boucle de code et de la boucle de porteuse, et un NCO code 16 et un NCO porteuse 17 intégrant ces corrections de vitesse et les corrections permettant de prendre en compte l'effet Doppler et la dérive d'horloge, pour fournir des estimées de phase (ou de temps) de code et de porteuse.

Le récepteur non GNSS 20 comprend notamment comme dans le cas de la figure 9, une PLL qui comporte de manière classique un démodulateur 22 de porteuse, un dispositif d'asservissement 24 incluant un filtre de boucle, un NCO intégrateur 27 et un générateur 21 de porteuse locale de la forme e' ^{ω t} . La sortie du NCO intégrateur de phase 27 fournit la mesure de la phase du signal reçu φ' _re çue estimée = φ'i∞aie-

Cet élément de calcul de la dérive d'horloge (le dispositif 24 sur la figure) est relié au correcteur (le NCO 27 sur la figure) de la porteuse du récepteur non GNSS, dans le cas où il est intéressant de corriger la porteuse locale avec cette dérive, comme représenté sur les figures. Il peut être intégré dans le récepteur non GNSS comme dans l'exemple de la figure 10 ou dans le récepteur GNSS. Le récepteur non GNSS est par exemple un dispositif de télécommunication mobile comportant éventuellement un émetteur.

Le signal dit GNSS peut provenir d'un système de positionnement par satellites et/ou pseudolites.