FR 3000806 A1 20140711 FR 1350211 A 20130110 Domaine technique [0001] La présente invention concerne le domaine de la radionavigation, plus précisément un dispositif de poursuite d'un signal de radionavigation (ou signal radioélectrique de navigation). Le champ d'application de l'invention est notamment la réception de signaux de radionavigation émis par des émetteurs d'un système de positionnement par satellite (abrégé GNSS, pour « Global Navigation Satellite System »), p. ex. GPS (de l'anglais « Global Positioning System), Galileo, Glonass, QZSS, Compass, IRNSS, etc. Etat de la technique [0002] De manière générale, les signaux de radionavigation émis par des satellites (ou pseudolites) d'un système de positionnement 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 pseudo-alé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). Accessoirement, certains signaux de positionnement par satellites peuvent aussi transporter 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. Dans ce qui suit, on fera abstraction de cette charge de données utiles. [0003] Dans le cas de GPS, les signaux de radionavigation sont transmis dans les bandes de fréquences L1, centrée sur 1575,42 MHz, L2, centrée sur 1227,6 MHz, et L5, centrée sur 1176,45 MHz. Les satellites du GNSS européen (aussi connu sous le nom de Galileo) transmettront dans les bandes E2-L1-E1 (la portion de bande médiane L1 étant la même que celle du GPS), E5a (qui, selon la nomenclature Galileo, représente la bande L5 prévue pour le GPS), E5b (centrée sur 1207,14 MHz) et E6 (centrée sur 1278,75 MHz). Dans ce qui suit, les bandes E5a et E5b seront traitées ensemble comme bande E5, avec 1191,795 MHz comme fréquence centrale. Pour le cas des signaux ouverts de Galileo, leur description complète peut être trouvée dans le « Galileo Open Service Signal-In-Space Interface Control Document », ou Galileo OS SIS ICD, disponible sur le site web http://ec.europa. eu/enterprise/policies/satnav/galileo/open-service/index_en.htm. On peut encore noter que les satellites de la constellation Compass transmettent ou transmettront dans la bande B1 (centrée sur 1561,098 MHz), B1-2 (centrée sur 1589,742 MHz), L1 (centrée sur 1575,42 MHz), B2 (centrée sur 1207,14 MHz), et B3 (centrée sur 1268,52 MHz). Le système GLONASS utilise les fréquences centrales 1602 MHz et 1246 MHz. Les fréquences centrales indiquées représentent les fréquences des porteuses des différents signaux. [0004] La réception d'un signal de radionavigation comprend normalement une multiplication du signal reçu par une réplique interne de la porteuse générée dans le récepteur par un oscillateur piloté par une boucle de poursuite de porteuse et une autre multiplication par une réplique interne de la forme d'onde d'étalement produite par un générateur de forme d'onde piloté par une boucle de poursuite de forme d'onde d'étalement (appelée aussi « boucle de poursuite de code »). Les signaux d'erreur ou d'asservissement des boucles de poursuite de porteuse et de forme d'onde d'étalement sont utilisés par le récepteur pour déterminer sa position. Le signal d'écart de phase entre la porteuse du signal reçu et la réplique interne de porteuse produit à chaque pas de temps par la boucle de poursuite de porteuse fournit une première mesure ou observable (la mesure ou l'observable de phase). Le signal de décalage temporel entre la forme d'onde d'étalement du signal reçu et la réplique interne de forme d'onde d'étalement produit à chaque pas de temps par la boucle de poursuite de forme d'onde d'étalement représente une deuxième mesure ou observable (la mesure ou l'observable de code). [0005] L'observable de phase n'est pas utilisée par tous les récepteurs. Surtout les récepteurs bon marché déterminent leur position uniquement sur la base des observations de code. Les mesures de phase sont mises en pratique, par exemple, dans les méthodes RTK (de l'anglais : « Real Time Kinematic ») et PPP (de l'anglais : « Precise Point Positioning »). [0006] Les mesures de code sont d'une précision métrique tandis que les mesures de phase sont d'une précision de quelques mm. Toutefois, les mesures de phase présentent le grand inconvénient qu'elles sont ambiguës dans le sens que le nombre de cycles entiers entre le satellite et le récepteur est inconnu au départ. Les mesures de phase sont modulo un cycle et ne livrent que la partie réelle de la différence de phase de la porteuse entre l'émission par le satellite et le récepteur. Afin de pouvoir bénéficier de la précision des mesures de phase, un récepteur doit résoudre les ambiguïtés dont elles sont entachées. [0007] La poursuite de la phase de la porteuse d'un signal de radionavigation est très sensible aux conditions de l'environnement. Le risque d'un décrochage est beaucoup plus élevé que pour la poursuite du code. De plus, la gestion des sauts de cycle est une tâche difficile. Dans un environnement difficile (par exemple en milieu urbain), la disponibilité des mesures de phase risque d'être très faible, ce qui réduit fortement l'intérêt des récepteurs capables d'effectuer et de traiter des observations de phase. Objet de l'invention [0008] Un objet de la présente invention est d'augmenter la robustesse des mesures de phase. Description générale de l'invention [0009] Selon l'invention, il est proposé un dispositif de poursuite de signaux de radionavigation, comprenant un premier étage de poursuite d'un premier signal de radionavigation contenu dans un signal entrant à appliquer au dispositif et un deuxième étage de poursuite d'un deuxième signal de radionavigation contenu dans le signal entrant. Le premier signal de radionavigation est supposé comprendre une première porteuse à une première fréquence modulée par une première forme d'onde d'étalement, tandis que le deuxième signal de radionavigation est supposé comprendre une deuxième porteuse à une deuxième fréquence, différente de la première fréquence, modulée par une deuxième forme d'onde d'étalement. Le premier étage de poursuite comprend une première boucle à verrouillage de phase de porteuse avec un mélangeur configuré pour multiplier le signal entrant avec une réplique locale de la première porteuse. La première boucle à verrouillage de phase de porteuse comprend un discriminateur de phase de la première porteuse configuré pour produire un premier signal d'erreur résultant d'un premier écart de phase entre la première porteuse et la réplique locale de la première porteuse. La première boucle à verrouillage de phase de porteuse est configurée pour régler la phase de la réplique locale de la première porteuse sur la base du premier signal d'erreur. Le deuxième étage de poursuite comprend une deuxième boucle à verrouillage de phase de porteuse avec un mélangeur configuré pour multiplier le signal entrant avec une réplique locale de la deuxième porteuse. La deuxième boucle à verrouillage de phase de porteuse comprend un discriminateur de la phase du battement entre la première porteuse et la deuxième porteuse configuré pour produire un deuxième signal d'erreur résultant d'une différence entre le premier écart de phase et un deuxième écart de phase entre la deuxième porteuse et la réplique locale de la deuxième porteuse. La deuxième boucle à verrouillage de phase est configurée pour régler la phase de la réplique locale de la deuxième porteuse sur la base des premier et deuxième signaux d'erreur. [0010] On note que le dispositif selon l'invention peut être employé dans le cadre d'un récepteur GNSS bi-fréquence ou multifréquence. Au lieu de poursuivre les premier et deuxième signaux de radionavigation de manière individuelle (c.-à-d. par des boucles de poursuite de phase séparées), le dispositif de poursuite de l'invention utilise des boucles à verrouillage de phase couplées entre elles. La phase du premier signal de radionavigation est poursuivie sur la base du premier signal d'erreur tandis que celle du deuxième signal de radionavigation est poursuivie grâce à l'utilisation du premier signal d'erreur et du deuxième signal d'erreur qui dépend de la phase du battement entre les deux porteuses. Dans ce document, le terme « battement » désigne la variation périodique de l'amplitude de l'oscillation résultant de la composition de deux porteuses. Le signal de battement a la fréquence Ifi-f21 où f1 et f2 représentent respectivement la première et la deuxième fréquence de porteuse. La phase du battement dépend de la phase de chacune des deux porteuses. Intuitivement, la connaissance de la phase d'une des porteuses et de la phase du battement permet de retrouver la phase de l'autre porteuse. Dans le dispositif de poursuite selon l'invention, la boucle de poursuite de phase de la deuxième porteuse règle la phase de la réplique locale de la deuxième porteuse en utilisant cette propriété. Un avantage de cette approche est que le dispositif selon l'invention ne poursuit directement que la phase de la première porteuse. La deuxième boucle à verrouillage de phase poursuit le battement, qui se trouve à une fréquence (Ifi-f21) nettement inférieure à la fréquence de la deuxième porteuse. Cette fréquence étant plus basse, le nombre de sauts de phase par unité de temps est considérablement réduit par rapport aux porteuses, ce qui entraîne une robustesse élevée de la mesure de phase de battement. La méthode utilisée par le dispositif selon l'invention permet ainsi une poursuite conjointe, plus robuste, de la phase de deux signaux de radionavigation en provenance du même émetteur (satellite ou pseudolite). [0011] De préférence, la première boucle à verrouillage de phase comprend un oscillateur à commande numérique (OCN), commandé par un premier filtre recevant le premier signal d'erreur en entrée. La deuxième boucle à verrouillage de phase comprend de préférence un oscillateur à commande numérique commandé par le premier filtre et un deuxième filtre (de la deuxième boucle de verrouillage de phase) recevant le deuxième signal d'erreur en entrée. Le deuxième filtre peut comprendre, p.ex. un filtre de Kalman, un filtre de Kalman étendu ou un filtre particulaire. [0012] De préférence, les premier et deuxième signaux de radionavigation sont choisis parmi les signaux Galileo El, E5 et E6 ou parmi les signaux GPS L5, L2C et Li ou encore parmi les signaux GLONASS L3, G2 et Gl. La méthode utilisée par le dispositif selon l'invention est particulièrement avantageuse si la poursuite de phase de porteuse de l'un des signaux est significativement plus robuste que la poursuite du ou des autres signaux. Dans une telle situation, le signal dont la poursuite de phase de porteuse est la plus robuste peut être choisi comme le premier signal de radionavigation. La phase de porteuse des autres signaux de radionavigation n'a pas besoin d'être poursuivie directement mais peut être trouvée grâce à la phase de battement. Selon un mode de réalisation avantageux de l'invention, le premier signal de radionavigation est un signal Galileo E5 et le deuxième signal de radionavigation est un signal Galileo El ou un signal Galileo E6. A cause de la modulation « AltBOC » du signal Galileo E5, la poursuite de la phase de porteuse de ce signal est plus robuste que celle de El (modulation CBOC) ou de E6. [0013] Le premier étage de poursuite comprend de préférence un premier corrélateur configuré pour corréler une réplique locale de la première forme d'onde d'étalement avec le signal entrant multiplié par la réplique locale de la première porteuse. Le discriminateur de phase de la première porteuse peut être configuré de sorte à déterminer le premier signal d'erreur sur la base du résultat de corrélation produit par le premier corrélateur. [0014] Le deuxième étage de poursuite comprend de préférence un deuxième corrélateur configuré pour corréler une réplique locale de la deuxième forme d'onde d'étalement avec le signal entrant multiplié par la réplique locale de la deuxième 3000 806 6 porteuse. Le discriminateur de phase du battement entre la première porteuse et la deuxième porteuse peut être configuré pour déterminer le deuxième signal d'erreur sur la base du résultat de corrélation produit par le premier corrélateur et du résultat de corrélation produit par le deuxième corrélateur. 5 [0015] De préférence, le discriminateur de phase du battement entre la première porteuse et la deuxième porteuse est configuré de sorte qu'il utilise soit le produit du résultat de corrélation fourni par le premier corrélateur et du conjugué complexe du résultat de corrélation fourni par le deuxième corrélateur, soit le produit du conjugué complexe du résultat de corrélation produit par le premier corrélateur et du résultat 10 de corrélation produit par le deuxième corrélateur pour produire le deuxième signal d'erreur. [0016] Le dispositif selon l'invention peut être configuré pour le traitement conjoint de trois signaux de radionavigation voire plus. Dans ce cas, un troisième étage de poursuite d'un troisième signal de radionavigation est à prévoir. Le troisième signal 15 de radionavigation est supposé comprendre une troisième porteuse à une troisième fréquence, différente des première et deuxième fréquences, modulée par une troisième forme d'onde d'étalement. Le troisième étage de poursuite est de préférence structuré de la même manière que le deuxième étage de poursuite. Il peut notamment comprendre une troisième boucle à verrouillage de phase de porteuse 20 avec un mélangeur configuré pour multiplier le signal entrant avec une réplique locale de la troisième porteuse. La troisième boucle à verrouillage de phase de porteuse peut comprendre un discriminateur de phase d'un battement entre la première porteuse et la troisième porteuse configuré pour produire un troisième signal d'erreur résultant d'une différence entre le premier écart de phase et un 25 troisième écart de phase entre la troisième porteuse et la réplique locale de la troisième porteuse. La troisième boucle à verrouillage de phase est alors configurée pour régler la phase de la réplique locale de la troisième porteuse sur la base des premier et troisième signaux d'erreur. De préférence, les premier, deuxième et troisième signaux de radionavigation sont choisis parmi les signaux Galileo El, E5 et 30 E6 ou parmi les signaux GPS L5, L2C et Li ou encore parmi les signaux GLONASS L3, G2 et Gl. [0017] La troisième boucle à verrouillage de phase comprend de préférence un oscillateur à commande numérique commandé par le premier filtre et un troisième filtre (de la troisième boucle de verrouillage de phase) recevant le troisième signal d'erreur en entrée. Autrement, l'oscillateur à commande numérique peut être commandé par le même filtre de Kalman, le filtre de Kalman étendu ou le filtre particulaire qui pilote l'oscillateur à commande numérique de la deuxième boucle. [0018] Le troisième étage de poursuite comprend de préférence un troisième corrélateur configuré pour corréler une réplique locale de la troisième forme d'onde d'étalement avec le signal entrant multiplié par la réplique locale de la troisième porteuse. Le discriminateur de phase du battement entre la première porteuse et la troisième porteuse peut alors déterminer le troisième signal d'erreur sur la base du résultat de corrélation produit par le premier corrélateur et du résultat de corrélation produit par le troisième corrélateur. [0019] Le discriminateur de phase du battement entre la première porteuse et la troisième porteuse est de préférence configuré de sorte à utiliser le produit du résultat de corrélation produit par le premier corrélateur et du conjugué complexe du résultat de corrélation produit par le troisième corrélateur ou le produit du conjugué complexe du résultat de corrélation produit par le premier corrélateur et du résultat de corrélation produit par le troisième corrélateur pour produire le troisième signal d'erreur. [0020] Il sera apprécié que le dispositif selon l'invention puisse être réalisé par un processeur de signaux numériques (DSP, de l'anglais « Digital Signal processor »). [0021] Un aspect de l'invention concerne un récepteur GNSS, comprenant un ou plusieurs dispositifs de poursuite de signaux de radionavigation tels que décrits. Brève description des dessins [0022] D'autres particularités et caractéristiques de l'invention ressortiront de la description détaillée de quelques modes de réalisation avantageux présentés ci- dessous, à titre d'illustration, avec référence aux dessins annexés. Ceux-ci montrent: Fig. 1: un schéma de principe d'un dispositif de poursuite selon un premier mode de réalisation avantageux de l'invention ; Fig. 2: un schéma de principe d'un dispositif de poursuite selon un deuxième mode de réalisation avantageux de l'invention ; Fig. 3: un schéma de principe d'une variante du dispositif de poursuite de la figure 2; Fig. 4: un schéma de principe d'un récepteur GNSS comprenant des dispositifs de poursuite selon l'invention. Description d'une exécution préférée [0023] Un récepteur de signaux de radionavigation peut réaliser des mesures de code (non ambiguës) et des mesures de phase (ambiguës d'un nombre entier de cycles) sur les signaux de radionavigation qu'il reçoit des différents satellites en visibilité (c.-à-d. au-dessus de l'horizon). Un récepteur multi-fréquence peut effectuer ces mesures sur au moins deux fréquences de porteuse f1 et f2 distinctes. En supposant un récepteur bi-fréquence, on peut donc avoir, pour chaque satellite en visibilité et à chaque pas de temps (tk), deux mesures de code, notées Pii(tk) et 132(tk), et deux mesures de phase, notées (1)ii(tk) et (I)2i(tk), sur les fréquences f1 et f2, où l'indice supérieur (j) indique l'émetteur (le satellite ou le pseudolite) dont le signal provient au récepteur. Comme les satellites orbitent autour de la Terre, seulement une partie d'entre eux est visible à un moment donné depuis le lieu du récepteur. Pour faciliter la notation, la dépendance du temps et l'indice de satellite ne seront pas toujours indiqués explicitement dans ce qui suit. [0024] On utilisera les notations suivantes : où c représente la vitesse de la lumière. [0025] Les mesures de code et de phase remplissent les équations (mesures à gauche, paramètres du modèle à droite) suivantes : P1=r+Sr+cAh P2 = r + (27 n + c Ah Àlg: ) = r + - + c Ah + ÀlAT (Eq. 1) 4(1)2 = r + À2W - + c Ah + 2N 2 où - r représente la distance de propagation entre le satellite et le récepteur, incluant l'élongation troposphérique, effets relativistes, etc. ; - Air et Sr° décrivent l'allongement ionosphérique qui varie approximativement en fonction du carré de la fréquence (d'où A17° sr° -f12 /122) ; - N1 et N2 représentent les ambiguïtés de phase (entières) des deux porteuses ; - W représente la contribution de l'effet « wind-up » en cycles ; et - Ah est la différence entre l'horloge du récepteur et l'horloge de l'émetteur. [0026] Dans le système d'équations ci-dessus, les mesures de code Pt P2 sont exprimées en unités de longueur, tandis que les mesures de phase (Pi, 02 sont exprimées en cycles. [0027] La combinaison widelane est définie comme suit : ', =4)1 -(1)2 ="1.-+Ni,' +LM - (147 - -1Air) (Eq. 2) 112 où - fw est la fréquence de battement, définie par L =11-12 et - Nw est l'ambiguïté widelane définie par Arw = NI -N2. [0028] On note =c/U1 - f2) = (1/ Ài -1/ À2)-1 la longueur d'onde du battement. La longueur d'onde du battement entre les fréquences El et E5 de Galileo s'élève, p.ex., à peu près à 78 cm. [0029] La figure 1 montre un dispositif de poursuite 10 selon un mode de réalisation préféré de l'invention. Le dispositif de poursuite 10 comprend un premier étage de poursuite 12 pour poursuivre un signal Galileo E5 contenu dans le signal entrant et un deuxième étage de poursuite 14 pour poursuivre un signal Galileo El (ou E6) contenu dans le signal entrant. Pour la suite, à titre d'exemple, on utilisera f1 et f2 pour désigner la fréquence de la porteuse de E5 et de El, respectivement. [0030] Le signal entrant comprend les signaux de radionavigation à la fréquence f1 et à la fréquence f2 transmis par les différents satellites en visibilité. Les signaux en provenance du même émetteur sont traités conjointement. Par contre, il n'est pas prévu de coupler la poursuite de signaux provenant d'émetteurs différents. Un récepteur GNSS comprendra donc un nombre de dispositifs de poursuite 10 correspondant au nombre de satellites (ou de pseudolites) que le récepteur doit être capable de poursuivre en parallèle. [0031] Le premier étage de poursuite 12 comprend un mélangeur 16 multipliant le signal entrant avec une réplique locale de la porteuse E5, générée par un oscillateur à commande numérique 18. En notation complexe, la réplique locale de la porteuse E5 s'écrit exp(27clIt+ 01), où est l'estimée de la fréquence (variable car soumise à l'effet Doppler) et (oi l'estimée de valeur réelle de la phase 91 de la porteuse (en radians). En sortie du mélangeur 16, on dispose d'une voie « en phase », notée I, et d'une voie « en quadrature de phase », notée Q, qui sont amenées à un étage de corrélation 24. [0032] A l'étage de corrélation 24 les voies I et Q sont chacune multipliées par une réplique locale de la forme d'onde d'étalement et intégrées. Les générateurs de réplique de forme d'onde d'étalement sont représentées sur la figure 1 comme des générateurs de réplique de code pseudo-aléatoire 26 afin de ne pas trop encombrer le dessin. Si le signal à recevoir contient une sous-porteuse (p.ex. une modulation BOC (de l'anglais « Binary Offset Carrier »), AltBOC (de l'anglais (( Alternative BOC ») ou autre), l'étage de corrélation 24 présente une plus grande complexité (voir p.ex. le brevet FR 2 941 535). La figure 1 est par ailleurs simplifiée en ce qu'elle ne montre que le corrélateur dit « prompt ». L'étage de corrélation 24 peut comprendre d'autres corrélateurs. L'architecture la plus performante comprend un corrélateur « prompt », un corrélateur « en avance » (en anglais : « early correlator ») et un corrélateur « en retard » (en anglais « late correlator »). [0033] La réplique du code pseudo-aléatoire qui est mélangée avec les voies I et Q par les mélangeurs 28 est notée Ci(t - fi) où fi est l'estimée de la valeur réelle du décalage de code Ti et C1 la fonction pouvant prendre les valeurs +1 et -1 décrivant le code pseudo-aléatoire. On a Pi = cf, en cas d'accrochage de la boucle de poursuite de code. [0034] La sortie des corrélateurs 30 est utilisée par une boucle à verrouillage de phase et de fréquence. Un discriminateur de phase 32 génère un signal d'erreur indiquant l'erreur de phase résiduelle, c.-à-d. la différence entre l'estimée (oi de la phase et de la valeur réelle 91 . Le signal d'erreur est noté (5(01 = - Le discriminateur de phase peut être un discriminateur « atan2 »; dans ce cas, (501(tk) = atan2(Qp1(tk),/pi(tk)), où Qpi(tk) est le résultat de la corrélation sur la voie prompte (index « P ») en quadrature de phase et Ipi(tk) est le résultat de la corrélation sur la voie prompte en phase. Un discriminateur de fréquence (non illustré sur la figure 1) produit un signal d'erreur indiquant l'erreur de fréquence résiduelle, c.-à-d. la différence (Sli Le discriminateur de fréquence peut être un discriminateur « cross-product » (produit croisé) calculant _ Pl(t k-1)QP1(t k) Pl(t k)QP1(t k-1) où T = tk-tk_i. Un filtre de boucle 34 reçoit les signaux d'erreur de phase et de fréquence et en déduit le signal de commande de l'oscillateur à commande numérique 18. [0035] La sortie des corrélateurs 30 est en outre utilisée par une boucle à verrouillage de code (ou boucle de poursuite de code). Un discriminateur dit « de code » 36 génère un signal d'erreur indiquant l'erreur de décalage de code résiduelle, c.-à-d. la différence (Sfi = - fi. Ce discriminateur peut être, p. ex., un discriminateur « early-late normalisé », qui calcule : El(t k)2 - QE1(t k)2 )1 1 Ll(t k)2 - QL1(t k)2 (5î-I -(ST où (ST -(1-A;L), avec AEL le )I/Ei(tk )2 ± - QE1(tk )2 ± ILl(t k)2 - Li@ k)2 décalage temporel entre les corrélateurs en avance et les corrélateurs en retard. lEt QE1 sont les sorties des corrélateurs en avance (« early ») et lu, Qu sont les sorties des corrélateurs en retard (« late »). Le signal d'erreur de décalage de code est amené à un filtre de boucle de code 38 qui en déduit le signal de commande des générateurs de la réplique de la forme d'onde d'étalement (ici représentés par les générateurs de code 26). [0036] Le deuxième étage de poursuite 14 comprend un mélangeur 40 multipliant le signal entrant avec une réplique locale de la porteuse El, générée par un oscillateur à commande numérique 42. En notation complexe, la réplique locale de la porteuse El s'écrit exp(27c12t + (02) , où 12 est l'estimée de la fréquence et (02 l'estimée de valeur réelle de la phase 92 de la porteuse El (en radians). La sortie « en phase », notée I, et la sortie « en quadrature de phase », notée Q, du mélangeur sont amenées à un étage de corrélation 44. [0037] A l'étage de corrélation 44, des mélangeurs 48 multiplient chacune des voies I et Q avec une réplique locale de la forme d'onde d'étalement et des corrélateurs 50 intègrent les voies I et Q. Les générateurs de réplique de forme d'onde d'étalement sont représentés sur la figure 1 comme des générateurs de réplique de code pseudoaléatoire 46. Si le signal à recevoir contient une sous-porteuse (p.ex. une modulation BOC (de l'anglais « Binary Offset Carrier »), CBOC (de l'anglais « Composite BOC ») ou autre), l'étage de corrélation 44 présente une plus grande complexité (voir p.ex. le brevet FR 2 941 535). La figure 1 est par ailleurs simplifiée en ce qu'elle ne montre que le corrélateur dit « prompt ». L'étage de corrélation 44 peut comprendre d'autres corrélateurs, notamment un corrélateur « en avance » et un corrélateur « en retard ». [0038] La sortie des corrélateurs 50 sert d'entrée à une boucle de poursuite de code. Un discriminateur de code 52 génère un signal d'erreur indiquant l'erreur de décalage de code résiduelle, c.-à-d. la différence (5f2 = T2 - /T2 . Ce discriminateur peut être, p. ex., un discriminateur « early-late normalisé ». Le signal d'erreur de décalage de code est amené à un filtre de boucle de code 54 qui en déduit le signal de commande des générateurs de la réplique de la forme d'onde d'étalement (ici représentés par les générateurs de code 46). Les boucles de poursuite de code des signaux E5 et El sont donc construites de manière analogue dans cet exemple. Il faut toutefois noter que le choix des différents discriminateurs doit être fait en fonction de la modulation des signaux de radionavigation et d'une éventuelle charge utile (message de navigation). Les différents critères pour choisir des discriminateurs appropriés sont bien connus dans le domaine technique et n'ont dont pas besoin d'être répétés dans ce contexte. [0039] La boucle à verrouillage de phase et de fréquence pour poursuivre le deuxième signal de radionavigation combine la sortie des corrélateurs prompts du premier étage de poursuite 12 avec la sortie des corrélateurs prompts du deuxième étage de poursuite 14. La sortie des corrélateurs prompts du premier étage de poursuite 12 s'écrit, en notation complexe, /pi (tk ) ./Qpi (tk ) ce qui peut être représenté par A1 exp( j(91 - 01)) = A1 exp( j - (SOI) . De manière analogue, la sortie des corrélateurs prompts du deuxième étage de poursuite 14 s'écrit, en notation complexe, /p ( i0 2 k ( ce qui peut être représenté par A2 exp(j(92 -02)) =A2 exp(j - (502) . Le module de conjugaison complexe 58 retourne exp(- j - (SOI) et le module de multiplication complexe 56 calcule donc A1A2 exp( j - 001 -(502). La grandeur (SOI -(502, qui sera notée (50,4 , correspond à la phase du battement entre la première porteuse et la deuxième porteuse. [0040] Le discriminateur de phase 60 génère un signal d'erreur résultant d'une différence entre le premier écart de phase (SOI et le deuxième écart de phase (502. Le discriminateur de phase 60, qui retourne (50 en tant que signal d'erreur, peut être un discriminateur de type « atan2 ». [0041] Un discriminateur de fréquence (non illustré sur la figure 1) produit un signal d'erreur indiquant l'erreur de la fréquence de battement résiduelle, c.-à-d. la différence (5'=f - I' . Le discriminateur de fréquence peut être un discriminateur « cross-product » (produit croisé). Un filtre de boucle 62 reçoit les signaux d'erreur de 15 phase et de fréquence de battement et en déduit un signal de correction de phase et de fréquence à combiner avec le signal de commande fourni par le filtre de boucle 34. Le signal qui commande l'oscillateur 42 est donc basé sur l'estimée de la deuxième fréquence 12, obtenue par 12 , et la phase (02, obtenue par 01+ 0' 20 [0042] Les sorties des discriminateurs 32, 36, 52 et 60 peuvent être utilisées pour la détermination de la position par un module de positionnement (non illustré sur la figure 1). [0043] La figure 2 montre un dispositif de poursuite 110 selon un deuxième mode de réalisation préféré de l'invention. Le dispositif de poursuite 110 comprend un premier 25 étage de poursuite 112 pour poursuivre un signal Galileo E5, un deuxième étage de poursuite 114 pour poursuivre un signal Galileo El et un troisième étage de poursuite 114' pour poursuivre un signal Galileo E6. Pour la suite, on notera f1, f2 et f3 les fréquences des porteuses E5, El et E6, respectivement. [0044] Le signal entrant comprend les signaux de radionavigation aux fréquences fl, 30 f2 et f3 transmis par les différents satellites en visibilité. Les signaux en provenance du même émetteur sont traités conjointement. Par contre, il n'est pas prévu de coupler la poursuite de signaux provenant d'émetteurs différents. Un récepteur GNSS comprendra donc un nombre de dispositifs de poursuite 110 correspondant au nombre de satellites (ou de pseudolites) que le récepteur doit être capable de poursuivre en parallèle. [0045] Le premier étage de poursuite 112 présente essentiellement la même architecture que le premier étage de poursuite du mode de réalisation de la figure 1. Dans un mélangeur 116, le signal entrant est multiplié avec une réplique locale de la porteuse E5, notée exp(27clIt + 01) , générée par un oscillateur à commande numérique 118. En sortie du mélangeur 116, on dispose d'une voie « en phase », notée I, et d'une voie « en quadrature de phase », notée Q, qui sont amenées au premier étage de corrélation 124. [0046] Dans le premier étage de corrélation, les voies I et Q sont mélangées avec des répliques de la forme d'onde d'étalement et intégrées. Les sorties des corrélateurs sont amenées à un discriminateur de code, qui génère un signal d'erreur indiquant l'erreur de décalage de code résiduelle, c.-à-d. la différence (S fi = - i. Un filtre de boucle de poursuite de code (non illustré) en déduit l'estimée fi du décalage de code et l'utilise pour commander un oscillateur à commande numérique, qui à son tour fournit le signal d'horloge qui pilote les générateurs de la réplique de forme d'onde d'étalement. [0047] La sortie des corrélateurs 130 (A1 exp( j(91 - 01)) =4 exp( j - (SOI) en notation complexe) est également utilisée par un discriminateur de phase 132, qui génère un signal d'erreur indiquant l'erreur de phase résiduelle, c.-à-d. la différence entre l'estimée (oi de la phase et de la valeur réelle 91 . Un discriminateur de fréquence (non illustré sur la figure 2) produit un signal d'erreur indiquant l'erreur de fréquence résiduelle, c.-à-d. la différence (S:fi = Un filtre de boucle 134 reçoit les signaux d'erreur de phase et de fréquence et en déduit le signal de commande de l'oscillateur à commande numérique 118. Optionnellement, un signal d'aide du navigateur (dépendant de la vitesse et ou de l'accélération du récepteur) peut être appliqué à la 30 boucle de poursuite de porteuse pour augmenter la robustesse de l'accrochage. [0048] Dans le deuxième étage de poursuite 114, le signal entrant est mélangé dans un mélangeur 140 avec une réplique locale de la porteuse El, notée exp(27c12t +02). Les voies I et Q ainsi obtenues sont ensuite corrélées avec une réplique de la forme d'onde d'étalement. La sortie des corrélateurs 150 (A2 exp( j(92 - 02)) = A2 exp( j - (502) en notation complexe) est combinée avec la sortie des corrélateurs 130 par un discriminateur 160 de la phase du battement entre la première et la deuxième porteuses. Le discriminateur de phase 160 génère un signal d'erreur résultant d'une différence entre le premier écart de phase (SOI et le deuxième écart de phase (502. Le discriminateur de phase 160 retourne (50,442 = (501 - (502 en tant que signal d'erreur. Le deuxième étage de poursuite 114 comprend en outre un discriminateur de fréquence (non illustré) produisant un signal d'erreur indiquant l'erreur de la fréquence de battement entre la première et la troisième porteuses, c.-à-d. la différence ew12 - L12 -L12 'avec f12 =f -f2. [0049] Le troisième étage de poursuite 114' est construit de manière analogue au deuxième étage de poursuite 114. Dans le troisième étage de poursuite 114', le signal entrant est mélangé dans un mélangeur 140' avec une réplique locale de la porteuse E6, notée exp(27c13t + 03). Les voies I et Q ainsi obtenues sont ensuite corrélées avec une réplique de la forme d'onde d'étalement. La sortie des corrélateurs 150' (A3 exp( j(93 -03)) = A exp( j - (503) en notation complexe) est combinée avec la sortie des corrélateurs 130 par un discriminateur 160' de la phase du battement entre la première et la troisième porteuses. Le discriminateur de phase 160' génère un signal d'erreur résultant d'une différence entre le premier écart de phase (SOI et le troisième écart de phase (503. Le discriminateur de phase 160' retourne (50,443 =(501 -(503 en tant que signal d'erreur. Le troisième étage de poursuite 114' comprend en outre un discriminateur de fréquence (non illustré) produisant un signal d'erreur indiquant l'erreur de la fréquence de battement entre la première et la troisième porteuses, c.-à-d. la différence ew13 L13 -j>14,13 avec fii3 = f3 . [0050] La sortie des corrélateurs 150' (A3 exp( j(93 -03)) = A3 exp( j - (503) en notation complexe) est en outre combinée avec la sortie des corrélateurs 150 par un discriminateur 160" de la phase du battement entre la deuxième et la troisième porteuses. Le discriminateur de phase 160" génère un signal d'erreur résultant d'une différence entre le deuxième écart de phase (502et le troisième écart de phase (503. Le discriminateur de phase 160" retourne donc (50,4,23 = (502 -5ç en tant que signal d'erreur. Un discriminateur de fréquence (non illustré) produisant un signal d'erreur indiquant l'erreur de la fréquence de battement entre la deuxième et la troisième porteuses, c.-à-d. la différence (5l,v23 =L23 -:f,v23, avec L23 = f2 - f3 , peut également être prévu. [0051] Les sorties des discriminateurs 136, 160, 160', 160" ainsi que des discriminateurs de fréquence (non illustrés) sont traitées dans un filtre de Kalman étendu 120. Le vecteur de mesures appliqué au filtre de Kalman étendu 120 peut donc s'écrire : z = [c(Sî-I, (5, 1442 / (5114,12 / (5014,13 / (5.71443 / (5014,23 / (5.714,23 ' Comme sortie du filtre de Kalman étendu, on a le vecteur d'état : x = où (Sb, désigne le changement estimé de la pseudodistance entre le satellite et le récepteur depuis le pas de temps précédent, (5,b9 le changement estimé, avec l'ambiguïté sur le nombre de cycles complets, de la pseudodistance entre le satellite et le récepteur depuis le pas de temps précédent, (5,^à le changement estimé de la vitesse relative entre le satellite et le récepteur depuis le pas de temps précédent et (51'à le changement estimé de l'accélération relative entre le satellite et le récepteur depuis le pas de temps précédent. En notant X = la règle de mise à jour peut être exprimée comme : Xk = FXk1 +x,,, où F désigne la matrice de transition (dépendant du modèle) et k est l'indice du pas de temps. On obtient l'estimée du retard de code -1"1 Par fi br I c, f2 par f2 c, f3 par f3 = 1 c, 1.12 par :11442 = ;(*).1- / 412 01442 Par 014,12- l'ô, 2r /À12, j>14,13 Par .^f.13 - Pr1413 et 01443 Par 01443 - Jço .271- 1413 ' Ces sorties sont utilisées pour la commande des OCN du dispositif de poursuite 110. [0052] La figure 3 montre une variante 110' du dispositif de poursuite de la figure 2. Les mêmes numéros de référence ont donc été utilisés pour désigner des composantes identiques. Selon la variante de la figure 3, le filtre de Kalman étendu sort les signaux de commande de toutes les boucles de poursuite de code et de porteuse. Le vecteur de mesures appliqué au filtre de Kalman étendu 120' peut s'écrire : z =[c(Si- , Ô, 1, (59,[, (501442/ (5^f14,12 / (501443 / (5114,13 r Le vecteur d'état est : où (Sb, désigne le changement estimé de la pseudodistance entre le satellite et le récepteur depuis le pas de temps précédent, (5,b9 le changement estimé, avec l'ambiguïté sur le nombre de cycles complets, de la pseudodistance entre le satellite et le récepteur depuis le pas de temps précédent, (5,^à le changement estimé de la vitesse relative entre le satellite et le récepteur depuis le pas de temps précédent et (5,*6 le changement estimé de l'accélération relative entre le satellite et le récepteur depuis le pas de temps précédent. En notant X = [b,,,b9,;b,1]T , la règle de mise à jour peut être exprimée comme : +x,,, où F désigne la matrice de transition (dépendant du modèle) et k est l'indice du pas de temps. On obtient l'estimée du retard de code fi par fi 2 par f 2 = bric, r3 par i:3 =bric, 1.12 par '7).r I ilw12 014,12 par 0,4,12 = ;09 - 27C / À,412 :11443 Par 114,13 '7).r 413 et 014,13 Par 014,13 = bço 413 ' Ces sorties sont utilisées pour la commande des OCN du dispositif de poursuite 110'. [0053] La figure 4 montre, de manière très schématique, un récepteur GNSS 64 comprenant une pluralité de dispositifs de poursuite 10 (ou 110 ou 110'). Chaque dispositif de poursuite 10 reçoit le signal entrant, transposé en fréquence (par exemple en bande de base) et numérisé par une tête HF 66 (en anglais « RF front end ») et poursuit les signaux d'un satellite (ou pseudolite) en visibilité depuis le récepteur. Les mesures faites par les dispositifs de poursuite 10 sont combinées par un module de positionnement 68, qui en déduit la position du récepteur, la vitesse, éventuellement l'accélération, et l'heure par résolution du système d'équations de navigation. Le récepteur 64 illustré comprend aussi un navigateur inertiel 70 dont les mesures sont également utilisées par le module de positionnement 68. L'information PVT (position, vitesse, temps) peut être affichée à l'utilisateur grâce à un dispositif d'affichage 72 (p.ex. sur une carte). Optionnellement, le module de positionnement 68 fournit également aux dispositifs de poursuite 10 un signal d'aide que le filtre de la boucle de la première porteuse peut utiliser pour une poursuite améliorée de la phase et de la fréquence de la première porteuse. REVENDICATIONS1. Dispositif de poursuite de signaux de radionavigation (10, 110, 110'), comprenant un premier étage de poursuite (12, 112) d'un premier signal de radionavigation contenu dans un signal entrant à appliquer au dispositif, le premier signal de radionavigation comprenant une première porteuse à une première fréquence modulée par une première forme d'onde d'étalement, ledit premier étage de poursuite comprenant une première boucle à verrouillage de phase de porteuse avec un mélangeur (16, 116) configuré pour multiplier le signal entrant avec une réplique locale de la première porteuse, la première boucle à verrouillage de phase de porteuse comprenant un discriminateur de phase de la première porteuse (32, 132) configuré pour produire un premier signal d'erreur résultant d'un premier écart de phase entre la première porteuse et la réplique locale de la première porteuse, la première boucle à verrouillage de phase de porteuse étant configurée pour régler la phase de la réplique locale de la première porteuse sur la base du premier signal d'erreur ; un deuxième étage de poursuite (14, 114) d'un deuxième signal de radionavigation contenu dans le signal entrant, le deuxième signal de radionavigation comprenant une deuxième porteuse à une deuxième fréquence, différente de la première fréquence, modulée par une deuxième forme d'onde d'étalement, le deuxième étage de poursuite comprenant une deuxième boucle à verrouillage de phase de porteuse avec un mélangeur (40, 140) configuré pour multiplier le signal entrant avec une réplique locale de la deuxième porteuse, caractérisé en ce que la deuxième boucle à verrouillage de phase de porteuse comprend un discriminateur de phase (60, 160) d'un battement entre la première porteuse et la deuxième porteuse configuré pour produire un deuxième signal d'erreur résultant d'une différence entre le premier écart de phase et un deuxième écart de phase entre la deuxième porteuse et la réplique locale de la deuxième porteuse ; et en ce que la deuxième boucle à verrouillage de phase est configurée pour régler la phase de la réplique locale de la deuxième porteuse sur la base des premier et deuxième signaux d'erreur. 2. Dispositif selon la revendication 1, dans lequel la première boucle à verrouillage de phase comprend un oscillateur à commande numérique (18, 118), commandé par un premier filtre (34, 134) recevant le premier signal d'erreur en entrée. 3. Dispositif selon la revendication 2, dans lequel la deuxième boucle à verrouillage de phase comprend un oscillateur à commande numérique (42) commandé par le premier filtre (34, 134) et un deuxième filtre (62, 120, 120') recevant le deuxième signal d'erreur en entrée. 4. Dispositif selon la revendication 3, dans lequel le deuxième filtre comprend un filtre de Kalman, un filtre de Kalman étendu (120, 120') ou un filtre particulaire. 5. Dispositif selon l'une quelconque des revendications 1 à 4, dans lequel les premier et deuxième signaux de radionavigation sont choisis parmi les signaux Galileo El, E5 et E6 ou parmi les signaux GPS L5, L2C et Li ou parmi les signaux GLONASS L3, G2 et Gl. 6. Dispositif selon l'une quelconque des revendications 1 à 4, dans lequel le premier signal de radionavigation est un signal Galileo E5 et dans lequel le deuxième signal de radionavigation est un signal Galileo El ou un signal Galileo E6. 7. Dispositif selon l'une quelconque des revendications 1 à 6, dans lequel le premier étage de poursuite comprend un premier corrélateur (30, 130) configuré pour corréler une réplique locale de la première forme d'onde d'étalement avec le signal entrant multiplié par la réplique locale de la première porteuse, et dans lequel le discriminateur de phase (32, 132) de la première porteuse détermine le premier signal d'erreur sur la base d'un résultat de corrélation produit par le premier corrélateur (30, 130). 8. Dispositif selon la revendication 7, dans lequel le deuxième étage de poursuite comprend un deuxième corrélateur (50, 150) configuré pour corréler une réplique locale de la deuxième forme d'onde d'étalement avec le signal entrant multiplié par la réplique locale de la deuxième porteuse, et dans lequel le discriminateur de phase (60, 160) du battement entre la première porteuse et la deuxième porteuses détermine le deuxième signal d'erreur sur la base du résultat de corrélation produit par le premier corrélateur et d'un résultat de corrélation produit par le deuxième corrélateur (50, 150). 9. Dispositif selon la revendication 8, dans lequel le discriminateur de phase (60, 160) du battement entre la première porteuse et la deuxième porteuse est configuré pour utiliser le produit du résultat de corrélation produit par le premier corrélateur (30, 130) et du conjugué complexe du résultat de corrélation produit par le deuxième corrélateur (50, 150) ou le produit du conjugué complexe du résultat de corrélation produit par le premier corrélateur (30, 130) et du résultat de corrélation produit par le deuxième corrélateur (50, 150) pour produire le deuxième signal d'erreur. 10. Dispositif selon l'une quelconque des revendications 1 à 9, comprenant un troisième étage de poursuite (114, 114') d'un troisième signal de radionavigation contenu dans le signal entrant, le troisième signal de radionavigation comprenant une troisième porteuse à une troisième fréquence, différente des première et deuxième fréquences, modulée par une troisième forme d'onde d'étalement, le troisième étage de poursuite comprenant une troisième boucle à verrouillage de phase de porteuse avec un mélangeur (140') configuré pour multiplier le signal entrant avec une réplique locale de la troisième porteuse, la troisième boucle à verrouillage de phase de porteuse comprenant un discriminateur de phase (160') d'un battement entre la première porteuse et la troisième porteuse configuré pour produire un troisième signal d'erreur résultant d'une différence entre le premier écart de phase et un troisième écart de phase entre la troisième porteuse et la réplique locale de la troisième porteuse ; et la troisième boucle à verrouillage de phase étant configurée pour régler la phase de la réplique locale de la troisième porteuse sur la base des premier et troisième signaux d'erreur. 11. Dispositif selon la revendication 10, dans lequel la troisième boucle à verrouillage de phase comprend un oscillateur à commande numérique commandé par le premier filtre (134) et un troisième filtre (120, 120') recevant le troisième signal d'erreur en entrée. 12. Dispositif selon la revendication 10 si dépendante de la revendication 4, dans lequel la troisième boucle à verrouillage de phase comprend un oscillateur à commande numérique commandé par le filtre de Kalman, le filtre de Kalman étendu (120, 120') ou le filtre particulaire. 13. Dispositif selon l'une quelconque des revendications 10 à 12, si dépendante de la revendication 7, dans lequel le troisième étage de poursuite comprend un troisième corrélateur (150') configuré pour corréler une réplique locale de la troisième forme d'onde d'étalement avec le signal entrant multiplié par la réplique locale de la troisième porteuse, et dans lequel le discriminateur de phase (160') du battement entre la première porteuse et la troisième porteuse détermine le troisième signal d'erreur sur la base du résultat de corrélation produit par le premier corrélateur (130) et d'un résultat de corrélation produit par le troisième corrélateur (150'). 14. Dispositif selon la revendication 13, dans lequel le discriminateur de phase (160') du battement entre la première porteuse et la troisième porteuse est configuré pour utiliser le produit du résultat de corrélation produit par le premier corrélateur (130) et du conjugué complexe du résultat de corrélation produit par le troisième corrélateur (150') ou le produit du conjugué complexe du résultat de corrélation produit par le premier corrélateur (130) et du résultat de corrélation produit par le troisième corrélateur (150') pour produire le troisième signal d'erreur. 15. Dispositif selon l'une quelconque des revendications 1 à 14, réalisé par un processeur de signaux numériques. 16. Récepteur GNSS (64), comprenant un ou plusieurs dispositifs de poursuite de signaux de radionavigation (10, 110, 110') selon l'une quelconque des revendications 1 à 15. 17. Dispositif de positionnement et/ou de datation incluant un récepteur GNSS (64) selon la revendication 16.