Brevet FR3002388A1 — publié le 20140822

MÉTHODE D'ACQUISITION D'UN SIGNAL GPS PAR DÉCODAGE ITÉRATIF DESCRIPTION DOMAINE TECHNIQUE La présente invention concerne de manière générale le domaine des systèmes de positionnement par satellites tels que le système GPS (Global Positioning System) ou le système Galileo. Elle s'applique plus particulièrement à la phase d'acquisition d'un signal satellitaire par un récepteur GPS. ÉTAT DE LA TECHNIQUE ANTÉRIEURE Un récepteur GPS détermine classiquement sa position en mesurant les temps de vol de signaux reçus d'une pluralité de satellites (en principe au moins quatre). Pour estimer le temps de vol d'un signal GPS émis par un satellite, le récepteur utilise un canal d'acquisition grossière dit canal C/A (Coarse Acquisition) sur lequel le satellite émet un message de navigation étalé spectralement au moyen d'une séquence pseudo-aléatoire (dite séquence d'étalement) possédant de bonnes propriétés d'auto-corrélation ainsi que d'inter-corrélation (avec les séquences PRN utilisées par les autres satellites). Les bonnes propriétés d'auto-corrélation des séquences d'étalement permettent d'optimiser la probabilité de détection d'un signal en présence de bruit alors les bonnes propriétés d'inter-corrélation de ces séquences permettent aux différents satellites de partager une même bande de transmission au moyen d'un accès multiple par répartition de codes (CDMA). La première tâche d'un récepteur GPS est de détecter les satellites en visibilité. Pour ce faire, le récepteur cherche à détecter le code d'étalement du satellite dont il souhaite tester la présence et, le cas échéant, se synchronise par rapport à ce signal. Cette première tâche est dénommée étape d'acquisition du signal GPS. De nombreuses méthodes d'acquisition du signal GPS sont connues de l'état de la technique, elles s'apparentent à la problématique de détection multi-utilisateur dans les systèmes CDMA. On peut notamment citer les articles de A. Polydoros et al. intitulés «A unified approach to serial search spread spectrum code acquisition » Parti : A general theory & Part II : a matched filter receiver, publiés dans IEEE Trans. on. Comm., Vol. 32, No. 5, 1984 et l'article de F. Principe et al. intitulé « Rapid acquisition of Gold codes and related sequences using iterative message passing on redundant graphical models » publié dans Proc. MILCOM '06, 2006. Ces méthodes d'acquisition sont généralement basées sur une corrélation du signal reçu par le récepteur avec le code d'étalement du satellite à tester. Toutefois, la méthode d'acquisition décrite dans l'article de F. Principe précité ne permet pas d'obtenir des niveaux de probabilité de détection ratée (missed detection probability) et de probabilité de fausse alarme suffisamment faible à bas rapport signal sur bruit (SNR). Par probabilité de détection ratée, on entend la probabilité de ne pas se synchroniser avec une séquence d'étalement d'un satellite alors que celle-ci est effectivement reçue par le récepteur. Par probabilité de fausse alarme, on entend la probabilité de se synchroniser avec une séquence d'étalement alors que celle-ci n'est pas reçue par le récepteur. A faible niveau de rapport signal sur bruit, la probabilité de fausse alarme est celle de détecter une séquence d'étalement dans du bruit. En outre, les méthodes d'acquisition par corrélation nécessitent de tester un à un les différents codes d'étalement possibles, ce qui se traduit, dans une configuration série, par un temps d'acquisition important et, dans une configuration parallèle, par une complexité élevée. Le but de la présente invention est de proposer une méthode d'acquisition d'un signal satellitaire, en particulier d'un signal GPS, permettant de détecter une séquence d'étalement avec une faible probabilité de détection ratée, même à bas niveau de rapport signal sur bruit. Un but subsidiaire de la présente invention est de proposer une méthode d'acquisition d'un signal satellitaire avec une faible probabilité de fausse alarme, même à bas niveau de rapport signal sur bruit. Un autre but de la présente invention est de proposer une méthode d'acquisition d'un signal satellitaire présentant un meilleur compromis temps d'acquisition/complexité que dans l'art antérieur. EXPOSÉ DE L'INVENTION La présente invention est définie par une méthode d'acquisition d'un signal satellitaire par un récepteur, ledit signal étant étalé par une séquence de Gold à une fréquence bribe, ladite séquence de Gold étant obtenue comme somme d'une première M-séquence générée au moyen d'un premier registre à décalage et d'une seconde M-séquence générée au moyen d'un second registre à décalage, ladite seconde M- séquence identifiant un satellite parmi une pluralité de satellites, ladite méthode comprenant les étapes suivantes : (a) le signal reçu en bande de base est échantillonné à la fréquence bribe à partir d'un instant déterminé pour fournir une séquence de premiers échantillons; (b) la séquence desdits premiers échantillons est multipliée, bribe à bribe, par ladite première M-séquence pour fournir une séquence de seconds échantillons; (c) la séquence de seconds échantillons est décodée au moyen d'un décodage itératif à passage de messages pour fournir une estimation du contenu du second registre à décalage au dit instant déterminé ; (d) le contenu du second registre à décalage ainsi estimé est comparé à une pluralité de contenus possibles du second registre à décalage pour les différents satellites de ladite pluralité, et en cas de succès de la comparaison pour un satellite identifié, le récepteur est synchronisé avec ledit satellite identifié à l'instant déterminé. Si le contenu du second registre à décalage ainsi estimé ne correspond à aucun desdits contenus possibles, l'instant déterminé est incrémenté d'une bribe et les étapes (a) à (d) sont itérées. Selon un premier mode de réalisation, l'estimation du contenu du second registre à décalage est avantageusement réalisée en prenant des décisions dures sur les r premières valeurs décodées fournies par le décodage itératif à passage de messages, OÙ r est le nombre de positions dans le premier/second registre à décalage. De préférence, le décodage itératif à passage de messages de l'étape (c) utilise la matrice de parité : gr - - - go 0 0 gr - - go 0 H= O 0 gr - - - go 0 0 gr go dans laquelle gr =1 etgo''gri sont les coefficients du polynôme générateur de la seconde M-séquence. La seconde M-séquence est avantageusement obtenue en sommant les contenus en deux positions a' du second registre à décalage, où i est un indice identifiant le satellite parmi ladite pluralité de satellites, les contenus possibles du second registre à décalage pour les différents satellites étant donnés par les vecteurs : A;nd = (Gya, + G pyi)Aly'ritt où A' est un vecteur donnant contenu du second registre à décalage à un instant de Y référence, et Gy la matrice de changement d'état du second registre à décalage. Selon un second mode de réalisation, la séquence de seconds échantillons est permutée à l'aide d'une relation de permutation ço(k)=- dk + h mod(N) où k est le rang de l'échantillon, d est un facteur de décimation, h est une valeur d'offset et N est la longueur de la première M-séquence, la permutation étant effectuée avant l'étape (c) ; l'étape (c) fournit une séquence de valeurs décodées ; la séquence de valeurs décodées est permutée à l'aide de la relation de permutation inverse à ladite relation de permutation pour fournir une séquence de valeurs permutées; l'estimation du contenu du second registre à décalage est réalisée en prenant des décisions dures sur les r premières valeurs de la séquence de valeurs permutées, où r est le nombre de positions dans le premier/second registre à décalage. Avantageusement, le décodage itératif à passage de messages de l'étape (c) utilise la matrice de parité : go 0 0 g --- go 0 0 H 0 -- - 0 g,. --- go 0 0 gr go dans laquelle gr =1 et sont les coefficients du polynôme générateur de la première M-séquence. Selon une variante des premier et second modes de réalisation, on vérifie que le récepteur est bien synchronisé avec ledit satellite identifié en : (e) multipliant la séquence desdits premiers échantillons, bribe à bribe, par la seconde M-séquence associée audit satellite identifié pour fournir une séquence de troisièmes échantillons ; (f) décodant la séquence de troisièmes échantillons au moyen d'un décodage itératif à passage de messages pour fournir une estimation du contenu du premier registre à décalage au dit instant déterminé. Si le contenu du premier registre à décalage ainsi estimé ne correspond au contenu initial du premier registre, l'instant déterminé est incrémenté d'une bribe et les étapes (a) à (f) sont itérées. De préférence, l'estimation du contenu du premier registre à décalage est réalisée en prenant des décisions dures sur les r premières valeurs décodées fournies par le décodage itératif à passage de messages, où r est le nombre de positions dans le premier/second registre à décalage. Avantageusement, le décodage itératif à passage de messages de l'étape (f) utilise la matrice de parité 0 - - - 0 gr - - - go 0 g, --- g0 0 0 H 0 - - - 0 g, --- go 0 0 O g- go dans laquelle gr =1 etgo''gri sont les coefficients du polynôme générateur de la première M-séquence. Selon un exemple avantageux de réalisation, le décodage itératif à passage de messages de l'étape (c) et le décodage itératif à passage de messages de l'étape (f) utilisent tous deux un algorithme de type « Min-Sum ». BRÈVE DESCRIPTION DES DESSINS D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de modes de réalisation préférentiels de l'invention faite en référence aux figures jointes parmi lesquelles : La Fig. 1 représente de manière schématique un générateur de séquence maximale selon une représentation de Fibonacci ; La Fig. 2 donne des paramètres de codes C/A pour différents satellites ; La Fig. 3 représente de manière schématique un générateur de code C/A dans un système GPS ; La Fig. 4 représente schématiquement une méthode d'acquisition d'un signal satellitaire selon un premier mode de réalisation de l'invention ; d'un signal La Fig. 5 représente schématiquement une méthode d'acquisition satellitaire selon une variante du premier mode de réalisation de l'invention ; d'un signal La Fig. 6 représente schématiquement une méthode d'acquisition satellitaire selon un second mode de réalisation de l'invention ; d'un signal La Fig. 7 représente schématiquement une méthode d'acquisition satellitaire selon une variante du second mode de réalisation de l'invention. EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS Nous rappellerons dans un premier temps comment le signal de synchronisation grossière C/A est émis par un satellite d'un système de positionnement tel que le système GPS ou le système Galileo. Le signal émis par un satellite sur le canal C/A peut s'exprimer sous la forme : k L LN _ x1(k)=d q(k) (1) où i est l'indice du satellite, d(e) est le Cr- symbole BPSK du message de navigation, k est un indice temporel, q(k) est la kième bribe (chip) du code d'étalement associé au satellite i, où N est la longueur du code d'étalement (en pratique N =1023 pour le système GPS) et L est le nombre de symboles BPSK du message de navigation. Pour le système GPS, chaque symbole ci(e) dure 20 ms, soit 20 périodes de répétition du code d'étalement, la fréquence bribe (ou rythme chip) étant de 1,023 MHz. Le signal satellitaire reçu par le récepteur sur le canal C/A peut alors s'écrire, en bande de base : p1x,(k)e27rj(kAF,T+6;) w(k) (2) 1,1 où / est le nombre maximum de satellites (I =32 dans le système GPS), p, est le coefficient d'atténuation du signal sur la liaison entre le satellite i et le récepteur, supposé constant sur la durée du message, AF, est l'offset de la fréquence porteuse (principalement dû au décalage Doppler) et 0, est la phase du signal en provenance du satellite i (supposée constante), et où enfin w(k) est un échantillon de bruit supposé blanc gaussien (AWGN) de variance o-2. Les séquences d'étalement c des différents satellites i =1,...,/ sont des séquences de Gold, chaque séquence de Gold étant construite à partir d'un somme de deux séquences de longueur maximale, encore dénommées M-séquences. On rappelle qu'une M-séquence est une suite périodique de valeurs produite par un registre à décalage à rétroaction linéaire (LFSR) qui explore toutes les valeurs pouvant être produites par le registre à décalage. Toute combinaison linéaire de deux M-séquences de même degré r (c'est-à-dire générées par des polynômes générateurs de même degré) ne donne pas nécessairement une séquence de Gold : une paire de M-séquences dont la somme donne une séquence de Gold est appelée paire de séquences préférentielles. Les séquences d'étalement c, du canal C/A sont chacune construites à partir d'une paire de M-séquences préférentielles, x et y. Plus précisément : ci(k)= x(k)+ y(k- ri) (3) où r, est un retard, exprimé en nombre de périodes bribes, spécifique au satellite j. Les M-séquences x et y sont respectivement générées à partir des polynômes générateurs : gx (D) -= Di° + +1 gy(D)=D10+E9-F-D8+D3+D2+1 (4) On rappelle qu'une séquence LFSR générée par un polynôme générateur r-1 g(D)=1-FIgpDP de degré r est obtenue par un registre à décalage à r positions p=.1 mémoire (dénommées ci-après plus simplement positions) et r prises, rebouclé sur lui-même, comme représenté en Fig. 1, selon la représentation dite de Fibonacci. Les coefficients multiplicateurs sur les différentes prises sont les coefficients du polynôme générateur. Une M-séquence est définie de manière univoque par son polynôme générateur, le contenu du registre à décalage à l'instant initial donnant sa valeur à l'origine. Une M-séquence générée par un polynôme générateur (primitif) de degré r est périodique et de longueur N = 2' -1 . Les M-séquences x et Y sont générées en initialisant les positions de leurs registres à décalage respectifs à 1, soit x(0) = x(1)=...= x(9)=1 et y(0) = y(1)=...= y(9)=1. La séquence retardée y définie par yi(k). y(k- r ,) est avantageusement obtenue en utilisant la propriété « add and shift » des M-séquences. Plus précisément, les délais ri sont choisis de sorte à ce que la séquence y puisse être générée à partir de la somme des contenus en deux positions a, et du registre à décalage générant la séquence y soit : y(k)=a(k)+a(k) (5) où aa,(k), a), (k) sont respectivement les contenus du registre à décalage aux positions a, et pi, à l'instant k, autrement dit aai (k) = y(k+ a,) et afl, (k) = y(k+ )6,). On a représenté en Fig. 2 une table donnant les paramètres a,, J6, et r, pour les différents satellites. Plus précisément, la colonne intitulée « code phase selection » donne les positions a,,,61 et la colonne intitulée « code delay chips » donne le retard r, . Enfin, la dernière colonne fournit le mot d'initialisation, c'est-à-dire les 10 premières valeurs de la séquence y, en représentation octale. La Fig. 3 représente de manière schématique un générateur de code C/A dans un système GPS. On distingue dans la partie supérieure le générateur, 310, de la M-séquence X et dans la partie inférieure le générateur, 320, de la M-séquence y. Les deux générateurs 310 et 320 utilisent des registres à décalage rebouclés sur eux-mêmes. Pour chaque générateur, les prises du registre à décalage correspondent aux monômes du polynôme générateur correspondant. Le sélecteur de phase, 330, vient sommer sélectivement certaines positions du registre à décalage de 320. Le choix de ces positions a, et est effectué selon la table de la Fig. 2, la somme des contenus du registre en ces deux positions fournissant la séquence retardée y. Cette implémentation est avantageuse dans la mesure où elle ne nécessite pas l'utilisation d'une importante ligne à retard, de longueur r,. La somme des M-séquences x et y en 350 donne la séquence de Gold c' autrement dit la séquence d'étalement relative au satellite I. La Fig. 4 représente schématiquement une méthode d'acquisition d'un signal satellitaire selon un premier mode de réalisation de l'invention. On considérera dans un premier temps la synchronisation du récepteur par rapport à un signal GPS d'un satellite I. Le récepteur échantillonne à la fréquence bribe le signal satellitaire reçu, après démodulation en bande de base. A chaque instant q, le récepteur forme en 410 un vecteur constitué de M échantillons successifs du signal ainsi reçu, soit rq(k)-= r(k+ q), k -1, où M = N+1 est la période de la M-séquence ( M = 1024 dans le système GPS). Le récepteur teste alors s'il est bien synchronisé avec la M-séquence X. Pour ce faire, il fait l'hypothèse qu'il est synchronisé avec cette séquence et effectue une multiplication bribe à bribe du signal reçu avec la M-séquence x, à l'étape 420, soit : q(k) = rq (k) X(k) (6) où X(k)= 1- 2 x(k) est la représentation bipolaire (BPSK) de la M-séquence X. Cette opération permet d'éliminer la M-séquence x du signal reçu : q(k) = elii(k)+ tiv(k) (7) où y(k) 1 - 2y (k) est la représentation bipolaire de la M-séquence y, et où vit(k) regroupe toutes les sources de bruit. On déduit de l'expression (7) que le vecteur r~ q(k), est une observation bruitée de M bribes consécutives de la séquence y. On effectue à l'étape 430 un décodage itératif à passage de messages (encore dénommé à propagation de croyance ou belief propagation) du mot rq(k) k = M -1 En effet, l'opération qui associe le contenu (à l'instant q) du registre à décalage d'un générateur d'une M-séquence à la M-séquence elle-même peut être considérée comme un codage linéaire cyclique de rendement Or -1 (le mot du registre est de taille r alors que la M-séquence est de taille N =2' -1). De manière générale, si l'on considère une M-séquence z définie par un polynôme générateur gz(D)= gi,DP , on a la contrainte de parité : P=0 gr_pz(pt+ k)- 0 (8) p=0 dans laquelle on a défini gr = I. La matrice de parité du code est la matrice de taille (M-r)x M: gr - go 0 - - - 0 gr - go 0 H= (9) 0 0 gr - - - g, 0 0 - - - 0 g. - - - g, Le décodage s'effectue alors comme pour un décodage par passage de messages entre noeuds de variable (les variables y,(k)) et noeuds de contrainte (les contraintes de parité données par les lignes de la matrice de parité H où les coefficients sont ceux du polynôme gy(D)). Une description de la méthode de décodage par passage de messages a été donnée dans l'article de F.R. Kschichang et al. intitulé « Factor graphs and the SumProduct algorithm » publié dans IEEE Trans. on Information Theory, Vol. 47, No. 2, Feb. 2001. Une version plus simple mais sous-optimale de cet algorithme a été proposée par V. Savin dans l'article intitulé « Self-corrected Min-Sum decoding of LDPC codes » publié dans Proc. of Int'l Symposium on Information Theory (ISIT 2008), Toronto, Juillet 2008. Quelle que soit la version utilisée, le décodage du mot, r q(k), k=0,...,M -1 (ou d'une partie de ce mot seulement puisque celui-ci présente un redondance très élevée) fournit un vecteur 1-",q(k), k= M - r - 1 (l'indice q a été ajouté pour souligner que ce vecteur a été obtenu en prenant pour hypothèse que la séquence était synchronisée à l'instant q) duquel on ne retient que les r premières valeurs donnant, au moyen de décisions dures, le contenu du registre à l'instant initial q. Les valeurs ?,q(k) sont des logarithmes de rapports de vraisemblance (LLR) et donc le contenu du registre à décalage à l'instant q est donné par :25 A,q(k), sgn(S'fq,), k= 0,...,r -1 (1o) où sgn0 est la fonction signe associant à toute valeur positive la valeur +1 et à toute valeur négative la valeur -1. A l'étape 440, on compare le vecteur A,9(k), k= r -1 avec les / = 32 mots possibles d'initialisation du registre à décalage de la séquence y, notés Ces mots d'initialisation se déduisent du mot d'initialisation Alyn't =( 1 .- - 1 )T du registre à décalage de la séquence y. En effet, on rappelle que, d'après l'expression (5) : y(k)- y(k+ G)+ y(k+ A) (11) ce qui se traduit matriciellement par : (Gya' +Gye')Ay (12) où Ai et Ay sont respectivement le mot contenu dans le registre à décalage de la séquence y et le mot contenu dans le registre à décalage de la séquence y, Gy est la matrice de transposition d'état, c'est-à-dire la matrice reliant les vecteurs de contenu du registre à décalage en des instants consécutifs, soit : Ay(k + 1) = GyAy(k) (13) La matrice de transposition d'états Gy est donnée par : 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 = (14) 0 0 0 0 0 0 1 1 g1 1 Y g7 -3 où les coefficients e..., g, sont ceux du polynôme générateur de la séquence y. Si l'on applique la relation (12) à l'instant d'initialisation, on a, en particulier : (G yai + G ) Ai nit (15) Les vecteurs A,n't peuvent être calculés à l'avance de manière à ne pas ralentir l'acquisition. S'il existe un indice io tel que : A=Ainit (16) 15 on considère que le récepteur est synchronisé avec le signal GPS du satellite io et la méthode d'acquisition se termine en 490. A défaut, la méthode d'acquisition se poursuit en incrémentant q en 445, c'est-à-dire en partant de l'échantillon suivant, et en retournant à l'étape 410. 20 La Fig. 5 représente schématiquement une variante de la méthode d'acquisition de la Fig. 4. Les étapes 510 à 545 sont identiques aux étapes 420 à 445 et leur description ne sera donc omise. 10 La variante de la Fig. 5 se distingue de celle de la Fig. 4 au sens où, une fois que le récepteur a déterminé l'indice io du satellite avec lequel il considère être synchronisé, il procède à une phase de vérification. Cette phase de vérification permet de réduire considérablement la probabilité de fausse alarme. Elle est particulièrement avantageuse lorsque la méthode de décodage par passage de messages est implémentée dans sa version simplifiée « Min-Sum » mentionnée plus haut. Cette phase de vérification est basée sur le principe suivant : si la séquence yo qui a été détectée est correcte, on peut alors l'utiliser pour retrouver la M-séquence X à partir du signal reçu rq(k), k= M -1. Plus précisément, à l'étape 550, on multiplie le signal reçu rq(k), r(k q), k M 1 avec la séquence Y,o(k) = 1 - 2y,0 (k), k 0 M - 1, représentation bipolaire de la M-séquence y, . On obtient la séquence définie par : fq(k) = rq(k)yo(k) (17) Cette opération élimine la séquence yo du signal reçu, si celle-ci y est bien présente : fq(k)= X(k)+ 1-47(k) (18) où X(k) = 1 - 2 x(k), k = M -1, est la représentation bipolaire de la M-séquence X et où W(k) regroupe toutes les sources de bruit. On comprendra que l'étape 550 est en quelque sorte la « symétrique » de l'étape 520. Le vecteur (k) k= M -1, étant une observation bruitée de M bribes q consécutives de la M-séquence X, on en effectue à l'étape 560 un décodage itératif à l'aide d'un algorithme à passage de messages, de type « Sum-Product» ou « Min-Sum ». 25 Par exemple, si un algorithme Min-Sum a été utilisé pour le décodage de rq(k), k= M-1 à l'étape 530, le type algorithme sera avantageusement utilisé pour le décodage de rq(k), k= 0, M -1. Les contraintes sont données par les lignes de la matrice H dans laquelle les coefficients sont ici celui du polynôme générateur gx(D). Quel que soit le type de l'algorithme, le décodage du mot, tq(k), k= M -1 (ou d'une partie de ce mot seulement puisque celui-ci présente un redondance très élevée) fournit un vecteur )-(q(k),k=0,...,M -r -1 (l'indice q a été ajouté pour souligner que ce vecteur a été obtenu en prenant pour hypothèse que la séquence était synchronisée à l'instant q) duquel on ne retient que les r premières valeurs donnant, au moyen de décisions dures, le contenu du registre à l'instant initial q. Les valeurs 5'q(k) sont des logarithmes de rapports de vraisemblance (LLR) et donc le contenu du registre à décalage à l'instant q est donné par : Aq(k)= sgn(Xq(k)), k=0,...,r -1 (19) A l'étape 570, on vérifie si le vecteur Aq est égal au vecteur d'initialisation de la M-séquence X, autrement dit si Aq(k)= 1, k=0,...,r -1. Si c'est bien le cas, c'est-à-dire si la vérification est positive, on conclut que le récepteur est bien synchronisé sur le signal GPS du satellite io. La méthode d'acquisition s'achève alors en 590. En revanche, si ce n'est pas le cas, la synchronisation avec le signal GPS du satellite io est infirmée. La méthode d'acquisition se poursuit en incrémentant q en 575, c'est-à-dire en partant de l'échantillon suivant, et en retournant à l'étape 510. La probabilité de fausse alarme dans le cadre de cette variante est égale à Pf2, où Pja est la probabilité de fausse alarme obtenue avec la méthode d'acquisition (sans vérification) de la Fig. 4. Ainsi, par exemple, si la probabilité de fausse alarme sans vérification est de iO3, celle avec vérification sera de l'ordre de 10-6. La Fig. 6 représente schématiquement une méthode d'acquisition d'un signal GPS selon un second mode de réalisation de l'invention. Ce mode de réalisation tire parti du fait que deux M-séquences sur le même corps de Galois GF( 2r) sont liées par une relation de décimation. Ainsi, si u et y sont deux M-séquences générées par des polynômes primitifs de degré r, on a: u(k) v(dk + h) (20) où d est un facteur de décimation et h est une constante dépendant du contenu initial des registres des M-séquences u et y. Pour une paire de séquences préférentielles x,y, le facteur de décimation vaut d = 2e +1 avec la condition ged(2e+1,2r -1)=1. Dans le cas des polynômes générateurs gx et g, donnés en (4), r =10, et l'on peut alors montrer que e = 6, d'où d = 65 et, compte tenu des mots d'initialisation des registres pour les séquences x, y : y(k) x (65k + 27300) (21) Si l'on définit la séquence x par x(k) = x (k -65r,) , on a : y,(k) = x. (65k +27300) (22) La relation entre la séquence y, et la séquence x, est indépendante de l'indice 1, autrement dit celle-ci reste valable quel que soit le satellite. On considérera là encore dans un premier temps la synchronisation du récepteur par rapport à un signal satellitaire. Le récepteur échantillonne le signal reçu à la fréquence bribe, après démodulation en bande de base et forme en 610, à chaque instant q , le vecteur rq(k), r(k + q), k= 0,..,M -1, constitué de M échantillons successifs du signal reçu. A l'étape 620, on effectue une multiplication bribe à bribe du signal reçu avec la M-séquence x pour obtenir les échantillons retraités Fq(k), k= 0,..,M -1 comme définis en relation (6). A l'étape 630, on effectue une permutation des M échantillons retraités au moyen de la permutation : q)(k) = dk + h mod(N) (23) c'est-à-dire dans le cas des polynômes générateurs g Quel que soit le type de l'algorithme, le décodage fournit un vecteur k=0,...,M -1, duquel on ne retient que les r premières composantes donnant, au moyen de décisions dures, le contenu du registre à l'instant initial q. soit Ag le vecteur constitué par les r premières composantes en question. A l'étape 785, on vérifie si le vecteur Ag est égal au vecteur d'initialisation de la M-séquence x, autrement dit si AI (k) ,k = 0 , ...,r -1. Si c'est bien le cas, c'est-à-dire si la vérification est positive, on conclut que le récepteur est bien synchronisé sur le signal GPS du satellite io. La méthode d'acquisition s'achève alors en 790. En revanche, si ce n'est pas le cas, la synchronisation avec le signal GPS du satellite io est infirmée. La méthode d'acquisition se poursuit en incrémentant q en 787, c'est-à-dire en partant de l'échantillon suivant, et en retournant à l'étape 710. Dans les modes de réalisation précédents, le récepteur détermine s'il est synchronisé à un instant donné par rapport au signal GPS du canal C/A pour un satellite io donné. En pratique, le canal C/A comprend des signaux d'une pluralité de satellites d'indice io,...,i où Q est le nombre de satellite visibles du récepteur. La méthode d'acquisition exposée plus haut fournit le satellite io présentant le rapport signal sur bruit le plus élevé au niveau de récepteur. Une fois ce premier signal acquis, on peut reconstruire le signal reçu du satellite io, en bande de base : r,o(k) = (,5 ,oebr ).,c,o(k)e2ff jk AF,'T (28) que l'on peut ensuite soustraire du signal reçu en bande de base (expression (2)). " DrAo est une estimation du coefficient d'atténuation complexe sur le canal entre le satellite et le récepteur (obtenu par exemple au moyen de symboles pilotes) et AFia est l'offset de fréquence estimé par le récepteur. x,o(k) est un message de synchronisation connu du récepteur. On poursuit ensuite avec l'acquisition du second satellite à partir du signal r(k)-17,(k) et ainsi de suite, selon un schéma SIC (Serial Interference Cancellation) connu en soi. REVENDICATIONS1. Méthode d'acquisition d'un signal satellitaire par un récepteur, ledit signal étant étalé par une séquence de Gold (c,) à une fréquence bribe, ladite séquence de Gold étant obtenue comme somme d'une première M-séquence (x ) générée au moyen d'un premier registre à décalage et d'une seconde M-séquence (y1) générée au moyen d'un second registre à décalage, ladite seconde M-séquence identifiant un satellite parmi une pluralité de satellites, caractérisée en ce que : (a) le signal reçu en bande de base est échantillonné (410,510,610,710) à la fréquence bribe à partir d'un instant déterminé (g) pour fournir une séquence de premiers échantillons (r,7(k),k = 0, M - 1 ); (b) la séquence desdits premiers échantillons est multipliée (420,520,620,720), bribe à bribe, par ladite première M-séquence pour fournir une séquence de seconds échantillons (( (k), k = 0 , , M - 1); (c) la séquence de seconds échantillons est décodée au moyen d'un décodage itératif à passage de messages (430,530,640,740) pour fournir une estimation du contenu du second registre à décalage au dit instant déterminé ; (d) le contenu du second registre à décalage ainsi estimé est comparé (440,540,660,760) à une pluralité de contenus possibles du second registre à décalage pour les différents satellites de ladite pluralité, et en cas de succès de la comparaison pour un satellite identifié (i0), le récepteur est synchronisé avec ledit satellite identifié à l'instant déterminé. 2. Méthode d'acquisition d'un signal satellitaire selon la revendication 1, caractérisée en ce que si le contenu du second registre à décalage ainsi estimé ne correspond à aucun desdits contenus possibles, l'instant déterminé est incrémenté d'une bribe et les étapes (a) à (d) sont itérées. 3. Méthode d'acquisition d'un signal satellitaire selon la revendication 1, caractérisée en ce que l'estimation du contenu du second registre à décalage est réalisée en prenant des décisions dures sur les r premières valeurs (17:1(k)) décodées fournies par le décodage itératif à passage de messages, où r est le nombre de positions dans le premier/second registre à décalage. 4. Méthode d'acquisition d'un signal satellitaire selon la revendication 3, caractérisée en ce que le décodage itératif à passage de messages de l'étape (c) utilise la matrice de parité : go 0 --- 0 gr 0 g. - - - go 0 0 H 0 0 g --- go 0 0 --- O gr -" go dans laquelle g,. =1 et sont les coefficients du polynôme générateur de la seconde M-séquence. 5. Méthode d'acquisition d'un signal satellitaire selon la revendication 1, caractérisée en ce que, la seconde M-séquence étant obtenue en sommant les contenus en deux positions a , g du second registre à décalage, où i est un indice identifiant le satellite parmi ladite pluralité de satellites, les contenus possibles du second registre à décalage pour les différents satellites sont donnés par les vecteurs : A;nif (G + G'35) )Alyn't où A.7 est un vecteur donnant contenu du second registre à décalage à un instant de référence, et Gy la matrice de changement d'état du second registre à décalage. 6. Méthode d'acquisition d'un signal satellitaire selon la revendication 1, caractérisée en ce que : la séquence de seconds échantillons est permutée (630,730) à l'aide d'une relation de permutation (p(k) = dk + h mod(N) où k est le rang de l'échantillon, d est un facteur de décimation, h est une valeur d'offset et N est la longueur de la première M-séquence, la permutation étant effectuée avant l'étape (c) ; l'étape (c) fournit une séquence de valeurs décodées ('(k)) ; la séquence de valeurs décodées est permutée (650,750) à l'aide de la relation de permutation inverse à ladite relation de permutation pour fournir une séquence de valeurs permutées ( È:9 (k) ); l'estimation du contenu du second registre à décalage (660,760) est réalisée en prenant des décisions dures sur les r premières valeurs de la séquence de valeurs permutées, où r est le nombre de positions dans le premier/second registre à décalage. 7. Méthode d'acquisition d'un signal satellitaire selon la revendication 6, caractérisée en ce que le décodage itératif à passage de messages de l'étape (c) utilise la matrice de parité : go 0 - - - O 0 gr - - - go 0 o H= 0 - - - 0 gr - - - go 0 0 - - - - - 0 g, - - - go dans laquelle gr =1 etgo'---'gri sont les coefficients du polynôme générateur de la première M-séquence. 8. Méthode d'acquisition d'un signal satellitaire selon l'une des revendications précédentes, caractérisée en ce que l'on vérifie que le récepteur est bien synchronisé avec ledit satellite identifié en : (e) multipliant la séquence desdits premiers échantillons (550,780), bribe à bribe, par la seconde M-séquence associée audit satellite identifié pour fournir une séquence de troisièmes échantillons ((i-, (k), k = 0,.., M -1); (f) décodant la séquence de troisièmes échantillons au moyen d'un décodage itératif à passage de messages (560, 783) pour fournir une estimation du contenu du premier registre à décalage au dit instant déterminé. 9. Méthode d'acquisition d'un signal satellitaire selon la revendication 8, caractérisée en ce que caractérisée en ce que si le contenu du premier registre à décalage ainsi estimé ne correspond au contenu initial du premier registre, l'instant déterminé est incrémenté d'une bribe et les étapes (a) à (f) sont itérées. 10. Méthode d'acquisition d'un signal satellitaire selon la revendication 9, caractérisée en ce que l'estimation du contenu du premier registre à décalage est réalisée en prenant des décisions dures sur les r premières valeurs ( ,(k),k 0, r -1) décodées fournies par le décodage itératif à passage de messages, où r est le nombre de positions dans le premier/second registre à décalage. 11. Méthode d'acquisition d'un signal satellitaire selon l'une des revendications 8 à 10, caractérisée en ce que le décodage itératif à passage de messages de l'étape (f) utilise la matrice de parité gr go O 0 0 g. - - go 0 0 H= O 0 g- - - go 0 O 0 g, godans laquelle gr =1 et go,...,g sont les coefficients du polynôme générateur de la première M-séquence. 12. Méthode d'acquisition d'un signal satellitaire selon la revendication 8, caractérisée en ce que le décodage itératif à passage de messages de l'étape (c) et le décodage itératif à passage de messages de l'étape (f) utilisent tous deux un algorithme de type « Min-Sum ».