Brevet FR2892202A1 — publié le 20070420

RECEPTEUR GNSS A PRECISION AMELIOREE UTILISANT DEUX PORTEUSES DE SIGNAL L'invention concerne les récepteurs de navigation par satellite utilisant les signaux émis par des satellites comportant chacun deux porteuses proches en fréquence modulées par des codes d'étalement de manière synchrone. Les systèmes de positionnement par satellite mettent en oeuvre plusieurs satellites transmettant leurs positions par des signaux radioélectriques à un récepteur placé à la position à localiser estimant les distances, dites pseudo-distances, qui le séparent des satellites à partir des temps de propagation des signaux des satellites captés. La localisation du ~o récepteur est effectuée par une technique voisine de la triangulation et cette localisation est d'autant plus précise que les positions des satellites sont connues avec précision du récepteur et que les mesures des pseudo-distances faites par le récepteur sont précises. Les positions des satellites sont déterminées à partir d'un réseau 15 de stations sol de poursuite indépendant des récepteurs de positionnement. Elles sont communiquées aux récepteurs de positionnement par les satellites eux-mêmes par les données modulant les signaux émis. Les pseudo-distances sont déduites par les récepteurs de positionnement des retards apparents présentés par les signaux reçus : Ces 20 signaux sont datés à l'émission par les horloges des satellites, qui sont synchronisée sur le temps système, et à la réception par l'horloge du récepteur, pouvant présenter un biais par rapport au temps du système. Les distances déduites des temps de propagation ainsi mesurés sont appelées pseudo-distances car elles sont entachées d'une erreur commune pouvant 25 être importante, due au biais de l'horloge du récepteur. Ce biais est éliminé lors de la résolution du point à condition que les signaux d'au moins quatre satellites soient reçus. Si la précision de la connaissance des positions des satellites du système de positionnement est indépendante des performances d'un 30 récepteur de positionnement, ce n'est pas le cas de celle des mesures de pseudo-distance qui dépend de la précision des mesures de temps de propagation des signaux, au niveau du récepteur. Les signaux radioélectriques émis par des satellites parcourant de grandes distances, et étant émis avec des puissances limitées, parviennent avec de très faibles puissances aux récepteurs, noyés dans un bruit radioélectrique. Pour améliorer leurs réceptions on a cherché à les rendre les moins sensibles possible aux parasites à bande étroite, en augmentant leurs largeurs de bande au moyen de la technique de la bande étalée. Les systèmes actuels et prévus dans un futur proche, de positionnement par satellites utilisent, pour les signaux radioélectriques émis par leurs satellites, la technique d'étalement de bande par modulation à l'aide de séquences binaires pseudo-aléatoires, technique dite DSSS (sigle tiré de l'anglo-saxon : Direct Sequence Spread Spectrum"). Cette modulation DSSS consiste après avoir mis les données à transmettre sous la forme d'une suite d'éléments binaires à débit régulier, à effectuer le produit de chaque élément binaire d'information avec une séquence binaire pseudo-aléatoire connue à l'avance de débit nettement plus rapide. L'étalement de bande obtenu est proportionnel au débit de la suite binaire pseudo-aléatoire d'étalement. On parle de code d'étalement. Les données à transmettre par les satellites, une fois mises sous forme d'une suite de données binaires étalée en fréquence par un code d'étalement selon la technique DSSS, sont transposées dans la plage de fréquence d'émission par modulation avec une porteuse d'émission. En réception, les données contenues dans un signal radioélectrique de satellite d'un système de positionnement sont extraites par deux démodulations successives : une première démodulation à l'aide d'une porteuse engendrée localement par un oscillateur piloté par une boucle de poursuite en fréquence et en phase dite PLL (sigle tiré de l'anglo-saxon : "Phase Lock Loop") - ou boucle de phase de porteuse - permettant de transposer le signal reçu en bande de base et une deuxième démodulation à l'aide de code d'étalement engendrés localement par un générateur de code piloté par une boucle de poursuite en retard dite DLL (sigle tiré de l'anglo-saxon : Delay Lock Loop) - ou boucle de code - permettant de désétaler la suite de données binaires présente dans le signal reçu. Les codes d'étalement engendrés localement sont identiques à ceux émis par le satellite, au retard prés dû à la propagation dans l'espace. Lorsque qu'une boucle de code est accrochée le code engendré localement est en phase avec le code contenu dans le signal reçu du satellite Ainsi les retards affectant les code d'étalement reçus sont accessibles au niveau des signaux d'asservissement des boucles de code. Les retards constatés par ces boucles permettent des mesures non ambiguës ou faiblement ambiguës, des temps de propagation des signaux car les nombres de séquences pseudo-aléatoires entières répétées s'écoulant pendant les trajets des signaux sont relativement petits. On parle de mesures de code. Par exemple, pour le système de positionnement par satellites GPS (sigle tiré de l'anglo-saxon :Global Positioning System), la séquence binaire pseudo-aléatoire répétée la plus courte, celle utilisée pour l'étalement des signaux de satellites de type C/A (sigle provenant de I'anglosaxon :"Coarse/Acquisition Code ou Clear/Acquisition Code), est composée de 1023 éléments binaires avec un débit de 1,023 MHz et une durée d'une milliseconde. Sa durée totale correspond à un trajet de 300 Km pour une onde radioélectrique et permet des mesures de distance modulo de 300 Km. La durée de 1 microseconde de chacun de ses éléments binaires autorise une précision de l'ordre de 0,1 microsecondes. L'ambiguïté des mesures de pseudo-distance obtenues à partir de la séquence binaire pseudo-aléatoire d'un code C/A due au fait que l'on affaire à des mesures modulo 300 Km est facile à lever dès que le récepteur reçoit plus de quatre satellites car il peut alors faire un contrôle de cohérence et ne retenir que la solution cohérente. En l'absence d'une telle possibilité, l'ambiguïté peut être aussi levée à l'aide d'une connaissance préalable très grossière de la position. Une telle ambiguïté de mesure ne se pose pas avec les signaux de satellite du type P du système GPS qui utilisent pour leur étalement une séquence binaire pseudo-aléatoire cryptée non répétée mais ces signaux ne sont pas à la libre disposition des utilisateurs. Les retards apparents des porteuses d'émission sont accessibles, modulo les périodes de ces porteuses, par les phases locales produites par les boucles de phase de porteuse PLL et pilotant les générateurs de porteuses locales. On parle de mesures de phase. Ces mesures sont très précises mais très ambiguës. Dans le cas du système GPS, les signaux accessible au public utilisent une porteuse à 1,575 GHz permettant des mesures de pseudo-distance modulo 0.19 m donc très ambiguës, puisque la distance aux satellites est de l'ordre de 20.000 Km. Certains systèmes de navigation par satellite utilisent des paires de signaux de fréquence porteuse proches de façon à pouvoir les utiliser de manière cohérente et ainsi améliorer la précision des mesures de position, du fait d'un spectre plus étendu. La figure 1 représente les spectres d'une paire de signaux émis par un satellite. Chaque satellite émet deux signaux a et b constitués chacun d'une porteuse modulée par un code pseudo-aléatoire qui étale le spectre. Les deux signaux, dans cet exemple E5a (signal a) et E5b (signal b) ayant des fréquences de porteuses différentes mais proches, sont émis de façon synchrone, les relations entre les phases des deux porteuses et des deux codes d'étalement à l'émission étant connues à tout instant. Par la suite, il sera convenu des notations et définitions suivantes : La fréquence de porteuse centrale est par définition la fréquence située au milieu des deux fréquences porteuses (par exemple E5a et E5b). La fréquence de sous-porteuse est par définition la distance entre la fréquence de porteuse du signal a ou b et la fréquence de porteuse centrale (où Fb ù Fp = Fp ù Fa ) Notations : Fp Fréquence de porteuse centrale Fp = ( Fa + Fb )/2 Fsp Fréquence de sous-porteuse Fsp = ( Fb - Fa )/2 Fa Fréquence de porteuse E5a Fa = Fp - Fsp Fb Fréquence de porteuse E5b Fb = Fp + Fsp Fc Fréquence du code (Op = 2nt . Fp Pulsation de porteuse centrale Wsp= 27t . Fsp Pulsation de sous-porteuse Wa = 2Tt . Fa Pulsation de porteuse E5a Wb = 2n . Fb Pulsation de porteuse E5b (op =(OJa+W) )/2 Wsp=( Wb - Wa)/2 Wa = (Op - Wsp Wb=Wp+Wsp X, = gît / Fp Longueur d'onde de la porteuse centrale 2 = 2n / Far, Longueur d'onde de la sous-porteuse Xa = 27c / Fa Longueur d'onde de la porteuse E5a Xb = 271 / Fb Longueur d'onde de la porteuse E5b Xcode = 1/Fc Longueur d'onde du code Dans la suite du docurnent on écrira expi(0) pour e'e où j2 = -1 Dans le cas de services bi-fréquences, les deux signaux dans le to récepteur peuvent être poursuivis indépendamment sur chaque bande par des boucles de poursuite PLL et DLL. Les mesures de pseudo-distance, qui reposent sur les retards constatés par les boucles de code, sont limitées en précision par la largeur de chaque bande disponible. Du fait de l'indépendance des boucles de poursuite, toute utilisation conjointe des 15 mesures issues de signaux de largeurs de bande différentes, conduira à une précision limitée par la moins bonne des deux poursuites. Il est cependant possible d'améliorer la précision des mesures de pseudo-distance dans le cadre d'un service bi-fréquence, à signaux synchrones de navigation, en combinant les signaux reçus sur les deux 20 fréquences pour bénéficier d'une largeur de bande de fréquence équivalente plus importante. La figure 2 montre les étages analogiques d'entrée d'un récepteur de position par satellite utilisant une paire de porteuse E5a, E5b telles que représentées à la figure 1. Les voies analogiques Va, Vb entre l'antenne et 25 les convertisseurs numériques ADC 12, 14 en fréquence intermédiaire (FI) ont des éléments distincts entre les deux signaux reçus S1 et S2, ce qui entraîne des différentiels de retard de propagation et de phase. La figure 3 montre un des deux étages numériques identiques d'un récepteur de positionnement par satellite de l'état de l'art attaqué par 30 une des deux composantes numérisées E5a ou E5b en fréquence intermédiaire (FI) en sortie de l'étage analogique de la figure 2. Chaque étage numérique comporte un circuit de corrélation de porteuse 20 suivi d'un circuit de corrélation de code 22. Pour chaque signal satellite (E5a ou E5b) on réalise un traitement 35 suivant l'état de l'art avec une boucle de phase de porteuse et une boucle de code. Pour cela on génère une porteuse locale PI et des codes locaux Ca, Cp, Cr pour démoduler le signal par corrélation. On parle de canal de poursuite. On asservit la phase de la porteuse locale PI et la position des codes locaux en phase avec le signal reçu grâce aux boucles de poursuite de chaque canal. Chaque boucle de poursuite comporte, pour la boucle de porteuse un discriminateur de porteuse DSP, un correcteur de porteuse CRP contrôlant un oscillateur NCOp 24 de porteuse générant une phase locale de porteuse pilotant un générateur de porteuse 26 fournissant la porteuse locale pour le circuit de corrélation de porteuse et, pour la boucle de code, un discriminateur de code DSC, un correcteur de code CRC contrôlant un oscillateur NCOc 28 de code générant une phase locale de code pilotant le générateur de code 30 générant les codes locaux Ca, Cp, Cr (avance, ponctuel et retard) pour le circuit de corrélation de code. Les discriminateurs de code et de phase de porteuse permettent de mesurer en sortie des corrélateurs après intégration par les intégrateurs INT les écarts de phase de porteuse et de code entre le signal reçu et le signal local, pour rétro-action dans les boucles de poursuite. La boucle de phase de porteuse aide la boucle de code pilotant l'oscillateur local de code NCOc afin de réduire le traînage de la boucle dû à la dynamique, ce qui permet de réduire la bande de boucle de code et donc le bruit sur la phase du code mesurée. On réalise ce traitement en parallèle pour chaque signal d'un satellite (soit un canal par signal satellite). Les codes différents et décorrélés d'un satellite à l'autre permettent de dissocier les signaux entre les satellites et donc d'allouer un canal à un satellite. Les mesures sont les phases de la porteuse locale Ppb (rad) et du code local (pcb (s) pour chaque canal. Les oscillateurs contrôlés numériquement NCOp ou NCOc ont le rôle de produire les phases des signaux locaux à haute cadence (> 10 MHz) à partir des commandes en vitesse de phase élaborées par le logiciel de traitement de signal à basse cadence ( 10 MHz). On désigne par démodulateur de code le multiplicateur entre le signal reçu démodulé en porteuse, complexe, et un code local. Le produit 5 résultant, complexe, est le signal reçu démodulé désétalé. On désigne par corrélateur complexe l'ensemble formé par un démodulateur de code et un intégrateur (INT) avec remise à zéro périodique (Integrate & Dump). Ici nous avons trois corrélateurs complexes produisant ZA, Zp, AR. Le traitement des deux signaux S1 et S2 (par exemple E5a et 10 E5b) au niveau des étages d'entrée du récepteur dans des voies analogiques distinctes introduisent des temps de propagations et des retards de phases différents sur les deux porteuses (qu'on nommera par la suite différentiels de retard et différentiels de phase) ce qui fait perdre la cohérence entre les deux composantes et a pour effet de réduire la précision 15 des mesures du récepteur (les signaux ne peuvent plus se sommer de manière constructive) et de biaiser les mesures. L'idée est de compenser les défauts introduits sur les signaux reçus par les voies analogiques au niveau du traitement numérique du signal, par une calibration effectuée grâce aux signaux reçus des satellites 20 eux-mêmes, en moyennant et en filtrant les mesures pour réduire les erreurs de calibration. A cet effet, l'invention propose un récepteur de positionnement par satellites comportant au moins un canal de réception (R1, R2,...Ri,...RN), chaque canal de réception étant destiné à effectuer un traitement combiné 25 d'un premier S1 et d'un second S2 signaux radioélectriques émis par un même satellite et séparés en fréquence, les signaux étant reçus par des voies analogiques du récepteur puis numérisés pour être traités dans des voies de réception numériques, caractérisé en ce que chaque canal de réception comporte : 30 - un premier et un second circuits de corrélation recevant respectivement le premier et le second signaux reçus numérisés, chaque circuit de corrélation ayant, une voie de corrélation de porteuse entre le respectif signal reçu et une respective porteuse locale complexe (Lp1, Lp2), des voies de corrélation de code entre le signal reçu corrélé en sortie des 35 voies de corrélation de porteuse et des respectifs codes locaux, un intégrateur par voie de corrélation de code, le récepteur comportant en outre - des discriminateurs de phase de porteuse et de code, les discriminateurs étant communs aux deux circuits de corrélation et fournissant, à partir des signaux en sortie des intégrateurs des deux circuits de corrélation, d'une part, après un correcteur de boucle de code, un signal de vitesse de sous-porteuse et de code (Vc) et, d'autre part, après un correcteur de boucle de porteuse, un signal de vitesse de porteuse centrale (Vp). - un générateur de codes locaux par voie de corrélation de code, chaque générateur de code étant piloté par un signal de phase locale de code en sortie d'au moins un oscillateur à commande numérique NCOc de code commandé par le signal de vitesse de sous-porteuse (ou vitesse de code), les codes locaux fournis par chaque générateur de codes étant décalés d'une correction de retard différentiel positive +Atca, pour un des circuits de corrélation et, d'une correction de retard différentiel négative -Dira, pour l'autre, ces corrections tendant à compenser le différentiel de retard entre les deux signaux reçus dans les voies analogiques du canal de réception ; - un générateur de porteuse locale complexe par voie de corrélation de porteuse, chaque générateur de porteuse locale étant piloté par un signal de phase locale de porteuse en sortie d'au moins un oscillateur à commande numérique NCOp de porteuse commandé par un signal de vitesse de porteuse centrale. Dans l'invention, on construit une seule boucle de porteuse et une seule boucle de code pour traiter les deux signaux reçus S1, S2. Dans une amélioration du récepteur selon l'invention, le récepteur comporte, en outre, un dispositif de correction de différentiel de phase A pua, tendant à compenser le différentiel de phase entre les deux signaux reçus se produisant dans les voies analogiques du canal de réception ; Dans une première réalisation de principe du récepteur selon l'invention, chaque canal de réception comporte en outre : - une voie de corrélation de sous-porteuse par signal reçu avec des sous-porteuses locales complexes générées par un respectif générateur de sous-porteuse piloté par un signal de phase locale de sous-porteuse en sortie de oscillateur à commande numérique NCOc de code, chaque voie de corrélation de sous-porteuse recevant les signaux corrélés de sa respective voie de corrélation de porteuse et fournissant des signaux de corrélation de sous-porteuse à sa respective voie de corrélation de code, le signal de phase locale de sous-porteuse étant décalé de la phase différentielle A(pca, tendant à compenser le différentiel de phase entre les deux signaux reçus se produisant dans les voies analogiques du canal de réception. La boucle de porteuse asservit la phase de la porteuse centrale locale, image de la moyenne des phases des deux porteuses locales, avec la phase de la porteuse centrale des signaux reçus, image de la moyenne des phases de porteuse des deux signaux reçues. La boucle de code (ou boucle de sous-porteuse) asservit simultanément la phase de la sous-porteuse locale, image de l'écart de phase entre les porteuses des deux porteuses locales, avec la phase de la sous-porteuse reçue, image de l'écart de phase entre les deux signaux reçus, et la position des codes locaux avec celle les codes reçus. En effet la sous-porteuse reçue reste en phase avec le code reçue au pied de l'antenne, car ils subissent le même retard de groupe. Le discriminateur de la boucle de code et sous-porteuse utilise bien la sortie des intégrateurs des voies de corrélation avance, ponctuel et retard. Le décalage induit par les voies analogiques distinctes du canal de réception doit être corrigé et à cet effet le récepteur, selon l'invention, comporte en outre un dispositif de calibration fournissant, à partir des valeurs des différentiels de phase et de retard (différentiels dus aux voies analogiques distinctes entre les deux signaux) estimées par les différents canaux du récepteur, les valeurs des corrections de différentiel de phase A(Pca, et de différentiel de retard A'Lca, à appliquer en sortie des oscillateurs à commande numérique NCO de porteuse (ou sous-porteuse) et de code, communes à tous les canaux. Dans une autre réalisation pratique du récepteur selon l'invention comportant N canaux de réception (R1, R2, ...RN), chaque canal de réception comporte un NCOc de code pour chaque composante du signal S1 et S2 et les respectives corrections de retard différentiels, positive +Dira, et négative -Atca,, sont effectuées par application aux entrées des NCOc de code de respectifs signaux de correction +dAtca,/dt et -dAtca,/dt dérivée du signal de retard différentiel +Atca, et -Dira, par rapport au temps t, et en ce qu'on ajoute en sortie des oscillateurs à commande numérique NCOc de code une respective correction +Ato; et -A'ro; propre à chaque canal i qui dépend de l'état initial des oscillateurs NCOc de code au moment où le canal de réception commence la poursuite en mode cohérent du signal. Le récepteur peut comporter en outre un dispositif de correction de différentiel de phase A(pca, recevant d'une part, un signal de correction de phase différentielle AÇ)ca, , d'autre part les signaux en sortie des intégrateurs des deux circuits de corrélation de code et fournissant au discriminateurs des signaux comportant la compensation de différentiel de phase Ocpcai . Cette réalisation du récepteur selon l'invention permet de conserver l'architecture matérielle utilisée dans le cas d'une poursuite indépendante des signaux S1 et S2 (suivant l'état de l'art), ainsi il est facile d'initialiser la poursuite en commençant par une phase de convergence suivant ce mode. L'invention sera mieux comprise à l'aide d'exemples de réalisations du récepteur selon l'invention, en référence aux dessins annexés, dans lesquels : - la figure 1, déjà décrite, représente les spectres d'une paire de 20 signaux émis par un satellite ; - la figure 2 déjà décrite, montre les étages analogiques d'entrée d'un récepteur de position par satellite ; - la figure 3 montre un des deux étages numériques identiques d'un récepteur de positionnement par satellite suivant l'état de l'art ; 25 - la figure 4 représente un schéma de principe d'un récepteur selon l'invention effectuant un traitement combiné d'un premier et d'un second signaux radioélectriques reçus ; - la figure 5 montre un dispositif de détermination des différentiels de phase et de retard du récepteur selon l'invention ; 30 -les figures 6, 7 et 8 montrent des fonctions d'autocorrélation du récepteur selon l'invention ; - la figure 9 représente une réalisation pratique d'un canal de réception du récepteur selon l'invention. La figure 4 représente un schéma de principe d'un canal de réception du récepteur selon l'invention effectuant un traitement combiné d'un premier S1 et d'un second S2 signaux radioélectriques reçus. Le canal de réception de la figure 4 comporte le premier circuit de corrélation 40 attaqué par un premier signal reçu S1 et le second circuit de corrélation 42 attaqué par le second signal reçu S2, chaque circuit de corrélation comportant une voie de corrélation de porteuse centrale 100, 102 entre son respectif signal reçu S1, S2 et deux respectives porteuses locales complexes Lpl pour la première voie de corrélation 100 et Lp2 pour la 1 o deuxième voie 102. Les porteuses locales complexes sont générées par un respectif générateur de porteuse locale 43, 45 piloté par un signal de phase locale cpp en sortie d'un oscillateur commun à commande numérique NCOp 104 de porteuse centrale contrôlé par le signal de vitesse de porteuse Vp en sortie 15 d'un circuit de discriminateurs (106) de phase de porteuse centrale DSP à travers un correcteur de phase de porteuse centrale CRP (107). Les signaux en sortie des voies de corrélation de porteuse centrale sont appliqués à des voies de corrélation de sous-porteuse 108, 110 avec des respectives sous-porteuses locales Lsl , pour la première voie de 20 corrélation de sous-porteuse et Ls2 pour l'autre voie de corrélation, puis aux voies de corrélation de code 112, 114. Les sous-porteuse locales Lsl et Ls2 sont générées par un par un respectif générateur de sous-porteuse locale 111, 113 piloté par un signal de phase locale sous-porteuse q)sp en sortie d'un oscillateur commun à 25 commande numérique NCOc de locale de sous-porteuse et de code 118. L'oscillateur de code NCOc 118 est commandé par un signal de vitesse de code Vc en sortie d'un discriminateur de phase de code DSC du circuit des discriminateurs 106 à travers un correcteur de phase de code CRC 120. 30 Un sornmateur de phase de sous-porteuse 122 reçoit à une entrée ep1 la phase de sous-porteuse en sortie de l'oscillateur de code NCOc 118 et à une autre entrée ep2 une correction de différentielle de phase Ocpca, tendant à compenser le différentiel de phase entre les deux signaux reçus se produisant dans les voie analogiques du récepteur, le sommateur de phase 35 de sous porteuse 122 fournissant, à sa sortie, un signal de phase locale des sous-porteuses décalées en phase pour attaquer le générateur de sous-porteuse des voie de corrélation de sous-porteuse. Les signaux en sortie des voies de corrélation de sous-porteuse 108, 110 sont ensuite corrélés dans les respectives voies de corrélation de code 112, 114 avec (les respectifs codes locaux fournis par les respectifs générateurs de code local 119, 121 pour la voie de signal S1 et la voie de signal S2. Les générateurs de codes locaux 119, 121 sont pilotés par un respectif signal de phase de code (PCA, (PCB en sortie de l'oscillateur de code commun NCOc 118, commandé par le signal de vitesse de sous-porteuse (ou de code) Vc, à travers, un premier sommateur de code local 124 recevant à une entrée ec1 la phase du code local en sortie de l'oscillateur de code local et à une autre entrée ec2 une correction de retard différentiel -ATca, tendant à compenser le retard différentiel, le premier sommateur 124 fournissant à sa sortie une phase de code local décalée pour piloter le générateur de codes local 119 de la voie de signal S1 et, un second sommateur de code local 126 recevant à une entrée ec1 la phase du code local en sortie de l'oscillateur de code local NCOc et à une autre entrée ec2 une correction de retard différentiel +Atca, tendant à compenser le retard différentiel, le second sommateur 126 de code local fournissant à sa sortie une phase de code local décalée pour piloter le générateur de code local 118 de la voie de signal S2. Comme dans le récepteur de la figure 3, les signaux en sortie des voies de corrélation de code sont ensuite intégrés par des respectifs intégrateurs de code INT 130, 132 pour fournir des signaux attaquant les discriminateurs 106 de porteuse DSP et de code DSC du récepteur. La correction du différentiel de phase A(Pca, et du différentiel de retard Dira, est commune à tous les canaux utilisés pour la réception des signaux émis par les satellites eux-mêmes La figure 5 montre un dispositif de détermination de retard différentiel de propagation Atca, et du retard différentiel de phase Acpca, du récepteur de positionnement par satellite, selon l'invention, à partir des satellites reçus. Les paires de fréquences reçues des satellites, par exemple les porteuses E5a et E5b, sont appliquées à N canaux de réception bifréquences R1, R2,...Ri .. RN de même type que celui décrit à la figure 4. Chaque canal récepteur bi-fréquences fournit à un filtre de calibration multicanaux 140 les valeurs estimées de l'erreur de correction (ou calibration) du différentiel de phase EL(pi et de l'erreur de correction (ou calibration) du différentiel de retard Eoti, qui calcule les corrections de différentiel de phase A(Pca, et de différentiel de retard Dira, qui seront appliqués aux sommateurs de phase de sous-porteuse et de code local en sortie des NCO, pour tous les canaux. Chaque canal applique en sortie des corrélateurs des discriminateurs qui permettent d'identifier séparément les défauts d'asservissement des boucles de poursuite et les défauts de calibration. Ces défauts de calibration estimés par chaque canal Eo(pi et EDii, (i étant le rang du canal) entachés d'erreurs propres aux canaux,alimentent un filtre de calibration multicanaux 140 dont le rôle est de moyenner entre les canaux et de filtrer dans le temps ces mesures pour remettre à jour les corrections A pua, et Atca,, en minimisant l'impact des erreurs de mesure sur la précision de la calibration. Nous allons par la suite définir l'expression des signaux aux différents étages du récepteur selon l'invention et à cet effet on défini les signaux reçus au niveau de la numérisation des signaux par les expressions suivantes : 1 - Signaux reçus : 1.a - Forme physique des signaux reçus : Sa reçu = A . cos( (Pa reçu(t) + &Pa ) . Ca( t reçu(t) + bta ) . Da( t reçu(t) ) Sb reçu = A . cos( (Pb reçu(t) + &pb) . Cb( t reçu(t) + Stb) . Db( t reçu(t) ) 1.b - Temps reçu : treçu(t) : Temps auquel a été émis par le satellite le signal reçu par le 35 récepteur à l'instant t (au pied de l'antenne) tpropagation(t) : Temps de propagation entre le satellite et l'antenne du récepteur treçu(t) = t ù tpropagation(t) 1.0 - Phases de porteuse reçues : (Pa reçu(t) Phase du signal reçu à l'instant t (au pied de l'antenne) sur E5a (Pb reçu(t) Phase du signal reçu à l'instant t (au pied de l'antenne) sur E5b &pa : Retard de phase sur la voie E5a dû au récepteur &p,: Retard de phase sur la voie E5b dû au récepteur (Pp reçu(t) = (Op.treçu(t) + (Pdiv(t) (Pa reçu(t) = COa.treçu(t) + (Pdiv(t) (Pa reçu(t) = (Op .treçu(t) - (Osp•treçu(t) + (Pdiv(t) (Pa reçu(t) = (Pp reçu(t) -(osp•treçu(t) (Pb reçu(t) = (Ob.treçu(t) + (Pdiv(t) (Pb reçu(t) = (Op .treçu(t) + Wsp•treçu(t) + (Pdiv(t) (Pb reçu(t) = (Pp reçu(t) + wsp•treçu(t) (Pdiv(t) : divergence entre la phase de la porteuse et le retard de groupe dû à l'ionosphère, qui justifie l'utilisation de boucles de code et de porteuse distinctes. 1.d - Codes d'étalement : Ca(t) : Code d'étalement émis par le satellite à l'instant t sur E5a 30 Cb(t) : Code d'étalement émis par le satellite à l'instant t sur E5b Ca( t reçu(t) ) : Code d'étalement sur E5a reçu par le récepteur à l'instant t (au pied de l'antenne) Cb( t reçu(t) ) : Code d'étalement sur E5b reçu par le récepteur à l'instant t 35 (au pied de l'antenne) 25 Sta : Retard de groupe sur la voie E5a dû au récepteur Stb : Retard de groupe sur la voie E5b dû au récepteur (*) Traitement analogique entre le pied de l'antenne et le convertisseur analogique / numérique 1.e - Données : Da( treçu(t)) : Données reçues à l'instant t sur le signal E5a Db( treçu(t)) : Données reçues à l'instant t sur le signal E5b Dans le cas où le signal dispose d'une voie pilote (c.a.d. sans données) on utilise cette voie (distinguée de la voie de données par le code d'étalement décorrélé). 1.f - Nouvelle expression des signaux reçus : Sa reçu = A . cos( (pp reçu(t) - (osp•treçu(t) + S(Pa ) . Ca( treçu(t) + ôta ) . D( t reçu(t) ) Sb reçu = A . cos( (Pp reçu(t) + Wsp•treçu(t) + S(Pb ) . Cb( treçu(t) + Stb ) . D( t reçu(t) ) 2 -Signaux locaux : On génère des signaux locaux que l'on corrèle avec le signal reçu. La corrélation se fait séparément pour la porteuse locale et le code locale. Cependant cela est mathématiquement équivalent à corréler le signal reçu avec le produit de la porteuse locale et du code local (produit appelé signal local). 2.a - Expression de la forme physique des signaux locaux : Sa local avance = A expi( - (pp(t) ) . expi( + (Psp(t) ) . Ca( (Pca(t) + d ) Sa local ponctuel = A expi( - c)p(t) ) . expi( + (psp(t) ) . Ca( (Pca(t) ) Sa local retard = A expi( - (pp(t) ) . expi( + (psp(t) ) . Ca( (Pca(t) ù d )35 Sb local avance = A . expi( - (pp(t) ) expi( - (psp(t) ) Cb( (Pcb(t) + d ) Sb local ponctuel = A . expi( - (pp(t) ) expi( - (psp(t) ) Cb( (Pcb(t) ) Sb local retard = A expi( - (pp(t) ) . expi( - (Psp(t) ) . Cb( (Pcb(t) ù d ) expi( - (pp(t) ) : Porteuse centrale, complexe (conjuguée) expi( - (psp(t)) : Sous-porteuse locale, complexe (conjuguée) Ca((pca(t)+E.d) : Code locaux avance, ponctuel, retard (E =-1,0,1) voie a, réels Cb((pcb(t)+E.d) : Code locaux avance, ponctuel, retard (E =-1,0,1) voie b, réels Notation : expi(0) = ele où j2 = -1 (pp(t) : Phase de la porteuse centrale locale (rad) (psp(t) : Phase de la sous-porteuse locale (rad) (Pca(t) : Phase du code local sur la voie a (s) (pcb(t) : Phase du code local sur la voie b (s) 1.b - Relation entre ternps reçu estimé et phase locales des codes : treçu estimé(t)Temps reçu estimé en sortie du NCO code (psp(t) = ()sp• treçu estimé(t) + A(Pcal (Pca(t) = treçu estimé(t) - ATcal (Pcb(t) = treçu estimé(t) + Atcal 1.c - Relation entre temps reçu et temps reçu estimé : cp : Déphasage entre la phase locale et la phase reçu de la porteuse centrale i : Retard du temps reçu estimé par rapport au temps reçu vrai (pp(t) _ (Pp reçu(t) - (p treçu estimé(t) = treçu(t) - 'C 1.d - Nouvelle expression des signaux locaux : Sa local ponctuel = A . expi( - (pp reçu(t) + + (Osp.treçu(t) - ()sp.t + A(Pcal ) . Ca( treçu(t) - ti - Otical ) Sb local ponctuel = A . expi( - (Pp reçu(t) + (p - (Osp.treçu(t) + (Osp.t - A(Pcal ) Cb( treçu(t) - ti + Atcal ) 3 - Sortie des corrélateurs : (blocs Integrate and Dump ) ZP = 1/T f [nT, (n+1)T) Sreçu(t) • Slocal ponctuel(t) dt ZA = 1/T ! [nT, (n+1)T] Sreçu(t) • Slocal avance(t) dt ZR = 1/T f [nT, (n+1)T] Sreçu(t) Slocal retard(t) dt (Voie ponctuelle) (Voie avance) (Voie retard) Zp a = "A . expi( () - (Osp.T + A(pcal + ô(Pa ) . R( T + ATcal + ôta ) Zp b ='hA . expi( () + (Osp.T -A(pcal + 8(pb ) . R( T - ATcal + ôtb ) ZA a =1hA . expi( (p - ()sp.T + A()cal + ôcpa ) . R( T + ATcal + ôta - d ) 15 ZA b = 1hA . expi( () + (Osp.T - A(pcal + ô(Pb ) . R( T - Aical + ôtb -d ) ZR a = 1/2A . expi( (p - Wsp.ti + A(pcal + 8()a) . R(T + ATcal + ôta + d ) ZR b "A . expi( (p + (1Jsp.T - A(pcal + &pb) . R(T - ATcal + ôtb + d ) ( Notation : (x+jy)* = (x-jy)* conjugué au sens des nombres complexes ) R : fonction d'autocorrélation du code d'étalement : R( 2 ) = 1/T J [nT, (n+1)T] Ca(t) . Ca(t-ti) dt = 11T f [nT, (-+1)T] Cb(t) . Cb(t-ti) dt R est symétrique : R( T) = R( -T ) 4 - Asservissement idéal : On se place dans un premier temps au point de fonctionnement idéal stabilisé en régime permanent, pour lequel l'asservissement maximise la puissance en sortie de la voie de corrélation ponctuelle, en phase. 35 En régime permanent on veut : (P - wsp•ti - A(pcal + 6()a = 0 Pour redresser Zp a (partie imaginaire nulle) (p + wsp.T + A()cal + &pb = 0 Pour redresser Zp b (partie imaginaire nulle) T - ATcal + ôta = 0 Pour maximiser R( T - ATcal + ôta ) T + ATcal + ôtb = 0 Pour maximiser R( T - ATcal + ôtb ) 20 25 30 (R(i) est maximale pour t = 0) On a quatre équations et quatre inconnues : y), 'r, A(Pcal et Aical Soit 'ro, (Po, Atcal o et A(Pcal 0 les solutions du problème : To = (- 8ta - 8tb )/2 Atcal 0 = ( 8ta - 8tb )/2 (Po = (-&Pa - &Pb )/2 A(Pcal 0 = ( &Pa - 8(Pb )/2 + COsp. ( 8ta + 8tb )/2 Zpa = "A. expi{((p-y)o) wsp•(ti-20) - (A(Pcal-A(Pcal 0 )}. R{(t -to) (Aical-Atical o)} Zpb ='/2A. expi{((ï'-(po) + Cosp.(T-2 o) + (O9cal-09cal 0 )}. R{(2 -To) + (ATcal-O2cal o)} ZAa = ~/2A. expi{(cp-90) COsp.('C-ti0) - (Oc)cal-OcPcal o )}. R{(t+d -2 o) - (02ca1-Otcal o)} ZAb = "A. expi{(y)-90) + COsp.(T-20) + (A9cal-Ay)ca10 )}. R{(i+d -20) + (A2cal-AicalO)} ZRa = ~/2A. expi{(y)-90) COsp.(T-20) -(0(Pcal-O9cal o )}. R{(2-d -'LO) - (Otical-Otical o)} ZRb ='/2A. expi{(y)-90) + wsp.(T-tio) + (Oc)cal-OcQcal o )}. R{(t-d -tio) + (OTcal-Atical o)} 20 5 - Changement de variables : (P-(Po devient (P i-io devient ~ A(Pcal-Ay)cal0 devient Acp 25 Atcal-Atica10 devient Ai Zpa=~/2A.expi{9 -cOsp.i -Ocp}.R{i - Di} Zpb=ihA.expi{yp+OJsp.T+Ay)}. R{ 'r +AT } 3o Zaa=lhA.expi{y) -wsp.T -AcP}.R{i+d - A•r} Zab=1/2A.expi{y)+wsp.i+Acp}. R{i+d+Ai} ZRa=1/2A.expi{cP -oJsp.T -Ay)}.R{T-d-Ai} ZRb=1/2A. expi{y)+wsp.i+Acp} . R{'[-d+Ai} 35 - Voie a + b : On cherche à sommer de manière cohérente les sorties des corrélateurs de voies a et b de manière à avoir le maximum de sensibilité vis-à-vis des erreurs de synchronisation des boucles de poursuite de code et de phase de porteuse. 19 Soit : 10 Zp = Zp a + ZPb ZA = expi( - wsp.d ) . ZA a + expi( +wsp.d ) . ZA b ZR = expi( +wsp.d ) . ZR a + expi( - wsp.d ) . ZR b NB : d est l'espacement entre les corrélateurs avance, ponctuel et retard. 15 Si Ar=0: Zp = A . expi((p) . cos{ wsp.T + Ocp } . R(t) 20 Si dT= 0 et Av.= 0 : Zp = A . expi((p) . cos{ wsp.(T ) } . R(i ) ='hA . expi((p) . Rcos(i ) ZA = A . expi((p) . cos{ wsp.(i+d) } . R(i+d) = Y2A . expi((p) . Rcas(i+d) ZR = A . expi((p) . cos{ (iosp.(i -d) } . R(T -d) ='hA . expi(cp) . Rcos(i -d) Avec Rcos(~) = cos( wsp.i) . R(r ) La figure 6 rnontre la fonction d'autocorrélation Rcos('r ). On a un pic de la fonction d'autocorrélation beaucoup plus 30 prononcé sur Ruos que sur R, ce qui améliore sensiblement la précision de la mesure. En outre on récupère le maximum d'amplitude sur la voie Zp ce qui rend plus précis et plus robuste l'asservissement de la phase de la porteuse. 25 Cela justifie l'objectif de maintenir Ai et A(p à zéro, par le choix de discriminateurs dédiés détectant des valeurs non nulles sur At et Acp et permettant une correction rétroactive par un filtre de calibration. 7 - Discriminateurs (pour chaque satellite reçu) : Phase de porteuse centrale : s(p (c), i, Acp, At) Estimation de cp Retard de code : s(p (cp, i, A(p, Ai) Estimation de 'r Différentiel de phase : so(p (cp, i, Acp, At) Estimation de Acp Différentiel de groupe soi (c), i, A(p, Ai) Estimation de At Le but des ces discriminateurs est d'estimer, chacun par une combinaison des sorties des 6 corrélateurs, les quantités cp, i, Acp, Ati (après changement de repère) qui indiquent ce qu'il manque respectivement sur (p, i, A(pca,, A-rca, (avant changement de repère) pour atteindre l'asservissement idéal donnant le maximum de précision et de robustesse sur les mesures de phase de porteuse centrale et de temps reçu. Cela permet par rétroaction de converger vers l'asservissement idéal. L'objectif est que : s(p (cf), ti, 0, 0) a (pi au premier ordre (c), 0, 0) a 'r au premier ordre Et que : so(p (cf), i, Acp, 0 ) -_ Acp au premier ordre soi (c), i, 0 , Ai) _= Ai au premier ordre 7.a - Discriminateur de phase : s(p = Arg( ZP ) = Arg( ZP a+ZP b )= arctan( lP a+lP b, QP a+QP b) sq) ((II, T, AT, Ai) = Arg[ expi((p) . Ra+b(t, Acp, Ai) l s(p ((p, T, A(p, Ai) = (p + Arg[ Ra+b(ti, Acp, Ai) ] avec Ra+b(t, Act), At) = expi( - Wsp.T - A(p ) . R(i-Ai) + expi( COsp.ti + A(p ) . R(i+Ai) Si AT= 0 on a bien : Ew ((P, i, A(p, 0) = (p car Ra+b(ti, A(p, 0) = cos( wsp.'r + A(p) . R(i) est une fonction réelle Cependant si A(p ≠ 0 Ra+b(ti, A(p, 0) 1 SiAT≠ O ona: Ra+b(2, A(p, A`t) = cos(-wsp.t -Acp) { R(i-AT) + R(ti+Ati) } + i.sin(-wsp.2 -Acp) { -R(ti-A'r) + R('r+Ai) } Ra+b(T, A(p, Ai) n'est plus une fonction réelle La mesure dépend de i donc on a un couplage entre la mesure de code et de phase. Ra+b(0, A(p, Ati) = cos(-A(p) { R(-Ai) + R(+Ai) } + i.sin(-A(p) {-R(-Ati) + R(+Ai) } Ra+b(0, A(p, AT) = 2 cos(-A(p) . R(At) si R est symétrique On a pas de biais pour i = 0 25 7.b - Discriminateur de code : le discriminateur de code utilise la sortie des voies de corrélation avance, ponctuelle et retard sur S1 et S2 (après calibration). E = Re[ (ZA-ZR) . Zp* ] / (Re : partie réelle, * : conjugué) 30 : facteur de normalisation pour avoir une pente unitaire à l'origine Si AT = 0 et A(p = 0 on a bien : ZA - ZR = A . expi(T) . [ Rcos('L+d) - Rcos(ti -d) ] 35 Zp* = A . expi(-T). Rcos(ti) 21 15 20 Et(c), i, 0, 0) = A?/ . [ Rcos(T+d) -Rcos('r-d) ] . Roos(T) 6,((p, 0, 0, 0) = 0 aeT/at((p, 0, 0, 0) _= 2A2/ . Rcos'(d) • Rcos(0) Et (cf), 'r, 0, 0) - 2A2/ . Rcos (d) . Rcos(0) . i ((p, i, 0, 0) 'r pour g = 2A2 Roos'(d) . Roos(0) et 'r 1 La figure 7 rnontre les fonctions Rops(t) et s ((p,'r, 0, 0) , Tchip étant la durée d'un créneau du code d'étalement. Si At=O et Acp;4O on a: &t(cp, i, A(p, 0) = A2/ . [ cos{wsp.('r+d)+A(p}.R(i+d) ù cos{wsp.('r-d)+A p}.R(T-d) ]. cos(wsp.'r+A p).R('r) La figure 8 montre les fonction cos{wsp .i-FA p).R(i) et Et((p,'r, 0, 0) et((), 0, A(p, 0) = A2/ . [ cos(wsp.d+A(p) ù cos(-wsp.d+A p) ] . R(d) . cos(A(p) . R(0) (cp, 0, A(p, 0) A2/ . 2sin(wsp.d) . A(p . R(d) . cos(Acp) . R(0) Ct(c), 0, Ocp, 0) p . A q) pour p = A2/ . 2sin(wsp.d+A p) . R(d). cos(4).R(0) On voit apparaître un biais et en outre on perd de la précision : D'où la nécessité de maintenir Ocp nul. 30 En fait, le biais va se traduire dans la boucle de code par une erreur sur le temps reçu estimé, c'est-à-dire, un 'r non nul. SiDi≠ O et A(p ==0 ona: ~ 0 ~ 5 2 0 2 5 35 E,r((p, 'r, 0, At) = ('/2A)2/ . [ cos{(osp.(i+d)}.{ R(i+d+Di) + R(t+d-Di) } - cos{wsp.(ti-d)}.{ R(t-d+At) + R(T-d-Ai) } ]. cos{wsp.t}.{ R(i+AT) + R(i-AT) } ET((p, 0, 0, AT) = ('hA)2/ . cos(wsp.d).{ R(d+AT) + R(d-At) - R(-d+Ai) - R(-d-Ai) } ].2.R(At) E2((p, 0, 0, At) = 0 car R est symétrique On n'a pas de biais par contre on perd de la puissance. 10 7.c - Discriminateur de différentiel de retard : EAq) = arg[ Zp b . Zp a* ]/2 - wsp . ET 15 Zp b . Zp a* = ('/2A)2 . expi{ 2wsp.i + 2A(p } . R{t -At} . R{r+AT} arg[ Zp b. Zp a* ]/2 = wsp.T + A(p Si At = O on a: aT(c), 'L, A(), 0) E t + [.A%) Eo(v((p, T, A(P, 0) = wsp.T + A(p - wsp . ( T + p.A(p ) E (1 - wsp.p)•A(p pour A(p 1 Variante : 25 On peut prendre : Eoq, = arg[ Zp b . Zp a* ]/2 - wsp . (Ecode b+Ecode a) / r1 (Ecode b+Ecode a)/rl = ~ pour i 1 et A'C 1 EaN((p, T, A(p, 0) = wsp.T + A(p - wsp . (Ecode b+Ecode a)/rl = A(p 7.d - Discriminateur de différentiel de retard : EAT = ( Ecode b - Ecode a ) / r1 1 Ecode a = Re[ (ZA a - ZR a) . Zp a* J 20 30 ZPa=~/2A.expi{cp -wsP.i -4cp}.R{i - Di} ZA a ='/2A . expi{ cp - wsp.i - Ocp } . R{ i+d - AT } ZRa=1hA. expi{(p -wsp.T -Ocp}. R{i -d -AT} Ecode a(c), ti, 0(p, OT) = (1/2A)2 . [ R(i+d - OT) R(i -d - Di) ] . R(i - Di) Ecode a((p, ti, 0(p, OT) = 2(1/2A)2 R'(d) . R(0) . (T-OT) Ecode b = Re[ (ZA b - ZR b) . ZP b* ] / ~1 1 o ZP b = "A . expi{ qp + wsp.i + 4(p } . R{ T + OT } ZAb=1/2A. expi{(p+wsp.i+0(p}. R{i+d+Di} ZRb=1/2A. expi{(p+wsp.i+0(p}. R{2-d+OT} Ecode b((p, 2 , A(p, AT) = ("A)2. [ R(i+d+ Ai) - R(i -d+ 4T) ] . R(T+ Oti) 2("A)2 15 R'(d) . R(0) . (ti+Di) EoT((P, T, A(p, Di) __ 2(1/2A)2/rl . R'(d) . R(0) . OT EAT((p, T, A p, AT) AT 20 pour '1 = 2(1A)2.R'(d).R(0) , T 1 et AT 1 8 - Estimation du différentiel phase pour la calibration : On recale Acpcal à partir des estimations de l'erreur de calibration 25 (0(p = A(pcal - A pcal o ) effectuées sur chaque satellite reçu. Cela permet de moyenner les bruits de mesure entre les satellites, en pondérant les mesures en fonction de leur qualité, qui dépend du rapport signal sur bruit et de l'inclinaison du satellite : (cf. figure 5) 30 0(pcal (k) = 0(pcal (k -1) -• Ei=1à N wi(k) . Eocp i (k) wi(k) = gain de recalage dépendant du satellite dans le canal i considéré (on a N canaux) 35 Par exemple : wi = C/NO, / (Em=1àN C/NOm) . e-dT/T avec C/NOi rapport signal sur bruit estimé sur le canal i Filtre du premier ordre : Xi=1à N wi(k) = e -dT/t période en k avec 'r constante de temps du filtre de calibration et dT EA(f) est une estimée de l'erreur de calibration du différentiel de phase : ESN j (k) = A(Pcal idéal - AcPcal (k -1) + bruit canal i 9 - Estimation du différentiel de retard pour la calibration : Idem. (cf. figure 5) Atcal (k) = ATcal (k-1) - ~i=1à N wi(k) . 4t i (k) wi = gain de recalage dépendant du satellite considéré 10 û Réalisation pratique : La figure 9 représente une réalisation pratique d'un canal de 20 réception du récepteur selon l'invention. L'idée est de conserver l'architecture matérielle d'un canal de réception, qui inclut : - les oscillateurs contrôlés numériquement de code NCOc et de porteuse NCOp, 25 - les générateurs de codes locaux et de porteuse locale et les corrélateurs (multiplicateurs + intégrateur), identique à celle d'un traitement dissocié entre les deux signaux d'entrée (E5a et E5b) conformément à l'état de l'art. On ne modifie que le traitement réalisé à basse cadence en logiciel pour obtenir un récepteur selon l'invention. A cet effet on applique la 30 compensation du différentiel de phase en sortie des corrélateurs, avant les discriminateurs. Le canal de réception de la figure 9 comporte une partie matérielle HD entourée en pointillé dans la figure 9, comportant les deux circuits de corrélation 160, 162 pour les deux signaux S1 et S2, dans cet exemple les15 signaux E5a et E5b. Chacun des circuits de corrélation 160, 162 comporte un intégrateur INT 164, 166, chacun fournissant des signaux : pour le corrélateur 160 de porteuse E5a : ZAa = lAa + IQAa Zpa = 1pa + IQPa ZRa = 'Ra + IQRa pour le corrélateur 162 de porteuse E5b : ZAb = lAb + IQAb Zpb = 'Pb + IQPb ZRb = lRb + IQRb Correspondant au code avance A, ponctuel P et retard R ; Une compensation du différentiel de phase en sortie des corrélateurs des deux circuits est appliquée avant les discriminateurs 181 par un module de compensation 182. Les discriminateurs fournissent un signal d'erreur de code > et un signal d'erreur de phase de porteuse E,p pour l'asservissement des boucles, ainsi qu'un signal d'erreur de calibration du différentiel de phase ,,gyp et un signal d'erreur de calibration du différentiel de retard eA, pour le filtre de calibration multi-canaux. On applique la compensation du différentiel de retard en entrée des NCOc de code des circuits de corrélation 160, 162, en dérivant la correction Atcal par rapport au temps soit dAtcal/dt, de sorte qu'en sortie des oscillateurs NCOc de code on obtienne la correction désirée. A cet effet le canal de réception de la figure 9 comporte : -un sommateur Spi recevant par une entré el un signal de vitesse de code après le correcteur de boucle de code CRC 180 et, par une entré e2 le signal de correction dérivé -dAtcal/dt, et fournissant à une sortie le signal somme des signaux appliquées à ses deux entrées el et e2 pour fournir la commande de l'oscillateur de code NCOc 190 du circuit de corrélation recevant la porteuse E5a et, - un sommateur Sp2 recevant par une entré el le signal de vitesse de code après le correcteur de boucle de code CRC 180 et, par une entrée e2 le signal de correction dérivé dAtca,/dt, et fournissant à une sortie le signal somme des signaux appliquées à ses deux entrées pour fournir la commande de l'oscillateur de code NCOc 192 de l'autre circuit de corrélation recevant la porteuse E5b. En outre on ajoute aux sorties des deux oscillateurs NCOc de code 190, 192, attaquant les respectifs générateurs de code 191, 193, une correction Aloi propre au canal i qui dépend de l'état initial des NCOc au moment où on commence la poursuite en mode cohérent. A cet effet, le canal de réception comporte aussi : - un sommateur Sp3 recevant par une entré el le signal en sortie de l'oscillateur de code NCOc 190 et, par une entrée e2 une correction -ATO; propre à chaque canal et, fournissant à une sortie le signal somme des signaux appliquées à ses deux entrées el et e2 pour fournir la phase des générateurs de code (191) du circuit de corrélation recevant la porteuse E5a et ; - un sommateur Sp4 recevant par une entrée el le signal en sortie de l'oscillateur de code NCOc 192 et, par une entré e2 une correction +Aio; propre à chaque canal i et fournissant à une sortie le signal somme des signaux appliquées à ses deux entrées el et e2 pour fournir la phase des générateurs de code (193) de l'autre circuit de corrélation recevant la porteuse E5b. Le canal de réception comporte en outre, par circuit de corrélation, un oscillateur NCOp 194, 195 de porteuse centrale piloté à partir du signal de vitesse de porteuse centrale Vp fournie par un correcteur de boucle de porteuse CRP 196 alimenté par les sorties des discriminateurs 181. La différence de réalisation entre la boucle de phase de porteuse (compensation après la corrélation) et la boucle de code (compensation avant la corrélation sur les phases des codes locaux) tient au fait qu'une correction de 8 sur la phase de la porteuse locale équivaut à une multiplication par eje de la porteuse locale complexe ou d'une multiplication par el de la sortie des corrélateurs complexes. Cela n'est pas généralisable au code car une modification de la phase du code local se traduit par une variation d'amplitude (suivant une loi non linéaire en fonction du retard) sur les corrélateurs (et non pas une rotation complexe). Pour chaque satellite on prend la moyenne des mesures : (Pp i = ((Ppa + (Ppb )/2 treçu estimé i = ((Pea '+ (Pc) )/235 10.a - Initialisation de la poursuite cohérente : On commence par acquérir et poursuivre chacun des deux signaux indépendamment tel que dans le récepteur représenté à figure 2, conformément à l'état de l'art. Une fois que chacun des asservissements (PLL et DLL) a convergé sur chaque satellite, on passe en mode de poursuite cohérent selon le récepteur représenté à la figure 4 avec calibration, en réorganisant le calcul des commandes en vitesse appliquées au NCO, suivant le schéma du récepteur de la figure 9. On ne change pas les valeurs des phases en sortie des NCO de porteuse, mais on les corrige en sortie des corrélateurs. Par contre, il faut corriger les phases en sortie des oscillateurs NCOc de code pour se conformer au schéma de la figure 4 dans lequel on a un seul NCOc code. Pour cela on applique une correction à chacun des NCO code, de signe contraire, dont le but est de supprimer les écarts de code et d'initialiser la correction de calibration (pilotée ensuite en vitesse de variation dAtcal / dt). 10.b - Initialisation de la calibration (pour l'ensemble des satellites) : A l'issue de la phase de convergence des boucles de poursuite séparées, on observe différentiels de retard Atcal(0) et de phase A(pcal(0) entre les deux voies analogiques E5a et E5b directement sur les phases des porteuses et codes locaux, moyennées sur tous les satellites : 25 O(Pcal(0) = Ei=1à N Wi(0) • [ ( (Ppb 0 i -qppa 0 i )/2 - wsp•( (Nb 0 i + (Pca 0 i )/2 ] &cal (0) = 1i=1à N Wi(0) . [ ( (Qcb 0 i - (Na 0 i )/2 ] 30 Avec : Ei=1à N Wi(k) = 1 (On retrouve une analogie avec les discriminateurs de calibration du différentiel de phase et de retard) Le filtre de calibration 140 est ensuite entretenu pour assurer un 35 suivi des différentiels de phase A(pcal et de retard &cal entre les voies E5a et E5b, qui peuvent dériver avec le temps et ou évoluer avec les variations de température. 5 10.c Compensation du différentiel de phase : (En sortie des corrélateurs, par logiciel) Soit i = indice du canal On corrige les sorties des corrélateurs par une rotation dépendant du déphasage observé sur le canal (à la fin de la phase de poursuite séparée sur E5a et E5b) et de la correction de calibration AcPcal issue de l'ensemble des canaux. 10 jB,cp ZA a i compensé = Ai . e + cal(0) • ZA a i jB,cp ZP a i compensé = Ai . e + cal(0) Z • P a i ZRaicompensé =Ai . e + cal(0) ZRai ZA b i compensé = Bi . e jA(p cal(0) . ZA b i jAcp ZP b i compensé = Bi . e cal(0) . ZP b i ZR b i compensé = Bi . e jA(p cal(0) . ZR b i 15 Ai = expi( - [ ( (PO 0 i - (Ppa 0 i )/2 - wsp.( (Pcb 0 i + (Pca 0 i )/2 ] ) Bi=expl(-I-[ ((ppb0i-(Ppa0i)/2-wsp•((pcb0i+(Pca0i)/2]) Ainsi juste après la correction : ZP ai compensé = Ai . e +jAcp cal(0) "A . expi( - (ppa o i + (Pa reçu i + &Pa) • R( • • • ) ZP b i compensé = Bi . e -1A'p cal(0) "A . expi(- (Ppb 0 i + (Pb reçu i + &Pb) • R( • • • ) ZP a i compensé = 1/2A. explL((ppb 0 i+(ppa 0 i)/2 + OJsp. ((Pcb 0 i+(Pca 0 i)/2 + jA(Pcal(0) +(Pa reçu i +S(Pa ].R(...) ZP b i compensé = 25 12A. expl[((Ppb o i+(Ppa 0 i)/2 - wsp•((Pcb 0 i+(Pca 0 i)/2 - j0cpcal(0) + (Pb reçu i +S(Pb]•R(•••) Ce qui est équivalent à : ZP ai compensé = "A . expi[ - (Ppa virtuel i + (Pa reçu i + &Pa ] • R( ) ZP b i compensé = "A . expi[ - (Ppb virtuel i + (Pb reçu i + 8(Pb ] R( . ) Avec : (Ppa virtuel i=(Pp 0 i + wsp •(Pc 0 i + OcPcal(0) (Ppb virtuel i=(pp 0 i wsp •(Pc 0 i - 0cpcal(0) 20 30 35 (Pp o i = ((Ppb 0 i+(Ppa 0 i)/2 (Pc o i= ((Pcb o i+(Pc:a o i)/2 (Pp 0 et cpsp 0 dépendent du canal (donc du satellite) alors que A(Pcal(0) est commune à tous les canaux Conclusion : On retrouve une expression conforme au schéma de la figure 4. 10.4 - Compensation du différentiel de retard : (En sortie des NCOc code) Initialisation : Correction initiale : Or, 0 i = - ( (Pcb 0 i -(Pca 0 i )/2 + Aical (0) 15 Cette correction, propre à chaque canal est appliquée en sortie des NCO de code, conformément au schéma de la figure 4 (elle pourrait aussi être induite en jouant sur les commandes NCO, par glissement, à I'initialisation) à I'initialisation. 20 Ainsi juste après la correction : (Na i=(Pca 0 i- Di o i=((Pcb 0 i+(Pca 0 i )/2 - Atical(0) qPcb i=(Pcb 0 i+ Ai o i= ( (Pcb 0 i+(Pca 0 i )/2 + Otical(0) Ce qui est équivalent à : (Pca i = (Pc 0 i - Atca1(0) (Pcb i = (Pc 0 i + Atcal(0).. DTD: cpc o i dépend du canal (donc du satellite) alors que Atcal(0) est commune à tous les canaux 30 Conclusion : On retrouve une expression conforme au schéma de la figure 3. Remarque : 35 Dans le cas idéal où pour chaque canal on a : 30 10 25 ( (Ppb i - (Ppa i )/2 - 'cosp•( (Pcb i + (Pca i )/2 = A(Pcal (0) ( (Pcb i - (Pca i )12 = Atcai (0) (autrement dit tous les termes de gauches seraient identiques, sur tous les canaux) On obtiendrait : Ai . e - jA p cal(o) = 1 Bi . e +jAcp cal(0) = 1I Di p i =0 Donc la calibration initiale n'apporterait rien ! Ce résultat n'est en fait pas surprenant dans la mesure où à l'issue de la phase de convergence initiale on aurait dans ce cas déjà atteint 15 l'objectif de la calibration : chaque canal s'est positionné au maximum d'énergie par les boucles de poursuite initiales sur E5a et sur E5b, et les valeurs Atcal(0) i et A(pcal(0) i du différentiel de retard et du différentiel de phase (respectivement) observées en sortie des NCO sont identiques sur tous les canaux. 20 Entretien de la compensation : Après l'initialisation on entretien la correction en ajoutant un complément aux commandes en vitesse des deux NCO code, conformément au schéma de la 25 figure 5 : dATcal/dt = [ Atcal (k+1) - ATcal (k) ] / dT avec dT : pas d'échantillonnage en k 30 10.5 - Calibration de l'amplitude : Les gains des deux voies analogiques peuvent être différents. Il est important de compenser ces différences lorsque l'on construit les discriminateurs, pour éviter les biais dus au déséquilibre entre les voies. Pour 35 cela on identifie le rapport entre les puissances en sortie des corrélateurs ponctuels des voies a et b : 10 a•(0) _ zi=1à N w,(0) ( I IZP bI I2 i I IZP aI I2 )i X(k) = X(k-1) -li=1à N wi(k) • [ (I IZP bI I2 / IIZP aI I2)i - %,(k-1) ] On applique la compensation en sortie des corrélateurs, avant les discriminateurs : Zp = %,. Zp a + Zp b ZA = a, .e+iwd ZA a + e- iwd Z Ab ZR = %, .e-iwd ZR a + e+ i id Z Rb On suppose que les erreurs de boucle et de calibration sont résorbées. En retour la compensation de la différence d'amplitude réduit 15 l'impact du déséquilibre d'amplitude sur les discriminateurs de boucle et de calibration. 10.6 - Traitement des données 20 Cas où le signal ne dispose pas d'une voie pilote (sans données) : Dans ce cas on doit redresser les bits de données avant d'appliquer les discriminateurs : 25 ZA a redressé = ZA a • signe( lm[ Zp a ] ) ZA b redressé = ZA b • signe( lm[ ZP b ] ) ZP a redressé = ZP a • signe( lm[ Zp a 1 ) ZP b redressé = ZP b • signe( lm[ ZP b ] ) ZR a redressé = ZR a . signe( lm[ Zp a 1 ) ZR b redressé = ZR b • signe( Im[ Zp b ] ) ( Im : partie réelle) 30 Cas où on souhaite utiliser les voies pilote et donnée : On applique les discriminateurs de sur la voiepilote uniquement, sauf le discriminateur de code qui utilise les deux voies : i0 35 = { Re[ (ZA pilote -ZR pilote) . ZP pilote* ] + Re[ (ZA données -ZR données) • ZP données* ] } / 2 / Cela permet de réduire le bruit thermique sur la mesure de code, à 5 condition que le rapport signal sur bruit soit largement positif dans la bande de prédétection (Bprédetection = 1/T). Le récepteur selon l'invention permet une meilleure précision de mesure du positionnement par une moindre sensibilité de la mesure de code vis-à-vis du bruit thermique et des multi-trajets. Il présente en outre une 10 meilleure robustesse (moins de saut de cycles et décrochage de la boucle de phase de porteuse du aux interférences). Cette technique de correction de phase et délai différentiel se produisant dans les étages analogiques de réception pour deux porteuses s'applique aussi lorsque le filtre analogique est commun aux signaux a et b 15 car il peut y avoir des déphasages et retards entre les deux composantes du signal, dus aux distorsions des fonctions de transfert des filtres analogiques. REVENDICATIONS 1. Récepteur de positionnement par satellites comportant au moins un canal de réception (R1, R2,...Ri,...RN), chaque canal de réception étant destiné à effectuer un traitement combiné d'un premier S1 et d'un second S2 signaux radioélectriques émis par un même satellite et séparés i o en fréquence, les signaux étant reçus par des voies analogiques du récepteur puis numérisés pour être traités dans des voies de réception numériques, caractérisé en ce que chaque canal de réception comporte : - un premier (40, 160) et un second (42, 162) circuits de corrélation recevant respectivement le premier et le second signaux reçus 15 numérisés, chaque corrélateur ayant, une voie de corrélation de porteuse (100, 102) entre le respectif signal reçu et une respective porteuse locale complexe (Lpl, Lp2), des voies de corrélation de code (112, 114) entre le signal reçu corrélé en sortie des voies de corrélation de porteuse et des respectifs codes locaux, un intégrateur (130, 132, 164, 166) par voie de 20 corrélation de code, le récepteur comportant en outre : - des discrirninateurs de phase de porteuse et de code (106, 181), les discriminateurs étant communs aux deux circuits de corrélation et fournissant, à partir des signaux en sortie des intégrateurs des deux circuits de corrélation, d'une part, après un correcteur de boucle de code, un signal 25 de vitesse de sous-porteuse et de code (Vc) et, d'autre part, après un correcteur de boucle de porteuse, un signal de vitesse de porteuse centrale (Vp) ; - un générateur de codes locaux (119, 121, 191, 193) par voie de corrélation de code, chaque générateur de code étant piloté par un signal de 30 phase locale de code en sortie d'au moins un oscillateur à commande numérique NCOc de code (118, 190, 192) commandé par le signal de vitesse de sous-porteuse (ou vitesse de code), les codes locaux fournis par chaque générateur de codes étant décalés d'une correction de retard différentiel positive +Atca, pour un des circuits de corrélation et, d'une 35 correction de retard différentiel négative -Amal pour l'autre, ces correctionstendant à compenser le différentiel de retard entre les deux signaux reçus dans les voies analogiques du canal de réception ; - un générateur de porteuse locale complexe (43, 45) par voie de corrélation de porteuse, chaque générateur de porteuse locale étant piloté par un signal de phase locale de porteuse en sortie d'au moins un oscillateur à commande numérique NCOp de porteuse (104, 194, 195) commandé par un signal de vitesse de porteuse centrale. 2. Récepteur de positionnement par satellites selon la revendication 1, caractérisé en ce qu'il comporte, en outre, un dispositif de correction de différentiel de phase A(Pca, tendant à compenser l'écart différentiel de phase entre les deux signaux reçus S1, S2 se produisant dans les voies analogiques du canal de réception. 3. Récepteur de positionnement par satellites selon l'une des revendications 1 ou 2, caractérisé en ce que chaque canal de réception comporte en outre : - une voie de corrélation de sous-porteuse par signal reçu S1, S2 avec des sous-porteuses locales complexes générées par un respectif générateur de sous-porteuse (111, 113) piloté par un signal de phase locale de sous-porteuse en sortie de oscillateur à commande numérique NCOc de code (118), chaque voie de corrélation de sous-porteuse recevant les signaux corrélés de sa respective voie de corrélation de porteuse (100, 102) et fournissant des signaux de corrélation de sous-porteuse à sa respective voie de corrélation de code (112, 114), le signal de phase locale de sous-porteuse étant décalé de la phase différentielle A(pca, tendant à compenser le différentiel de phase entre les deux signaux reçus se produisant dans les voies analogiques du canal de réception. 4. Récepteur de positionnement par satellites selon l'une des revendications 1 à 3, caractérisé en ce qu'il comporte en outre un dispositif de calibration fournissant, à partir des valeurs des différentiels de phase et de retard estimées par les différents canaux (R1, R2,...RN) du récepteur, les valeurs des corrections de différentiel de phase A(pca, et de différentiel de retard &cal à appliquer en sortie des oscillateurs à commande numériqueNCO de porteuse (ou sous-porteuse) et de code communes à tous les canaux. 5. Récepteur de localisation selon l'une des revendications 1 à 4, caractérisé en ce que chaque canal comporte le premier circuit de corrélation (40) attaqué par un premier signal reçu S1 et le second circuit de corrélation (42) attaqué par le second signal reçu S2, chaque circuit de corrélation comportant un voie de corrélation de porteuse centrale (100, 102) entre son respectif signal reçu S1, S2 et deux respectives porteuses locales complexes Lpl pour la première voie de corrélation (100) et Lp2 pour la deuxième voie (102), ces porteuses locales étant générées par un respectif générateur de porteuse locale (43, 45) piloté par un signal de phase locale cpp en sortie d'un oscillateur comrnun à commande numérique NCOp (104) de porteuse centrale contrôlé par le signal de vitesse de porteuse Vp en sortie d'un circuit de discriminateurs (106) de phase de porteuse centrale DSP à travers un correcteur de phase de porteuse centrale CRP (107), les signaux en sortie des voies de corrélation de porteuse centrale étant appliqués à des voies de corrélation de sous-porteuse (108, 110) avec des respectives sous-porteuses locales Lsl;, pour la première voie de corrélation de sous-porteuse et Ls2 pour l'autre voie de corrélation, puis aux voies de corrélation de code (112, 114). 6. Récepteur de localisation selon la revendication 5, caractérisé en ce que les sous-porteuse locales Lsl et Ls2 sont générées par un respectif générateur de sous-porteuse locale (111, 113) piloté par un signal de phase locale sous-porteuse (psp en sortie d'un oscillateur commun à commande numérique NCOc de locale de sous-porteuse et de code (118) commandé par un signal de vitesse de code Vc en sortie d'un discriminateur de phase de code DSC du circuit des discriminateurs (106) à travers un correcteur de phase de code CRC (120) et un sommateur de phase de sous-porteuse (122) recevant, à une entrée ep1, la phase de sous-porteuse en sortie de l'oscillateur NCOc de code (118) et, à une autre entrée ep2, une correction de différentiel de phase A pua, tendant à compenser le différentiel de phase entre les deux signaux reçus se produisant dans les voies analogiques du récepteur, le sommateur de phase de sous porteuse (122)fournissant, à sa sortie, un signal de phase locale des sous-porteuses décalées en phase pour attaquer le générateur de sous-porteuse des voies de corrélation de sous-porteuse. 7. Récepteur de localisation selon la revendication 6, caractérisé en ce que les signaux en sortie des voies de corrélation de sous-porteuse (108, 110) sont ensuite corrélés dans les respectives voies de corrélation de code (112, 114) avec les respectifs codes locaux fournis par les respectifs générateurs de code local (119, 121) pour la voie de signal S1 et la voie de signal S2, les générateurs de code local (119, 121) étant pilotés par un respectif signal de phase de code (PCA, (PCB en sortie de l'oscillateur de code commun NCOc (118) commandé par le signal de vitesse de sous-porteuse (ou de code) Vc, à travers, un premier sommateur de code local (124) recevant à une entrée ec1 la phase du code local en sortie de l'oscillateur de code local et à une autre entrée ec2 une correction de retard différentiel -Dira, tendant à compenser le retard différentiel, le premier sommateur (124) fournissant à une sortie une phase de code local décalée pour piloter le générateur de codes local (119) de la voie de signal S1 et, un second sommateur de code local (126) recevant à une entrée ec1 la phase du code local en sortie de l'oscillateur de code local NCOc et à une autre entrée ec2 une correction de retard différentiel +Atca, tendant à compenser le retard différentiel, le second sommateur (126) de code local fournissant à une sortie une phase de code local décalée pour piloter le générateur de code local (118) de la voie de signal S2. 8. Récepteur de positionnement par satellites selon l'une revendications 1 ou 2, caractérisé en ce que chaque canal de réception comporte un NCOc de code (190, 192) pour chaque composante du signal S1 et S2 et en ce que les respectives corrections de retard différentiels, positive +Dira, et négative -A'tca,, sont effectuées par application aux entrées des NCOc de code de respectifs signaux de correction +dA'rca,/dt et -dAtca,/dt dérivée du signal de retard différentiel +Dira, et -Dira, par rapport au temps t, et en ce qu'on ajoute en sortie des oscillateurs à commande numérique NCOc de code (190, 192) une respective correction +A'o; et -ATo; propre à chaque canal i qui dépend de l'état initial des oscillateurs NCOc de code aumoment où le canal de réception commence la poursuite en mode cohérent du signal. 9. Récepteur de positionnement par satellites selon l'une des revendications 5 à 8, caractérisé en ce que chacun des circuits de corrélation (160, 162) du canal de réception comporte un intégrateur INT (164, 166), chacun fournissant des signaux : pour le corrélateur (160) de porteuse E5a : ZAa = lAa + IQAa Zpa = Ipa + IQPa ZRa=IRa+IQRa, pour le corrélateur (162) de porteuse E5b : ZAb = lAb + IQAb Zpb = Ipb + IQPb ZRb = IRb + iQRb correspondant aux codes avance A, ponctuel P et retard R ; 10. Récepteur de positionnement par satellites selon l'une des revendications 5 à 10, caractérisé en ce qu'une compensation du différentiel de phase en sortie des corrélateurs des deux circuits est appliquée avant les discriminateurs (181) par un module de compensation (182), les discriminateurs fournissant un signal d'erreur de code c,p et un signal d'erreur de vitesse de phase s,P pour l'asservissement des boucles. 11. Récepteur de localisation selon l'une des revendications 5 à 10, caractérisé en ce que chaque canal de réception comporte : -un sommateur Spi recevant par une entré el un signal de vitesse de code après une correction de vitesse de code CRC (180) et, par une entré e2 le signal de correction dérivé -dAtca,/dt, et fournissant à une 3o sortie le signal somme des signaux appliquées à ses deux entrées el et e2 pour attaquer la comrnande de oscillateur de code NCOc (190) du circuit de corrélation recevant la porteuse E5a et, -un sommateur Sp2 recevant par une entré el le signal de vitesse de code Et après la correction de vitesse de code CRC (180) et, par une 35 entrée e2 le signal de correction dérivé dAtcal/dt, et fournissant à une sortie lesignal somme des signaux appliquées à ses deux entrées pour attaquer la commande de oscillateur de code NCOc (192) de l'autre circuit de corrélation recevant la porteuse E5b ; -un sommateur Sp3 recevant par une entré el le signal en sortie de l'oscillateur de code NCOc (190) et, par une entrée e2 une correction -A'r propre à chaque canal et, fournissant à une sortie le signal somme des signaux appliquées à ses deux entrées el et e2 pour fournir la phase des générateurs de code (191) du circuit de corrélation recevant la porteuse E5a et ; -un sommateur Sp4 recevant, par une entrée el, le signal en sortie de l'oscillateur de code NCOc (192) et, par une entré e2, une correction +At0 propre au canal i et fournissant à une sortie le signal somme des signaux appliquées à ses deux entrées el et e2 pour fournir la phase des générateurs de code (193) de l'autre circuit de corrélation recevant la porteuse E5b. - par circuit de corrélation, un oscillateur NCOp (194, 195) de porteuse centrale piloté à partir du signal de vitesse de porteuse centrale Vp fournie par un correcteur de boucle de porteuse CRP (196) alimenté par les sorties des discrirninateurs (181). 12. Récepteur de positionnement par satellites selon l'une des revendications 1 à 11, caractérisé en ce que le discriminateur de code utilise la sortie des voies de corrélation avance, ponctuelle et retard sur S1 et S2 (après calibration).25