FR 2936669 A1 20100402 FR 0805320 A 20080926 Procédé pour optimiser une acquisition d'un signal à spectre étalé provenant d'un satellite par un récepteur mobile La présente invention concerne un procédé pour optimiser un processus d'acquisition d'un signal à spectre étalé, mis en oeuvre par un récepteur mobile. L'invention s'applique par exemple à tout signal de navigation provenant d'un système global de navigation par satellites de type GPS (Global Positioning System), Galileo. On entend ici par � récepteur mobile � aussi bien les dispositifs (ou récepteurs) dédiés exclusivement au positionnement par satellites portables ou embarqués dans un véhicule terrestre, maritime ou aérien que les terminaux de communication équipés d'un dispositif de positionnement par satellites, comme par exemple des téléphones mobiles, ou d'ordinateurs portables ou des assistants personnels numériques (ou PDA), éventuellement de type communicant. Dans un système de positionnement par satellites utilisant des récepteurs du type GNSS (Global Navigation Satellite System), les signaux de données permettant au récepteur de calculer son positionnement proviennent de différents satellites appartenant à une constellation de satellites de positionnement (quatre au minimum pour déterminer quatre inconnues correspondant aux coordonnées géographiques x, y, z et temporelles t du récepteur). Le positionnement par satellites consiste en un enchaînement de deux étapes. La première étape, dite d'acquisition, consiste à déterminer, au niveau du récepteur mobile concerné, les codes d'étalement pseudo-aléatoires qui modulent les signaux provenant des satellites appartenant à la constellation et rapportés à un temps de référence. Il s'agit en fait de � comparer � les signaux reçus des satellites à des répliques de signaux élaborés localement par le récepteur et résultant d'hypothèses sur le temps de référence et sur la fréquence de cadencement des satellites, afin d'en déduire les codes pseudo- aléatoires qui modulent lesdits signaux reçus ou en d'autres termes de synchroniser l'horloge de cadencement du récepteur et sa fréquence sur l'horloge et la fréquence de chaque satellite. Pour cela, une recherche en temps-fréquence de l'énergie du signal provenant du satellite est effectuée, cette recherche étant habituellement réalisée par des mesures de corrélation reposant sur des couples d'hypothèses temporelle et fréquentielle afin de déterminer la corrélation maximale entre le signal reçu et la réplique locale du récepteur. La seconde étape consiste à déterminer la position du récepteur mobile à partir des codes acquis et de données de navigation notamment contenues dans les signaux reçus. Cette seconde étape peut être plus précisément subdivisée en trois sous-étapes : une sous-étape de détermination à partir des codes pseudo-aléatoires acquis, des temps de propagation des signaux entre chacun des satellites et le récepteur, une sous-étape de détermination, à partir des données de navigation contenues dans les signaux et des temps de propagation, de pseudo-distances entre le récepteur et chacun des satellites, et une sous-étape de détermination de la position du récepteur à partir des pseudo-distances. La précision de chaque temps de propagation, et donc de chaque pseudo-distance, détermine directement la précision de la position. Or, la précision de chaque temps de propagation dépend de la qualité de l'acquisition des codes pseudo-aléatoires du signal reçu correspondant, laquelle est dépendante de la qualité dudit signal reçu. Par conséquent, lorsque l'un au moins des signaux reçus d'un satellite est de mauvaise qualité, ce qui est relativement fréquent, notamment dans les environnements accidentés ou encombrés tels que les milieux urbains, la position déterminée est généralement entachée d'erreur. On peut même se trouver momentanément dans l'impossibilité de déterminer la position du récepteur, alors même que les signaux provenant des autres satellites sont de bonne qualité. Le récepteur possède trois incertitudes fréquentielles auxquelles s'ajoute une inconnue sur la date le conduisant à effectuer la recherche en temps et en fréquence de l'énergie provenant d'un satellite. Ces trois incertitudes fréquentielles sont l'effet Doppler lié à la mobilité du satellite, l'incertitude lié à la précision de l'horloge du récepteur, et l'effet Doppler lié à la mobilité du récepteur. L'effet Doppler lié au mouvement du satellite peut être déterminé de manière connue en utilisant par exemple un serveur d'assistance tel que celui utilisé dans les techniques de localisation de type AGPS (Assisted GPS). Les oscillateurs locaux des horloges sont de plus en plus performants et deviennent de plus en plus stables. L'effet Doppler lié au mouvement du récepteur devient la source prépondérante des incertitudes sur la localisation du pic de corrélation et donc du récepteur. A notre connaissance, aucune des méthodes connues ne prend en compte l'effet Doppler lié au récepteur. Ceci peut avoir un effet dramatique dans le cas de l'acquisition de signaux de faible énergie qui nécessitent d'intégrer le signal reçu de manière cohérente sur une longue durée temporelle. En effet, la largeur des fenêtres d'intégration des hypothèses fréquentielles étant inversement proportionnelle à la durée d'intégration, plus le temps d'intégration de manière cohérente est long, plus la largeur des fenêtres d'intégration des hypothèses fréquentielles est petite et donc plus le nombre d'hypothèses fréquentielles est grand. Les balayages en temps et en fréquences impliquent alors une puissance de calcul et un temps de traitement des données très importants pour le récepteur et augmentent l'incertitude sur la localisation du récepteur. Par ailleurs, dans un schéma d'acquisition classique, l'intégration cohérente s'effectue sur une hypothèse de fréquence qui reste la même tout au long de la durée d'intégration, ce qui suppose que la fréquence effective du signal reçu est stable sur la durée d'intégration. Lors d'une intégration cohérente longue, le doppler lié au mouvement de l'utilisateur a pour effet de faire varier la fréquence effective reçue durant ladite intégration et ce qui rend inopérant l'intégration si celle-ci s'effectue sur une hypothèse de stabilité. La présente invention a pour but de remédier à ce problème en proposant un procédé pour optimiser une acquisition d'un signal de navigation provenant d'un satellite par un récepteur mobile prenant en compte l'effet Doppler lié au récepteur et permettant de réduire le nombre d'hypothèses fréquentielles à explorer lors de la recherche du pic de corrélation du signal et de diminuer le temps de calcul au sein du récepteur. A cet effet, l'invention a pour objet un procédé pour optimiser une acquisition d'un signal à spectre étalé provenant d'un satellite par un récepteur mobile, consistant : à recevoir un signal émis par au moins un émetteur, le signal étant modulé par un code pseudo-aléatoire, à élaborer une réplique locale du code pseudo-aléatoire du signal 10 émis, à émettre une pluralité d'hypothèses relatives à la fréquence du signal reçu à un instant initial to et à l'évolution de cette fréquence au cours du temps jusqu'à un instant final tri correspondant à un temps d'intégration, 15 à calculer les fonctions de corrélation de tous les couples d'hypothèses fréquence - temps, caractérisé en ce qu'il consiste en outre : à choisir une première hypothèse de fréquence initiale f, correspondant à l'instant initial to, 20 à partir de la même hypothèse de fréquence initiale à choisir successivement, les uns après les autres, tous les chemins possibles d'évolution de la fréquence fi entre les instants to et tn et à déterminer leur énergie totale respective, à choisir successivement toutes les autres hypothèses de 25 fréquences initiales possibles et pour chaque hypothèse de fréquence initiale à sélectionner le chemin d'évolution dont l'énergie totale est la plus forte, à décider que, parmi tous les chemins sélectionnés à chaque hypothèse de fréquence, le chemin qui a l'énergie la plus forte est 30 le chemin le plus vraisemblable et correspond au chemin réel suivi par le récepteur. Avantageusement, l'énergie totale de chaque chemin possible d'évolution est obtenue en sommant de manière cohérente les résultats de corrélation entre les instants initial to et final t,,, la fréquence prenant successivement toutes les valeurs de fréquences situées sur le chemin choisi. Préférentiellement, pour chaque hypothèse de fréquence initiale, les chemins possibles d'évolution de la fréquence sont déterminés en prenant en compte des informations délivrées par au moins un capteur inertiel localisé dans le récepteur . Avantageusement, les informations délivrées par le capteur inertiel sont prises en compte dans une enveloppe d'erreur établie à partir des caractéristiques du capteur. Optionnellement, lorsque la variation des informations délivrées par le capteur, entre deux instants de mesure, est inférieure à un seuil de variation prédéterminé, il est décidé que le récepteur n'a pas bougé entre les deux instants de mesure, et il est alors fait l'hypothèse d'une fréquence stable lors de l'intégration. D'autres particularités et avantages de l'invention apparaîtront clairement dans la suite de la description donnée à titre d'exemple purement illustratif et non limitatif, en référence aux dessins schématiques annexés qui représentent : figure 1 : un exemple d'architecture d'un récepteur d'un système global de navigation par satellite, selon l'invention ; figure 2 : un exemple de graphique fréquence-temps illustrant les hypothèses de fréquences à explorer lors d'une acquisition sur 30 un temps d'intégration donné; - figure 3 : un exemple des chemins possibles à explorer pour une hypothèse de fréquence initiale fi donnée, selon l'invention; figure 4: un exemple des chemins possibles pour une hypothèse de fréquence initiale fi donnée, en tenant compte de l'information délivrée par un capteur inertiel, selon l'invention. L'architecture représentée sur la figure 1 comporte des moyens de réception 10 des signaux radiofréquences émis par un satellite, des moyens d'acquisition 11 des signaux reçus pour déterminer les codes d'étalement pseudo-aléatoires qui modulent les signaux provenant du satellite, des moyens de poursuite 12 des signaux permettant de compenser la dynamique du récepteur, des moyens de calcul 13 de la position du récepteur. En outre, l'architecture comporte des capteurs inertiels 14 pour délivrer, après traitement par des moyens de traitement 15, des informations relatives à la vitesse et à l'attitude du récepteur par rapport au satellite. Le graphique de la figure 2 comporte un axe vertical correspondant aux hypothèses de fréquence à explorer et un axe horizontal correspondant au temps d'intégration. Les hypothèses de fréquences à explorer sont comprises entre une fréquence minimale liée à un effet Doppler minimal et une fréquence maximale liée à un effet Doppler maximal. L'axe des temps et l'axe des fréquences sont découpés respectivement en durées temporelles élémentaires et en bandes de fréquences élémentaires et chaque couple fréquence-temps correspond à une case particulière du quadrillage de la figure 2. A chaque hypothèse de fréquence correspond une disposition relative récepteur-satellites et un décalage horloge récepteur possibles à un instant donné. Entre deux instants donnés to et tn, la position du récepteur peut évoluer suivant une pluralité de chemins possibles. Pour connaître la position du récepteur à chaque instant, l'invention consiste donc à explorer chaque chemin possible et à calculer l'énergie du signal pour chacun de ces chemins. Le chemin qui possède l'énergie la plus forte est le chemin suivi par le récepteur. A chaque fréquence f; donnée et à un instant tk donné correspondant à une case élémentaire du graphique de la figure 2, l'énergie du signal est calculée. Pour cela, le récepteur corrèle le signal reçu avec une réplique du code pseudo-aléatoire du satellite qui a envoyé le signal. Le code pseudo-aléatoire du satellite est connu et ne change pas pendant toute la vie du satellite. Le résultat de la corrélation à la fréquence f; et à l'instant tk s'écrit : k+1)T s(f , tk , z-) =--T r(t) • c (t û kT û z-)--2- R~ • 2d t fC Où r(t) est le signal reçu par le récepteur, c(t) est le code d'étalement utilisé dans la modulation du signal reçu, T est le temps d'intégration 10 cohérent qui généralement correspond à une période du code d'étalement c(t), T est le décalage estimé entre le code d'étalement reçu et le code d'étalement émis, fc est la fréquence de la porteuse (typiquement 1575,42MHz pour GPS L1 C/A), Rc est le rythme de bribe du code d'étalement (typiquement 1.023Mbribes/s pour GPS L1 C/A). Par exemple, 15 dans le cas d'un signal GPS C/A, T est typiquement égal à 1 ms. Mais T pourrait avoir une toute autre valeur. Lorsque toutes les corrélations élémentaires correspondant à chaque hypothèse de fréquence ont été réalisées, le procédé selon l'invention consiste à choisir une première hypothèse de fréquence initiale fi 20 correspondant à un instant initial to, puis à choisir un chemin particulier d'évolution de la fréquence du signal entre to et tn parmi tous les chemins possibles et à déterminer l'énergie totale du signal correspondant à ce chemin par sommation de chaque élément s(fi, tk ,'r,) constituant le chemin. tn correspond à l'instant final de la recherche. La figure 3 illustre l'ensemble des 25 chemins possibles à partir d'une hypothèse de fréquence initiale fi à un instant to jusqu'à l'instant t,,. L'énergie totale correspondant à un chemin choisi est obtenue en sommant de manière cohérente les résultats de corrélation s(fi, tk, r), tk prenant successivement toutes les valeurs entre to et tn et fi prenant successivement toutes les valeurs de fréquences situées sur5 le chemin sélectionné. Par exemple, à l'instant th, la valeur de la fréquence peut avoir évoluée à la valeur f2 ou f3 ou être restée inchangée. A partir de la même hypothèse de fréquence, tous les chemins possibles sont ensuite sélectionnés successivement, les uns après les autres, et leur énergie respective est calculée de façon similaire. Parmi tous les chemins possibles, il n'existe qu'un unique chemin permettant une recombinaison des résultats de corrélation non destructive, illustré en trait épaissi 3 sur la figure 3. Ce chemin dépend entre autre de la dynamique relative entre le satellite et le récepteur. A ce chemin correspond donc une énergie du signal maximale. Le procédé selon l'invention consiste alors pour chaque hypothèse de fréquence initiale, à sélectionner le chemin dont l'énergie totale est la plus forte et à décider que, parmi tous les chemins sélectionnés à chaque hypothèse de fréquence, le chemin qui a l'énergie la plus forte est le chemin le plus vraisemblable et correspond au chemin réel suivi par le récepteur. Ce procédé nécessitant d'explorer tous les chemins possibles pour chaque hypothèse de fréquence initiale est très lourd et très long en temps de traitement. Avantageusement, pour réduire le nombre de chemins à explorer pour chaque hypothèse de fréquence initiale, entre les instants initial et final, le procédé comporte une étape intermédiaire consistant à utiliser au moins un capteur inertiel, par exemple de type MEMS (Micro-Electro- Mechanical Systems). Le capteur inertiel est placé par exemple sur le récepteur et délivre des informations relatives à la position du récepteur telles que par exemple l'accélération et/ou la position angulaire du récepteur, et/ou le cap par rapport au nord magnétique, ces informations étant délivrées pour les trois directions de l'espace. Pour cela, le capteur inertiel peut comporter un accéléromètre trois axes et/ou un gyroscope trois axes et/ou un magnétomètre trois axes. Les informations délivrées par le capteur inertiel sont ensuite traitées par les moyens de traitement pour obtenir l'attitude et/ou la position et/ou la vitesse du récepteur. Ces données sont alors transmises aux moyens d'acquisition pour lors de la phase d'acquisition tenir compte des positions successives du récepteur données par le capteur inertiel et limiter ainsi les hypothèses des évolutions de fréquence au cours du temps et donc le temps de traitement pendant la phase d'acquisition. En effet, une hypothèse fréquentielle correspond à une hypothèse sur l'effet Doppler et l'effet Doppler est proportionnel à la vitesse relative entre l'utilisateur et le satellite. Par conséquent, l'information donnée par le capteur inertiel est une information relative à l'évolution de l'effet Doppler pendant l'acquisition. Cette information permet de borner les possibilités de recombinaisons possibles entre les instants initial et final. La figure 4 est un exemple illustrant les chemins possibles pour une hypothèse de fréquence initiale fi donnée, en tenant compte de l'information délivrée par le capteur inertiel. Le chemin réel 3 correspondant à l'énergie maximale du signal est représenté en trait épaissi, les chemins 2 à explorer en raison de l'incertitude sur les informations issues du capteur inertiel sont représentés dans l'enveloppe d'erreur 1 en traits pointillés. En effet, les informations délivrées par un capteur inertiel présentent une marge d'erreur qui augmente, notamment, avec le temps de fonctionnement du capteur et la température. Les caractéristiques des capteurs étant données par le fabriquant, après traitement, la marge d'erreur sur les informations délivrées par le capteur peut être convertie en fréquence et une enveloppe d'erreur de la fréquence correspondant à un capteur donné peut être déterminée aisément à partir de ces caractéristiques. Ainsi, comme le montre la figure 4, le procédé consiste alors, lorsque toutes les corrélations élémentaires correspondant à chaque hypothèse de fréquence ont été réalisées, à sélectionner une hypothèse de fréquence initiale fi correspondant à un instant initial to, à suivre le chemin extrait des informations données par le capteur jusqu'à un instant tk correspondant au moment où les données délivrées par le capteur commencent à présenter une marge d'erreur, puis entre l'instant tk et l'instant final tn, à sélectionner un chemin parmi tous les chemins possibles compris dans l'enveloppe d'erreur de fréquence établie préalablement. Avantageusement, le procédé peut comporter aussi une option consistant à réaliser une discrimination de l'information donnée par le capteur lorsque celui-ci indique que le récepteur n'a pas bougé. La discrimination est réalisée en introduisant un seuil de déplacement au-dessous duquel il est décidé de suivre l'information indiquée par le capteur. Ainsi, dans le cas ou le capteur indique que, entre deux instants de mesure consécutifs, le récepteur a bougé d'une distance inférieure au seuil de déplacement, le procédé considère que le récepteur, et donc la fréquence, n'a pas bougé pendant ces deux instants de mesure. Le procédé présente les avantages de permettre la détection de toute l'énergie du signal reçu pendant tout le temps d'intégration sans pénaliser le temps de traitement et donc le temps de réponse du récepteur. En outre, le récepteur gagne en sensibilité dans sa capacité de détection des énergies faibles car le procédé permet la prise en compte d'un plus grand nombre hypothèses fréquence-temps. Bien que l'invention ait été décrite en relation avec un mode de réalisation particulier, il est bien évident qu'elle n'y est nullement limitée et qu'elle comprend tous les équivalents techniques des moyens décrits ainsi que leurs combinaisons si celles-ci entrent dans le cadre de l'invention.20 REVENDICATIONS1. procédé pour optimiser une acquisition d'un signal à spectre étalé provenant d'un satellite par un récepteur mobile, consistant : à recevoir (10) un signal émis par au moins un émetteur, le signal étant modulé par un code pseudo-aléatoire, à élaborer une réplique locale du code pseudo-aléatoire du signal émis, - à émettre une pluralité d'hypothèses relatives à la fréquence f; du signal reçu à un instant initial to et à l'évolution de cette fréquence au cours du temps jusqu'à un instant final tä correspondant à un temps d'intégration, à calculer les fonctions de corrélation s(f;, tk, T) de tous les couples d'hypothèses fréquence û temps (fi, tk), caractérisé en ce qu'il consiste en outre : à choisir une première hypothèse de fréquence initiale fi correspondant à l'instant initial to, à partir de la même hypothèse de fréquence initiale f1, à choisir successivement, les uns après les autres, tous les chemins possibles d'évolution de la fréquence fi entre les instants to et tn et à déterminer leur énergie totale respective, à choisir successivement toutes les autres hypothèses de fréquences initiales possibles et pour chaque hypothèse de fréquence initiale à sélectionner le chemin d'évolution dont l'énergie totale est la plus forte, à décider que, parmi tous les chemins sélectionnés à chaque hypothèse de fréquence, le chemin qui a l'énergie la plus forte (3) est le chemin le plus vraisemblable et correspond au chemin réel suivi par le récepteur. 2. Procédé selon la revendication 1, caractérisé en ce que l'énergie totale de chaque chemin possible d'évolution est obtenue en sommant de manière cohérente les résultats de corrélation entre les instants initial to et final tn, la fréquence prenant successivement toutes les valeurs de fréquences situées sur le chemin choisi. 3. Procédé selon l'une des revendications 1 ou 2, caractérisé en ce que pour chaque hypothèse de fréquence initiale, les chemins possibles d'évolution de la fréquence sont déterminés en prenant en compte des informations délivrées par au moins un capteur inertiel (14) localisé dans le récepteur. 4. Procédé selon la revendication 3, caractérisé en ce que les informations délivrées par le capteur inertiel (14) sont prises en compte dans une enveloppe d'erreur (1) établie à partir des caractéristiques du capteur. 5. Procédé selon la revendication 4, caractérisé en ce que lorsque la variation des informations délivrées par le capteur, entre deux instants de mesure, est inférieure à un seuil de variation prédéterminé, il est décidé que le récepteur n'a pas bougé entre les deux instants de mesure et il est alors fait l'hypothèse d'une fréquence stable lors de l'intégration.25