FR 2696851 A1 19940415 FR 9211960 A 19921008 PROCEDE DE CALCUL DE LA POSITION D'UN MOBILE PAR UN RECEPTEUR GPS La présente invention concerne un procédé de calcul de la position spatiale et temporelle d'un mobile par un récepteur GPS. Le système GPS (Global Positioning System) est un système bien connu mis au point pour déterminer, à partir de la position, connue à chaque instant, d'un certain nombre de10 satellites, la position spatiale et temporelle d'un mobile muni d'un récepteur adapté au système, dit récepteur GPS, et capable de recevoir et de traiter les données fournies par les satellites du système Le système GPS comporte dix-huit satellites disposés de sorte que la plupart des points de la15 Terre soient "visibles" par au moins quatre satellites. Actuellement, on utilise jusqu'à neuf satellites pour déterminer la position d'un mobile sur la Terre, mais on comprend bien que seuls quatre satellites sont nécessaires pour déterminer la position spatiale et temporelle qui20 dépend de quatre variables: x, y et z, coordonnées spatiales du mobile dans un référentiel prédéterminé, et t, coordonnée temporelle du mobile (c'est-à-dire date à laquelle le mobile se trouve à la position calculée, cette date étant donnée dans le système de temps universel). 25 Le récepteur GPS associé à un mobile comprend des moyens lui permettant: de sélectionner un ensemble de quatre satellites (ou plus) parmi tous les satellites en vue depuis le point o se trouve le mobile associé,30 de calculer, à partir des données fournies par les satellites sélectionnés, les positions spatiale et temporelle du mobile associé. Le système GPS classique garantit de manière sûre à % une erreur maximale de 100 mètres sur la position horizontale calculée et de 157 mètres sur la position verticale calculée En d'autres termes, il y a 95 % de chances que l'erreur sur la position horizontale soit inférieure ou égale à 100 mètres et que l'erreur sur la position verticale soit inférieure ou égale à 157 mètres On dit que le système GPS garantit une précision globale, cette 5 précision globale étant une constante du système. Habituellement, on utilise un critère de sélection des quatre satellites à utiliser par le récepteur GPS pour le calcul de position Ce critère permet de sélectionner, parmi toutes les constellations, c'est-à-dire tous les ensembles,10 de quatre satellites possibles à partir de tous les satellites en vue, celle qui permet d'obtenir, dans le calcul de position, la précision globale garantie par le système GPS Ce critère, noté C, est proportionnel à i Trl(AT A)-1 l, o15 la matrice A donnée par a 11 a 12 a 131 a 21 a 22 a 23 1 A = a 31 a 32 a 33 1 a 41 a 42 a 43 i aij, pour i variant de 1 à 4 et j variant de 1 à 3, étant le cosinus directeur de l'angle entre le j-ième axe de coordonnées et la droite joignant le mobile au i-ième20 satellite, les indices 1, 2 et 3 étant respectivement associés aux axes des abscisses, des ordonnées et des cotes, Tr symbolise la trace de la matrice (AT A)- l, c'est-à-dire la somme de ses coefficients diagonaux, AT est la transposée de la matrice A et A-1 l'inverse de la matrice A.25 Le coefficient de proportionnalité utilisé dans le calcul de C prend en compte toutes les marges d'imprécision du système GPS Ce coefficient, constant, a donc une valeur très importante, et ce pour des questions de sécurité que l'on comprend bien.30 Le critère C est associé à un seuil de tolérance, de sorte que s'il dépasse ce seuil, on suppose que la constellation correspondante ne permet pas de garantir la précision globale; le récepteur GPS sélectionne alors une autre constellation de quatre satellites, et ainsi de suite jusqu'à trouver le constellation correspondant à une valeur de C inférieure au seuil de tolérance. 5 Compte tenu de la valeur importante du coefficient de proportionnalité, le critère C est systématiquement pénalisé Ainsi, il peut arriver qu'une constellation soit rejetée parce qu'elle ne satisfait pas au critère global C alors qu'en pratique, la précision qu'elle permet d'obtenir10 sur la position spécifique calculée est tout de même acceptable. D'autre part, dans certaines applications, comme par exemple la construction de tunnels, l'installation d'antennes de satellites, les rallyes sportifs, la15 navigation de plaisance, entre autres, le récepteur GPS peut effectuer le calcul de position en intégrant sur des durées très importantes, de sorte que les différents bruits sont minimisés et que la position calculée a finalement une précision plus grande que la précision globale garantie par20 le système GPS. Ceci n'est pas possible dans d'autres applications, comme l'aviation civile notamment En effet, dans ce cas, il n'est pas possible d'attendre pour permettre au système d'intégrer sur une durée importante La précision obtenue25 est alors la précision globale garantie, ce qui n'est pas suffisant dans le domaine de l'aviation De plus, les 5 % d'incertitude sont intolérables, on le comprend bien, dans ce domaine. En conséquence, compte tenu des dangers qu'il fait encourir, le système GPS classique n'a actuellement qu'un usage restreint aux applications dans lesquelles il n'est pas très pénalisant de calculer la position en intégrant sur des durées très importantes. L'utilisation du critère C est donc d'une part trop contraignante et pénalisante, et d'autre part insuffisamment précise dans de nombreuses applications. Dans le domaine militaire, le système GPS utilise une fréquence de transmission différente de celle du système classique Dans ce cas, les données transmises sont codées et impossibles à déchiffrer par un utilisateur du domaine 5 civil qui ne possède pas le décodeur adapté Le système GPS garantit alors à 95 % une erreur maximale de 30 mètres dans les trois dimensions Le critère C utilisé est le même. La globalité de la précision garantie dans ce cas, et l'utilisation du critère C, posent les mêmes problèmes que l'utilisation du système GPS dans le domaine civil. Le but de la présente invention est par conséquent de mettre au point un procédé permettant de fournir à l'utilisateur du système GPS classique un critère moins contraignant et plus précis que le critère C.15 Le but de la présente invention est également de permettre une sélection plus affinée de la constellation de satellites à utiliser dans le calcul de la position d'un mobile. La présente invention propose à cet effet un procédé de calcul, par un récepteur GPS porté par un mobile, de la position spatiale et temporelle dudit mobile dans un référentiel prédéterminé, les coordonnées de ladite position étant notées x, y, z et t, o x, y et z sont respectivement l'abscisse, l'ordonnée et la cote dudit mobile dans ledit référentiel, et t la date, dans le système de temps universel, à laquelle ledit mobile se trouve dans la position (x, y, z), ladite position étant calculée à partir de données fournies par au moins quatre satellites du système GPS sélectionnés par ledit récepteur GPS pour former30 une constellation dite de calcul, ledit procédé étant caractérisé en ce qu'il comporte une étape d'estimation de l'erreur commise dans ledit calcul de position, ladite estimation de l'erreur étant représentée par un vecteur noté e dont les composantes sont35 respectivement l'estimée de l'erreur sur l'abscisse dx, l'estimée de l'erreur sur l'ordonnée dy, l'estimée de l'erreur sur la cote dz et l'estimée de l'erreur temporelle c.dt, c étant la célérité de la lumière dans le vide, et en ce que ladite estimation est effectuée par ledit récepteur GPS pour chaque position calculée, de sorte que ladite 5 estimation de l'erreur commise varie avec la position dudit mobile. Selon l'invention, on estime l'erreur commise lors de chaque calcul de position Par conséquent, la précision n'est plus globale, c'est-à-dire uniforme quelle que soit la10 position du mobile, mais varie en fonction de cette dernière Ceci permet à l'utilisateur d'accepter ou de rejeter la position calculée, selon que l'erreur commise dans le calcul de cette position et estimée au moyen du procédé de l'invention dépasse ou non un seuil de tolérance15 prédéterminé par exemple L'estimation de l'erreur est donc dynamique, c'est-à-dire variable en fonction de la position du mobile, et instantanée. De manière très avantageuse, la constellation comporte au plus tous les satellites du système GPS visibles depuis le mobile. Plus le nombre de satellites utilisés dans le calcul est important, plus la fiabilité des résultats obtenus avec le procédé selon l'invention est grande. En pratique, a est donné par la relation: Sr = K s o K est une matrice à N lignes et quatre colonnes, N étant le nombre de satellites que comporte ladite constellation, telle que a 11 a 12 a 13 i a 21 a 22 a 23 1 K a Nl a N 2 a N 3 i aij, pour i variant de 1 à N et j variant de 1 à 3, étant le cosinus directeur de l'angle entre le j-ième axe de coordonnées et la droite joignant le mobile au i-ième satellite de la constellation de calcul, les indices 1, 2 et 5 3 étant respectivement associés aux axes des abscisses, des ordonnées et des cotes, 8 r est un vecteur à N composantes notées URE 1 à UREN, ou UR Ei, pour i variant de 1 à N, est une donnée fournie par le système GPS et donnant l'écart radial moyen entre la10 position théorique du i-ième satellite et sa position réelle sur la droite joignant le i-ième satellite à un point de la zone géographique de la Terre vue par le i-ième satellite, la moyenne étant effectuée sur les points de cette zone géographique o se trouve une station de relais au sol du15 système GPS. Lorsque N est strictement supérieur à 4, la résolution de l'équation précédente par la méthode des moindres carrés conduit à la solution suivante donnant la meilleure estimée E de e:20 E = (KT K)> 1 KT 8 r. La solution apportée par l'invention provient de la constatation que la caractéristique URE (User Range Error) fournie par le système GPS régulièrement dans le temps, permet d'évaluer l'erreur commise lors de chaque calcul de25 position En pondérant le vecteur 8 r par la matrice K, on obtient alors la meilleure estimée de l'erreur commise dans le calcul de la position. Afin de fournir à l'utilisateur un critère facilement exploitable, le procédé selon l'invention peut comporter en outre les opérations suivantes: une valeur P, appelée précision, est calculée à partir des composantes de a ou de E, cette précision étant donnée par la formule: p = gdx 2 + dy 2 + dz 2 + (c dt)2, P est comparée à un seuil de tolérance préalablement fixé, si P dépasse le seuil précédent, le récepteur GPS peut soit sélectionner une autre constellation de calcul, soit rejeter la position calculée. Si l'on utilise E au lieu de ú, on remplace les composantes de a par celles de E dans le calcul de P. D'autres caractéristiques et avantages de la présente invention apparaîtront dans la description suivante donnée à titre d'exemple illustratif et nullement limitatif. Il est possible de déterminer de manière théorique une équation faisant intervenir le vecteur e représentant l'erreur commise sur la position calculée et dont les composantes sont respectivement l'erreur sur l'abscisse, l'erreur sur l'ordonnée, l'erreur sur la cote et l'erreur temporelle Cette équation est la suivante Èr = A e ( 1), O br est un vecteur à quatre composantes notées AR 1 à AR 4. A Ri, pour i variant de 1 à 4, est un biais théorique. L'équation ( 1) est utilisée classiquement pour prouver que le critère C est suffisant pour garantir la précision globale sur la position calculée.20 On pourrait penser à utiliser cette équation pour mettre en oeuvre le procédé de l'invention et déterminer l'erreur exacte commise sur le calcul de position. Toutefois, les A Ri sont des données théoriques qui n'ont pas de réalité physique, c'est-à-dire qui n'existent pas, au25 niveau du système GPS Ces données théoriques correspondent aux erreurs radiales moyennes sur la distance entre chaque point de la terre et le satellite considéré Il est bien évident qu'une telle donnée n'est pas calculable. Le mode de mise en oeuvre préféré de la présente invention repose sur le fait que, pour avoir une bonne estimation du vecteur e, on peut remplacer, dans l'équation ( 1), les A Ri par des données notées UR Ei dont on va à présent expliquer l'origine et l'intérêt. Dans le système GPS, chaque satellite est associé à une valeur d'URE (User Range Error) Pour déterminer l'URE relatif à un satellite donné, toutes les stations de relais au sol situées en vue du satellite considéré calculent l'écart radial, c'est-à-dire l'écart sur la droite station au sol-satellite, entre la position réelle du satellite et sa position théorique sur la trajectoire prévue avant le 5 lancement Cela est possible puisque chaque station au sol connaît exactement sa propre position Ainsi, chaque station au sol fait calculer sa position par les satellites du système et en déduit la position réelle de ces derniers à partir de la différence entre sa position calculée et sa10 position exacte L'URE est alors calculé comme la moyenne des écarts radiaux déterminés par toutes les stations au sol situées en vue du satellite considéré. L'URE est une donnée ayant une réalité physique au niveau du système GPS Elle est par conséquent accessible au récepteur GPS Elle varie au cours du temps et avec la zone géographique dans laquelle se trouve le mobile. De plus, les URE sont calculés et rechargés régulièrement au niveau des satellites par les stations au sol, de sorte qu'ils évoluent avec la dérive du satellite.20 Enfin, l'URE de chaque satellite est affecté d'un biais qui peut permettre au récepteur GPS d'estimer et de tenir compte de l'évolution de l'URE entre deux chargements depuis les stations au sol. Ainsi, le calcul effectué en remplaçant les A Ri par les UR Ei dans l'équation ( 1) permet d'estimer de manière rigoureuse le vecteur e On obtient alors une estimation de la précision de la position calculée, cette estimation n'étant plus globale, c'est-à-dire uniforme pour tous les points de la terre, mais associée à une zone géographique30 relativement petite, de sorte que l'erreur estimée peut être assimilée à l'erreur réelle sur la position calculée. Pour estimer e, le récepteur GPS commence par sélectionner une constellation de satellites de la même manière que celle décrite dans le préambule, c'est-à-dire à35 l'aide du critère C éventuellement généralisé pour pouvoir utiliser N satellites, o N est supérieur ou égal à 4 et inférieur ou égal au nombre de satellites visibles par le mobile. Si N est égal à 4, on peut calculer directement un vecteur 8, dont les composantes sont l'estimée de l'erreur sur l'abscisse, notée dx, l'estimée de l'erreur sur l'ordonnée, notée dy, l'estimée de l'erreur sur la cote, notée dz, et l'estimée de l'erreur temporelle notée c dt, par l'équation er = K ú, o: K est une matrice à N lignes et quatre colonnes telle que: a 11 a 12 a 13 1 a 21 a 22 a 23 1 K= a N 1 a N 2 a N 3 aij, pour i variant de 1 à N et j variant de 1 à 3, étant le cosinus directeur de l'angle entre le j-ième axe de coordonnées et la droite joignant le mobile au i-ième15 satellite, les indices 1, 2 et 3 étant respectivement associés aux axes des abscisses, des ordonnées et des cotes, ar est un vecteur à N composantes notées URE 1 à UREN, O UR Ei, pour i variant de 1 à N, est la valeur de l'URE associée au i-ième satellite.20 Lorsque N est strictement supérieur à 4, il y a plus d'équations que d'inconnues La meilleure estimée E de a est alors donnée, en utilisant la méthode des moindres carrés, par l'équation suivante: E = (KT K)-l KT -r, o KT est la transposée de la matrice K et o (KT K)-1 est l'inverse de la matrice (KT K). Dans certaines applications, on peut utiliser directement les composantes de E (notées Dx, Dy, Dz et c.Dt) Dans ce cas, par exemple, on aura défini au préalable30 un vecteur de référence Eréf, et la comparaison des composantes de même rang des vecteurs E et Eréf permet de conserver la position calculée ou de la rejeter si les écarts entre certaines des coordonnées de E et de Eréf sont trop importants. Il est toutefois préférable d'effectuer la comparaison sur les normes des vecteurs E et Eréf On calcule ainsi une valeur, appelée précision et notée P, définie par la formule: p = 4 Dx 2 + Dy 2 + Dz 2 + (c Dt)2, o Dx, Dy, Dz et c Dt sont les composantes de E. Les considérations précédentes concernant le calcul de P s'appliquent également en remplaçant E par a et Dx, Dy, Dz et c Dt respectivement par dx, dy,-dz et c dt. Si P dépasse un seuil Préf, déterminé au préalable en fonction de la précision recherchée, le récepteur GPS peut par exemple l'indiquer à l'utilisateur à l'aide d'une alarme visuelle ou sonore Le résultat de la comparaison entre P et Préf peut également servir à poursuivre le programme implanté dans le récepteur d'une manière différente, en20 générant par exemple un bit à zéro ou à un selon que P est inférieure ou supérieure à Préf. L'invention permet donc d'établir un critère de rejet simple et très fiable de la position calculée par le récepteur GPS S'il y a rejet, le récepteur peut25 sélectionner une autre constellation, et effectuer le calcul de position une seconde fois Si la précision n'est pas encore satisfaisante, le récepteur GPS sélectionne d'autres constellations jusqu'à ce que la précision souhaitée soit obtenue.30 Dans le cas de l'utilisation dans l'aviation civile, l'utilisateur peut, grâce à l'invention, soit accepter la position calculée par le récepteur GPS si la précision déterminée selon l'invention est suffisante, soit rejeter cette position et faire appel à la centrale inertielle de35 l'avion pour déterminer sa position. il Dans les applications o l'intégration sur de longues durées est acceptable, l'invention permet de raccourcir de manière significative les temps de calcul, ce qui est également très avantageux. Le procédé selon l'invention peut être implanté sans problème sur une carte à l'intérieur d'un récepteur GPS Les URE sont prélevés très facilement dans le message fourni par chaque satellite de la constellation choisie par le récepteur, et le calcul de E est rapide puisqu'il ne10 comprend qu'une inversion de matrice carrée d'ordre égal au maximum à 9 ou 10 et des produits de matrices. Bien entendu, l'invention n'est pas limitée au mode de réalisation qui vient d'être décrit et l'on pourra remplacer tout moyen par un moyen équivalent sans sortir du cadre de15 l'invention. REVENDICATIONS 1/ Procédé de calcul, par un récepteur GPS porté par un mobile, de la position spatiale et temporelle dudit mobile dans un référentiel prédéterminé, les coordonnées de ladite position étant notées x, y, z et t, o x, y et z sont respectivement l'abscisse, l'ordonnée et la cote dudit mobile dans ledit référentiel, et t la date, dans le système de temps universel, à laquelle ledit mobile se trouve dans10 la position (x, y, z), ladite position étant calculée à partir de données fournies par au moins quatre satellites du système GPS sélectionnés par ledit récepteur GPS pour former une constellation dite de calcul, ledit procédé étant caractérisé en ce qu'il comporte une15 étape d'estimation de l'erreur commise dans ledit calcul de position, ladite estimation de l'erreur étant représentée par un vecteur noté a dont les composantes sont respectivement l'estimée de l'erreur sur l'abscisse dx, l'estimée de l'erreur sur l'ordonnée dy, l'estimée de20 l'erreur sur la cote dz et l'estimée de l'erreur temporelle c.dt, c étant la célérité de la lumière dans le vide, et en ce que ladite estimation est effectuée par ledit récepteur GPS pour chaque position calculée, de sorte que ladite estimation de l'erreur commise varie avec la position dudit25 mobile. 2/ Procédé selon la revendication 1 caractérisé en ce que ladite constellation comporte au plus tous les satellites du système GPS visibles depuis ledit mobile. 3/ Procédé selon la revendication 2 caractérisé en ce que e est donné par la relation: v er = K e o: K est une matrice à N lignes et quatre colonnes, N étant le nombre de satellites que comporte ladite constellation, telle que: a 11 a 12 a 13 1 a 21 a 22 a 23 1 K = a N 1 a N 2 a N 3 1 aij, pour i variant de 1 à N et j variant de 1 à 3, étant le cosinus directeur de l'angle entre le j-ième axe de coordonnées et la droite joignant ledit mobile au i-ième 5 satellite de ladite constellation, les indices 1, 2 et 3 étant respectivement associés aux axes des abscisses, des ordonnées et des cotes, 6 r est un vecteur à N composantes notées URE 1 à UREN, o UR Ei, pour i variant de 1 à N, est une donnée fournie par ledit i-ième satellite et donnant l'écart radial moyen entre la position théorique du i-ième satellite et sa position réelle sur la droite joignant ledit i-ième satellite à un point de la zone géographique de la Terre vue par ledit i-ième satellite, la moyenne étant effectuée sur les points15 de ladite zone géographique o se trouve une station de relais au sol du système GPS. 4/ Procédé selon la revendication 3 caractérisé en ce que, lorsque N est strictement supérieur à 4, la meilleure estimée E de 8 est donnée par la formule suivante:20 E = (KT K)-i KT 8 r, et l'on utilise E à la place de 8 en remplaçant les composantes dx, dy, dz et c dt de e respectivement par les composantes notées Dx, Dy, Dz et c Dt de E. 5/ Procédé selon l'une des revendications 1 à 4 caractérisé en ce qu'il comporte en outre les opérations suivantes: une valeur P, appelée précision, est calculée à partir des composantes de a, ladite précision étant donnée par la formule: P = -dx 2 + dy 2 + dz 2 + (c dt)2, P est comparée à un seuil de tolérance préalablement fixé, ledit récepteur GPS sélectionne une autre constellation si P dépasse ledit seuil. 6/ Procédé selon l'une des revendications 1 à 4 caractérisé en ce qu'il comporte en outre les opérations suivantes: une valeur P, appelée précision, est calculée à partir des composantes de a, ladite précision étant donnée par la formule: p = x Jdx 2 + dy 2 + dz 2 + (c dt)2, P est comparée à un seuil de tolérance préalablement fixé, ledit récepteur GPS indique que ladite position calculée est à rejeter.