WO2003054576A1 - Procede d'amelioration de la determination de l'attitude d'un vehicule a l'aide de signaux de radionavigation par satellite - Google Patents

Procede d'amelioration de la determination de l'attitude d'un vehicule a l'aide de signaux de radionavigation par satellite Download PDF

Info

Publication number
WO2003054576A1
WO2003054576A1 PCT/FR2002/004397 FR0204397W WO03054576A1 WO 2003054576 A1 WO2003054576 A1 WO 2003054576A1 FR 0204397 W FR0204397 W FR 0204397W WO 03054576 A1 WO03054576 A1 WO 03054576A1
Authority
WO
WIPO (PCT)
Prior art keywords
signals
antennas
phase
attitude
intercorrelation
Prior art date
Application number
PCT/FR2002/004397
Other languages
English (en)
Inventor
Marc Chenus
Marc Revol
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
Priority to EP02805393A priority Critical patent/EP1459093A1/fr
Priority to CA002470550A priority patent/CA2470550A1/fr
Priority to US10/499,199 priority patent/US7454289B2/en
Publication of WO2003054576A1 publication Critical patent/WO2003054576A1/fr

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/30—Acquisition or tracking or demodulation of signals transmitted by the system code related
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42—Determining position
    • G01S19/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
    • G01S19/44—Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/53—Determining attitude
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/22—Multipath-related issues

Definitions

  • the present invention relates to a method for improving the determination of the attitude of a vehicle using radio navigation signals by satellite.
  • the determination of attitude using radio navigation signals consists in comparing the differences in satellite-carrier distances with respect to a known base of antennas. A very high measurement resolution is therefore necessary. To achieve this resolution, phase measurements of the GPS signal are used.
  • attitude determination by multi-antenna GPS thus carried out requires improvements to achieve the high level of integrity required by “safety-of-life” applications, that is to say with vital safety requirements.
  • attitude determination uses the phase measurements delivered by the loop.
  • carrier tracking Phase Locked Loop, PLL
  • PLL Phase Locked Loop
  • the object of the invention is therefore to improve the determination of the attitude of a vehicle from a radio navigation signal, in particular a GPS signal.
  • the attitude of a carrier here includes the angles of heading, pitch and roll in a local landmark.
  • the GPS signal is received simultaneously by a pair (or pairs) of antennas, making it possible to establish the phase difference of the signals emitted by the GPS satellites and measured respectively by each of the two antennas of the pair considered.
  • a GPS receiver tracks the GPS signal in code using a DLL (thus measuring the pseudo-distance according to the position of the code received) and in carrier using a PLL (measuring the pseudo-speed according to the phase speed received).
  • a reception channel is dedicated to each tracked satellite.
  • the DLL is therefore essential for the localization function in relation to the PLL. This can also be temporarily replaced by speed assistance (so-called “code only” mode).
  • speed assistance so-called “code only” mode.
  • phase measurements when receiving the GPS signal, known methods are rejection of interference sources by spectral or spatial processing applied to the input signal (before demodulation by the code or the carrier).
  • temporary replacement of the PLL with speed assistance (“code only” mode) allows the DLL to be maintained for positioning in the event of strong interference, but the PLL is no longer operational, the phase measurements are no longer available for attitude measurement.
  • Phase measurement and attitude determination are therefore not always available in the event of strong interference even if GPS positioning can be maintained.
  • the known methods do not make it possible to reject all the multiple paths which may specifically degrade the determination of attitude, in particular for multiple paths with low delay not detectable by known methods at the DLL level (resolution of the DLL insufficient for the need for attitude measurement).
  • Initializing attitude determination, ambiguity resolution (LA) conditions the reliability of the entire subsequent operating phase.
  • LA ambiguity resolution
  • the known methods of removing ambiguity proceed by maximum likelihood among a set of potential solutions preselected or not. Without a priori knowledge, these potential solutions describe a large domain corresponding to a space whose dimension corresponds to the number of satellites processed by pair of antenna, and having for half-axes the basic length counted in number of wavelengths, in all directions of the search space.
  • the size of the research area determines both the duration of the ambiguity and above all the risk of error on the initial attitude.
  • a known solution uses the intermodulation of two or three GPS carriers (L1 x L2, L1x L5, ...) which has the effect of increasing the apparent wavelength.
  • the details sought may require the use of a base between pairs of very long antennas, even greater than the wavelength.
  • an inertial unit (1RS class 2Nmi / h) can provide the initial attitude (heading and vertical),
  • - accelerometers or level sensors can provide the initial pitch and / or roll (ICS class 1 mg),
  • a magnetometer (or a flow valve) can provide an initial heading, ...
  • the assistance methods to reduce the search area for ambiguity removal no longer apply if the course is not available, the assistance then being incomplete or partial (search for an optimized navigation system in terms of cost, without IRS-2Nmi / h, nor magnetometer).
  • the integrity performance of the LA can be very insufficient: of the order of 1% when the need is 10 "7.
  • current method does not allow to know the quality of the ambiguity resolution obtained, which would be essential before selecting the GNSS signals as source of attitude and heading of the carrier.
  • the subject of the present invention is a method for improving the determination of the attitude of a vehicle using radio navigation signals, a method which simultaneously makes this determination with the best possible resolution, without being sensitive neither to interference nor to multi-routes.
  • the method according to the invention is characterized in that it consists in producing a carrier phase loop (PLL) by summing the two signals from '' a pair of antennas distant by half a wavelength, these signals being made coherent before summation, to carry out a measurement of phase difference on the signals coming from a pair of distant antennas of several lengths d wave, passing through the frequency domain, then to perform a reduction in the area of research of the initial ambiguity removal for determining the attitude of a vehicle by interferometric GPS measurement and implementing a statistical test for the selection of ambiguity
  • PLL carrier phase loop
  • FIG. 1 is a diagram illustrating the path difference of a GPS signal, seen by a pair of antennas,
  • FIG. 2 is a block diagram of a GPS reception channel used by the method of the invention
  • FIG. 3 is a diagram explaining the step of determining attitude according to the present invention.
  • FIG. 4 is a diagram illustrating the processing of robust PLL measurement in accordance with the invention
  • FIG. 5 is a diagram illustrating the processing of measurement of phase difference with frequency excision, in accordance with the invention
  • FIG. 6 is a diagram illustrating the processing for precise measurement of phase difference with a robust PLL, in accordance with the invention.
  • FIG. 7 is a diagram illustrating the first step in defining a research area to resolve ambiguity without heading assistance, in accordance with the invention.
  • the present invention is described below with reference to the determination of the attitude of an aircraft using radio navigation signals emitted by a constellation of GPS satellites, but it is understood that it is not limited to this single application, and that it can be implemented for other kinds of vehicles, and that it can have recourse to other radio navigation signals.
  • the path difference of the GPS signal seen by a pair of antenna is representative of the angle ⁇ between the antenna base b and the direction (LOS) of the satellite, as shown in Figure 1. This path difference is measured by the differential phase between the GPS antennas referenced 1 (master antenna) and 2 (slave antenna).
  • the GPS delivers two ambiguous phase measurements (fractional part compared at wavelength) respectively for each pair of antenna and for the satellite in view considered.
  • the various successive wave planes of the signal of the satellite in question have been shown in broken lines.
  • the path difference is defined geometrically by the projection, (or the scalar product) of the base antenna length on the direction of the satellite considered:
  • a receiver receives and processes the GPS signal using a code tracking loop (DLL) and a carrier tracking loop (PLL) as shown in Figure 2, which relates to a reception channel 3.
  • DLL code tracking loop
  • PLL carrier tracking loop
  • circuit 4 comprising a demodulation filter, with automatic gain control and analog-digital conversion of the received signal.
  • Circuit 4 attacks two demodulators 5 and 6.
  • a local digital oscillator 7 directly attacks demodulator 5, and via a 90 ° phase shifter 8 demodulator 6.
  • a reference clock 9 is connected to circuits 4 and 7.
  • Demodulators 5 and 6 attack correlators 11 and Q1 referenced 10 and 11, and correlators 12 and Q2 referenced 12 and 13.
  • Correlators 10 and 11 supply S1_0 to a circuit corrector DLL and to a circuit 15 corrector PLL with Arctg function, while the correlators 12 and 13 supply S1_ ⁇ to circuit 14.
  • Circuit 14 is connected to a local code generator 16 which is connected on the one hand via an advance-delay circuit 17 to correlators 12 and 13, on the other hand directly to correlators 10 and 11, and which provides the pseudo-distance information.
  • the corrector 15 is connected to the oscillator 7 and provides the pseudo-speed information. Ambiguous phase information is available at the output of oscillator 7. In the case of two antennas, the signals received demodulated by the code and the carrier are written (neglecting the delay on the data):
  • - s1 (t) and s2 (t) are the satellite signals received by each antenna respectively
  • n1 (t) and n2 (t) are the noises received by each antenna respectively - f is the residual carrier frequency
  • the determination of the attitude is obtained by calculation from the measurements of ambiguous phases, by adding the ambiguity previously initialized and maintained, then by reversing the projection of the baseline on the directions of the different satellites, according to the synoptic of the figure 3.
  • the satellite LOS are taken from the “Position-speed-time resolved” (PVT) set resolved by GPS positioning, and from the “Satellite position-calculated speed” (PVS) set calculated by the receiver at from the navigation data received (ephemeris or almanac) and the time resolved.
  • the basic vector (18) is determined according to the knowledge of the direction of the satellite (LOS) and of the ambiguous phases, after removal of ambiguity (19) operated (20) on the phases ambiguous and performed when the device was initialized, then maintained in real time, knowing the phase velocities. From the determination of the base vector, the attitude of the vehicle (20) and therefore the GPS attitude parameters of this vehicle are extracted.
  • the initial LA consists in selecting for each satellite the whole number of wavelengths corresponding to the whole part of the path difference of the GPS signal.
  • ⁇ carry out the carrier phase loop (PLL) by summing the two respective signals from a pair of antennas close, half a wavelength apart to reinforce interference resistance and in operation with a low signal-to-noise ratio.
  • the summation is implemented using a signal processing method making it possible to make the two signals coherent before summation (by normalized intercorrelation) independently of any bias of the RF chains.
  • This alignment gives the PLL a spatial rejection of multi-routes which are not in the preferred direction (moreover, and in a manner already known, this pair of close antennas can also be useful for a first direct and unambiguous initialization of attitude)
  • ⁇ perform the phase difference measurement using a new signal processing technique on two signals respectively from a pair of antennas sufficiently spaced apart (several wavelengths, to achieve good angular precision for attitude measurement).
  • the measurement of phase difference is obtained by passing through the frequency domain, which makes it possible to reject signals which are not coherent in frequency (multipath with doppler shifted by the preferred frequency, etc.).
  • the invention consists in slaving the two respective PLLs of the pair of antennas by the same command resulting from a summation of the two detected signals. This summation increases the signal-to-noise ratio seen by the PLL. This summation is previously made consistent by canceling the delay of one signal with respect to the other. This delay is first measured by standardized complex intercorrelation.
  • a corollary advantage to this consistency is the partial spatial rejection of multi-path signals whose directions seen by the two antennas are not identical to each other.
  • the invention consists in: - calculating the phase difference by interferometry (normalized intercorrelation 21) between homologous signals, that is to say coming from the same satellite (S1_0 and S1 b_0), issued respectively by the two antennas (22, 23) of the same pair of antennas, these antennas being less than ⁇ / 2 apart, then filtered, demodulated and despread (24, 25),
  • Restoring the coherence of the demodulated signals on the basis of the standardized intercorrelation output (representing a complex argument signal equal to the instantaneous phase difference, the noise remaining negligible) has the consequence of reducing the effect of the decorrelated noises. spatially, close to a “patch” type antenna with beam formation.
  • this summing makes it possible to gain 3 dB on the signal to noise ratio in the loop (therefore on the measurement accuracy) in the case of independent Gaussian noises, and thus to operate in the presence of even higher disturbances.
  • any interference signal not arriving in the direction of the signal will be attenuated, but with an efficiency dependent on the bandwidth of this interference with respect to the separation distance of the two antennas:
  • the antennas are brought closer to one another better than ⁇ / 2, and in this case the attenuation of a coherent wave is carried out without any other angular ambiguity on any type of interference (narrow band and broad band), Either the antennas are separated by several wavelengths, and in this case, the directivity becomes illusory on narrow band signals, due to spatial ambiguities (unambiguous spatial directivity can only be achieved in the case of noise of very broadband, not compatible with conventional receiver bandwidth). This provision is not retained.
  • the invention proposes to carry out an intercorrelation passing through the frequency domain.
  • the phase difference measurement is associated with processing in the frequency domain by means of a direct fast Fourier transformation (DFT) of the two signals received on the two antennas (antennas 22 and 23A - see FIG. 5). , these two signals being from the same satellite.
  • DFT direct fast Fourier transformation
  • This transformation is followed by a complex product in the frequency domain, a frequency excision processing to reject the signals of abnormal frequency (non-coherent doppler multipath ...) then an inverse DFT considered at the offset of zero time (the phase difference between the two carrier signals being measured by the value of the intercorrelation for a zero time offset).
  • This processing is carried out at the rate of the correlation of the code (typically of the order of 1 kHz) after Fourier transformation by FFT of the correlation “I” and “Q” outputs.
  • a method for carrying it out is for example:
  • the advantage of this method is to allow the calculation of phase differences in a relatively wide band compatible with significant dynamic changes (without having to integrate excessively), by benefiting from an elimination of the noises which may appear in the band.
  • the invention consists in:
  • noise suppression 34, 35
  • spectral normalization carried out by means of a sliding window serving as a noise calculation reference
  • FIG. 6, combining FIGS. 4 and 5, presents the complete processing in the receiver, for the measurement of precise phase difference with robust PLL.
  • the same elements as those of Figures 4 and 5 are assigned the same reference numerals.
  • the invention consists, for the preselection of the field of potential ambiguities, in:
  • the search area is reduced, even in the event of partial assistance, thus optimizing the initialization time and the robustness of the attitude determination by GPS.
  • Figure 7 illustrates the first step in defining the domain:
  • n M [N E ( ⁇ • b A , DE (Pitch A
  • ⁇ MAX E (û " • b A i D E ( Tan g a g e A i D E.Azimut SATELLITE ))
  • ⁇ MIN ⁇ MIN _ dn
  • MAX ⁇ MAX + n
  • the present invention provides a solution to the problem of transfer of ambiguity to the long base in the presence of propagation bias.
  • This process makes it possible to compensate for differences in the path of the wave between antennas with sufficient precision, to improve the signal-to-noise ratio towards the satellites, and (as natural fallout) reduce the effect of correlated sources of interference in the 'space.
  • An important advantage of this method based on the use of the three cited antennas is to allow the direct resolution of ambiguities on simple phase difference measurements independently of the phase biases of the RF channels.
  • the present invention provides an improvement in LA, thanks to a new statistical evaluation criterion. This new method also develops a criterion for evaluating resolved ambiguity.
  • phase residues are normalized by their estimated standard deviation
  • phase measurement for more robustness of the phase measurement, it is possible to use simultaneously: close antennas (or elements of a multi-antenna type CRPA) allowing better elimination of coherent noise sources spatially, and usable to initialize the removal of ambiguity, distant antennas making it possible to carry out precise measurements of phase difference, initialized from closely spaced antennas, or any other inertial attitude calculation means.
  • close antennas or elements of a multi-antenna type CRPA
  • the test can be applied to several recurrences of attitude measurements (which reduces the renewal rate).

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

Le procédé conforme à l'invention, mis en oeuvre à l'aide d'au moins une paire d'antennes, consiste à effectuer : une analyse spectrale normalisée avant intercorrélation des signaux démodulés I, Q, par le code et la porteuse, une remise en cohérence de phase des signaux avant recombinaison, une mesure d'écart de phase réalisée par une vraie intercorrélation des voies satellites homologues par interférométrie des deux signaux issus d'un même satellite et reçus respectivement par une paire d'antennes, puis à effectuer une réduction du domaine de recherche du lever d'ambiguïté initial pour la détermination d'attitude d'un véhicule par mesure GPS interférométrique et mettre en oeuvre un test statistique pour la sélection de l'ambiguïté.

Description

PROCEDE D'AMELIORATION DE LA DETERMINATION DE L'ATTITUDE D'UN VEHICULE A L'AIDE DE SIGNAUX DE RADIONAVIGATION PAR
SATELLITE
La présente invention se rapporte à un procédé d'amélioration de la détermination de l'attitude d'un véhicule à l'aide de signaux de radionavigation par satellite.
La détermination d'attitude à l'aide de signaux de radionavigation, en particulier de signaux GPS, consiste à comparer les différences de distances satellites-porteur par rapport à une base d'antennes connue. Une très grande résolution de mesure est donc nécessaire. Pour atteindre cette résolution, les mesures de phase du signal GPS sont utilisées.
Cependant, la détermination d'attitude par GPS multi-antenne ainsi réalisée nécessite des améliorations pour atteindre le haut niveau d'intégrité exigé par les applications « safety-of-life », c'est-à-dire à exigences de sécurité vitales
En effet contrairement au positionnement GPS qui exploite les pseudo-distances délivrées par la boucle de poursuite du code (Delay Locked loop, DLL) à partir d'une seule antenne, la détermination d'attitude exploite les mesures de phases délivrées par la boucle de poursuite de porteuse (Phase Locked Loop, PLL) et en mode différentiel entre paire d'antennes, cela se traduit par un besoin d'amélioration des aspects suivants : - amélioration du fonctionnement de la PLL en présence de perturbations (interférences, multitrajets de structure ou autres, sauts de phase PLL...)
- amélioration de l'initialisation de la détermination d'attitude (lever d'ambiguïté), en particulier pour la réduction du domaine de recherche, même en ayant recours à une aide partielle
(sans aide en cap) afin de réduire le temps d'initialisation et surtout d'atteindre la fiabilité nécessaire en réduisant le risque d'erreur,
- amélioration de l'intégrité de l'attitude GPS (détection des erreurs).
L'invention a donc pour objet l'amélioration de la détermination de l'attitude d'un véhicule à partir d'un signal de radionavigation, en particulier d'un signal GPS. L'attitude d'un porteur inclut ici les angles de cap, tangage et roulis dans un repère terrestre local.
Dans le cas de la détermination d'attitude, le signal GPS est reçu simultanément par une paire (ou des paires) d'antennes, permettant d'établir l'écart des phases des signaux émis par les satellites GPS et mesurées respectivement par chacune des deux antennes de la paire considérée.
Dans le cas plus classique d'un simple positionnement, hors mesure d'attitude, un récepteur GPS poursuit le signal GPS en code à l'aide d'une DLL (mesurant ainsi la pseudo-distance suivant la position du code reçu) et en porteuse à l'aide d'une PLL (mesurant la pseudo-vitesse suivant la vitesse de phase reçue). Un canal de réception est dédié à chaque satellite poursuivi. La DLL est donc primordiale pour la fonction localisation par rapport à la PLL. Celle-ci peut d'ailleurs être momentanément remplacée par une aide en vitesse (mode dit « code only »). Mais pour déterminer une attitude, si on voulait faire la différence entre pseudo-distances des DLL entre deux antennes, la résolution obtenue ne serait pas suffisante.
Par contre, la comparaison des phases reçues issues des PLL permet d'atteindre la résolution nécessaire. Cependant, cette mesure différentielle pose des problèmes spécifiques par rapport à l'utilisation classique du positionnement GPS , problèmes qui sont :
- PLL moins robuste aux interférences que la DLL,
- Mesure de phase modulo la longueur d'onde, d'où ce qui est communément appelé « l'ambiguïté » de la détermination d'attitude par GPS,
- Mesure de phase sensible aux transitoires PLL : bruit, saut de phase,
- Mesure différentielle entre antennes induisant une sensibilité particulière au multitrajet local (multitrajet de structure,...). II apparaît nécessaire, pour satisfaire le besoin de mesure d'attitude, en particulier dans le domaine aéronautique, avec les exigences de sécurité requises, d'améliorer la robustesse et la précision des mesures de phase du signal GPS, et d'améliorer et de fiabiliser le lever d'ambiguïté.
Pour améliorer la robustesse et la précision des mesures de phase, lors de la réception du signal GPS, on connaît des méthodes de réjection des sources d'interférences par traitement spectral ou spatial appliqué au signal d'entrée (avant démodulation par le code ou la porteuse). De plus le remplacement temporaire de la PLL par une aide en vitesse (mode « code only ») permet de maintenir la DLL pour le positionnement en cas de forte interférence, mais la PLL n'étant plus opérationnelle, les mesures de phase ne sont plus disponibles pour la mesure d'attitude.
En ce qui concerne la robustesse des mesures GPS vis-à-vis des trajets multiples, on connaît des méthodes au niveau de la DLL ; corrélateur étroit, ou méthodes inverses s'appuyant sur la connaissance a priori de la forme de la fonction d'intercorrélation du code local et du code reçu.
Les solutions actuelles de traitement des interférences assurent principalement le bon fonctionnement de la partie DLL et du positionnement GPS (traitements avant démodulation, mode dit « code only »).
La mesure de phase et la détermination d'attitude ne sont donc pas toujours disponibles en cas de forte interférence même si le positionnement GPS peut être maintenu.
De même, en ce qui concerne la robustesse des mesures GPS vis-à-vis des trajets multiples, les méthodes connues ne permettent pas de rejeter tous les trajets multiples pouvant dégrader spécifiquement la détermination d'attitude, en particulier pour des trajets multiples à faible retard non détectables par les méthodes connues au niveau de la DLL (résolution de la DLL insuffisante pour le besoin de mesure d'attitude).
Initialisant la détermination d'attitude, le lever d'ambiguïté (LA) conditionne la fiabilité de toute la phase d'exploitation ultérieure. Les procédés connus de lever d'ambiguïté procèdent par maximum de vraisemblance parmi un ensemble de solutions potentielles présélectionnées ou non. Sans connaissance a priori, ces solutions potentielles décrivent un large domaine correspondant à un espace dont la dimension correspond au nombre de satellites traités par paire d'antenne, et ayant pour demi-axes la longueur de base comptée en nombre de longueurs d'onde, dans toutes les directions de l'espace de recherche.
La taille du domaine de recherche détermine à la fois la durée du lever d'ambiguïté et surtout le risque d'erreur sur l'attitude initiale. Pour réduire le domaine de recherche, une solution connue utilise l'intermodulation de deux ou trois porteuses GPS (L1 x L2, L1x L5,...) ce qui a pour effet d'augmenter la longueur d'onde apparente. Cependant les précisions recherchées peuvent nécessiter l'utilisation d'une base entre paire d'antennes de grande longueur, encore supérieure à la longueur d'onde.
Des aides sont alors communément utilisées pour fournir des attitudes initiales grossières et réduire ainsi l'espace de recherche. Par exemple :
- une centrale inertielle (1RS de la classe 2Nmi/h) peut fournir l'attitude initiale (cap et verticale),
- des accéléromètres ou des capteurs de niveau peuvent fournir le tangage et/ou le roulis initiaux (ICS de la classe 1 mg),
- un magnétomètre (ou une vanne de flux) peut fournir un cap initial,...
Les méthodes d'aide pour réduire le domaine de recherche du lever d'ambiguïté ne s'appliquent plus si le cap n'est pas disponible, l'aide étant alors incomplète ou partielle (recherche d'un système de navigation optimisé en terme de coût, sans IRS-2Nmi/h, ni magnétomètre).
En l'absence de réduction du domaine de recherche, la performance d'intégrité du LA peut être très insuffisante : de l'ordre de 1 % lorsque le besoin est de 10"7. En plus du besoin de lever des ambiguïtés potentielles, la méthode actuelle ne permet pas de connaître la qualité du lever d'ambiguïté obtenu, ce qui serait pourtant essentiel avant de sélectionner les signaux GNSS comme source d'attitude et de cap du porteur.
La présente invention a pour objet un procédé d'amélioration de la détermination de l'attitude d'un véhicule à l'aide de signaux de radionavigation, procédé qui permet à la fois d'effectuer cette détermination avec la meilleure résolution possible, sans être sensible ni aux interférences ni aux multitrajets.
Le procédé conforme à l'invention, mis en œuvre à l'aide d'au moins une paire d'antennes, est caractérisé en ce qu'il consiste à réaliser une boucle de phase porteuse (PLL) en sommant les deux signaux issus d'une paire d'antennes distantes d'une demi-longueur d'onde, ces signaux étant rendus cohérents avant sommation, à réaliser une mesure d'écart de phase sur les signaux issus d'une paire d'antennes distantes de plusieurs longueurs d'onde, en passant par le domaine fréquentiel, puis à effectuer une réduction du domaine de recherche du lever d'ambiguïté initial pour la détermination d'attitude d'un véhicule par mesure GPS interférométrique et mettre en œuvre un test statistique pour la sélection de l'ambiguïté
La présente invention sera mieux comprise à la lecture de la description d'un mode de mise en œuvre, pris à titre d'exemple non limitatif et illustré par le dessin annexé, sur lequel :
- la figure 1 est un schéma illustrant la différence de marche d'un signal GPS, vue par une paire d'antennes,
- la figure 2 est un bloc-diagramme d'un canal de réception GPS utilisé par le procédé de l'invention,
- la figure 3 est un diagramme expliquant l'étape de détermination d'attitude selon là présente invention,
- la figure 4 est un diagramme illustrant le traitement de mesure PLL robuste conformément à l'invention, - la figure 5 est un diagramme illustrant le traitement de mesure d'écart de phase avec excision fréquentielle, conformément à l'invention,
- la figure 6 est un diagramme illustrant le traitement de mesure précise d'écart de phase avec PLL robuste, conformément à l'invention, et
- la figure 7 est un diagramme illustrant la première étape de définition de domaine de recherche de lever d'ambiguïté sans aide en cap, conformément à l'invention.
La présente invention est décrite ci-dessous en référence à la détermination de l'attitude d'un aéronef à l'aide de signaux de radionavigation émis par une constellation de satellites GPS, mais il est bien entendu qu'elle n'est pas limitée à cette seule application, et qu'elle peut être mise en œuvre pour d'autres sortes de véhicules, et qu'elle peut avoir recours à d'autres signaux de radionavigation. La différence de marche du signal GPS vue par une paire d'antenne est représentative de l'angle θ entre la base d'antenne b et la direction (LOS) du satellite, comme le montre la figure 1. Cette différence de marche est mesurée par la phase différentielle entre les antennes GPS référencées 1 (antenne maître) et 2 (antenne esclave). Pour cela, le GPS délivre deux mesures de phases ambiguës (partie fractionnaire par rapport à la longueur d'onde) respectivement pour chaque paire d'antenne et pour le satellite en vue considéré. On a représenté en traits interrompus les différents plans d'onde successifs du signal du satellite considéré.
La différence de marche est définie géométriquement par la projection, (ou le produit scalaire) de la longueur de base d'antenne sur la direction du satellite considéré :
ΔΦ = (n + Fract)* λ = û « b, avec,
ΔΦ, différence de marche , n, partie entière (ambigϋité),
Fract, partie fractionnaire, délivrée par le calculateur GPS λ, longueur d'onde
0, direction du satellite , b, longueur de base d'antenne
Pour délivrer une mesure de phase, un récepteur reçoit et traite le signal GPS à l'aide d'une boucle de poursuite en code (DLL) et d'une boucle de poursuite en porteuse (PLL) comme le montre la figure 2, qui se rapporte à un canal de réception 3.
Une des antennes, par exemple l'antenne 1 , est suivie d'un circuit 4 comportant un filtre de démodulation, avec commande automatique de gain et conversion analogique-numérique du signal reçu. Le circuit 4 attaque deux démodulateurs 5 et 6. Un oscillateur local numérique 7 attaque directement le démodulateur 5, et via un déphaseur 8 de 90° le démodulateur 6. Une horloge de référence 9 est reliée aux circuits 4 et 7. Les démodulateurs 5 et 6 attaquent des corrélateurs 11 et Q1 référencés 10 et 11 , et des corrélateurs 12 et Q2 référencés 12 et 13. Les corrélateurs 10 et 11 fournissent S1_0 à un circuit 14 correcteur DLL et à un circuit 15 correcteur PLL à fonction Arctg, tandis que les corrélateurs 12 et 13 fournissent S1_τ au circuit 14. Le circuit 14 est relié à un générateur de code local 16 qui est relié d'une part via un circuit 17 d'avance-retard aux corrélateurs 12 et 13, d'autre part directement aux corrélateurs 10 et 11 , et qui fournit l'information de pseudo-distance. Le correcteur 15 est relié à l'oscillateur 7 et fournit l'information de pseudo-vitesse. L'information de phase ambiguë est disponible à la sortie de l'oscillateur 7. Dans le cas de deux antennes, les signaux reçus démodulés par le code et la porteuse s'écrivent (en négligeant le retard sur les données) :
r, (t) = s-, (t) + n-, (t) = a.D(t). exp j(2πft + φ) + n-, (t) r2(t) = s2(t) + n2(t) = a.D(t).expj(2πft + 2πf0τ + φ) + n2(t) OÙ,
- s1 (t) et s2(t) sont respectivement les signaux satellites reçus par chaque antenne
- n1 (t) et n2(t) sont respectivement les bruits reçus par chaque antenne - f est la fréquence porteuse résiduelle,
- fO est la fréquence porteuse reçue
- φ est la phase initiale du signal (indéterminée)
- D(t) est le signe de la donnée (qui module éventuellement la porteuse) - a est l'amplitude du signal traité
- τ est le retard de propagation de l'onde entre les deux antennes (= — — ) , séparées de d et d'incidence θ. c
La détermination de l'attitude est obtenue par calcul à partir des mesures de phases ambiguës, en additionnant l'ambiguïté préalablement initialisée et entretenue, puis en inversant la projection de la ligne de base sur les directions des différents satellites, suivant le synoptique de la figure 3. Les LOS des satellites sont issues de l'ensemble « Position-vitesse-temps résolus » (PVT) résolu par le positionnement GPS, et de l'ensemble « Position-vitesse satellite calculées » (PVS) calculé par le récepteur à partir des données de navigation reçues (ephémérides ou almanach) et du temps résolu.
Sur le diagramme de la figure 3, le vecteur de base (18) est déterminé d'après la connaissance de la direction du satellite (LOS) et des phases ambiguës, après lever d'ambiguïté (19) opéré (20) sur les phases ambiguës et réalisé à l'initialisation du dispositif, puis entretenu en temps réel, en connaissant les vitesses de phase. A partir de la détermination du vecteur de base, on extrait l'attitude du véhicule (20) et donc les paramètres d'attitude GPS de ce véhicule. Le LA initial consiste à sélectionner pour chaque satellite le nombre entier de longueurs d'onde correspondant à la partie entière de la différence de marche du signal GPS.
Du fait de l'existence de biais de phase différents dans chaque canal RF du récepteur, on réalise pratiquement la double différence des écarts de phase pris sur des axes à vue différents ; ce procédé utilisé classiquement pour la résolution des biais différentiels entre récepteurs GPS permet ici la résolution d'attitude.
Cette sélection est prise dans un domaine de recherche incluant au départ toutes les valeurs possible d'ambiguïté. En absence de connaissance à priori, le domaine est borné pour chaque satellite par le nombre entier de longueur d'onde contenu dans la ligne de base.
L'utilisation d'une aide consiste à réduire cet intervalle de recherche. Pour cela la prédiction de cette valeur est donnée par l'utilisation de la formule suivante (cas d'une longueur de base définie par deux antennes, et ses angles géographiques tangage et cap) :
ΠDOMAINE = E(ϋ • bA|DE ) ± dn avec
Figure imgf000010_0001
ou,
Cos(tangageA,DE ).Cos(CapA,DE ) Cos(tangageA|DE ).Sin(CapA]DE ) BODY Sin(tangageA,DE)
Pi , paramètres d'aide, respectivement de tangage, de cap dPi , incertitude à priori sur la connaissance du paramètre βoDY étant la longueur de base du véhicule considéré Pour améliorer la robustesse et la précision des mesures d'écart de phase entre antennes d'une paire d'antennes, utiles à la détermination d'attitude, l'invention propose deux nouveaux traitements :
→ réaliser la boucle de phase porteuse (PLL) en sommant les deux signaux respectifs issus d'une paire d'antennes proches, distantes d'une demi-longueur d'onde pour renforcer la tenue aux interférences et en fonctionnement à faible rapport signal à bruit. La sommation est implémentée à l'aide d'une méthode de traitement de signal permettant de rendre cohérents les deux signaux avant sommation (par intercorrélation normalisée) indépendamment de tout biais des chaînes RF. Cette mise en cohérence confère à la PLL une réjection spatiale des multitrajets qui ne sont pas dans la direction privilégiée (par ailleurs, et de façon déjà connue, cette paire d'antennes rapprochées peut également être utile pour une première initialisation directe et non ambiguë de l'attitude),
→ réaliser la mesure d'écart de phase à l'aide d'une nouvelle technique de traitement de signal sur deux signaux issus respectivement d'une paire d'antennes suffisamment écartées (plusieurs longueurs d'onde, pour accéder à une bonne précision angulaire pour la mesure d'attitude). La mesure d'écart de phase est obtenue en passant par le domaine fréquentiel, ce qui permet de rejeter des signaux non cohérents en fréquence (multitrajet à doppler décalé de la fréquence privilégiée,...). Pour effectuer la sommation sur deux antennes rapprochées, l'invention consiste à asservir les deux PLL respectives de la paire d'antennes par la même commande issue d'une sommation des deux signaux détectés. Cette sommation permet d'augmenter le rapport signal à bruit vu par la PLL. Cette sommation est préalablement rendue cohérente en annulant le retard d'un signal par rapport à l'autre. Ce retard est d'abord mesuré par intercorrélation complexe normalisée.
Un avantage corollaire à cette mise en cohérence est le rejet spatial partiel des signaux de multitrajets dont les directions vues par les deux antennes ne sont pas identiques entre elles.
Dans le détail (voir figure 4), l'invention consiste à : - calculer l'écart de phase par interférométrie (intercorrélation normalisée 21 ) entre signaux homologues, c'est-à-dire issus d'un même satellite (S1_0 et S1 b_0), délivrés respectivement par les deux antennes (22, 23) d'une même paire d'antennes, ces antennes étant distantes de moins de λ/2, puis filtrés, démodulés et desétalés (24, 25),
- réaliser une remise en cohérence par multiplication (26) des signaux (S1_0 et S1 b_0) par la moitié du retard relatif (27) évalué par l'intercorrélation normalisée précédente. Les deux signaux sont préalablement retardés (28, 29) de façon à compenser le retard introduit par l'intercorrélation ( 20ms dans le présent exemple), - asservir (30) la boucle de phase PLL à partir de l'addition vectorielle (31 ) des deux signaux d'antennes remis en cohérence. Les PLL sont alors asservies sur la phase moyenne entre les deux signaux respectivement reçus par les deux antennes, avec un gain de 3dB sur le rapport signal à bruit.
La remise en cohérence des signaux démodulés sur la base de la sortie d'intercorrélation normalisée (représentant un signal complexe d'argument égal à l'écart de phase instantané, le bruit restant négligeable) a pour conséquence de réduire l'effet des bruits décorrélés spatialement, de manière proche d'une antenne de type « deux patches » à formation de faisceaux.
Après remise en cohérence, les signaux étant sommés avant d'asservir la PLL, cette sommation permet de gagner 3 dB sur le rapport signal à bruit dans la boucle (donc sur la précision de mesure) dans le cas de bruits indépendants gaussiens, et ainsi de fonctionner en présence de perturbations d'autant plus élevées.
De plus, tout signal d'interférence n'arrivant pas dans la direction du signal sera atténué, mais avec une efficacité dépendant de la largeur de bande de cette interférence vis à vis de la distance de séparation des deux antennes :
Soit, selon l'invention, les antennes sont rapprochées à mieux que λ/2, et dans ce cas l'atténuation d'une onde cohérente est réalisée sans autre ambiguïté angulaire sur tout type d'interférence (bande étroite et bande large), Soit les antennes sont écartées de plusieurs longueurs d'onde, et dans ce cas, la directivité devient illusoire sur les signaux à bande étroite, du fait des ambiguïtés spatiales (une directivité spatiale non ambiguë ne pouvant être atteinte que dans le cas de bruit de très large bande, non compatible avec la bande passante classique des récepteurs). Cette disposition n'est pas retenue.
Pour améliorer la robustesse des mesures d'écart de phase, l'invention propose de réaliser une intercorrélation en passant par le domaine fréquentiel. Ainsi on associe à la mesure d'écart de phase un traitement dans le domaine fréquentiel par l'intermédiaire d'une transformation de Fourier rapide directe (DFT) des deux signaux reçus sur les deux antennes (antennes 22 et 23A -voir figure 5), ces deux signaux étant issus d'un même satellite. Cette transformation étant suivie d'un produit complexe dans le domaine fréquentiel, d'un traitement d'excision de fréquence pour rejeter les signaux de fréquence anormale (multitrajet à doppler non cohérent...) puis d'une DFT inverse considérée au décalage de temps nul (l'écart de phase entre les deux signaux porteuses étant mesuré par la valeur de l'intercorrélation pour un décalage de temps nul).
Cette solution permet de "nettoyer" le spectre des signaux de porteuse démodulés, des raies ou bandes parasites ayant pu se superposer au signal de porteuse (interférences à bande étroite, trajets multiples décalés en Doppler...)
Ce traitement est réalisé au rythme de la corrélation du code (typiquement de l'ordre de 1 kHz) après transformation de Fourier par FFT des sorties « I » et « Q » de corrélation. Un procédé pour le réaliser est par exemple:
- une normalisation du spectre (centrage et normalisation de l'écart-type du bruit de fond, par fenêtrage glissant par exemple, permettant d'estimer le bruit de fond en dehors des niveaux fréqϋentiels dus au signal utile)
- une détection des pointes de bruits, et identification de la porteuse (a priori, de rapport signal à bruit le plus important et le plus stable ; typiquement 10 dB dans 1 KHZ en GPS)
- la suppression des pointes de bruit dans le domaine fréquentiel des deux signaux (excision fréquentielle par mise à zéro ou au niveau du bruit de fond des amplitudes spectrales aux fréquences détectées),
- la multiplication complexe des spectres conjugués filtrés (intercorrélation) - le retour en temporel par DFT inverse considérée à zéro, pour obtenir l'intercorrélation entre les 2 signaux et donc leur écart de phase relatif.
L'avantage de cette méthode est de permettre le calcul des écarts de phase dans une bande relativement large compatible avec des évolutions de dynamique importantes (sans devoir intégrer outre mesure), en bénéficiant d'une élimination des bruits pouvant apparaître dans la bande.
Dans le détail, comme représenté en figure 5, l'invention consiste à :
- réaliser les transformées de Fourier (32, 33), par exemple par FFT-128 points, des deux signaux (échantillonnés à environ
1 KHz)
- réaliser une suppression des bruits (34, 35) de chacun des deux spectres complexes obtenus, et incluant une normalisation spectrale menée au moyen d'une fenêtre glissante servant de référence de calcul de bruit (typiquement
2x10 points dans le cas d'une FFT à 128 points)
- une détection des bosses de bruit (36, 37) sur la base du spectre normalisé, et la détermination de la raie porteuse et de son pistage (sur un critère de niveau et stabilité du rapport signal/bruit et de sa position en fréquence), et élimination des bosses identifiées à du bruit ou des trajets multiples, de manière conjointe sur les deux spectres
- réaliser une intercorrélation par produit complexe (38) des spectres (le deuxième signal étant conjugué : 39) - réaliser le calcul de son argument à décalage nul (pour l'abscisse zéro) par une DFT inverse complexe (40), pour exprimer la phase différentielle (41 ) entre les deux signaux de départ. La date de validité de cette phase mesurée est synchrone du milieu de l'intervalle de temps d'analyse (20 ms par exemple). Pour un rapport signal à bruit typique de 40dBHz, la précision ainsi obtenue est de l'ordre de 0,5° à 5Hz (5° à 50Hz).
La figure 6, combinant les figures 4 et 5, présente le traitement complet dans le récepteur, pour la mesure d'écart de phase précise avec PLL robuste. Sur cette figure, les mêmes éléments que ceux des figures 4 et 5 sont affectés des mêmes références numériques.
On va maintenant décrire la partie de procédé de l'invention se rapportant au lever d'ambiguïté avec aide partielle.
Dans le cas d'une aide partielle (pas d'aide en cap) l'invention consiste, pour la présélection du domaine des ambiguïtés potentielles à:
- Définir un premier domaine de recherche d'ambiguïté par deux hypothèses particulières sur le cap d'aide, et en utilisant les autres aides connues.
- Puis élargir ce domaine avec les incertitudes sur les autres paramètres d'aide.
Ainsi, le domaine de recherche est réduit, même en cas d'aide partielle, optimisant ainsi le temps d'initialisation et la robustesse de la détermination d'attitude par GPS.
Ces hypothèses particulières de cap a priori sont les suivantes et correspondent à :
- cap d'aide égal à l'azimut du satellite (maximum du produit scalaire de projection de la base d'antenne sur cette direction)
- cap d'aide égal à l'azimut du satellite + 180° (minimum de projection). La figure 7 illustre la première étape de définition du domaine :
Sur cette figure, on a repéré les vecteurs de tangage d'aide et d'azimut du satellite considéré, ainsi que les vecteurs b et u définis précédemment.
On a alors : nM[N = E(ϋ • bA,DE (Tangage A|DE,AzimutSATELLlTE + π))
ΠMAX = E(û" • bAiDE(Tangage AiDE.AzimutSATELLITE)) L'élargissement du domaine est ensuite donné par : ΠMIN = ΠMIN _dn ΠMAX = ΠMAX + n avec
Figure imgf000016_0001
et Pi paramètres d'aides restant disponibles (tangage accélérométrique)
Dans le cas d'une base complémentaire courte, la présente invention apporte une solution au problème de transfert d'ambiguïté sur la base longue en présence de biais de propagation.
En fait, la démonstration suivante montre que l'utilisation de doubles différences au lieu de simples différences, appliquée à cette configuration de trois antennes (base longue+ base courte : voir en figure 6 les antennes 22, 23 et 23A), permet de lever directement l'ambiguïté sur la phase, ou de la réduire à quelques unités si le retard de propagation RF est plus important qu'une demi longueur d'onde. L'ambiguïté sur l'écart de phase sur la base longue étant résolue, il est alors possible de calculer l'attitude du porteur.
Ayant obtenu l'attitude du porteur, il est alors possible de reconstituer (par inversion de l'estimation) le biais de propagation RF sur les bases courtes et longues (par estimation de la différence entre les mesures d'écart de phase entre antennes et les mesures d'écart de phase théoriques reconstitués à partir de l'estimation de l'attitude).
Ce procédé permet de compenser les différences de marche de l'onde entre antennes avec suffisamment de précision, pour améliorer le rapport signal à bruit en direction des satellites, et (comme retombée naturelle) diminuer l'effet des sources d'interférences corrélées dans l'espace.
Le rapport entre les écarts de phases entre les antennes 1 et 2 (distantes de D12) et 1 et 3 (distantes de D13), s'écrit:
D ? • (1 ) Δφ12 = Δφ13 x— Lέ- u13 et, pour un satellite i: Δφ13i = 2njπ + δφ13i + φ13 ΔΦl2i = δφ12i + φ12
de même pour un satellite j: Δφi3j = 2njπ + δ(p13j + φ13
Figure imgf000017_0001
δφ13j : écart_de_phase_mesuré - sur__base_longue Δcp12j : écart_de_phase_mesuré - sur_base_courte φ13:biais_RF_base_longue φ12:biais_RF_base_courte nj;ambiguité_de_mesure_de_phase_base_longue
Δφ12i = δ(p12i + φ12 = Δφ13i x — ^- = (2niπ + δφ13i + φ13)x-^ u13 u
(1 )=> 13
D 12 D 12 _ D 12 δΦl2i + Φl2 ~ Φl3 x D = l 3j x ^ = (2niπ + δΦl3i)x
13 D 13 D 13
de même, pour un satellite j: δ<Pl2i + Φl2 ~ Φl3 x TT = (2njπ + δφl3i )x ^- U13 U13 d'où l'expression de la double différence:
Δ<Pi2i - Δφ12j
Figure imgf000017_0002
(n, - nj) = _ (Δφ12i - Δφ12j )x— ^- (δφ13i - δφ13j )
2π U12
Un avantage important de cette méthode basée sur l'utilisation des trois antennes citées est de permettre la résolution directe des ambiguïtés sur les mesures d'écart de phase simples indépendamment des biais de phase des chaînes RF.
Toute autre méthode basée sur deux antennes seulement nécessite de travailler sur les doubles différences de phases inter-satellites, qui réduit de 3 dB le rapport signal à bruit de mesure. Du fait des incertitudes de mesure sur la phase, la solution précédente est en fait encadrée par les intervalles d'incertitude de mesures
(grossis du fait de la projection des longueurs); cependant on pourra se limiter à seulement quelques ambiguïtés (+/-1 par exemple en bruit de phase nominal, environ 5mm).
Cette résolution directe des ambiguïtés, sans aides externes, permet le calcul de l'attitude, puis autorise la calibration des retards RF(φ13,
Une fois le domaine de recherche défini, il est nécessaire de disposer d'une méthode de LA fiable. La présente invention propose une amélioration du LA, grâce à un nouveau critère d'évaluation statistique. Cette nouvelle méthode élabore également un critère d'évaluation de l'ambiguïté résolue.
Le principe est le suivant : - les mesures d'écart de phase (non ambiguës) et l'attitude porteur sont reliées par transformation linéaire, H, fonction de la géométrie des axes à vue des satellites,
- le calcul de l'attitude sur la base de l'ensemble des combinaisons d'ambiguïté de phase conduit pour chaque hypothèse à un vecteur de résidus de mesures de phase
(écart entre l'écart de phase mesuré et l'écart de phase théorique reconstitué), ces résidus de phases sont normalisés par leur écart type estimé
(dépendant de l'imprécision des mesures d'écart de phase élémentaires et des écarts de phase reconstitués). Ils sont significatifs de la distance existant entre l'attitude estimée et
- l'attitude réelle (sous l'hypothèse qu'aucune mesure n'est erronée).
- sous l'hypothèse que les bruits de mesure de phase suivent une distribution Gaussienne, on peut alors appliquer un test de détection quadratique, selon une loi du Khi-2 (somme des carrés des écarts types pondérés sur les différents axes satellites, éventuellement intégré sur plusieurs récurrences temporelles) - les solutions concurrentes sont sélectionnées sur la base d'un seuillage fonction des probabilités de détection et de fausse alarme choisies, fonctions de la loi du Khi-2,
- dans le cas où plusieurs solutions restent en concurrence, un test complémentaire de résidus portant sur la longueur de base apparente de la solution d'attitude est réalisé,
- s'il reste malgré tout encore des solutions candidates (liées à la non observabilité des solutions), il convient de tester un autre sous-ensemble de satellites visibles (ce qui permet d'écarter une erreur d'un jeu de satellites).
Selon des variantes du procédé de l'invention, pour plus de robustesse de la mesure de phase, on peut utiliser simultanément : des antennes rapprochées (ou des éléments d'une multi-antenne type CRPA) permettant de mieux éliminer des sources de bruit cohérentes spatialement, et utilisable pour initialiser le lever d'ambiguïté, des antennes éloignées permettant de réaliser des mesures précises d'écart de phase, initialisées à partir des antennes rapprochées, ou de tout autre moyen inertiel de calcul d'attitude.
Pour plus de robustesse du test de sélection des ambiguïtés, on peut élargir le domaine de recherche de quelques ambiguïtés supplémentaires :
- une ou deux ambiguïtés pour compenser les troncatures dues aux calculs par parties entières
- une ambiguïté pour le bruit de mesure sur la partie fractionnaire de la phase interférométrique
Pour améliorer le contrôle d'intégrité, le test peut être appliqué sur plusieurs récurrences de mesures d'attitude ( ce qui réduit la cadence de renouvellement).

Claims

REVENDICATIONS
1. Procédé d'amélioration de la détermination de l'attitude d'un véhicule à l'aide de signaux de radionavigation, mis en œuvre à l'aide d'au moins une paire d'antennes (22, 23, 23A), caractérisé en ce qu'il consiste à réaliser une boucle de phase porteuse (PLL) en sommant les deux signaux issus d'une paire d'antennes distantes d'une demi-longueur d'onde, ces signaux étant rendus cohérents avant sommation, à réaliser une mesure d'écart de phase sur les signaux issus d'une paire d'antennes distantes de plusieurs longueurs d'onde, en passant par le domaine fréquentiel, puis à effectuer une réduction du domaine de recherche du lever d'ambiguïté initial pour la détermination d'attitude d'un véhicule par mesure GPS interférométrique et mettre en œuvre un test statistique pour la sélection de l'ambiguïté.
2. Procédé selon la revendication 1 , caractérisé en ce que pour réaliser une boucle de phase porteuse, on calcule l'écart de phase par interférométrie (intercorrélation normalisée 21 ) entre signaux homologues, c'est-à-dire issus d'un même satellite (S1_0 et S1 b_0) délivrés respectivement par les deux antennes (22, 23) d'une même paire d'antennes, ces antennes étant distantes de moins de λ/2, puis filtrés, démodulés et desétalés (24, 25), on réalise une remise en cohérence par multiplication (26) d'un des signaux (S1_0) par le retard relatif (27) évalué par l'intercorrélation normalisée précédente, les deux signaux étant préalablement retardés (28, 29) de façon à compenser le retard introduit par l'intercorrélation ( 20ms dans le présent exemple), et que l'on asservit (30) la boucle de phase PLL à partir de l'addition vectorielle (31 ) des deux signaux d'antennes remis en cohérence.
3. Procédé selon la revendication 1 ou 2, caractérisé en ce que pour réaliser une mesure d'écart de phase, on réalise une intercorrélation en passant par le domaine fréquentiel.
4. Procédé selon la revendication 3, caractérisé en ce que pour réaliser le traitement d'intercorrelation, après transformation de Fourier par
FFT des sorties « I » et « Q » de corrélation, on effectue une normalisation du spectre par centrage et normalisation de l'écart-type du bruit de fond, par fenêtrage glissant , puis une détection des pointes de bruits, et identification de la porteuse, puis la suppression des pointes de bruit dans le domaine fréquentiel des deux signaux par excision fréquentielle par mise à zéro ou au niveau du bruit de fond des amplitudes spectrales aux fréquences détectées, puis la multiplication complexe des spectres conjugués filtrés (intercorrélation) et le retour en temporel par DFT inverse considérée à zéro, pour obtenir l'intercorrélation entre les 2 signaux et donc leur écart de phase relatif.
5. Procédé selon l'une des revendications précédentes, caractérisé en ce que pour réduire le domaine de recherche, avec une aide partielle, pour la présélection du domaine des ambiguïtés potentielles, on définit un premier domaine de recherche d'ambiguïté par deux hypothèses particulières sur le cap d'aide, et en utilisant les autres aides connues, puis on élargit ce domaine avec les incertitudes sur les autres paramètres d'aide.
6. Procédé selon la revendication 5, caractérisé en ce que, si l'on ne dispose pas d'aide en cap, les hypothèses particulières sont : cap d'aide égal à l'azimut du satellite c'est-à-dire maximum du produit scalaire de projection de la base d'antenne sur cette direction, et cap d'aide égal à l'azimut du satellite + 180°, c'est-à-dire minimum de projection.
7. Procédé selon l'une des revendications précédentes, caractérisé en ce que pour améliorer la mesure d'écart de phase, on utilise simultanément des antennes rapprochées et des antennes éloignées permettant de réaliser des mesures précises d'écart de phase simple, par opposition aux doubles différences, non ambiguës et non affectées des biais de phase RF, initialisées à partir des antennes rapprochées, ou de tout autre moyen inertiel de calcul d'attitude.
8. Procédé selon l'une des revendications précédentes, caractérisé en ce que pour augmenter la robustesse du test de sélection des ambiguïtés, on élargit le domaine de recherche de quelques ambiguïtés supplémentaires.
9. Procédé selon la revendication 8, caractérisé en ce que les ambiguïtés supplémentaires sont : une ou deux ambiguïtés pour compenser les troncatures dues aux calculs par parties entières et une ambiguïté pour le bruit de mesure sur la partie fractionnaire de la phase interférométrique.
10. Procédé selon l'une des revendications précédentes, caractérisé en ce que pour améliorer le contrôle d'intégrité, le test est appliqué sur plusieurs récurrences de mesures d'attitude.
PCT/FR2002/004397 2001-12-20 2002-12-17 Procede d'amelioration de la determination de l'attitude d'un vehicule a l'aide de signaux de radionavigation par satellite WO2003054576A1 (fr)

Priority Applications (3)

Application Number Priority Date Filing Date Title
EP02805393A EP1459093A1 (fr) 2001-12-20 2002-12-17 Procede d'amelioration de la determination de l'attitude d'un vehicule a l'aide de signaux de radionavigation par satellite
CA002470550A CA2470550A1 (fr) 2001-12-20 2002-12-17 Procede d'amelioration de la determination de l'attitude d'un vehicule a l'aide de signaux de radionavigation par satellite
US10/499,199 US7454289B2 (en) 2001-12-20 2002-12-17 Method of improving the determination of the attitude of a vehicle with the aid of satellite radionavigation signals

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR01/16561 2001-12-20
FR0116561A FR2834069B1 (fr) 2001-12-20 2001-12-20 Procede d'amelioration de la determination de l'attitude d'un vehicule a l'aide de signaux de radionavigation par satellite

Publications (1)

Publication Number Publication Date
WO2003054576A1 true WO2003054576A1 (fr) 2003-07-03

Family

ID=8870758

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/FR2002/004397 WO2003054576A1 (fr) 2001-12-20 2002-12-17 Procede d'amelioration de la determination de l'attitude d'un vehicule a l'aide de signaux de radionavigation par satellite

Country Status (5)

Country Link
US (1) US7454289B2 (fr)
EP (1) EP1459093A1 (fr)
CA (1) CA2470550A1 (fr)
FR (1) FR2834069B1 (fr)
WO (1) WO2003054576A1 (fr)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7693211B2 (en) * 2006-04-20 2010-04-06 Sirf Technology Holdings, Inc. Fast fourier transform based phase locked loop for navigational receivers
EP2083282A1 (fr) * 2008-01-28 2009-07-29 Technische Universiteit Delft Système d'émetteur-récepteur
EP2748634B1 (fr) * 2011-08-25 2018-11-21 BAE SYSTEMS Information and Electronic Systems Integration Inc. Détection de vecteur d'orientation gps avec un seul récepteur
CN107248891B (zh) * 2017-06-13 2023-08-04 施浒立 一种用于移动通信天线指向监测的测向测姿装置
JP6821091B2 (ja) * 2018-05-07 2021-01-27 三菱電機株式会社 到来波数推定装置及び到来波数到来方向推定装置
CN109283555A (zh) * 2018-09-28 2019-01-29 西南电子技术研究所（中国电子科技集团公司第十研究所） 卫导波束形成抑制干扰方法
FR3088489B1 (fr) * 2018-11-09 2020-10-23 Office National Detudes Rech Aerospatiales Determination de dephasages d'emission pour un radar a plusieurs voies d'emission juxtaposees

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6061631A (en) * 1997-07-03 2000-05-09 Trimble Navigation, Ltd. Hybrid approach for antenna baseline self-survey and line bias calibration using GPS carrier phase
US6166683A (en) * 1998-02-19 2000-12-26 Rockwell International Corporation System and method for high-integrity detection and correction of cycle slip in a carrier phase-related system
US6259398B1 (en) * 2000-05-19 2001-07-10 Sri International Multi-valued variable ambiguity resolution for satellite navigation signal carrier wave path length determination

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5185610A (en) * 1990-08-20 1993-02-09 Texas Instruments Incorporated GPS system and method for deriving pointing or attitude from a single GPS receiver
US5659318A (en) * 1996-05-31 1997-08-19 California Institute Of Technology Interferometric SAR processor for elevation
US6005514A (en) * 1997-09-15 1999-12-21 United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Method for attitude determination using GPS carrier phase measurements from nonaligned antennas
US6313789B1 (en) * 1998-06-10 2001-11-06 Topcon Positioning Systems, Inc. Joint tracking of the carrier phases of the signals received from different satellites
FR2783929B1 (fr) 1998-09-25 2000-12-08 Sextant Avionique Procede et dispositif de traitement en reception d'un signal l2 de satellite gps
US6516021B1 (en) * 1999-09-14 2003-02-04 The Aerospace Corporation Global positioning systems and inertial measuring unit ultratight coupling method
US7162367B2 (en) * 1999-11-29 2007-01-09 American Gnc Corporation Self-contained/interruption-free positioning method and system thereof

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6061631A (en) * 1997-07-03 2000-05-09 Trimble Navigation, Ltd. Hybrid approach for antenna baseline self-survey and line bias calibration using GPS carrier phase
US6166683A (en) * 1998-02-19 2000-12-26 Rockwell International Corporation System and method for high-integrity detection and correction of cycle slip in a carrier phase-related system
US6259398B1 (en) * 2000-05-19 2001-07-10 Sri International Multi-valued variable ambiguity resolution for satellite navigation signal carrier wave path length determination

Also Published As

Publication number Publication date
US20050043887A1 (en) 2005-02-24
FR2834069B1 (fr) 2004-02-27
EP1459093A1 (fr) 2004-09-22
FR2834069A1 (fr) 2003-06-27
US7454289B2 (en) 2008-11-18
CA2470550A1 (fr) 2003-07-03

Similar Documents

Publication Publication Date Title
EP2520949B1 (fr) Dispositif de réception d&#39;un système de positionnement par satellite comprenant une fonction de détection de faux accrochages
EP3709058B1 (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é.
US20180074207A1 (en) Cepstrum-based multipath mitigation of a spread spectrum radiocommunication signal
EP3081962B1 (fr) Dispositif de poursuite intègre de signaux gnss
US10191158B2 (en) GNSS receiver calculating a non-ambiguous discriminator to resolve subcarrier tracking ambiguities
CA2486079C (fr) Recepteur de positionnement par satellite utilisant deux porteuses de signal
EP2587691B1 (fr) Procédé de traitement coordonné de signaux émis par des balises
EP2095150B1 (fr) Procédé et dispositif de réception d&#39;un signal de radionavigation à modulation boc
EP2609447A1 (fr) Dispositif spatio temporel multi-antennes multi-correlateurs pour la rejection des multi-trajets des systemes de navigation
EP3543745B1 (fr) Dispositif multi-antennes pour la rejection de multi-trajets dans un systeme de navigation par satellite et procede associe
EP2511734A2 (fr) Récepteur bi-fréquences de positionnement par satellites et procédé de réception associé
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
EP1459093A1 (fr) Procede d&#39;amelioration de la determination de l&#39;attitude d&#39;un vehicule a l&#39;aide de signaux de radionavigation par satellite
KR20200099958A (ko) 글로벌 네비게이션 위성 시스템 수신기에서 다중 경로 완화를 제공하는 시스템 및 방법
FR2918765A1 (fr) Procede de determination d&#39;une erreur d&#39;asservissement dans une boucle de poursuite d&#39;un code pseudo-aleatoire.
EP2784546B1 (fr) Procédé de formation de faisceau signaux d&#39;un récepteur de signaux d&#39;un système de navigation par satellites pour améliorer la résistance au brouillage
FR2913114A1 (fr) Procede et systeme de surveillance de boucles de poursuite de code dans un recepteur de positionnement par satellites
WO2024008883A1 (fr) Procédé de détection d&#39;un signal interférant d&#39;un récepteur gnss et dispositif de détection associé
EP4359825A1 (fr) Procédé de détection d&#39;une pluralité de signaux de géolocalisation
EP4343379A1 (fr) Procédé et dispositif de détection d obstacles proches par radar passif multistatique gnss pour plateformes mobiles

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): CA US

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR IE IT LU MC NL PT SE SI SK TR

121 Ep: the epo has been informed by wipo that ep was designated in this application
DFPE Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101)
WWE Wipo information: entry into national phase

Ref document number: 2002805393

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 2470550

Country of ref document: CA

WWE Wipo information: entry into national phase

Ref document number: 10499199

Country of ref document: US

WWP Wipo information: published in national office

Ref document number: 2002805393

Country of ref document: EP