FR 2779234 A1 19991203 FR 9806767 A 19980528 i POSITIONNEMENT PAR CALCUL DE PSEUDO-VITESSES DANS UN SYSTEME DE NAVIGATION PAR SATELLITE La présente invention a pour objet un récepteur pour un système de navigation comprenant une constellation de satellites émettant des signaux. L'invention concerne aussi un procédé de localisation d'un récepteur dans un système de navigation comprenant une constellation de satellites émettant des signaux. Elle s'applique notamment aux satellites en orbite basse. L'invention concerne le domaine des systèmes de navigation par satellites; ces systèmes, à partir d'une constellation de satellites émettant des signaux connus, permettent à un récepteur de calculer sa position dans un repère géocentrique. De tels systèmes sont utilisés pour la navigation des véhicules de toutes natures, et en particulier, pour les bateaux et les aéronefs. Les systèmes existant actuellement sont les systèmes GPS et GLONASS, qui sont constitués de constellations de satellites situés à 20 000 km d'altitude. Ces systèmes sont fondés sur la mesure de pseudo- distances entre les satellites et le récepteur. Cette pseudo-distance est évaluée à partir du temps de transmission d'un signal émis à un instant déterminé par le satellite, et reçu par le récepteur. Pour cela, le récepteur mesure un temps d'arrivée du signal codé émis par un satellite de la constellation, et il lit dans le message reçu le temps d'émission du signal par le satellite.Le récepteur évalue le temps de trajet en déterminant l'instant de réception du message, ou en calculant la phase du signal. La connaissance de la position des satellites permet de déterminer la position du récepteur, par triangulation à partir des pseudo-distances. Les systèmes existant dans l'art antérieur utilisent quatre satellites pour le calcul de la position du récepteur; on obtient ainsi quatre mesures de pseudo-distances, qui permettent de déterminer le temps et la position de la station du récepteur. L'une des caractéristiques principales de ces systèmes est que le nombre d'utilisateurs est illimité, i. e. que le système est non saturable; en effet les satellites sont uniquement émetteurs, et le calcul de la position s'effectue au niveau de chaque récepteur. Celui ci est donc passif, en mode "écoute". Ceci implique aussi que les mesures de pseudo distances se font sur un trajet simple satellite-récepteur. D'autres détails sur la structure et le fonctionnement de ces systèmes de navigation se trouvent dans l'art antérieur maintenant cité. US-A-5 420 592 décrit une utilisation d'un système de navigation GPS pour le positionnement de ballons météorologiques; il est proposé de séparer le récepteur GPS entre une partie mobile portée par le ballon et une partie fixe au sol. Un filtre de Kalman est utilisé pour prévoir la position du ballon et permettre un positionnement même avec une visibilité réduite de la constellation de satellites. Le système corrige les décalages induits par l'effet Doppler sur les signaux des satellites, pour permettre une correction du signal du satellite pour le calcul rapide de la pseudodistance. EP-A-0 460 862 décrit un récepteur de système de positionnement global, du type GPS, susceptible d'être implémenté dans un circuit VLSI. Le récepteur comprend des moyens pour évaluer la pseudo distance, et la variation de pseudo- distance, de sorte à supprimer l'effet Doppler sur la porteuse provenant des satellites. US-A-5 436 632 propose d'utiliser l'effet Doppler pour calculer des corrections à apporter au taux de variation des pseudo-distances, dans un système de surveillance de l'intégrité des signaux d'un système de type GPS. US-A-5 343 209 propose d'évaluer les décalages Doppler en utilisant une estimation de la vitesse du récepteur et des variations de ligne de vision (LOS pour l'anglais "line of sight") d'un système de positionnement type GPS. Cette évaluation d'obtenir permet une estimation de la phase de la porteuse, de la fréquence de la porteuse et de la modulation de phase de chacun des signaux reçus qui est meilleure que celle que pourrait fournir chaque signal indépendamment. EP-A-0 518 146 décrit un système de positionnement multi-point, pour un navire et des bouées tractées. Ce système utilise plusieurs récepteurs GPS, et des filtres de Kalman pour obtenir une évaluation de la qualité de chaque information de position. Il est encore proposé de corriger le décalage Doppler pour l'évaluation des pseudo-distances. US-A-5 594 453 décrit un récepteur GPS avec un mode d'attente, dans lequel la consommation de puissance est réduite; pour permettre une acquisition rapide de la position en sortie du mode d'attente, ce document propose de corriger l'influence de la température sur la fréquence d'horloge. L'effet Doppler est pris en compte pour calculer l'estimation de la fréquence du satellite dans la nouvelle position, en sortie du mode d'attente. US-A-5 557 284 propose un système de détection d'émission de signaux factices dans un système de positionnement global ("spoofing" en langue anglaise); il est proposé d'utiliser les différences d'effet Doppler entre deux récepteurs placés sur un aéronef ou sur le véhicule pour détecter les signaux factices. Dans tous ces documents, comme dans l'art antérieur en général, l'effet Doppler sur les signaux des satellites est perçu comme un inconvénient; il est proposé de corriger les signaux reçus pour supprimer l'effet Doppler, ou de tenir compte de l'effet Doppler pour une prévision de la position du récepteur. Aucun de ces documents ne propose d'utiliser les mesures de pseudo-vitesses par effet Doppler pour améliorer la précision de la position déterminée. En outre, ces documents concernent tous des constellations de satellites en orbite moyenne (MEO ou medium earth orbit en langue anglaise), telles que le GPS ou le GLONASS. L'invention propose au contraire d'utiliser les informations de pseudovitesses fournies par l'évaluation de l'effet Doppler pour le calcul de la position du récepteur. Elle s'applique tout particulièrement à une constellation de satellites en orbite basse (LEO ou low earth orbit en langue anglaise), pour lesquels l'effet Doppler et la variabilité géométrique des satellites visibles sont plus importants. Plus précisément, l'invention propose un récepteur pour un système de navigation comprenant une constellation de satellites émettant des signaux, le récepteur comprenant des moyens de mesure de pseudo-distances entre le récepteur et des satellites à partir des signaux reçus depuis les satellites, des moyens de mesure de pseudo-vitesses entre le récepteur et des satellites à partir des signaux reçus depuis les satellites, et des moyens de calcul de la position instantanée du récepteur à partir de la combinaison des mesures de pseudo-distances et de pseudo-vitesses. Avantageusement, les moyens de mesure des pseudo-vitesses mesurent les pseudo-vitesses à partir du décalage Doppler des signaux reçus depuis les satellites. Dans un mode de réalisation, le système comprend en outre des moyens de calcul de la vitesse du récepteur à partir des pseudo-vitesses mesurées. De préférence, les moyens de calcul calculent la position ou la vitesse du récepteur par une méthode de moindres carrés. Dans un mode de réalisation, les moyens de calcul calculent la position ou la vitesse du récepteur en minimisant la matrice de covariance des pseudo-distances et des pseudo-vitesses mesurées. De préférence, les moyens de calcul comprennent en outre des moyens de prédiction de la position du récepteur, par exemple à l'aide de filtres de Kalman. Avantageusement, la constellation de satellites est une constellation de satellites en orbite basse. L'invention propose en outre un procédé de localisation d'un récepteur dans un système de navigation comprenant une constellation de satellites émettant des signaux, comprenant les étapes de - mesure de pseudodistances entre le récepteur et des satellites à partir des signaux reçus depuis les satellites; - mesure de pseudo-vitesses entre le récepteur et des satellites à partir des signaux reçus depuis les satellites, et calcul de la position du récepteur à partir de la combinaison des mesures de pseudo-distances et de pseudo-vitesses. Avantageusement, l'étape de mesure des pseudo-vitesses comprend la mesure du décalage Doppler des signaux reçus depuis les satellites, et le calcul des pseudo- vitesses à partir du décalage mesuré. Dans un mode de mise en oeuvre, l'étape de calcul de la position du récepteur comprend aussi une étape de calcul de la vitesse du récepteur à partir des pseudo- distances et des pseudo-vitesses mesurées. De préférence, l'étape de calcul comprend le calcul de la position ou de la vitesse du récepteur par une méthode de moindres carrés. Avantageusement, l'étape de calcul comprend le calcul de la position ou de la vitesse du récepteur en minimisant la matrice de covariance des pseudo-distances et des pseudo-vitesses mesurées. L'étape de calcul peut aussi comprendre une étape de prédiction de la position du récepteur, par exemple à l'aide de filtres de Kalman. L'invention concerne enfin l'application d'un tel procédé à la localisation d'un récepteur à l'aide des signaux d'une constellation de satellites en orbite basse. D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de la description qui suit de modes de réalisation de l'invention, donnée à titre d'exemple et en référence aux dessins annexés qui montrent: - figure 1 des courbes de dilution de précision horizontale et verticale dans un système selon l'invention; - figure 2 une courbe d'erreur de navigation horizontale, dans un système selon l'invention; - figure 3 une courbe d'erreur de navigation verticale, dans un système selon l'invention. L'invention propose un système de positionnement global ou système de navigation, dans lequel on utilise les mesures de pseudo-distances, mais aussi des mesures de pseudo-vitesses évaluées par effet Doppler pour calculer la position d'un récepteur. On appelle pseudo-vitesse la vitesse radiale du récepteur par rapport à un satellite. Dans la suite de la présente description, l'invention est décrite en référence à un système de positionnement global utilisant une constellation de satellites en orbite basse (LEO); elle s'applique aussi à des constellations de satellites à des altitudes différentes. Selon l'invention, les satellites de la constellation sont disposés à une altitude moyenne de 1450 km; par rapport aux systèmes existants, l'altitude plus faible induit une plus grande variabilité géométrique. Cette variabilité est a priori moins favorable pour les applications de navigation, dans lesquelles on utilise classiquement quatre satellites visibles en permanence avec une bonne configuration géométrique. L'invention permet, même dans une telle configuration, d'améliorer une performance de positionnement. Dans un premier mode de réalisation, l'invention propose d'utiliser les mesures de pseudo-vitesses en instantané, pour améliorer la résolution du positionnement. On peut alors linéariser les équations reliant la position à estimer et les pseudo-distances à mesurer. En outre, on linéarise aussi les équations reliant la position à estimer et les pseudo-vitesses. On peut alors combiner les mesures de pseudo-distance et de pseudo-vitesses dans une même équation linéaire; cette équation peut être résolue au sens des moindres carrés pour obtenir la position du système. De la sorte. les moyens de calcul utilisent dès l'origine les mesures de pseudo-distances et de pseudo- vitesses; les mesures des pseudo-vitesses sont dès l'origine une variable de la fonction de calcul des positions instantanées, au même titre que les mesures de pseudo-vitesses, à l'inverse de l'art antérieur, o elles n'interviennent que pour la correction de la position calculée auparavant à partir des seules mesures de pseudo- vitesses. On utilise les notations suivantes: Sk: satellite numéro k de la constellation; N: nombre de satellites visibles dans la constellation; Rk: pseudo-distance entre le récepteur et le satellite Sk; c: vitesse de la lumière; pk: distance géométrique entre le récepteur et le satellite Sk; b: biais de l'horloge du récepteur par rapport au temps de référence du système de postionnement global; Ek: bruit de mesure sur la pseudo-distance Rk; (xsk XSk XSk) coordonnées du satellite Sk; (X1 2 3 (Xl, X2, X3): coordonnées du récepteur (point M); (X 1,0, X2,0, X3,0): coordonnées initiales du récepteur (point Mo); 6X = (XI - XI0 X2 - X2,0 X3 - X3.)0): coordonnées du vecteur du point M0 au point M; u0k: vecteur unitaire du point Mo au satellite Sk. Pour chaque satellite, on recherche dans le récepteur le maximum de corrélation entre le code de bruit pseudo-aléatoire PN du signal émis par le satellite à un instant donné et le paradigme PN du satellite, stocké dans le récepteur. On peut de façon connue en soi utiliser pour cela une boucle à verrouillage de phase. Cette recherche du maximum de corrélation fournit une différence de temps AT et permet d'évaluer la pseudo-distance, par simple produit par la vitesse de la lumière; on peut relier la pseudo-distance à la distance entre le récepteur et le satellite par l'équation: Rk = c.AT = pk + c.b + Ek (1) On peut linéariser cette équation par rapport aux quatre inconnues que sont la position et le biais de l'horloge, en partant de la position initiale connue M0 du satellite; au premier ordre, on trouve: 6R/ = Rk - R = pk _po +c.b+Ek R ' ( pk P k + c b + k apk X2 p X3 _ k _( p k k a k 8Rk = __P _P k '8X'+c.b+ s' ag ax=, a Le calcul des dérivées partielles fournit: apk xSk -X pk d'o: Lapk P) pk_k' ax, ax, aDX3 En considérant tous les N satellites visibles. et en notant 8R le vecteur composé des mesures, en différence par rapport à la postion M0 connue du récepteur aR = (R'... R) - i! -u0 1 A=...... -U Ni 1 -- 0 et Y le vecteur des inconnues, L,=L "C b, on obtient 8R=A-Y+ On utilise un minimum de quatre satellites visibles pour déterminer quatre pseudo-vitesses et ainsi déterminer les trois inconnues de position et l'inconnu de biais d'horloge. L'équation peut être résolue suivant une méthode des moindres carrés: y (A- A) A'.(R Toutefois, l'invention propose d'utiliser non seulement cette simple méthode de positionnement à partir de la mesure des pseudo-distances, mais en outre une mesure des pseudo-vitesses à partir de l'effet Doppler, pour déterminer la position du récepteur. L'observation de l'effet Doppler introduit l'équation suivante: - + f c ou -X.Af =-D= p +c.b +l avec les notations suivantes: vitesse radiale vraie; b À dérive de l'horloge du récepteur; D mesure Doppler convertie en distance; TD temps d'acquisition de la mesure de Doppler; D' D*TD; l À bruit résiduel de la mesure Doppler; uk: vecteur unitaire du récepteur vers le satellite d'indice k; VS: vitesse du satellite Sk; V: vitesse du récepteur. L'invention propose de linéariser le système en introduisant une position approchée M0 et une vitesse approchée Vo du récepteur. On peut alors écrire, pour un satellite Sk donné: 1 5 D - =(VV u). u + c.b +rlk -D =(V -6+(V 6V)) uo + ( +c-b+n avec ôV = V- VO aUk aUk aUl i k axi axi axi C Uk) CôUk au k aUk _ U 5u I- 2 3 ax ÈX =X2 5X3 = X DX2 DX2 aX3 aX3 aX3, Or, YUk a a MS k a À XSA __ x axX, ax, J) ax. J pk auk x x', -Z X ap k ô xjj; X- x j x, aX j ( k) pk (pk)2 p k aUx, 6 k U OXi p k pk d'o (D y X) k -, - 3) En développant cette équation et en négligeant le terme du second ordre AV'.X', il vient: 8DkS =_k___"- _Vo' uk 8Dk = -D (' -() -. U 8& = -ôo V + (V j-o)(k*( I3)*A ÈD' probables des récepteur. L'équation précédente devient alors: VS' r 1.'l -13)' D =....... ..À +.; P N 0- p AZ (?L U 13) -U0 lJ On combine les mesures de pseudo- distances et de pseudo-vitesses dans la même équation linéaire -uo 0 1 0 N,.... ... ... Y' c. + X 1 _R OU' 1 0 I =- V''1 /, (U / -uI 0 1 c + 1 8D, V I1 L LI1 _ O 1 C r =D."... ....C-I. V N U*' -I3) -U 0 1 _] -U Pour assurer l'homogénéité des inconnues, on peut multiplier certains termes par le temps Td d'acquisition de la mesure des pseudo-vitesses: I t 0 10 ... _... 8..DTD: 86R -u0 N 0 1 0 V.TD N =_ V SVTD/P (uI u") - o IL cb. It i t' Vs'.TD/P'. (U.U -I3) -U O 1 vS,'.v/pN.$. uNtI3) _u 0 1) On obtient une équation du type Z = A.Y + s avec: 58 Y = 8 TI) Y= et Z= /c L Cette équation comprend 8 inconnues pour 2N mesures, et elle peut être résolue en utilisant quatre satellites visibles. Il est clair que si plus de quatre satellites sont visible, on peut utiliser certains pour la mesure de pseudo-distance et d'autres pour la mesure de pseudovitesses; toute combinaison des satellites permettant d'obtenir un nombre de mesures suffisant pour résoudre l'équation peut être utilisée. L'équation est avantageusement résolue au sens des moindres carrés; selon l'invention, pour assurer une solution optimale, on recherche une solution qui minimise la matrice de covariance des mesures de pseudo-vitesses et de pseudo- distances. En notant: Y = cov(Z) la matrice de covariance des mesures, on recherche à minimiser min Z- A - La solution est donnée par Y('.'A) A' A '.Z De la sorte, l'invention ne propose pas simplement une juxtaposition des calculs de pseudo-vitesses et de pseudo-distances, mais bien une utilisation combinée des mesures de pseudo-distance et des mesures de pseudo-vitesses pour calculer la position du récepteur. L'invention permet d'obtenir une meilleure précision dans la détermination de la position du récepteur. La précision de la localisation peut s'évaluer par le calcul de la covariance de la solution obtenue, qui donne une indication sur la position, la vitesse, et le temps. Cette covariance s'écrit: Cov(Y)= (A'. -' A) Les tableaux suivants montrent les résultats des mesures de pseudo-distances, dans le cas d'un système de navigation utilisant une constellation de satellites en orbite basse, et pour un récepteur disposé à bord d'un avion. On a utilisé pour les tableaux le calcul de la précision à l'aide de la covariance de la solution. suivant l'équation précédente. La première colonne donne les résultats d'un positionnement utilisant simplement les mesures de pseudo-vitesses comme dans les systèmes de l'art antérieur. La deuxième colonne donne les résultats d'un positionnement selon l'invention. Les modes correspondent respectivement à - mode 0: mode simple ou absolu, monofréquence sans correction différentielle; - mode 1: mode différentiel local; - mode 2: mode bifréquence; - mode 3: mode différentiel local amélioré. Dans le mode 3, on utilise un corrélateur étroit dans le récepteur; le code est filtré par la phase dans le récepteur de la station différentielle; on utilise aussi avantageusement une correction différentielle sur les pseudo-distances; une correction différentielle sur les pseudo-vitesses pourrait aussi être mise en oeuvre, mais elle n'est pas appliquée pour l'obtention des résultats des tableaux. Erreur de navigation horizontale (mètres) Mode pseudo- distances invention mode 0 20 15 mode 1 4,5 4 mode 2 18 11 mode 3 1,7 1,6 La limite supérieure de l'erreur de navigation ou NSE (acronyme de l'anglais "navigation sensor error") pour les système de catégorie 1 est de 18,7 m; l'invention permet d'atteindre cette limite, pour l'ensemble des modes. Erreur de navigation verticale (mètres) Mode pseudo- distances invention mode 0 46 15 mode 1 11 7,3 mode 2 41 9,7 mode 3 4,2 3,8 La limite supérieure de l'erreur de navigation ou NSE pour les système de catégorie 1 est de 5,4 m; l'invention permet d'atteindre cette limite, au moins pour le mode 3. l1 Dans tous les cas, l'invention permet d'améliorer significativement les résultats d'un positionnement classique obtenu par les seules mesures de pseudodistances. Ces résultats apparaissent aussi dans les figures I à 3. qui sont des représentations graphiques de simulations de l'évolution dans le temps des NSE, pour le mode 0, et pour une constellation de satellites en orbite basse. La figure 1 sans dimension en ordonnées, et en heures en abscisses. Elle montre en trait pleins le nombre de satellites visibles, en trait interrompus longs la dilution de précision horizontale, et en traits interrompus courts la dilution de précision verticale. La dilution de précision horizontale présente une valeur moyenne de 0,92, tandis que la dilution de précision verticale présente une valeur moyenne de 2. La figure 2 montre les valeurs de l'erreur de navigation, pour les mêmes simulations que celles des tableaux. Les ordonnées sont graduées en mètres, tandis que les abscisses sont encore en heures. Le trait interrompu horizontal à la valeur de 18,7 mètres correspond à la limite de la catégorie 1. Le trait plein montre l'erreur de navigation dans un système classique, o la position est évaluée uniquement par un calcul de pseudo-distance; comme indiqué dans le tableau ci-dessus, la moyenne de l'erreur de navigation est de l'ordre de 20 mètres, et est supérieure à la limite de la catégorie 1. Le trait pointillés montre l'erreur de navigation dans le récepteur de l'invention, qui présente une moyenne de 15 mètres. On constate en outre que l'erreur de navigation selon l'invention présente une excursion relativement faible, et reste toujours inférieure au seuil limite de la catégorie 1. La figure 3 montre avec les mêmes notations les résultats correspondant pour l'erreur verticale. La limite de la catégorie 1 est de 5,4 mètres, et les valeurs moyennes selon l'invention et selon la méthode classique sont respectivement de 46 et 15 mètres. On constate encore sur les figures que l'invention permet d'obtenir de meilleurs résultats, non seulement en valeur moyenne, mais à tout instant. Dans un autre mode de réalisation, l'invention propose en outre d'utiliser un modèle dynamique pour la détermination des coordonnées du récepteur, à partir des mesures de position et de vitesses. On peut notamment utiliser à cet effet un filtrage de Kalman pour prédire la position du récepteur. Le filtre de Kalman utilisé fonctionne avec un vecteur d'état En: En = (X, V, A, u)t o X est la position du récepteur, V sa vitesse, A l'accélération, et u le bruit à l'instant n. On note P0 la matrice de covariance de E0 à l'instant initial, ou dans la position M0, avec la vitesse Vo et l'accélération A0, et Pn la matrice de covariance de En. L'équation du modèle dynamique est la suivante: En+l = F. En + w o n est l'indice de l'instant. wl est un terme d'erreur du modèle, dépendant du bruit, qui est supposé être un bruit blanc gaussien, et F est la matrice de propagation associée au filtre de Kalman. On note Q la covariance de wl. Q = cov(wl)= c-2 13 Dans cette expression, C est le facteur d'innovation, qui caractérise la confiance que l'on place dans le modèle de propagation, et qui permet un réglage du filtre. Une valeur de l'ordre de ?? pour le facteur d'innovation s'avère fournir de bons résultats. La matrice de propagation F s'écrit: At2 At3 1 3 At3213 6 13 At2 A1 3 At3 2 13 F= 13 At- 13 L'équation de mesure s'écrit Y= H. E + w2 o, comme dans le mode de réalisation précédent, Y est la matrice des inconnues, E est le vecteur d'état estimé par le filtrage de Kalman, H est la matrice de mesure définie plus bas, et w2 est le bruit de mesure, soit l'erreur de positionnement. La matrice de bruit R, qui est égale à la covariance de ce bruit de mesure est donnée par l'étape de positionnement décrit dans le mode de réalisation précédent. Ainsi, dans ce mode de réalisation comme dans le précédent, on utilise la matrice de covariance des positions et des vitesses, et on utilise cette matrice de covariance pour le calcul de la position approchée servant d'entrée au modèle dynamique. A chaque instant de mesure, le filtrage de Kalman se compose de deux étapes, une étape de prédiction - passage de En à une valeur estimée de En+l - et une étape de correction de cette valeur extimée. On peut utiliser le filtrage suivant: En+l = F.En et Pn+l = F.Pn.Ft + Q pour l'étape de prédiction: et K = Pn+I.Ht.-yI, qui traduit le gain de Kalman, Z = R + H.P.Ht. H étant définie plus bas (Question: indice sur P ???) En+l = En+l + K.(Yn+l - H.En+l) Pn+l = (I - K.H).Pn+1 pour l'étape de correction de la prédiction. P0 s'écrit Ro Rv Po = R, R0 étant la matrice des covariances des mesures de position et de vitesse, Ra est choisie comme 0.l2 Ro = 0.12 en, (m/s2)2. 0.12 On peut prendre pour l'état initial E0 la position et la vitesse calculées, et une valeur approchée de l'accélération. La matrice de mesure H s'écrit 13 03 03 03 03 13 03 03 On procède ainsi, à chaque instant, à une mesure des pseudo-distances et des pseudo-vitesses, puis ensuite à un calcul de la position et de la vitesse approchées à partir de ces mesures. On peut utiliser comme dans le premier mode de réalisation une méthode de moindre carrés itératifs, en minimisant les erreurs. La position et la vitesse approchés sont utilisée en entrée du modèle dynamique pour obtenir une prédiction du vecteur d'état du récepteur. Cette prédiction est corrigée pour fournir une estimation de la vitesse et de la position du récepteur. L'invention permet encore dans ce cas d'améliorer la précision dans le positionnement, comme le montrent les tableaux suivants, toujours dans le cas d'une constellation de satellites en orbite basse. Les notations sont les mêmes que dans le mode de réalisation précédent; on a en outre, pour le cas du positionnement par calcul uniquement des pseudo-distances, mentionné les résultats avec ou sans filtrage de Kalman; pour le positionnement selon l'invention, on a aussi indiqué les résultats avec ou sans utilisation du modèle dynamique. Erreur de navigation horizontale (mètres) Mode pseudo- distances invention sans avec sans avec filtrage filtrage filtrage filtrage mode 0 24 24 9%8 2.8 mode 1 4,4 3,2 3.8 1,9 mode3 1,8 1,4 1,7 1,3 Erreur de navigation verticale (mètres) Mode pseudo-distances invention sans avec sans avec filtrage filtrage filtrage filtrage mode 0 32 32 11,4 4.2 mode 1 12 7,4 7,8 5.6 mode 3 4,6 3,6 4 2,4 De nouveau, les tableaux montrent que l'invention permet d'obtenir des résultats plus précis que l'art antérieur utilisant simplement des pseudo-distances. Le mode de réalisation de l'invention utilisant un modèle dynamique pour prédire la position du récepteur permet encore d'améliorer les résultats obtenus. Par rapport aux limites de la catégorie 1 de systèmes de navigation, l'invention permet d'atteindre dans les modes 0, 1 et 3 la catégorie 1, en erreur horizontale; elle permet encore, dans ce deuxième mode de réalisation, d'atteindre cette catégorie pour l'erreur verticale, très facilement pour le mode 3, et de façon marginale pour les modes 0 et 1. Ce deuxième mode de réalisation permet donc d'obtenir un positionnement plus précis, toujours en utilisant à la fois les mesures de pseudodistances et les mesures de pseudo-vitesses pour la détermination de la position. Bien entendu, la présente invention n'est pas limitée aux exemples et modes de réalisation décrits et représentés, mais elle est susceptible de nombreuses variantes accessibles à l'homme de l'art. On pourrait ainsi utiliser des satellites différents pour les différentes mesures, si plus de quatre satellites sont visibles. Il est aussi clair que l'invention n'est pas limitée au mode de réalisation préféré d'une constellation en orbite basse, même si c'est pour ce type d'orbite que les variations de pseudo-vitesses mesurées par le décalage Doppler sont les plus importantes. REVENDICATIONS 1.- Récepteur pour un système de navigation comprenant une constellation de satellites émettant des signaux, le récepteur comprenant des moyens de mesure de pseudo-distances entre le récepteur et des satellites à partir des signaux reçus depuis les satellites, des moyens de mesure de pseudo-vitesses entre le récepteur et des satellites à partir des signaux reçus depuis les satellites, et des moyens de calcul de la position instantanée du récepteur à partir de la combinaison des mesures de pseudo- distances et de pseudo-vitesses. 2.- Système selon la revendication 1, caractérisé en ce que les moyens de mesure des pseudo-vitesses mesurent les pseudo-vitesses à partir du décalage Doppler des signaux reçus depuis les satellites. 3.- Système selon la revendication 1 ou 2, caractérisé en ce qu'il comprend en outre des moyens de calcul de la vitesse du récepteur à partir des pseudo-vitesses mesurées. 4.- Récepteur selon la revendication 1, 2 ou 3, caractérisé en ce que les moyens de calcul calculent la position ou la vitesse du récepteur par une méthode de moindres carrés. 5.- Récepteur selon l'une des revendications 1 à 4, caractérisé en ce que les moyens de calcul calculent la position ou la vitesse du récepteur en minimisant la matrice de covariance des pseudo-distances et des pseudo-vitesses mesurées. 6.- Récepteur selon l'une des revendications 1 à 5, caractérisé en ce que les moyens de calcul comprennent en outre des moyens de prédiction de la position du récepteur, par exemple à l'aide de filtres de Kalman. 7.- Récepteur selon l'une des revendications 1 à 6, caractérisé en ce que la constellation de satellites est une constellation de satellites en orbite basse. 8.- Procédé de localisation d'un récepteur dans un système de navigation comprenant une constellation de satellites émettant des signaux, comprenant les étapes de - mesure de pseudo-distances entre le récepteur et des satellites à partir des signaux reçus depuis les satellites; mesure de pseudo-vitesses entre le récepteur et des satellites à partir des signaux reçus depuis les satellites, et - calcul de la position du récepteur à partir de la combinaison des mesures de pseudo-distances et de pseudo-vitesses. 9.- Procédé selon la revendication 8, caractérisé en ce que l'étape de mesure des pseudo-vitesses comprend la mesure du décalage Doppler des signaux reçus depuis les satellites, et le calcul des pseudo-vitesses à partir du décalage mesuré. 10.- Procédé selon la revendication 8 ou 9, caractérisé en ce que l'étape de calcul de la position du récepteur comprend aussi une étape de calcul de la vitesse du récepteur à partir des pseudo-distances et des pseudo-vitesses mesurées. 11.- Procédé selon la revendication 8, 9 ou 10, caractérisé en ce que l'étape de calcul comprend le calcul de la position ou de la vitesse du récepteur par une méthode de moindres carrés. 12.- Procédé selon l'une des revendications 8 à 11, caractérisé en ce que l'étape de calcul comprend le calcul de la position ou de la vitesse du récepteur en minimisant la matrice de covariance des pseudo-distances et des pseudo-vitesses mesurées. 13.- Procédé selon l'une des revendications 8 à 12, caractérisé en ce que l'étape de calcul comprend une étape de prédiction de la position du récepteur, par exemple à l'aide de filtres de Kalman. 14.- Application du procédé selon l'une des revendications 8 à 13 à la localisation d'un récepteur à l'aide des signaux d'une constellation de satellites en orbite basse.