WO2006048394A1 - Procede et dispositif de reception d'un signal de radionavigation degrade - Google Patents

Procede et dispositif de reception d'un signal de radionavigation degrade Download PDF

Info

Publication number
WO2006048394A1
WO2006048394A1 PCT/EP2005/055544 EP2005055544W WO2006048394A1 WO 2006048394 A1 WO2006048394 A1 WO 2006048394A1 EP 2005055544 W EP2005055544 W EP 2005055544W WO 2006048394 A1 WO2006048394 A1 WO 2006048394A1
Authority
WO
WIPO (PCT)
Prior art keywords
signal
code
duration
correlation
correlator
Prior art date
Application number
PCT/EP2005/055544
Other languages
English (en)
Inventor
Jean-Pierre Raffegeau
Original Assignee
Thales
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Thales filed Critical Thales
Publication of WO2006048394A1 publication Critical patent/WO2006048394A1/fr

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/69Spread spectrum techniques
    • H04B1/707Spread spectrum techniques using direct sequence modulation
    • H04B1/7073Synchronisation aspects
    • H04B1/7075Synchronisation aspects with code phase acquisition
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B2201/00Indexing scheme relating to details of transmission systems not covered by a single group of H04B3/00 - H04B13/00
    • H04B2201/69Orthogonal indexing scheme relating to spread spectrum techniques in general
    • H04B2201/707Orthogonal indexing scheme relating to spread spectrum techniques in general relating to direct sequence modulation
    • H04B2201/70715Orthogonal indexing scheme relating to spread spectrum techniques in general relating to direct sequence modulation with application-specific features

Definitions

  • the invention relates to the reception of a radionavigation signal from a satellite positioning system such as the GPS system (acronym for the English expression "Global Positioning System”).
  • a satellite positioning system such as the GPS system (acronym for the English expression "Global Positioning System”).
  • the GPS system consists of a constellation of 28 satellites and a terrestrial network of ground-based reference stations. Each satellite gravitates 20 000 km from the earth with a period of revolution of 12 hours. They emit two signals, one at 1575.452 MHz for civil applications and the other at 1227.6 MHz for reserved access applications. In the following, we only consider the civil frequency.
  • the signal emitted by a satellite consists of a 1575.452 MHz frequency carrier, modulated by a known spreading code and possibly by unknown data also known as databits.
  • the satellites all transmit on the same frequency and the emitted signals are differentiated by their code.
  • T for example 1 ms and typically consist of 1023 chips.
  • the positioning of the receiver is obtained by measuring the distance between a satellite and the receiver from the duration of propagation of the signal between this satellite and the receiver.
  • a replica of the transmitted code is generated locally; the offset between the received signal and the local signal (i.e., the replica) corresponds to the desired propagation time.
  • This offset is measured by phasing the received signal and the local signal; the phasing criterion corresponds to the maximization of the correlation function of the two signals, that is to say to the search for a peak of correlation. Correlation calculations are performed from the real and imaginary components of the received signal.
  • this computation is carried out by integrating over several periods of code, the correlations of the same component (real or imaginary) and the same offset between the signal and the code.
  • the coherent integration time is less than 100ms and more typically set at 20ms, that is, the period of the databits.
  • the reception sensitivity for a coherent duration of 20 ms is located around -144 dBm according to the probabilities of absence and false detections of the radionavigation signal.
  • a so-called non-coherent integration performed over a non-coherent duration Tnc is also generally performed, starting from the coherent integrations for the two components I and Q calculated over a single code period T.
  • the non-coherent integration can be achieved by means of a dedicated integrated circuit or a microprocessor.
  • the non-coherent integration realized by means of a dedicated integrated circuit makes it possible to add only coherent integrations computed over a single period of code T.
  • the non-coherent integration of coherent integrations computed over a period greater than the code period T can only be performed by successive readings of the coherent integrations and an addition; the signal is then stored on 1 to 20 ms depending on the number of readings required, which requires a lot of resources of the microprocessor.
  • An important object of the invention is therefore to be able to increase the integration time and consequently the sensitivity of the receiver, especially when the reception is degraded.
  • the invention proposes a method for receiving a radionavigation signal, which comprises, for a first determined duration, a correlation calculation step of the sampled received signal with a locally generated period code T, the signal and the code being shifted by n relative to each other by a correlation to the next correlation.
  • the correlation calculation step comprises a step of storing the signal in a memory during a determined second duration and in that the signal is written in memory according to a write frequency fe and read according to a reading frequency such that, for the correlation calculation step, the sampled signal comprising a real component I and an imaginary component Q, it comprises the steps of adding the results of the correlation calculations of the same component and the same offset n, and obtained over a duration Tc equal to k c code periods, k c being an integer greater than 1, so as to obtain for each component, a coherent integration respectively designated
  • the coherent integration is calculated by integrating over several periods of code, the correlations of the same component (real or imaginary) and the same offset n between the signal and the code.
  • the correlation calculation step is performed for several codes and the memorized signal is read for each code. Similarly, the correlation calculation step is performed for several frequency hypotheses considered to compensate for the Doppler effect and the memorized signal is read in memory for each of these hypotheses.
  • the subject of the invention is also a receiver of a radionavigation signal, characterized in that it comprises a microprocessor capable of implementing the method as described.
  • Such a receiver makes it possible to mix non-coherent and coherent integrations, the latter being calculated over a period greater than the period of the code T and to achieve this mixing rapidly.
  • Tc 10 ms
  • the configuration flexibility between the coherent duration and the number of coherent integrations makes it possible to optimize the detection power of the weak signals.
  • it comprises connected to the correlator a third memory capable of storing the received signal at a write frequency fe and to restore it to the correlator at a read frequency f1 such that f1> fe.
  • the correlator comprises a non-rotating shift register for the code and at least one non-rotating shift register for the signal.
  • FIG. 1 schematically represents an example of a known correlator with a rotary shift register for the code, and a non-rotational shift register for the signal
  • Figure 2 schematically shows an example of a correlator with a non-rotating shift register for the code, and a slip shift register for the signal
  • FIG. 3 diagrammatically represents an example of improvement of the correlator of FIG. 2
  • FIG. 4 schematically represents an example of a receiver according to the invention able to implement a mix of coherent and non-coherent integrations with post-processing.
  • the coherent integration is calculated by integrating over several periods of code, the correlations of the same component (real or imaginary) and the same offset n between the signal and the code.
  • the modification of the offset n is obtained by shifting the code, the signal remaining fixed.
  • This technique proceeds in two phases using for each component of the signal, a structure such as that of Figure 1.
  • N of the shift register represents the number of code chips considered over the period of the code, ie 1 ms for example for the GPS. Then the signal remains fixed in its shift register 1 and the code performs N rotations on itself in the rotating shift register 2. A new correlation is calculated for each rotation, by means of N multipliers 3 and an adder 4 at N entries.
  • the first phase loads the signal over a period of code.
  • the second phase calculates a correlation with each rotation of the code on itself.
  • the 2 real and imaginary components of the signal are processed using a shift register and can for example be expressed as follows:
  • SRr (i, t k ) are the samples of the real component of the signal at the instant ix 1 ms / 2046 + t k ,
  • Code (i, n) is the index code (i + n) modulo 2046, i the index of the shift register varying from 0 to 2045, n the delay index between the signal and the code, varying from 0 to 2045.
  • V [Code (i, n) * SRr (i, t k )] Code (0, n) * SRr (O, t k ) + Code (l, n) * SRr (I, t k ) + .. .
  • V [Code (i, n) * SRi (i, t k )] Code (0, n) * SRi (O, t k ) + Code (l, n) * SRi (I, t k ) + .. .
  • the code is now fixed and only the signal is shifted: it slides in its shift register at each clock cycle.
  • the slice of signal is always equal to a period of code but this slice is slippery.
  • This embodiment is implemented using for each component of the signal, a structure such as that of FIG. 2.
  • the code is entered in a non-rotating shift register 2 'and the component of the compensated signal of its Doppler is continuously fed into a sliding shift register 1.
  • the length N of the shift register represents the number of chips of the code considered on the period of the code, ie 1 ms for example for the GPS.
  • the two structures for the real and imaginary components are combined as shown in the correlator example of FIG. correlator 20 is provided with a non-rotating shift register 2 'for the code, common to both components, of length N.
  • the other elements were already present in the example of FIG. 3, namely a shift register 1 of length N for each real and imaginary component of the signal (respectively designated signal I and signal Q in the figure), 2N multipliers 3 (N for each component) and 2 adders 4 to N inputs each (1 adder for each component).
  • SRr i, t k , n
  • SRr are the samples of the real component of the signal at the instant ix 1 ms / 2046 + t k , n ,
  • tk, n represents the (k + 1) th code period, n varying from 0 to 2045 and represents the delay index between the signal and the code.
  • V [Code (i) * SRi (i, t kn )] Code (0) * SRi (O, t kn ) + Code (l) * SRi (I, t kn ) + ...
  • the process according to the invention proceeds in 3 steps.
  • the first step calculates a coherent integration over a duration Tc with its 2 real and imaginary components.
  • the second step calculates an amplitude or a power from the coherent integrations on the 2 components.
  • the third step calculates over a duration Tnc the sum of the amplitudes or powers having the same offset.
  • the coherent integration (1 st step) is performed for a duration Tc, multiple of the period of the code, in the multiple occurrence of 1 ms; it is performed separately on both the real and imaginary components.
  • Code (i, n) is the index code (i + n) modulo 2046 as previously indicated.
  • This coherent integration can be obtained by means of a microprocessor capable of implementing the preceding calculations or by means of an adder 6 associated with a memory 7 with 2N boxes able to be read and then written such as a RAM memory. In known manner, the memory 7 is read and written as the results provided by the correlator 20.
  • the coherent integration requires taking into account the phase of the data bits by means of a known device 5: the result of each correlation is multiplied by ⁇ 1 according to the value of the databit on the code period during which the correlations are made .
  • This second step is obtained in known manner by a microprocessor or an integrated circuit 8.
  • the non-coherent integration consists of calculating over the period of the code, the amplitude or the power of the correlations on the two components, and then summing them.
  • This non-coherent integration can be obtained by means of a microprocessor capable of implementing the preceding calculations or by means of an adder 9 associated with a memory 10 with N boxes able to be read and written such as a RAM memory. In known manner, the memory 10 is read and then written as the results provided by the memory 7.
  • the signal is sampled and processed for a satellite and a frequency assumption (considered to compensate for the Doppler effect) at a time.
  • the samples are directly used for the processing of the signal which lasts Tnc.
  • the signal is sampled and stored in memory 11 preferably during the non-coherent duration Tnc which can reach 16 s. Then the samples are read into memory for each satellite and each Doppler frequency hypothesis; that is, the I and Q signals from the memory 11 are converted into baseband and compensated for a frequency Doppler by the element 12 before being inserted into the shift register of the correlator 20. The samples are read at a frequency f 1 much higher than their memory write frequency.
  • This post-processing mode makes it possible to reduce the overall search time of the satellites as it appears on the following comparison.
  • NBSAT the number of satellites
  • BPDOPPLER the frequency band to be considered to compensate the Doppler effect for the search for Tnc satellites: the duration of non-coherent integration
  • T tr the total real-time processing time
  • Tpt the overall processing time post-processing
  • Ttr NBSAT x BPDOPPLER x Tc x Tnc
  • Tpt NBSAT x BPDOPPLER x Tc x Tnc / K + Tnc,
  • Tpt Ttr / K + Tnc
  • the received signal may be a satellite and / or pseudolite radionavigation signal.
  • GPS Global Positioning System
  • GALILEO Global System for Mobile Communications

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

L'invention concerne un procédé de réception d'un signal de radionavigation, qui comprend sur une première durée déterminée une étape de calculs de corrélation du signal reçu échantillonné avec un code de période T généré localement, le signal et le code étant décalés de n l'un par rapport à l'autre d'une corrélation à la corrélation suivante. Préalablement à l'étape de calculs de corrélation, il comporte une étape de mémorisation du signal dans une mémoire pendant une deuxième durée déterminée et le signal est écrit en mémoire selon une fréquence d'écriture fe et lu selon une fréquence de lecture fl telle que fl>fe en vue de l'étape de calculs de corrélation. Le signal échantillonné comportant une composante réelle I et une composante imaginaire Q, il comprend les étapes consistant à additionner les résultats des calculs de corrélation de même composante et de même décalage n, et obtenus sur une durée Tc égale à k<SUB>c</SUB> périodes de code T, k<SUB>c</SUB> étant un entier supérieur à 1, de manière à obtenir pour chaque composante, une intégration cohérente respectivement désignée I<SUB>TC</SUB>(n) et Q<SUB>RC</SUB>(n), calculer la somme STc(n)=I<SUB>TC</SUB>(n)<SUP>2</SUP> + Q<SUB>RC</SUB>(n)<SUP>2</SUP>, et additionner les sommes vS<SUB>TC</SUB>(n) de même décalage n, obtenues sur une durée Tnc, Tnc = K x Tc, K étant un entier supérieur à 1, de manière à obtenir pour chaque décalage n, une intégration non cohérente C(n).

Description

PROCEDE ET DISPOSITIF DE RECEPTION D'UN SIGNAL DE RADIONAVIGATION DEGRADE
L'invention concerne la réception d'un signal de radionavigation provenant d'un système de positionnement par satellites tel que le système GPS (acronyme de l'expression anglo-saxonne « Global Positioning System »).
Pour bien fonctionner, les récepteurs actuels nécessitent généralement une réception à vue directe des satellites. Le positionnement se détériore rapidement notamment en termes de précision et de durée d'acquisition, lorsque la réception est perturbée comme c'est le cas à l'intérieur d'un bâtiment ou plus généralement en milieu dégradé.
On rappelle succinctement le fonctionnement du système GPS. Il est constitué d'une constellation de 28 satellites et d'un réseau terrestre de stations de référence à terre. Chaque satellite gravite à 20 000 km de la terre avec une période de révolution de 12h. Ils émettent deux signaux, l'un à 1575.452 MHz pour les applications civiles et l'autre à 1227.6 MHz pour les applications à accès réservés. Dans la suite, on ne considère que la fréquence civile. Le signal émis par un satellite est constitué d'une porteuse de fréquence 1575.452 MHz, modulée par un code d'étalement connu et éventuellement par des données inconnues aussi désignées databits. Les satellites émettent tous sur la même fréquence et les signaux émis se différencient par leur code.
Ces codes présentent une période T, par exemple 1 ms et sont typiquement constitués de 1023 chips.
Le positionnement du récepteur est obtenu en mesurant la distance entre un satellite et le récepteur à partir de la durée de propagation du signal entre ce satellite et le récepteur. Dans le récepteur, une réplique du code émis est générée localement ; le décalage entre le signal reçu et le signal local (c'est-à-dire la réplique) correspond à la durée de propagation recherchée. Ce décalage est mesuré en mettant en phase le signal reçu et le signal local ; le critère de mise en phase correspond à la maximisation de la fonction de corrélation des deux signaux, c'est-à-dire à la recherche d'un pic de corrélation. Les calculs de corrélation sont effectués à partir des composantes réelle et imaginaire du signal reçu. Lorsque la réception est dégradée, le calcul d'un pic de corrélation sur une seule période de code T n'est pas possible : ce calcul est effectué en intégrant sur plusieurs périodes de code, les corrélations de même composante (réelle ou imaginaire) et de même décalage entre le signal et le code.
Ce calcul de corrélation est généralement effectué de demi-chip de code en demi-chip de code sur un intervalle d'intégration que l'on peut faire varier. Pour un intervalle d'intégration de 1 ms, on obtient une durée de calcul de corrélation d'environ 2s. (2 x1023 x 1 ms = 2s). Ce calcul est multiplié par un facteur P fonction de la dérive de l'horloge locale (ou pilote) du récepteur et du nombre d'hypothèses sur la fréquence du signal à considérer pour compenser l'effet Doppler. Pour une incertitude d'horloge d'environ ±10 kHz, P=21 ; on obtient alors de calcul d'environ 2s x 21 , soit 42s par code c'est-à-dire par satellite. En prenant un intervalle de 10ms (c'est-à-dire 10 périodes de code), on gagne environ 10 dB, mais le calcul dure alors environ 20s multiplié par un facteur P ≡201 , soit environ 4020 s par satellite. On n'a finalement gagné que 10 dB pour une durée de positionnement qui a considérablement augmenté. Une solution pour réduire cette durée de calcul est d'effectuer les
2 x 1023 multiplications en parallèle au lieu de les calculer en série comme précédemment. On obtient alors une durée de calcul sur un intervalle d'intégration de 10 ms d'environ 2s (201 x 10 ms=2s).
Dans la suite on se place dans le cadre de cette solution. Les intégrations sont effectuées à partir des composantes réelle et imaginaire du signal échantillonné. On peut effectuer une intégration dite cohérente sur une durée cohérente Tc éventuellement supérieure à T.
Une durée cohérente important permet d'avoir une bande passante étroite et de faire ressortir le signal du bruit. Par contre, il n'est pas possible d'augmenter fortement la durée cohérente. Celle-ci se trouve très vite limitée par : le nombre d'hypothèses de fréquences à considérer pour compenser l'effet Doppler, l'accélération potentielle du mobile, les caractéristiques et le prix de l'oscillateur pilote dans le récepteur.
En raison de ces limitations, la durée de l'intégration cohérente est inférieure à 100ms et plus généralement fixée à 20ms, c'est-à-dire à la période des databits. La sensibilité de réception pour une durée cohérente de 20ms est située autour de -144dBm suivant les probabilités d'absence et de fausses détections du signal de radionavigation.
Pour repousser ces limites, on effectue aussi en général une intégration dite non cohérente réalisée sur une durée non cohérente Tnc, à partir des intégrations cohérentes pour les deux composantes I et Q calculées sur une seule période de code T.
L'intégration non cohérente peut être réalisée au moyen d'un circuit intégré dédié ou d'un microprocesseur.
L'intégration non cohérente réalisée au moyen d'un circuit intégré dédié, ne permet d'additionner que des intégrations cohérentes calculées sur une seule période de code T.
En utilisant un microprocesseur, l'intégration non cohérente d'intégrations cohérentes calculées sur une durée supérieure à la période de code T ne peut être effectuée que par des lectures successives des intégrations cohérentes et une addition ; le signal est alors mémorisé sur 1 à 20 ms selon le nombre de lectures nécessaires, ce qui nécessite beaucoup de ressources du microprocesseur.
Un but important de l'invention est donc de pouvoir accroître la durée d'intégration et par conséquent la sensibilité du récepteur, notamment lorsque la réception est dégradée.
Pour atteindre ce but, l'invention propose un procédé de réception d'un signal de radionavigation, qui comprend sur une première durée déterminée une étape de calculs de corrélation du signal reçu échantillonné avec un code de période T généré localement, le signal et le code étant décalés de n l'un par rapport à l'autre d'une corrélation à la corrélation suivante.
Il est principalement caractérisé en ce que préalablement à l'étape de calculs de corrélation, il comporte une étape de mémorisation du signal dans une mémoire pendant une deuxième durée déterminée et en ce que le signal est écrit en mémoire selon une fréquence d'écriture fe et lu selon une fréquence de lecture fl telle que fl>fe en vue de l'étape de calculs de corrélation, le signal échantillonné comportant une composante réelle I et une composante imaginaire Q, il comprend les étapes consistant à additionner les résultats des calculs de corrélation de même composante et de même décalage n, et obtenus sur une durée Tc égale à kc périodes de code, kc étant un entier supérieur à 1 , de manière à obtenir pour chaque composante, une intégration cohérente respectivement désignée
Iτc(n) et Qτc(n), calculer la somme Sτc(n)=lτc(n)2 + Qτc(n)2 , et additionner les sommes Sτc(n) ou sJSτc (n) de même décalage n, obtenues sur une durée Tnc, Tnc = K x Tc, K étant un entier supérieur à 1 , de manière à obtenir pour chaque décalage n, une intégration non cohérente C(n).
Ainsi, on mixe une intégration cohérente sur une durée cohérente Tc égale à kc périodes de code T, kc étant un entier supérieur à 1 et une intégration non cohérente sur une durée Tnc, avec Tnc = K x Tc .
L'intégration cohérente est calculée en intégrant sur plusieurs périodes de code, les corrélations de même composante (réelle ou imaginaire) et de même décalage n entre le signal et le code. L'intégration non cohérente est calculée en effectuant la somme sur une durée Tnc égale à un multiple de Tc (Tnc = K x Tc), des amplitudes ou puissances sur les deux composantes des intégrations cohérentes ayant le même décalage n.
Selon une caractéristique de l'invention, l'étape de calculs de corrélation est réalisée pour plusieurs codes et le signal mémorisé est lu pour chaque code. De même, l'étape de calculs de corrélation est réalisée pour plusieurs hypothèses de fréquence considérées pour compenser l'effet Doppler et le signal mémorisé est lu en mémoire pour chacune de ces hypothèses.
L'invention a aussi pour objet, un récepteur d'un signal de radionavigation, caractérisé en ce qu'il comprend un microprocesseur apte à mettre en œuvre le procédé tel que décrit. Enfin, l'invention concerne également un récepteur d'un signal de radionavigation pourvu d'un corrélateur apte à corréler le signal reçu qui comporte une composante réelle et une composante imaginaire, avec un code de longueur N et de période T généré localement, caractérisé en ce qu'il comporte relié au corrélateur, une première mémoire et un premier additionneur aptes à mémoriser 2N résultats issus du corrélateur pour les deux composantes et à les incrémenter respectivement au fur et à mesure des résultats fournis par le corrélateur, pendant une durée Tc égale à kc périodes de code, kc étant un entier supérieur à 1 ,pour réaliser 2N intégrations cohérentes respectivement désignées pour chaque composante lτc et QTC , des moyens pour obtenir N amplitudes ou puissances de sommes Sτc(n)=lτc(n)2 + Qτc(n)2, une deuxième mémoire et un deuxième additionneur aptes à mémoriser les N amplitudes ou puissances des sommes STc(n), et à les incrémenter respectivement pendant une durée Tnc, Tnc = K x Tc, K étant un entier supérieur à 1 , au fur et à mesure que les amplitudes ou les puissances des sommes Sτc(n) sont obtenues, pour réaliser N intégrations non cohérentes.
Un tel récepteur permet de mixer des intégrations non cohérentes et cohérentes, ces dernières étant calculées sur une durée supérieure à la période du code T et de réaliser ce mixage rapidement. Avec Tc = 10 ms, on gagne plus de 30% en durée de traitement du signal. La souplesse de configuration entre la durée cohérente et le nombre d'intégrations cohérentes permet d'optimiser le pouvoir de détection des faibles signaux.
Avantageusement, il comprend relié au corrélateur une troisième mémoire apte à mémoriser le signal reçu selon une fréquence d'écriture fe et à le restituer vers le corrélateur selon une fréquence de lecture fl telle que fl>fe.
Selon une caractéristique de l'invention, le corrélateur comprend un registre à décalages non tournant pour le code et au moins un registre à décalages non tournant pour le signal. D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de la description détaillée qui suit, faite à titre d'exemple non limitatif et en référence aux dessins annexés dans lesquels : la figure 1 représente schématiquement un exemple de corrélateur connu avec un registre à décalages tournant pour le code, et un registre à décalages non tournant pour le signal, la figure 2 représente schématiquement un exemple de corrélateur avec un registre à décalages non tournant pour le code, et un registre à décalages glissant pour le signal, la figure 3 représente schématiquement un exemple de perfectionnement du corrélateur de la figure 2, la figure 4 représente schématiquement un exemple de récepteur selon l'invention apte à mettre en œuvre un mixage d'intégrations cohérentes et non cohérentes avec post-traitement.
Selon un mode de réalisation de l'invention, on mixe une intégration cohérente sur une durée cohérente Tc une durée Tc égale à kc périodes de code T, kc étant un entier supérieur à 1 , et une intégration non cohérente sur une durée Tnc, avec Tnc = K x Tc . L'intégration cohérente est calculée en intégrant sur plusieurs périodes de code, les corrélations de même composante (réelle ou imaginaire) et de même décalage n entre le signal et le code.
On va dans un premier temps détailler ces corrélations sur un exemple en prenant N = 2046 = 2 x 1023 car le calcul de corrélation est effectué de demi-chip de code en demi-chip de code. Les chips de code peuvent être divisés par 2 comme dans cet exemple ou par 3, 4, ..., voire même 3.3, etc.
Selon un premier mode de réalisation, d'une corrélation à l'autre, la modification du décalage n est obtenue en décalant le code, le signal restant fixe. Cette technique procède en 2 phases en utilisant pour chaque composante du signal, une structure comme par exemple celle de la figure 1.
Le signal reçu dont la fréquence est compensée pour prendre en compte l'effet Doppler, est chargé dans un registre à décalages non tournant
1 et le code est chargé dans un registre à décalages tournant 2. La longueur N du registre à décalages représente le nombre de chips du code considérés sur la période du code, soit 1 ms par exemple pour le GPS. Puis le signal reste fixe dans son registre à décalages 1 et le code effectue N rotations sur lui-même dans le registre à décalages tournant 2. Une nouvelle corrélation est calculée pour chaque rotation, au moyen de N multiplieurs 3 et d'un additionneur 4 à N entrées.
La première phase charge le signal sur une période du code. La seconde phase calcule une corrélation à chaque rotation du code sur lui- même.
Les 2 composantes réelle et imaginaire du signal sont traitées à l'aide d'un registre à décalages et peuvent être par exemple exprimées de la façon suivante :
SRr(i, tk) sont les échantillons de la composante réelle du signal à l'instant i x 1 ms / 2046 + tk,
SRi(i, tk) sont les échantillons de la composante imaginaire du signal à l'instant i x 1 ms / 2046 + tk, i l'indice du registre à décalages variant de 0 à 2045, tk = to + k x 1 ms. tk représente la (k+1 )ième période de code. Code(i, n) est le code à l'indice (i + n) modulo 2046, i l'indice du registre à décalages variant de 0 à 2045, n l'indice de retard entre le signal et le code, variant de 0 à 2045.
A la k+1 ième période de code, le résultat de la corrélation sur la composante réelle, fourni par l'additionneur pour un décalage n, est :
2045
V [Code(i,n) * SRr(i, tk )] = Code(0, n) * SRr(O, tk ) + Code(l , n) * SRr(I , tk ) + ...
1=0
+ Code(2044, n) * SRr(2044, tk ) + Code(2045, n) * SRr(2045, tk )
De même, à la k+1 ième période de code, le résultat de la corrélation sur la composante imaginaire, fourni par l'additionneur pour un décalage n, est :
2045
V [Code(i, n) * SRi(i, tk )] = Code(0, n) * SRi(O, tk ) + Code(l , n) * SRi(I, tk ) + ...
1=0
+ Code(2044, n) * SRi(2044, tk ) + Code(2045, n) * SRi(2045, tk ) Selon un autre mode de réalisation, d'une corrélation à l'autre, le code est maintenant fixe et seul le signal est décalé : il glisse dans son registre à décalages à chaque cycle d'horloge. La tranche de signal est toujours égale à une période de code mais cette tranche est glissante.
Ce mode de réalisation est mis en œuvre en utilisant pour chaque composante du signal, une structure comme par exemple celle de la figure 2.
Le code est entré dans un registre à décalage non tournant 2' et la composante du signal compensé de son Doppler est rentrée en continu dans un registre à décalages glissant 1. La longueur N du registre à décalages représente le nombre de chips du code considérés sur la période du code, soit 1 ms par exemple pour le GPS.
L'un des avantages de ce corrélateur réside dans la mise en parallèle des deux phases des solutions conventionnelles (chargement du signal dans le registre à décalages et calcul des corrélations). La durée de traitement du signal, donc de recherche des satellites est ainsi divisée par 2. Selon un autre mode de réalisation, les deux structures pour les composantes réelle et imaginaire sont combinées comme représenté sur l'exemple de corrélateur de la figure 3. Ce corrélateur 20 est pourvu d'un registre à décalages non tournant 2' pour le code, commun aux deux composantes, de longueur N. Les autres éléments étaient déjà présents dans l'exemple de la figure 3, à savoir un registre à décalage 1 de longueur N pour chaque composante réelle et imaginaire du signal (respectivement désignées signal I et signal Q sur la figure), 2N multiplieurs 3 (N pour chaque composante) et 2 additionneurs 4 à N entrées chacun (1 additionneur pour chaque composante).
Les 2 composantes réelle et imaginaire du signal peuvent être par exemple exprimées de la façon suivante : SRr(i, tk,n) sont les échantillons de la composante réelle du signal à l'instant i x 1 ms / 2046 + tk,n ,
SRi(i, tk,n) sont les échantillons de la composante imaginaire du signal à l'instant i x 1 ms / 2046 + tk,n, i l'indice du registre à décalages varie de 0 à 2045, tk,n = to + k x 1ms + n x 1ms / 2046. tk,n représente la (k+1 )ième période de code, n variant de 0 à 2045 et représente l'indice de retard entre le signal et le code.
A la k+1 ième période de code, le résultat de la corrélation sur la composante réelle, fourni par l'additionneur pour un décalage n, est :
2045 y [Code(i) * SRr(i, tk n )] = Code(0) * SRr(O, tk n ) + Code(l ) * SRr(I , tk n )
L^ L + ...
I=O
+ Code(2044) * SRr(2044, tk n ) + Code(2045) * SRr(2045, tk n )
De même, à la k+1 ième période de code, le résultat de la corrélation sur la composante imaginaire, fourni par l'additionneur pour un décalage n, est :
2045
V [Code(i) * SRi(i, tk n )] = Code(0) * SRi(O, tk n ) + Code(l) * SRi(I, tk n ) + ...
1=0
+ Code(2044) * SRi(2044, tk n ) + Code(2045) * SRi(2045, tk n )
Le procédé selon l'invention procède en 3 étapes. La première étape calcule une intégration cohérente sur une durée Tc avec ses 2 composantes réelle et imaginaire. La deuxième étape calcule une amplitude ou une puissance à partir des intégrations cohérentes sur les 2 composantes. La troisième étape calcule sur une durée Tnc la somme des amplitudes ou des puissances ayant le même décalage. Les 3 étapes sont répétées N fois. Dans notre exemple N = 2046. Ce mixage est mis en œuvre par exemple par des moyens représentés figure 4.
L'intégration cohérente (1ere étape) est réalisée pendant une durée Tc, multiple de la période du code, en l'occurrence multiple de 1 ms ; elle est réalisée de manière séparée sur les deux composantes réelle et imaginaire.
Lorsque les corrélations ont été calculées à partir d'un code tournant et d'un signal fixe, on obtient les résultats suivants. Pour un décalage n entre le signal et le code, il résulte de l'intégration cohérente, l(n, Tc) (ou Iτc(n)) pour la composante réelle et Q(n, Tc) (ou Qτc(n)) pour la composante imaginaire : eo 5(n,Tc) (Code(i,n) * SRr(Mk))
Figure imgf000011_0001
Figure imgf000011_0002
avec SRr(i, tk) et SRi(i, tk) à l'instant i x 1 ms / 2046 + tk, tk = t0 + k x 1 ms.
Code(i, n) est le code à l'indice (i + n) modulo 2046 comme indiqué précédemment.
Lorsque les corrélations ont été calculées à partir d'un code fixe et d'un signal glissant, on obtient les résultats suivants. Pour un décalage n entre le signal et le code, il résulte de l'intégration cohérente, l(n, Tc) (ou Iτc(n)) pour la composante réelle et Q(n, Tc) (ou Qτc(n)) pour la composante imaginaire :
*SRr(Mk,j)
Figure imgf000011_0003
QS'(n,Tc) (Code(i)* SRi(Mkij)
Figure imgf000011_0004
avec SRr(i, tk,n) et SRi(i, tk,n) à l'instant i x 1 ms / 2046 + tk,n, tk,n = to + k x 1 ms + n x 1 ms/2046.
Cette intégration cohérente peut être obtenue au moyen d'un microprocesseur apte à mettre en œuvre les calculs précédents ou au moyen d'un additionneur 6 associé à une mémoire 7 à 2N cases apte à être lue puis écrite telle qu'une mémoire RAM. De manière connue, la mémoire 7 est lue et écrite au fur et à mesure des résultats fournis par le corrélateur 20.
Si une inversion de phase due à une modification de data bit intervient durant la durée cohérente, il peut être difficile voir impossible de détecter le pic de corrélation. L'intégration cohérente nécessite la prise en compte de la phase des data bits au moyen d'un dispositif connu 5 : le résultat de chaque corrélation est multiplié par ±1 selon la valeur du databit sur la période de code pendant laquelle les corrélations sont effectuées.
Une amplitude ou une puissance est ensuite calculée à partir des 2 composantes des intégrations cohérentes (2è étape) :
Figure imgf000012_0001
ou S(n,Tc) =[l(n, Tc)2 + Q(n, Tc)2 ]
On désigne aussi S(n,Tc) par STc(n).
Cette 2ème étape est obtenue de manière connue par un microprocesseur ou un circuit intégré 8.
Dans l'art antérieur, l'intégration non cohérente consiste à calculer sur la période du code, l'amplitude ou la puissance des corrélations sur les deux composantes, puis à les sommer. Selon l'invention, l'intégration non cohérente (3è étape) est calculée en effectuant la somme sur une durée Tnc multiple d'un nombe entier de durée cohérente Tc (Tnc = K x Tc), des amplitudes ou puissances d'intégrations cohérentes ayant le même décalage n. On obtient :
[l(n, kTc)2 + Q(n, kTc)2]
Figure imgf000012_0002
Cette intégration non cohérente peut être obtenue au moyen d'un microprocesseur apte à mettre en œuvre les calculs précédents ou au moyen d'un additionneur 9 associé à une mémoire 10 à N cases apte à être lue et écrite telle qu'une mémoire RAM. De manière connue, la mémoire 10 est lue puis écrite au fur et à mesure des résultats fournis par la mémoire 7.
Ces trois étapes sont répétées pour chaque décalage n, c'est-à- dire N fois. Dans notre exemple N = 2046.
Selon un mode de réalisation de l'invention, désigné mode temps réel, le signal est échantillonné et traité pour un satellite et une hypothèse de fréquence (considérée pour compenser l'effet Doppler) à la fois. Il n'y a pas de mémorisation du signal : les échantillons sont directement utilisés pour le traitement du signal qui dure Tnc.
Selon un mode préférentiel désigné mode post-traitement, le signal est échantillonné et sauvegardé en mémoire 11 de préférence pendant la durée non cohérente Tnc qui peut atteindre 16 s. Puis les échantillons sont lus en mémoire pour chaque satellite et chaque hypothèse de fréquence Doppler ; c'est-à-dire que les signaux I et Q issus de la mémoire 11 sont convertis en bande de base et compensés d'une fréquence Doppler par l'élément 12 avant d'être insérés dans le registre à décalages du corrélateur 20. Les échantillons sont lus à une fréquence fl bien supérieure à leur fréquence fe d'écriture en mémoire. Ce mode post¬ traitement permet de réduire la durée de recherche globale des satellites comme il apparaît sur le comparatif suivant.
On désigne :
NBSAT : le nombre de satellites
BPDOPPLER : la bande de fréquences à considérer pour compenser l'effet Doppler pour la recherche des satellites Tnc : la durée d'intégration non cohérente
Tc : la durée d'intégration cohérente
K : le rapport des fréquences d'écriture et de lecture en post traitement T tr : la durée de traitement global temps réel Tpt : la durée de traitement globale post-traitement
On a :
Ttr = NBSAT x BPDOPPLER x Tc x Tnc Tpt = NBSAT x BPDOPPLER x Tc x Tnc / K + Tnc ,
Tpt = Ttr / K + Tnc
Tnc est négligeable devant Ttr / K ; on a typiquement K = 25 ou plus.
Dans l'exemple décrit, on a considéré des signaux à deux états ( ±1 ) ; bien sûr, l'invention s'applique de la même façon à des signaux à plus de deux états.
Jusqu'à présent on a considéré qu'un seul code était généré localement. Lorsque l'on souhaite considérer d'autres codes, ceux-ci peuvent être considérés en série en reprogrammant le code local et les fréquences Doppler associées, au fur et à mesure que l'on change de code. Ils peuvent être considérés en parallèle : il suffit alors, selon une méthode connue de l'homme du métier, de pourvoir le récepteur de structures parallèles, c'est-à-dire de multiplier la structure comportant le corrélateur (ainsi que la mémoire pour le post-traitement si nécessaire, et/ou les intégrateurs si nécessaire), par le nombre de codes locaux à considérer. Le signal reçu peut être un signal de radionavigation par satellites et/ou par pseudolites.
Il s'agit par exemple d'un signal de positionnement de type GPS (acronyme de l'expression anglo-saxonne « Global Positioning System »), ou GALILEO ou tout autre système de positionnement.

Claims

REVENDICATIONS
1. Procédé de réception d'un signal de radionavigation, qui comprend sur une première durée déterminée une étape de calculs de corrélation du signal reçu échantillonné avec un code de période T généré localement, le signal et le code étant décalés de n l'un par rapport à l'autre d'une corrélation à la corrélation suivante, caractérisé en ce que, préalablement à l'étape de calculs de corrélation, il comporte une étape de mémorisation du signal dans une mémoire pendant une deuxième durée déterminée et en ce que le signal est écrit en mémoire selon une fréquence d'écriture fe et lu selon une fréquence de lecture fl telle que fl>fe en vue de l'étape de calculs de corrélation, le signal échantillonné comportant une composante réelle I et une composante imaginaire Q, il comprend les étapes consistant à additionner les résultats des calculs de corrélation de même composante et de même décalage n, et obtenus sur une durée Tc égale à ko périodes de code T, ko étant un entier supérieur à 1 ,de manière à obtenir pour chaque composante, une intégration cohérente respectivement désignée lTc(n) et Qτc(n), calculer la somme STc(n)=lτc(n)2 + Qτc(n)2 , et additionner les sommes
Figure imgf000015_0001
de même décalage n, obtenues sur une durée Tnc, Tnc = K x Tc, K étant un entier supérieur à 1 , de manière à obtenir pour chaque décalage n, une intégration non cohérente C(n).
2. Procédé selon la revendication précédente, caractérisé en ce que l'étape de calculs de corrélation est réalisée pour plusieurs codes et en ce que le signal mémorisé est lu pour chaque code.
3. Procédé selon l'une des revendications précédentes, caractérisé en ce que l'étape de calculs de corrélation est réalisée pour plusieurs hypothèses de fréquence considérées pour compenser l'effet Doppler et en ce que le signal mémorisé est lu en mémoire pour chacune de ces hypothèses.
4. Procédé selon l'une des revendications précédentes, caractérisé en ce que la deuxième durée déterminée est égale à Tnc.
5. Procédé selon l'une des revendications précédentes, caractérisé en ce que sur la première durée déterminée, le signal est décalé d'une corrélation à la suivante, le code étant fixe.
6. Procédé selon l'une des revendications précédentes, caractérisé en ce que la première durée déterminée est égale à la période du code T.
7. Procédé selon l'une des revendications précédentes, caractérisé en ce que le signal reçu est un signal de radionavigation par satellites et/ou par pseudolites.
8. Récepteur d'un signal de radionavigation, caractérisé en ce qu'il comprend un microprocesseur apte à mettre en œuvre le procédé selon l'une des revendications 1 à 7.
9. Récepteur d'un signal de radionavigation pourvu d'un corrélateur (20) apte à corréler le signal reçu qui comporte une composante réelle et une composante imaginaire, avec un code de longueur N et de période T généré localement, caractérisé en ce qu'il comporte relié au corrélateur (20), une première mémoire (7) et un premier additionneur (6) aptes à mémoriser 2N résultats issus du corrélateur pour les deux composantes et à les incrémenter respectivement au fur et à mesure des résultats fournis par le corrélateur (20), pendant une durée Tc égale à kc périodes de code T, kc étant un entier supérieur à 1 ,pour réaliser 2N intégrations cohérentes respectivement désignées pour chaque composante lτc et QTC , des moyens (8) pour obtenir N amplitudes ou puissances de sommes STc(n)=lτc(n)2 + Qτc(n)2, une deuxième mémoire (10) et un deuxième additionneur (9) aptes à mémoriser les N amplitudes ou puissances des sommes STc(n), et à les incrémenter respectivement pendant une durée Tnc, Tnc = K x Tc, K étant un entier supérieur à 1 , au fur et à mesure que les amplitudes ou les puissances des sommes Sτc(n) sont obtenues, pour réaliser N intégrations non cohérentes.
10. Récepteur selon la revendication précédente, caractérisé en ce qu'il comprend relié au corrélateur (20) une troisième mémoire (11 ) apte à mémoriser le signal reçu selon une fréquence d'écriture fe et à le restituer vers le corrélateur selon une fréquence de lecture fl telle que fl>fe.
11. Récepteur selon l'une des revendications 9 à 10, caractérisé en ce que le signal comportant des databits présentant une phase, il comprend entre le corrélateur (20) et le premier additionneur (6), un dispositif (5) apte à prendre en compte la phase des databits.
12. Récepteur selon l'une des revendications 10 à 11 , caractérisé en ce que le corrélateur comprend un registre à décalages non tournant (2') pour le code et au moins un registre à décalages non tournant (1) pour le signal.
PCT/EP2005/055544 2004-11-05 2005-10-25 Procede et dispositif de reception d'un signal de radionavigation degrade WO2006048394A1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR0411812 2004-11-05
FR0411812A FR2877803B1 (fr) 2004-11-05 2004-11-05 Procede et dispositif de reception d'un signal de radionavigation degrade

Publications (1)

Publication Number Publication Date
WO2006048394A1 true WO2006048394A1 (fr) 2006-05-11

Family

ID=34954070

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/EP2005/055544 WO2006048394A1 (fr) 2004-11-05 2005-10-25 Procede et dispositif de reception d'un signal de radionavigation degrade

Country Status (2)

Country Link
FR (1) FR2877803B1 (fr)
WO (1) WO2006048394A1 (fr)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5768306A (en) * 1993-09-06 1998-06-16 Ntt Mobile Communications Network, Inc. Sliding correlator used in CDMA systems to establish initial synchronization
US6118808A (en) * 1996-07-12 2000-09-12 General Electric Company GPS receiver with efficient signal acquisition
US6363108B1 (en) * 1999-03-31 2002-03-26 Qualcomm Inc. Programmable matched filter searcher

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5768306A (en) * 1993-09-06 1998-06-16 Ntt Mobile Communications Network, Inc. Sliding correlator used in CDMA systems to establish initial synchronization
US6118808A (en) * 1996-07-12 2000-09-12 General Electric Company GPS receiver with efficient signal acquisition
US6363108B1 (en) * 1999-03-31 2002-03-26 Qualcomm Inc. Programmable matched filter searcher

Also Published As

Publication number Publication date
FR2877803A1 (fr) 2006-05-12
FR2877803B1 (fr) 2007-04-20

Similar Documents

Publication Publication Date Title
EP2140285B1 (fr) Procédé de traitement de signaux de radionavigation
EP2520949B1 (fr) Dispositif de réception d&#39;un système de positionnement par satellite comprenant une fonction de détection de faux accrochages
EP2588884B1 (fr) Procédé de réduction de l&#39;éblouissement d&#39;un récepteur recevant des signaux depuis des émetteurs
EP1514131B1 (fr) Recepteur de positionnement par satellite utilisant deux proteuses de signal
EP3709058A1 (fr) Procédé de contrôle d&#39;intégrité d&#39;un signal de radionavigation par satellite et recepteur apte à implementer ce procédé
EP2095150B1 (fr) Procédé et dispositif de réception d&#39;un signal de radionavigation à modulation boc
EP1613976A1 (fr) Procede d&#39;acquisition d&#39;un signal de radionavigation par satellite
EP1628139B1 (fr) Procédé de détermination de position d&#39;un récepteur de signaux radiofréquences sans connaissance de la position initiale, et récepteur pour sa mise en oeuvre
EP2487507B1 (fr) Procédé d&#39;acquisition avec multi-corrélation fréquentielle améliorée
EP3537181B1 (fr) Système de suivi de signal de système de navigation satellite mondial à taux minimal de faible puissance
EP2221628B1 (fr) Procédé de lever d&#39;ambiguïté, procédé de localisation d&#39;un récepteur par radionavigation comprenant une étape de lever d&#39;ambiguïté et récepteur de localisation
EP3168648B1 (fr) Procede de detection d&#39;une fausse synchronisation d&#39;un recepteur avec un satellite, module et produit programme d&#39;ordinateur associes
EP1170597B1 (fr) Récepteur GPS basse puissance avec redistribution de tâches de synchronisation
WO2006048394A1 (fr) Procede et dispositif de reception d&#39;un signal de radionavigation degrade
EP1807711A2 (fr) Procede et dispositif de reception d&#39;un signal de radionavigation degrade
WO2009071468A1 (fr) Procede de determination d&#39;une position d&#39;un recepteur de signaux radiofrequences basse puissance
EP1265371B1 (fr) Récepteur et procédé de reception de signaux radiofréquences avec des moyens de correction des effets des signaux à trajets multiplex
WO2021255386A2 (fr) Récepteur de signaux de radionavigation comprenant un calculateur d&#39;indicateur de puissance de corrélation
EP1255123B1 (fr) Récepteur de signaux radiofréquences ayant des moyens pour améliorer la dynamique de réception desdits signaux
FR3137763A1 (fr) Procédé de détection d&#39;un signal interférant d&#39;un récepteur GNSS et dispositif de détection associé
EP2724179B1 (fr) Optimisation d&#39;un systeme de localisation comprenant des stations de base emettrices terrestres

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BW BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE EG ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KM KP KR KZ LC LK LR LS LT LU LV LY MA MD MG MK MN MW MX MZ NA NG NI NO NZ OM PG PH PL PT RO RU SC SD SE SG SK SL SM SY TJ TM TN TR TT TZ UA UG US UZ VC VN YU ZA ZM ZW

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): BW GH GM KE LS MW MZ NA SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IS IT LT LU LV MC NL PL PT RO SE SI SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG

DPE1 Request for preliminary examination filed after expiration of 19th month from priority date (pct application filed from 20040101)
121 Ep: the epo has been informed by wipo that ep was designated in this application
NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 05801507

Country of ref document: EP

Kind code of ref document: A1