Brevet FR3111435A1 — publié le 20211217

Récepteur de signaux de radionavigation (10), comprenant un étage radiofréquence (E1); un générateur (E2) d’un signal de réplique de porteuse (Srep(t)) ; un mélangeur (47) du signal de réplique de porteuse (Srep(t)) avec un signal GNSS numérisé (Sr(t)) ; une boucle d’asservissement de position de code (E4) comprenant un générateur de signaux de réplique de code (E4a); une pluralité de corrélateurs (E4c) aptes à calculer une fonction de corrélation entre le signal de réplique de code principal (SC R(L) (t)) et le signal en bande de base Sbb(t) afin de produire une corrélation principale (50); un calculateur d’indicateur de puissance de corrélation (E4d) configuré pour déterminer des puissances de corrélation et déterminer un maximum de puissance de corrélation; un calculateur de correction (E4e) configuré pour déterminer une fréquence de code ajustée; la boucle d’asservissement de position de code (E4) comprenant une sortie délivrant au moins la corrélation principale (50). Figure 1 Récepteur de signaux de radionavigation comprenant un calculateur d’indicateur de puissance de corrélation La présente invention concerne un récepteur de signaux de radionavigation. ETAT DE LA TECHNIQUE ANTERIEURE Par radionavigation il faut comprendre l’utilisation de signaux radioélectriques émis par des satellites GNSS pour le positionnement ou la navigation. Dans le texte, les termes GNSS de l’anglais « Global Navigation Satellite System » sont les termes génériques qui englobent tous les systèmes de navigation par satellite. Ils se composent de plusieurs systèmes distincts dont par exemple le GPS « Global Positionning System » ou encore Galileo exploités par différents pays ou régions. De manière générale, plusieurs bandes fréquentielles sont allouées pour l’émission des signaux GNSS. Les signaux GNSS sont construits sur le même modèle de signal c’est à dire une ou plusieurs fréquences porteuses sinusoïdales sur lesquelles sont véhiculés différents signaux binaires modulés en phase et contenant les informations nécessaires à la navigation. Il peut y avoir plusieurs types de signaux binaires véhiculés par les porteuses des signaux GNSS : un message de navigation contenant plusieurs paramètres nécessaires au récepteur pour calculer sa position (heure d’envoi du message, éphémérides du satellite, état de la constellation, etc.), un code pseudo-aléatoire permettant l’identification du satellite émetteur et la réalisation de mesures de pseudo-distance avec, éventuellement, une sous-porteuse rectangulaire pour certains signaux ayant une modulation BOC (de la dénomination anglaise Binary Offset Carrier). Les codes pseudo-aléatoires utilisés dans les signaux GNSS sont des codes PRN (de l’anglais « Pseudo Random Noise ») composés d’éléments binaires appelés "bribes". A chaque satellite est associé un ou plusieurs codes PRN - appelé code dans la suite du texte - qui lui sont propre. De manière générale, les codes pseudo-aléatoires ou codes sont des codes binaires quasi orthogonaux entre eux qui permettent aux satellites d’émettre simultanément dans la même bande de fréquence. Dans le domaine de la métrologie, des signaux GNSS dégradés peuvent être utilisés. Notamment, le sondage atmosphérique est étudié par radio-occultation GNSS. Dans cette application, il est parfois difficile de rester synchronisé avec les signaux GNSS dégradés et donc de réaliser des mesures avec eux. Une solution connue classiquement consiste à utiliser un récepteur GNSS muni de boucles de poursuite fermées. Des canaux de réception du signal sont alors munis de boucles de code et de boucles de porteuse imbriquées de manière à réaliser des mesures de pseudo-distance sur le code et sur la porteuse du signal. Dans la plupart des récepteurs GNSS, la boucle de code est aidée par la boucle de porteuse. Cette solution n’est cependant pas satisfaisante pour de tels signaux dégradés car les boucles peuvent décrocher, en particulier la boucle de porteuse, la plus sensible Une autre solution, lorsque le mouvement du récepteur est connu, est de prévoir une aide en vitesse Doppler pour la boucle de porteuse, qui elle-même aide la boucle de code. C’est le cas notamment pour un récepteur en orbite spatiale comme cela est décrit dans le brevet « Procédé de réduction autonome des seuils d’acquisition et de poursuite des porteuses reçues en orbite » (FR2746232B1). De plus, ce document décrit la possibilité d’ouvrir totalement la boucle de porteuse si l’aide en vitesse Doppler est de qualité suffisante. On se trouve alors en mode dit « boucle de porteuse ouverte », seule la boucle de code reste fermée et seule la mesure de pseudo-distance de code est réalisée. Dans le cas de la radio-occultation, la mesure d’intérêt est la mesure de phase de porteuse. Ainsi, la dernière architecture décrite ci-dessus est utilisée en y ajoutant la récupération des sorties des corrélateurs de la voie de code synchrone du code reçue (dite « Prompt » dans la littérature anglo-saxonne, dite « principale » dans la suite du brevet) pour en extraire l’erreur entre la porteuse reçue et la réplique de porteuse générée d’après le modèle de vitesse Doppler injectée, ce qui permet ensuite de reconstruire la mesure de phase de porteuse. Cette reconstruction peut être réalisée au cours d’un post-traitement éventuellement effectué au sol en temps peu différé. C’est le cas par exemple sur les missions METOP ou COSMIC. Une condition nécessaire à ce que la mesure de phase de porteuse puisse se faire est que la boucle de code reste correctement accrochée sur le pic de corrélation du signal ce qui est difficile à réaliser en présence de signaux dont le retard et/ou l’amplitude peuvent varier fortement en cas de fort retard troposphérique. La présente invention a pour but de fournir un récepteur de signaux de radionavigation, permettant de répondre à tout ou partie des problèmes présentés ci-avant. Notamment, un but est de fournir une solution répondant à l’objectif suivant : - permettre une poursuite du code du signal GNSS dégradé de façon robuste et obtenir une mesure de la phase de porteuse signal GNSS dégradé sur une plage étendue. Ce but peut être atteint grâce à un récepteur de signaux de radionavigation, comprenant : - un étage radiofréquence pour produire un signal GNSS Sr(t), à partir d’un signal GNSS incident; - un générateur d’un signal de réplique de porteuse Srep(t), le signal de réplique de porteuse Srep(t) étant une estimée de la porteuse du signal GNSS Sr(t); - un mélangeur du signal de réplique de porteuse Srep(t) avec le signal GNSS numérisé Sr(t), configuré pour obtenir un signal en bande de base Sbb(t) ; - une boucle d’asservissement de position de code comprenant: - un générateur de signaux de réplique de code comprenant un intégrateur de code configuré pour délivrer, en entrée d’un générateur de code, une position de code estimée, le générateur de code comprenant d’une part un générateur de signal de réplique de code principal SC R(L) (t) délivrant en sortie un signal de réplique de code principal SC R(L) (t) et d’autre part une pluralité de générateurs de signaux de répliques de code en avance et en retard SC R(k) (t) présentant chacun un retard ou une avance de code décalés par rapport au signal de réplique de code principal SC R(L) (t) ; chaque réplique de la pluralité étant indexée par un indice entier compris entre 1 et n ; - une pluralité de corrélateurs aptes à calculer une fonction de corrélation : - entre le signal de réplique de code principal SC R(L) (t) et le signal en bande de base Sbb(t) afin de produire une corrélation principale sur un temps d’intégration cohérent Tic; - entre chacun des signaux de répliques de code en avance et en retard SC R(k) (t) et le signal en bande de base Sbb(t) afin de produire des corrélations d’avance et de retard sur le temps d’intégration cohérent Tic ; chaque corrélateur de la pluralité étant indexé par un indice ; - un calculateur d’indicateur de puissance de corrélation configuré pour : - déterminer des puissances de corrélation à partir de la corrélation principale (50) et à partir des corrélations d’avance et de retard ; et - déterminer un maximum de puissance de corrélation à partir des puissances de corrélation de sorte à obtenir un indice primaire kmax correspondant à l’indice du corrélateur montrant le maximum de puissance de corrélation ; - un calculateur de correction configuré pour déterminer une fréquence de code ajustée, en fonction de l’indice primaire kmax, en fonction d’un indice L correspondant au corrélateur traitant le signal de réplique de code principal SC R(L) (t), en fonction d’une période de fonctionnement T fb et en fonction d’un gain ajustable ; la commande de fréquence de bribe de code étant la somme entre la fréquence de code ajustée issue du calculateur de correction et un signal d’aide de porteuse issu du modèle physique multiplié par un coefficient d’échelle, la boucle d’asservissement de position de code comprenant une sortie délivrant au moins la corrélation principale. Certains aspects préférés mais non limitatifs sont les suivants. Dans une mise en œuvre du récepteur de signaux de radionavigation, les corrélateurs d’indices adjacents traitent des signaux de répliques de code SC R(L) (t), SC R(k) (t) séparés entre eux par une valeur constante d d’une bribe du code du signal GNSS. Dans une mise en œuvre du récepteur de signaux de radionavigation, la pluralité de corrélateurs comprend un nombre impair de corrélateurs, le nombre de signaux de répliques de code en avance par rapport au signal de réplique de code principal SCRP(t) et le nombre de signaux de répliques de code en retard SC R(k) (t) par rapport au signal de réplique de code principal SCRP(t) étant égaux. Dans une mise en œuvre du récepteur de signaux de radionavigation, le signal GNSS incident étant doté d’une fréquence porteuse émise et d’un rythme de code, récepteur dans lequel le coefficient d’échelle est égal au rapport entre le rythme de code et de la fréquence de porteuse. Dans une mise en œuvre du récepteur de signaux de radionavigation, la position de code délivrée par l’intégrateur de code correspond à la position de code de la réplique principale. Dans une mise en œuvre du récepteur de signaux de radionavigation, la position de code délivrée par l’intégrateur de code correspond à la position de code de la réplique de code principale translatée d’un écart constant. Dans une mise en œuvre du récepteur de signaux de radionavigation, la position de code délivrée par l’intégrateur de code correspond à une position de code d’une réplique de code la plus en avance. Dans une mise en œuvre du récepteur de signaux de radionavigation, le signal de réplique de porteuse Srep(t), le signal en bande de base Sbb(t), la corrélation principale et les corrélations d’avance et de retard sont complexes ; le calculateur d’indicateur de puissance de corrélation est configuré pour déterminer des puissances de corrélation également à partir de corrélations complexes issues de corrélateurs secondaires ; un calculateur de complexe conjugué , apte à calculer un complexe conjugué à partir de la corrélation complexe principale, est agencé en sortie de la corrélation complexe principale de sorte à délivrer le complexe conjugué de la corrélation principale ; et les corrélateurs secondaires sont agencés de sorte que chacune de la corrélation complexe principale et des corrélations complexes en avance et en retard soit corrélée au complexe conjugué de la corrélation principale sur la période de fonctionnement T fb ; la sortie des corrélateurs secondaires constituant l’entrée du calculateur d’indicateur de puissance de corrélation complexe. Dans une mise en œuvre du récepteur de signaux de radionavigation, la période de fonctionnement T fb correspond à un temps d’intégration non-cohérent T nc . Dans une mise en œuvre du récepteur de signaux de radionavigation, le temps d’intégration non-cohérent T nc est égal à un multiple du temps d’intégration cohérent Tic de la corrélation principale. Dans une mise en œuvre du récepteur de signaux de radionavigation, la période de fonctionnement T fb correspond à un temps d’intégration cohérent. Dans une mise en œuvre, le récepteur de signaux de radionavigation la pluralité de corrélateurs comprend un nombre supérieur ou égal à trois de corrélateurs. Dans une mise en œuvre du récepteur de signaux, le signal GNSS numérisé Sr(t) comporte une fréquence porteuse intermédiaire f i , la commande de fréquence en entrée du générateur de signal de réplique de porteuse étant obtenue à partir de la fréquence Doppler du signal Sr(t) estimée par le modèle externe et de la fréquence f i . Dans une mise en œuvre du récepteur de signaux, la numérisation du signal Sr(t) est réelle en fréquence intermédiaire f i, et le mélangeur du signal de réplique de porteuse Srep(t) avec le signal GNSS numérisé Sr(t) est apte à traiter la numérisation réelle du signal Sr(t). Dans une mise en œuvre du récepteur de signaux, la numérisation du signal Sr(t) est complexe en fréquence intermédiaire f i, et le mélangeur du signal de réplique de porteuse avec le signal GNSS numérisé Sr(t) est apte à traiter la numérisation complexe du signal Sr(t). Dans une mise en œuvre du récepteur de signaux, la boucle d’asservissement de position de code comprend : une pluralité de corrélateurs aptes à calculer une fonction de corrélation complexe : - entre le signal de réplique de code principal et le signal en bande de base afin de produire une corrélation complexe principale sur un temps d’intégration cohérent ; - entre chacun des signaux de répliques de code en avance et en retard et le signal en bande de base afin de produire des corrélations complexes d’avance et de retard sur le temps d’intégration cohérent ; chaque corrélateur de la pluralité étant indexé par un indice ; - un calculateur d’indicateur de puissance de corrélation configuré pour : - déterminer des puissances de corrélation à partir de la corrélation complexe principale et à partir des corrélations complexes d’avance et de retard; et - déterminer un maximum de puissance de corrélation complexe à partir des puissances de corrélation de sorte à obtenir un indice primaire correspondant à l’indice du corrélateur montrant le maximum de puissance de corrélation complexe ; - un calculateur de correction configuré pour déterminer une fréquence de code ajustée, en fonction de l’indice primaire (kmax), en fonction d’un indice correspondant au corrélateur traitant le signal de réplique de code principal, en fonction d’une période de fonctionnement et en fonction d’un gain ajustable ; la commande de fréquence de bribe de code étant la somme entre la fréquence de code ajustée issue du calculateur de correction et un signal d’aide de porteuse issu du modèle physique multiplié par un coefficient d’échelle, la boucle d’asservissement de position de code comprenant une sortie délivrant au moins la corrélation complexe principale. Dans une mise en œuvre du récepteur de signaux, le signal en bande de base est complexe. D’autres aspects, buts, avantages et caractéristiques de l’invention apparaîtront mieux à la lecture de la description détaillée suivante de modes de réalisation préférés de celle-ci, donnée à titre d’exemple non limitatif, et faite en référence aux dessins annexés sur lesquels : illustre une vue schématique d’un exemple de récepteur de signaux de radionavigation selon l’invention; illustre une vue schématique d’un exemple supplémentaire d’un récepteur de signaux de radionavigation selon l’invention où un calculateur de complexe conjugué est agencé en sortie de la corrélation complexe principale et où une pluralité de corrélateurs secondaires est agencée de sorte que chacune de la corrélation complexe principale et des corrélations complexes en avance et en retard soit corrélée au complexe conjugué de la corrélation principale, la sortie des corrélateurs secondaires constituant l’entrée du calculateur d’indicateur de puissance de corrélation complexe. EXPOSE DETAILLE DE MODES DE REALISATION PARTICULIERS Sur les figures 1 et 2 annexées et dans la suite de la description, des éléments identiques ou similaires en termes fonctionnels sont repérés par les mêmes références. De plus, les différents éléments ne sont pas représentés à l’échelle de manière à privilégier la clarté des figures pour en faciliter la compréhension. Par ailleurs, les différents modes de réalisation et variantes ne sont pas exclusifs les uns des autres et peuvent, au contraire, être combinés entre eux. Dans la suite de la description, sauf indication contraire, les termes « sensiblement », « environ », « globalement » et « de l'ordre de » signifient « à 10 % près ». De façon générale et dans la suite du texte, un signal ayant une valeur complexe ou dit complexe, est un signal ayant une composante en phase qui correspond à la partie réelle et une composante en quadrature qui correspond à la partie imaginaire. Dans l’ensemble du texte, le terme « période » est homogène à un temps et pourra être utilisé de façon équivalente. Dans l’ensemble des modes de réalisations de l’invention et comme illustré sur les figures 1 et 2, le récepteur de signaux de radionavigation 10 comprend un étage radiofréquence E1. L’étage radiofréquence E1 peut comprendre par exemple une ou plusieurs antennes, des éléments de filtrage, d’amplification, de translation de fréquence, de numérisation, et un contrôle automatique de gain. L’étage radiofréquence E1 produit un signal GNSS numérisé ou analogique Sr(t), dépendant du temps t, à la fréquence d’échantillonnage Fe qui vaut typiquement entre 1 et 50 MHz, à partir du signal GNSS incident. Le signal GNSS numérisé Sr(t) comprend ainsi une bande de fréquence de signal numérisé. Le signal GNSS numérisé Sr(t) peut être réel ou complexe. La fréquence de la porteuse du signal numérisé est généralement abaissée vers une fréquence intermédiaire fi. fi peut être proche de zéro voire égale à zéro, dans ce cas on dit que la numérisation se fait en bande de base. Dans la suite du texte les signaux indiqués avec (t) sont dépendants du temps t. Entre la réception du signal et les mesures effectuées sur le signal numérisé, le récepteur de signaux de radionavigation 10 doit effectuer plusieurs étapes pour extraire les informations nécessaires. En particulier, des étapes d’acquisition et de poursuite sont réalisées qui ont pour but de synchroniser le récepteur avec les signaux de navigation. Dans le texte, les termes « modèle physique » correspondent à une aide Doppler ou une aide de vitesse fine. Le récepteur de signaux de radionavigation 10 comprend également un générateur E2 de signal de réplique de porteuse Srep(t) pour former une réplique locale de la porteuse. Ce signal de réplique de porteuse Srep(t) est un signal complexe, représenté par une valeur complexe et généré à la cadence Fe. Le signal de réplique de porteuse Srep(t) consiste en une estimée de la porteuse du signal d’entrée. Ce signal est défini par : Srep(t) = e j.ϕ(t) où ϕ(t) est la phase de ce signal de réplique de porteuse Srep(t). La phase du signal de réplique de porteuse, est obtenue en intégrant une commande de fréquence en entrée d’un intégrateur de porteuse 40. L’intégrateur de porteuse 40 est un intégrateur pur de type oscillateur à commande numérique (ou NCO / « Numerically controlled oscillator » selon la terminologie anglo-saxonne). Dans un exemple, le générateur E2 d’un signal de réplique de porteuse Srep(t) complexe à la même cadence d’échantillonnage Fe, dont une phase ϕ(t) est obtenue en intégrant, avec un intégrateur de porteuse 40, une commande de fréquence en entrée du générateur E2 de signal de réplique de porteuse, la commande de fréquence en entrée du générateur E2 de signal de réplique de porteuse étant obtenue à partir d’au moins une fréquence Doppler du signal Sr(t) estimée par un modèle physique 41, le générateur E2 comprenant un dispositif 40a apte à générer le signal de réplique de porteuse Srep(t) complexe à partir de la sortie de l’intégrateur de porteuse 40. Dans un exemple de mise en œuvre du récepteur de signaux de radionavigation 10, le générateur E2 comprend un dispositif 40a apte à générer le signal de réplique de porteuse Srep(t) complexe à partir de la sortie de l’intégrateur de porteuse 40. Le récepteur de signaux de radionavigation 10 comprend en outre un mélangeur 47 du signal de réplique de porteuse Srep(t) avec le signal GNSS numérisé Sr(t). Le signal de sortie du mélangeur 47 est un signal en bande de base Sbb(t), échantillonné à la cadence Fe. Le signal en bande de base Sbb(t) et le signal de réplique de porteuse Srep(t) sont complexes. Dans une mise en œuvre du récepteur de signaux de radionavigation 10, la numérisation du signal Sr(t) est complexe. Dans une autre mise en œuvre du récepteur de signaux de radionavigation 10, la numérisation du signal Sr(t) est réelle. Le mélangeur 47 du signal de réplique de porteuse Srep(t) avec le signal GNSS numérisé Sr(t) est prévu pour être apte à traiter la numérisation réelle ou complexe du signal Sr(t) selon l’équation Math. 1 suivante : Sbb(t) = Sr(t) x Srep(t) Afin de réaliser la poursuite de code, le récepteur de signaux de radionavigation 10 comprend une boucle d’asservissement de position de code E4. La boucle d’asservissement de position de code E4 peut également être appelée poursuite de code. La boucle d’asservissement de position de code E4 comprend tout d’abord un générateur de signaux de réplique de code E4a fonctionnant à la cadence Fe. Le générateur de signaux de réplique de code E4a comprend notamment un intégrateur de code 45 et un générateur de code 43. L’intégrateur de code 45 est un intégrateur pur de type oscillateur à commande numérique (ou NCO / « Numerically controlled oscillator » selon la terminologie anglo-saxonne) qui délivre par exemple une position de code dite principale à la cadence du signal échantillonné Fe. Cet intégrateur de code 45 est configuré pour délivrer, en entrée d’un générateur de code 43, une position de code estimée, comme la position τ(t) du code principal ou, selon les implémentations, de tout autre réplique de code, typiquement la réplique la plus en avance. Pour cela il reçoit en entrée une commande de fréquence des bribes de code constituée de la somme entre la fréquence issue de la boucle de code, un biais 80 de rythme de code constant et valant généralement le produit RC GNSS et un signal d’aide de porteuse issu du modèle de Doppler 41 multiplié par un coefficient d’échelle 44 appelé K = RC GNSS /f GNSS où RC GNSS est le rythme de code transmis par le GNSS et f GNSS est la fréquence de porteuse. Le générateur de code 43 comprend un générateur 43a de signal de réplique de code principal SC RP (t) délivrant en sortie un signal réel de réplique de code principal SC RP (t). Le générateur de code 43 comprend également une pluralité de générateurs 43b de signaux réels de répliques de code en avance et en retard. Une autre notation des signaux réplique de code est : SC R(k) (t) où l’indice k est le rang de la réplique. Chacun des signaux de répliques de code en avance SC R(k) (t) présente une avance de code différente et échelonnée, autrement dit décalée, de façon quelconque ou discrète selon la mise en œuvre choisie, par rapport au signal de réplique de code principal SC RP (t). Chacun des signaux de répliques de code en retard SC R(k) (t) présente une position de code en retard, autrement dit présente un retard de code différent et échelonné de façon quelconque ou discrète par rapport au signal de réplique de code principal SC RP (t). On définit L comme étant le rang de la réplique principale de code, c’est à dire : , Dans tout le texte, les termes SC RP (t) et SC R(L) (t) sont équivalents. Les positions (avance si valeur positive, retard si valeur négative) des signaux répliques par rapport au compteur issu de l’intégrateur de code 45 sont notés Δτ k . Une réplique de code SC R(k) (t) est en avance si : Δτ k > Δτ L. Une réplique de code SC R(k) (t) est en retard si : Δτ k L . Dans un autre exemple de mise en œuvre du récepteur de signaux de radionavigation 10, chacun des signaux de répliques de code principal SC RP (t) et en avance et en retard SC R(k) (t) sont séparés entre eux par une valeur constante d exprimée en bribes du code du signal GNSS. En d’autres termes, l’intervalle entre deux répliques de code adjacentes est identique. Dans le cas où l’intégrateur de code 45 correspond à la réplique la plus en avance, les écarts de positions Δτ k entre chaque réplique de code SC R(k) (t) et la sortie de l’intégrateur de code 45 valent : {Δτ 1 ; Δτ 2 ;…; Δτ L, ; … ; Δτ n-1 ; Δτ n } = {-(n-1)d ; -(n-2)d ; … ; -(n-L)d ; … ; d ; 0 } Dans un autre exemple, le nombre de signaux de répliques de code en retard SC R(k) (t) est égal au nombre de signaux de répliques de code en avance SC R(k) (t). Dans cet exemple, l’indice L de la réplique principale vaut : L=(n+1)/2, le nombre total de répliques étant impair dans ce cas. Si ces deux derniers cas sont combinés, alors les écarts de position des répliques de code dans le cas où l’intégrateur de code 45 correspond à la réplique principale valent : { Δτ 1 ; … ; Δτ L-1 ; , Δτ L, ; Δτ L+1 ; … ; Δτ n } = {-(n-1)d/2 ; … ; -d ; 0 ; d ; … ; (n-1)d/2}. La boucle d’asservissement de position de code E4 comprend en outre la même pluralité n de corrélateurs E4c. Selon un mode de réalisation, un corrélateur parmi la pluralité de corrélateur E4c comprend un mélangeur et un sommateur avec remise à zéro. Les mélangeurs et les sommateurs constituent ensemble des corrélateurs E4c aptes à réaliser des fonctions de corrélations complexes connues de l’homme du métier. Par corrélation complexe il faut comprendre que les parties en phase et en quadrature des signaux sont traitées lors de la corrélation. Les parties en phase de la corrélation sont notées avec un indice I : C I(k) (t) et les parties en quadrature sont notées avec un indice Q : C Q(k) (t). Chacun des signaux de répliques de code principal SC R(L) (t) et en avance et en retard SC R(k) (t) subit une opération de corrélation complexe sur une période appelée temps d’intégration cohérent égale à Tic avec le signal en bande de base Sbb(t). Tic présente une valeur comprise typiquement entre 1 et quelques dizaines de millisecondes. Le résultat de la corrélation complexe C (k) (t) avec la réplique d’ordre k, à l’instant t est donné par l’équation Math. 3 suivante: T e représente l’intervalle entre 2 échantillons numérisés: T e = 1/Fe . La corrélation complexe principale 50 aussi notée C (L) (t j ) est aussi appelée point de corrélation « Prompt » ou encore appelée corrélation « synchrone ». Des corrélateurs sont également agencés entre chacun des signaux de répliques de code en avance SC R(k) (t) et le signal de bande de base Sbb(t) afin de produire des corrélations en avance 51, aussi notées C (k) (T j ) et appelées point de corrélation « en avance », sur un temps d’intégration cohérent Tic. Des corrélateurs sont également agencés entre chacun des signaux de répliques de code en retard SC R(k) (t) et le signal en bande de base Sbb(t) afin de produire des corrélations de retard 51, aussi notées C (k) (t j ) et appelées point de corrélation « En retard », suivant le temps d’intégration cohérent Tic. Dans un autre exemple, la pluralité E4c de mélangeurs et de corrélateurs comprend un nombre impair de mélangeurs et le même nombre impair de corrélateurs E4c. Dans un exemple, de façon avantageuse, le nombre de signaux de répliques de code SC R(k) (t) en avance par rapport au signal de réplique de code principal SC R(L) (t) est égal au nombre de signaux de répliques de code SC R(k) (t) en retard SC R(k) (t) par rapport au signal de réplique de code principal SC R(L) (t). La boucle d’asservissement de position de code E4 comprend également un calculateur d’indicateur de puissance de corrélation E4d. Celui-ci est configuré pour déterminer un indicateur représentant la puissance ou une grandeur liée à la puissance comme une amplitude de corrélation par exemple. Ici, celui-ci est configuré pour déterminer des puissances de corrélation à partir de la corrélation complexe principale 50 et à partir des différentes corrélations complexes d’avance et de retard 51 ou à partir de corrélations complexes issues de corrélateurs secondaires 71. La détermination d’une puissance de corrélation par le calculateur d’indicateur de puissance de corrélation E4d consiste à calculer le carré du module de chacune des corrélations complexes. Ici cela consiste à calculer le carré du module d’une corrélation complexe liée à un corrélateur d’indice k donné. L’équation Math. 4 décrit le calcul d’une puissance de corrélation appelée ici Power_Corr(k) où k correspond à un indice de corrélateur. Le calculateur d’indicateur de puissance de corrélation E4d est également configuré pour déterminer un maximum de puissance de corrélation entre chacune des puissances de corrélation calculées comme décrit dans le paragraphe précédent. Ainsi il est obtenu un indice primaire correspondant à l’indice du corrélateur montrant le maximum de puissance et qui, pour l’exemple qui suit, s’appelle kmax. La boucle d’asservissement de position de code E4 comprend également un calculateur de correction E4e. Le calculateur de correction E4e est configuré pour déterminer une correction de fréquence de code calculée à chaque itération de la boucle. La fréquence de code est ainsi ajustée avec une correction établie en fonction de l’écart de position de la réplique de l’indice primaire : Δτ kmax (t) et de l’indice L : Δτ L (t) désignant la réplique du code principal SC R(L) (t). La correction de fréquence δ dépend également d’une période de fonctionnement T fb de la boucle et d’un gain, ici appelé G pour l’exemple, qui peut être modulable pour optimiser le fonctionnement de la boucle. Une valeur nominale de G est égale à 1. La correction de fréquence δ est alors exprimée dans l’équation Math. 5 suivante : Dans le cas où les répliques de codes sont espacées d’un intervalle régulier d exprimé en bribe de code, on obtient la position du signal réplique ayant la meilleure corrélation par rapport au signal de réplique de code principal SC R(L) (t) en effectuant (kmax-L)xd, et la correction δ de fréquence de l’intégrateur de code 45 devient : Les corrections de fréquence δ sont appliquées à l’intégrateur de code 45 pendant la période de fonctionnement T fb qui suit jusqu’à ce qu’une nouvelle correction de fréquence soit calculée d’après les sorties de corrélateur à la période de fonctionnement T fb . T fb correspond, dans une mise en œuvre, à un temps cohérent d’intégration. Avantageusement, la période de fonctionnement de la boucle sera égale au temps d’intégration cohérent du corrélateur suivant : T fb = T ic . Une fois réinjectée vers l’intégrateur de code 45, une telle correction de fréquence issue du calculateur d’indicateur de puissance de corrélation E4d et du calculateur de correction E4e permet avantageusement de rendre plus robuste la corrélation du code du signal reçu de manière à réaliser une mesure de phase de porteuse sur les sorties de la corrélation complexe principale 50. En effet, les techniques classiques en boucle fermée ne permettent pas de poursuivre des signaux fortement perturbés par exemple pour l’utilisation en radio-occultation ou en navigation en environnement sévère. Ainsi la boucle d’asservissement de position de code E4 comprend une sortie délivrant la corrélation complexe principale 50 à chaque itération de boucle. Cette sortie est, comme illustré sur les figures 1 et 2, placée avant la détermination de puissance de corrélation. La corrélation principale 50 et les paramètres des signaux dont elle est issue sont ainsi transmis pour des traitements ultérieurs soit en temps réel, soit en temps différé par exemple dans des stations terrestres de traitement des données de radio-occultation. La commande de fréquence en entrée du générateur E2 est la somme de la fréquence Doppler du signal Sr(t) estimée par un modèle externe 41 et de la fréquence f i . Ainsi la poursuite reste stable et l’extraction de la phase de porteuse est assurée. Ceci est essentiel dans des applications par exemple de radio occultation. Dans un exemple supplémentaire de réalisation illustré sur la figure 2, le récepteur de signaux de radionavigation 10 comprend un calculateur de complexe conjugué 70. Le calculateur de complexe conjugué 70 est apte à calculer un complexe conjugué à partir de la corrélation complexe principale 50 ou à partir du signal de réplique de code principal SC R(L) (t). Le calculateur de complexe conjugué 70 est agencé en sortie de la corrélation complexe principale 50 de sorte à délivrer le complexe conjugué de la corrélation principale 50 ou du signal de réplique de code principal SC R(L) (t). Les corrélateurs secondaires 71 sont agencés de sorte que chacune de la corrélation complexe principale 50 et des corrélations complexes en avance et en retard 51 soit corrélée au complexe conjugué 70 de la corrélation principale 50 ou au complexe conjugué 70 du signal de réplique de code principal SC R(L) t. Ces corrélateurs secondaires 71 réalisent alors une intégration non cohérente dont la période t nc est multiple entier du temps d’intégration cohérent T ic des corrélateurs principaux 50 tel que : t nc = N x t ic . La sortie des corrélateurs secondaires 71 constitue ensuite l’entrée du calculateur d’indicateur de puissance de corrélation complexe E4d. Dans une implémentation, la période de fonctionnement t fb de la boucle de code vaudra avantageusement le temps d’intégration cohérent : T fb = T nc. Cette architecture est avantageuse car elle permet d’étendre le temps d’intégration total et d’augmenter le rapport signal sur bruit. Récepteur de signaux de radionavigation (10), comprenant : - un étage radiofréquence (E1) pour produire un signal GNSS (Sr(t)), à partir d’un signal GNSS incident; - un générateur (E2) d’un signal de réplique de porteuse (Srep(t)) , le signal de réplique de porteuse (Srep(t)) étant une estimée de la porteuse du signal GNSS (Sr(t)) ; - un mélangeur (47) du signal de réplique de porteuse (Srep(t)) avec le signal GNSS numérisé (Sr(t)), configuré pour obtenir un signal en bande de base Sbb(t) ; - une boucle d’asservissement de position de code (E4) comprenant: - un générateur de signaux de réplique de code (E4a) comprenant un intégrateur de code (45) configuré pour délivrer, en entrée d’un générateur de code (43), une position de code estimée, le générateur de code (43) comprenant d’une part un générateur (43a) de signal de réplique de code principal (SC R(L) (t)) délivrant en sortie un signal de réplique de code principal (SC R(L) (t)) et d’autre part une pluralité de générateurs (43b) de signaux de répliques de code en avance et en retard (SC R(k) (t)) présentant chacun un retard ou une avance de code décalés par rapport au signal de réplique de code principal (SC R(L) (t)) ; chaque réplique de la pluralité étant indexée par un indice entier compris entre 1 et n ; - une pluralité de corrélateurs (E4c) aptes à calculer une fonction de corrélation : - entre le signal de réplique de code principal (SC R(L) (t)) et le signal en bande de base Sbb(t) afin de produire une corrélation principale (50) sur un temps d’intégration cohérent (Tic); - entre chacun des signaux de répliques de code en avance et en retard (SC R(k) (t)) et le signal en bande de base (Sbb(t)) afin de produire des corrélations d’avance et de retard (51) sur le temps d’intégration cohérent (Tic) ; chaque corrélateur de la pluralité étant indexé par un indice ; - un calculateur d’indicateur de puissance de corrélation (E4d) configuré pour : - déterminer des puissances de corrélation à partir de la corrélation principale (50) et à partir des corrélations d’avance et de retard (51); et - déterminer un maximum de puissance de corrélation à partir des puissances de corrélation de sorte à obtenir un indice primaire (kmax) correspondant à l’indice du corrélateur montrant le maximum de puissance de corrélation ; - un calculateur de correction (E4e) configuré pour déterminer une fréquence de code ajustée, en fonction de l’indice primaire (kmax), en fonction d’un indice (L) correspondant au corrélateur traitant le signal de réplique de code principal (SC R(L) (t)), en fonction d’une période de fonctionnement (T fb ) et en fonction d’un gain ajustable ; la commande de fréquence de bribe de code étant la somme entre la fréquence de code ajustée issue du calculateur de correction (E4e) et un signal d’aide de porteuse issu du modèle physique (41) multiplié par un coefficient d’échelle (44), la boucle d’asservissement de position de code (E4) comprenant une sortie délivrant au moins la corrélation principale (50). Récepteur de signaux de radionavigation (10) selon la revendication 1 dans lequel les corrélateurs d’indices adjacents traitent des signaux de répliques de code (SC R(L) (t)), (SC R(k) (t)) séparés entre eux par une valeur constante d d’une bribe du code du signal GNSS. Récepteur de signaux de radionavigation (10) selon l’une quelconque des revendications 1 ou 2 dans lequel la pluralité de corrélateurs (E4c) comprend un nombre impair de corrélateurs, le nombre de signaux de répliques de code en avance par rapport au signal de réplique de code principal (SCRP(t)) et le nombre de signaux de répliques de code en retard (SC R(k) (t)) par rapport au signal de réplique de code principal (SCRP(t)) étant égaux. Récepteur de signaux de radionavigation (10) selon l’une des revendications 1 à 3, le signal GNSS incident étant doté d’une fréquence porteuse émise et d’un rythme de code, récepteur dans lequel le coefficient d’échelle (44) est égal au rapport entre le rythme de code et de la fréquence de porteuse. Récepteur de signaux de radionavigation (10) selon l’une des revendications 1 à 4, dans lequel la position de code délivrée par l’intégrateur de code (45) correspond à la position de code de la réplique principale. Récepteur de signaux de radionavigation (10) selon l’une des revendications 1 à 5, dans lequel la position de code délivrée par l’intégrateur de code (45) correspond à la position de code de la réplique de code principale translatée d’un écart constant. Récepteur de signaux de radionavigation (10) selon la revendication 6, dans lequel la position de code délivrée par l’intégrateur de code (45) correspond à une position de code d’une réplique de code la plus en avance. Récepteur de signaux de radionavigation (10) selon l’une quelconque des revendications 1 à 7, dans lequel le signal de réplique de porteuse (Srep(t)), le signal en bande de base Sbb(t), la corrélation principale (50) et les corrélations d’avance et de retard (51) sont complexes ; dans lequel le calculateur d’indicateur de puissance de corrélation (E4d) est configuré pour déterminer des puissances de corrélation également à partir de corrélations complexes issues de corrélateurs secondaires (71) ; dans lequel un calculateur de complexe conjugué (70), apte à calculer un complexe conjugué à partir de la corrélation complexe principale (50), est agencé en sortie de la corrélation complexe principale (50) de sorte à délivrer le complexe conjugué de la corrélation principale (50) ; et dans lequel les corrélateurs secondaires (71) sont agencés de sorte que chacune de la corrélation complexe principale (50) et des corrélations complexes en avance et en retard (51) soit corrélée au complexe conjugué (70) de la corrélation principale (50) sur la période de fonctionnement (T fb ) ; la sortie des corrélateurs secondaires (71) constituant l’entrée du calculateur d’indicateur de puissance de corrélation complexe (E4d). Récepteur de signaux de radionavigation (10) selon la revendication 8, dans lequel la période de fonctionnement (T fb ) correspond à un temps d’intégration non-cohérent (T nc ). Récepteur de signaux de radionavigation (10) selon la revendication 9, dans lequel le temps d’intégration non-cohérent (T nc ) est égal à un multiple du temps d’intégration cohérent (Tic) de la corrélation principale. Récepteur de signaux de radionavigation (10) selon l’une des revendications 1 à 8, dans lequel la période de fonctionnement (T fb ) correspond à un temps d’intégration cohérent.