FR 3013851 A1 20150529 FR 1302697 A 20131122 Système et procédé de détermination de l'erreur de position d'un récepteur de localisation satellitaire L'invention concerne le domaine de la localisation satellitaire. En particulier il concerne la détermination d'une erreur de positionnement que réalise un récepteur en utilisant un système de localisation satellitaire. Ces systèmes sont par exemple implémentés dans les systèmes d'amélioration de la précision de la localisation fournie par le système GPS (acronyme anglais de « Global Positionning System » ou système de localisation mondial en français). Ces systèmes sont aussi nommés systèmes de renforcement de la localisation. Ils sont aussi connus sous l'acronyme anglais de SBAS pour « Satellite-Based Augmentation Systems ». Il est également connu un système connu sous l'acronyme anglais EGNOS pour « European Geostationary Navigation Overlay Service » ou Service Européen de Navigation par Surcouche Géostationnaire comme système de renforcement du GPS. La connaissance de cette erreur de positionnement du récepteur, permet de déterminer un volume pour lequel la probabilité de présence du récepteur est supérieure à un seuil fixé par un standard. Ainsi la connaissance de ce volume permet par exemple de déterminer la distance minimale à laquelle deux aéronefs ne doivent pas se rapprocher. La connaissance de cette information est en particulier nécessaire pour les services d'intégrité. On appelle service d'intégrité la capacité d'un système à fournir une alarme au pilote lorsque le système de navigation ne peut plus être utilisé avec la performance requise en termes de risques pour l'usager. Ces systèmes de détermination d'erreur sont utilisables par 30 exemple dans des aéronefs mais également dans des véhicules terrestres ou des navires. Il est connu dans l'état de la technique des systèmes dans lesquelles des informations représentatives des positions des satellites, et 35 des temps de traversée de la ionosphère ainsi que des informations représentatives de l'erreur sur ces positions et ces temps de traversée sont envoyés aux différents récepteurs. La connaissance de ces informations permet de déterminer l'erreur de position du récepteur, on parle aussi d'intégrité de la position de l'objet localisé par le récepteur. Il est connu dans l'état de la technique que les informations représentatives des erreurs soient des écarts types margés (ai) des distributions des erreurs commises. Ainsi la distribution des erreurs est modélisée par une loi gaussienne centrée de la forme N(0, o-i2). Cependant la modélisation de la distribution des erreurs sous forme d'une loi gaussienne centrée est trop grossière, et amène à devoir prendre une marge de sureté, qui dans certains cas est trop importante. C'est en particulier le cas lorsque la combinaison des écarts types margés (6i), de l'erreur de position de chaque satellite utilisé, doit être faible pour permettre une précision suffisante à la réalisation des manoeuvres à exécuter. Dans de tels cas, la marge, prise pour couvrir la non adaptation d'un modèle basé sur des gaussiennes centrées, rend le service trop souvent indisponible. On a de plus pu montrer que l'utilisation de gaussiennes centrées n'est correcte du point de vue mathématique, que si l'on peut faire l'hypothèse que les erreurs sont distribuées de façon uni-modale et symétrique, que rien ne garantit de façon générale. La présente invention vise à remédier à ces problèmes en proposant un système de détermination de l'erreur de position d'un récepteur dont la précision est accrue par rapport à la précision des systèmes basés sur l'échange des écarts types margés de la distribution de l'erreur sur la position des satellites . Il est proposé, selon un aspect de l'invention un système de détermination d'une distribution d'une erreur de position d'un récepteur de signaux de localisation. Les signaux étant envoyés par au moins un satellite. Le système comporte ; - le récepteur dont une position est dite première position et est affectée par une erreur, dite première erreur, présentant une distribution, dite première distribution; - un premier dispositif de détermination d'au moins une positions du ou de chaque ou d'au moins un satellite, dite deuxième position, - un dispositif de transmission de la ou de chaque ou d'au moins une deuxième position du premier dispositif de détermination vers le récepteur. Le système est caractérisé en ce que - la première distribution est définie par au moins un premier cumulant, d'ordre supérieur à 2, - le premier dispositif est en outre adapté pour la détermination d'au moins un deuxième cumulant, d'ordre supérieur à 2, représentatif d'une deuxième distribution d'une deuxième erreur sur la deuxième position, - le dispositif de transmission est en outre adapté pour la transmission du ou de chaque ou d'au moins un deuxième cumulant du premier dispositif de détermination vers le récepteur, - et en ce que le récepteur comporte: o un deuxième dispositif de détermination du ou de chaque ou d'au moins un premier cumulant, à partir de la deuxième position, du deuxième cumulant et d'un modèle de détermination de la première position du récepteur à partir de distances entre le récepteur et le ou les satellites et o un troisième dispositif de détermination de la première distribution, à partir dudit premier cumulant. Les premiers et deuxièmes cumulants d'ordre supérieur à deux 25 peuvent aussi être nommés sous l'expression jeu de cumulants ou ensemble de cumulants. Les cumulants Kn de la variable aléatoire X sont définis par la fonction génératrice des cumulants g(t) : 30 CO (el g (t) = ln((eiXt) x) = 1K n! n n=0 Dans cette équation ( )x représente l'espérance mathématique de la variable aléatoire X et Kn le cumulant d'ordre n. Le système permet donc la transmission aux récepteurs de signaux satellitaires d'une information représentative de la distribution de l'erreur de position des satellites. De plus l'utilisation de cette information permet de déterminer l'erreur de positionnement que réalise le récepteur en utilisant le système de localisation satellitaire De plus, l'utilisation des cumulants, permet de transmettre aux usagers une information représentative de la distribution de l'erreur de la position du satellite qui est plus fidèle que dans le cas des systèmes connus dans l'état de la technique. L'utilisation de ces cumulants évite alors d'avoir à prendre autant de marges dans la position et ce pour deux raisons : - Elle permet de décrire plus finement les distributions d'erreur de positionnement/synchronisation des sources de signal de positionnement, et de ce fait permet d'améliorer la précision de l'information de localisation calculée par le récepteur. La zone de tolérance au risque d'intégrité peut être plus petite que dans les solutions actuelles. - Elle permet de ne pas avoir à faire une hypothèse dont la base mathématique est fausse, et qui consiste modéliser les erreurs commises pour la position de chaque satellite par une combinaison de gaussiennes centrées N(0,0-12) alors que la distribution n'est pas unimodale ou symétrique. En d'autres termes, la forme de la distribution permise par l'invention est beaucoup plus précise qu'une distribution Gaussienne centrée. Cependant cette forme générique peut être décrite par quelques paramètres seulement, les cumulants. On est donc en mesure de transmettre une information beaucoup plus fine sur la distribution de probabilité des deuxièmes erreurs de positionnement des satellites. Cette information précise permet alors de modéliser les premières erreurs de position du récepteur de façon aussi fidèle que possible. Ainsi, il 35 n'est pas nécessaire de prendre des marges de sécurité importantes pour assurer la sureté, et il est possible d'utiliser le système de localisation pour la réalisation de manoeuvres dans lesquelles l'exigence de précision est importante. Selon une caractéristique technique le deuxième dispositif de détermination est adapté pour appliquer le modèle au ou aux deuxièmes cumulants. Selon une caractéristique technique le premier dispositif de détermination est en outre adapté pour déterminer les deuxièmes cumulants d'ordre 1 à 5, et 10 dans lequel le deuxième dispositif de détermination est en outre adapté pour déterminer les premiers cumulants d'ordre 1 à 5. Plus il y a de cumulants plus la modélisation va être précise, mais ils sont plus compliqués à mesurer expérimentalement. 15 Selon une caractéristique technique le troisième dispositif sont adaptés en outre pour utiliser un développement d'Edgeworth. L'invention concerne aussi un procédé de détermination d'une distribution 20 d'une erreur de position d'un récepteur de signaux de localisation, les signaux étant envoyés par au moins un satellite. Le procédé comporte ; - une étape de réception, par un récepteur, des signaux de localisation satellitaire, une position du récepteur étant dite première position et est affectée par une erreur, dite première erreur, présentant une distribution, 25 dite première distribution, - une première étape de détermination, par un premier dispositif de détermination, d'au moins une position du ou de chaque ou d'au moins un satellite, dites deuxièmes positions, - une étape de transmission, par un dispositif de transmission, de la ou de 30 chaque ou d'au moins une deuxième position du premier dispositif de détermination vers le récepteur, Le procédé est caractérisé en ce que - la première distribution est définie par au moins un premier cumulant, d'ordre supérieur à 2, - la première étape de détermination est en outre adaptée pour la détermination d'au moins un deuxième cumulant, d'ordre supérieur à 2, représentatif d'une deuxième distribution d'une deuxième erreur sur la deuxième position, - l'étape de transmission est en outre adaptée pour la transmission du ou de chaque ou d'au moins un deuxième cumulant, associé à ladite deuxième position du premier dispositif de détermination vers le récepteur, - de plus le procédé comporte : o une deuxième étape de détermination, par un deuxième dispositif de détermination dudit récepteur, du ou de chaque ou d'au moins un premier cumulant, à partir de la deuxième position, du deuxième cumulant et d'un modèle de détermination de la première position du récepteur à partir de distances entre le récepteur et le ou les satellites et o une troisième étape de détermination par un troisième dispositif de détermination du récepteur, de la première distribution, à partir du premier cumulant Selon une caractéristique technique la deuxième étape de détermination est 20 adaptée pour appliquer le modèle au ou aux deuxièmes cumulants. Selon une caractéristique technique la première étape de détermination de détermination est en outre adaptée pour déterminer les deuxièmes cumulants d'ordre 1 à 5, et dans lequel la deuxième étape de détermination 25 est en outre adaptée pour déterminer les premiers cumulants d'ordre 1 à 5. Selon une caractéristique technique la troisième étape de détermination est adaptée en outre pour utiliser un développement d'Edgeworth. 30 L'invention sera mieux comprise et d'autres avantages apparaîtront à la lecture de la description détaillée, faite à titre d'exemple non limitatif. Cette description détaillée est réalisée à l'aide des figures suivantes : la figure 1 présente un premier mode de réalisation du système présenté dans cette invention 35 3013 851 7 La figure 1 présente le système comportant un satellite 101 et un récepteur de signaux satellitaires 102. Le système permet au récepteur de déterminer la distribution d'une première erreur associée à une première position du récepteur. Cette première distribution de la première erreur est 5 modélisée par au moins un premier ensemble de cumulants d'ordre supérieur à deux. Le système permet, afin de réaliser cette détermination, la transmission (via un dispositif de transmission 103) de deuxièmes cumulants 10 qui représentent une deuxième distribution représentant une deuxième erreur associée à la deuxième position d'un satellite. La détermination de ces éléments est réalisée par un premier dispositif de détermination 104. Les premiers cumulants sont déterminés par le récepteur par 15 l'utilisation d'un deuxième dispositif de détermination 105. Enfin un troisième dispositif 106 permet la détermination de la première distribution, à partir des premiers cumulants Cette modélisation est basée sur l'utilisation du développement d'Edgeworth de la densité de probabilité de l'erreur associée à la position d'un satellite. Les cumulants d'une variable aléatoire X distribuée selon une densité de probabilité f (on note X f), sont déterminés en introduisant la fonction cp(t) = (eixt)x. - e représente la fonction exponentielle - ( )x représente la moyenne sur les valeurs de X - i étant l'unité imaginaire (i2 = -1). On remarque que le développement de cette fonction, en fonction des puissances de l'exposant, est une série qui fait intervenir les moments d'ordre n de f : pin = (Xn).35 CO e (p, (t) = 1 (ln!l pin n=0 Il est aussi possible d'effectuer le développement de la fonction ln((eixt)x), on obtient alors un ensemble de coefficients Kn, que l'on définit de la manière suivante : CO 111((eiXt)X) = 1 (intr Kn n=0 Chaque Kn ainsi défini est le cumulant d'ordre n de la distribution f. Les deux premiers cumulants sont la moyenne et la variance de la distribution. De plus si X et Y sont deux variables aléatoires distribuées selon f et g respectivement, et dont les cumulants d'ordre n sont respectivement Kn [f] et Kit [g], alors les cumulants d'ordre n de la distribution h associée à la variable aléatoire Z = pX + qY , sont donnés par: Kn[h] = pnKn[f] + qnKn[g] De plus on sait que toute distribution, qui résulte de la combinaison de m variables aléatoires peut-être représentée par un développement, connu sous le nom de développement d'Edgeworth et ayant la forme suivante : Fn(x) = 1 +I-- P (Ki . - - Ki lif (X) [1 CO = 12 1 i 1 j -, ) , P X) Dans cette équation les variables sont les suivantes : - P(x) est une fonction de référence au choix (par exemple Gaussienne) - Kp est le cumulant d'ordre p de la distribution des erreurs d'orbite et/ou d'horloge - Pi est un polynôme d'ordre 3j en x, qui fait intervenir les j premiers Kp dans ses coefficients, et dont l'expression dépend du choix de (x) - n représente le nombre de variables combinées pour obtenir x De plus on sait que ce développement converge lorsque n tend vers 00. En s'appuyant sur les concepts mathématiques ci-dessus, l'invention propose la détermination de la première erreur de position du récepteur de la manière suivante : - Un dispositif de transmission 103 fournit les informations sur la distribution des erreurs de position et de synchronisation des satellites sous la forme de cumulants d'ordre supérieur à deux de cette distribution. Cette transmission est réalisée pour chacune des NS sources de signal de positionnement (par exemple des satellites émettant un signal respectant la norme GPS) - Le récepteur détermine sa position et un temps de référence en utilisant une combinaison linéaire des mesures de (pseudo-)distances pi faites entre son antenne et les NS sources de signal utilisés pour le positionnement. - Le récepteur détermine les m premiers cumulants (Kn) de la première distribution de l'erreur associée à sa position, à partir des cumulants Kij transmis, en utilisant la relation suivante: NS Kn,p = (111pi r Kn j=1 Avec n = 1, --- ,m représentant l'ordre du cumulant, j le satellite, j le coefficient p,j de la matrice qui permet de déterminer la deuxième position du récepteur à partir des distances entre le récepteur et les satellites, p représente la direction (x, y ou z) pour laquelle le cumulant est déterminée. Dans un mode de réalisation il est possible d'utiliser la méthode des moindres carrés pour déterminer la matrice Mni. Dans ce mode de réalisation le vecteur des distances entre le récepteur et les satellites se modélise de la manière suivante suivante : p= ÛX+ E Dans cette équation p = [Pi PNS] est le vecteur des distances entre le récepteur et les satellites, E = [El ENS] est le vecteur des erreurs des distances entre le récepteur et les satellites et X = [x,y,z, Atusr] est le vecteur de la deuxième position et du décalage d'horloge du récepteur et =F exei eyei ezei c exeNs e3,eNs ezeNsc 10 est la matrice d'observation du problème. Dans la matrice Ô, c représente la vitesse de la lumière et exej est le cosinus de l'angle entre le vecteur dans la direction x et le vecteur vers le satellite j. En utilisant la méthode des moindres carrés, la relation entre la deuxième position du récepteur et les distances entre le récepteur et les satellites peut 15 s'écrire Xest = cte-ip. Dans cette relation C = est la matrice de corrélation des erreurs. On peut alors déterminer dans ce mode de réalisation que M = (cte-i c)-i cte-1. - En utilisant le développement d'Edgeworth tronqué à l'ordre m, qui fait 20 intervenir les premiers cumulants Kn, on détermine la première distribution F qui approche la distribution de la première erreur associée au positionnement du récepteur. Enfin cette distribution de la première erreur peut-être utilisée pour 25 déterminer la dimension d'une zone de confiance (c'est-à-dire une zone dans laquelle la probabilité de trouver le récepteur est supérieure ou égale à un seuil déterminé). Il est possible de trouver cette zone de confiance est résolvant l'équation suivante. -Rp 00 fF(x)dx f F(x)dx = PHAil -00 Rp Cette détermination doit être effectuée pour chaque direction d'espace (verticale, horizontale). P - HMI représente la probabilité tolérée de non-intégrité, ceci pour s'assurer que ces Rp sont inférieurs aux dimensions de la zone de tolérance (par exemple les rayons d'alerte utilisés dans la navigation civile). Il aussi possible de trouver directement le risque de se trouver à l'extérieur de la zone de tolérance requise (Ra), pour ceci on peut utiliser l'équation suivante : -Ra co fF(x)dx f F(x)dx La diffusion des cumulants vers les récepteurs est réalisé avec un dispositif d'alerte et/ou par la remise à jour des valeurs des cumulants au cas où elles se révéleraient être inadaptées aux besoins d'intégrité suite à un changement de l'état du système. La distribution des premières erreurs est réalisé en utilisant les cumulants diffusés et modélisant la fonction de référence IF par une gaussienne centrée sur le premier cumulant d'ordre 1 et de largeur le premier cumulant d'ordre 2 Ensuite on peut évaluer la disponibilité du service et donc la zone 10 dans laquelle probabilité de présence du récepteur dépasse un seuil en utilisant les équations précédentes. Le premier dispositif de détermination 104 des positions du ou des satellites peut être localisé au sol ou dans un des satellites. 15 Les différents dispositifs de détermination décrits dans cette invention peuvent être des ordinateurs ou des processeurs programmés de manière à réaliser les différentes opérations réalisées par les dispositifs. Il est aussi possible d'utiliser des composants dédiés, des circuits de logique 20 programmable, des réseaux logiques programmables (connus aussi sous l'acronyme anglais FPGA pour Field-Programmable Gate Array) ou des circuits intégrés propre à une application (connus aussi sous l'acronyme anglais ASIC pour Application-Specific Integrated Circuit) programmés de manière à réaliser les différentes opérations réalisées par les dispositifs. 25 La présente invention peut également s'implémenter à partir d'éléments matériels et logiciels. Elle peut être disponible en tant que produit programme d'ordinateur sur un support lisible par ordinateur. Le support peut être électronique, magnétique, optique, électromagnétique ou être un support 30 de diffusion de type infrarouge. De tels supports sont par exemple, des mémoires à semi-conducteur (Random Access Memory RAM, Read-Only Memory ROM), des bandes, des disquettes ou disques magnétiques ou optiques (Compact Disk - Read Only Memory (CD-ROM), Compact Disk - Read/Write (CD-R/V1/) and DVD). 35 REVENDICATIONS1. Système de détermination d'une distribution d'une erreur de position d'un récepteur (102) de signaux de localisation, lesdits signaux étant envoyés par au moins un satellite, ledit système comportant ; - ledit récepteur (102) dont une position est dite première position et est affectée par une erreur, dite première erreur, présentant une distribution, dite première distribution; - un premier dispositif de détermination (104) d'au moins une position dudit ou d'au moins un dit ou de chaque dit satellite, dite deuxième position, - un dispositif de transmission (103) de ladite ou d'au moins une dite ou de chaque dite deuxième position dudit premier dispositif de détermination vers ledit récepteur, ledit système étant caractérisé en ce que - ladite première distribution est définie par au moins un premier cumulant, d'ordre supérieur à 2, - ledit premier dispositif est en outre adapté pour la détermination d'au moins un deuxième cumulant, d'ordre supérieur à 2, représentatif d'une deuxième distribution d'une deuxième erreur sur ladite deuxième position, - ledit dispositif de transmission (103) est en outre adapté pour la transmission dudit ou d'au moins un dit ou de chaque dit deuxième cumulant dudit premier dispositif de détermination vers ledit récepteur, - et en ce que ledit récepteur comporte: o un deuxième dispositif de détermination (105) dudit ou d'au moins un dit ou de chaque dit premier cumulant, à partir de ladite deuxième position, dudit deuxième cumulant et d'un modèle de détermination de ladite première position du récepteur à partir de distances entre ledit récepteur et ledit ou lesdits satellites et o un troisième dispositif de détermination (106) de ladite première distribution, à partir dudit ou d'au moins un dit ou de chaque dit premier cumulant. 2. Système selon la revendication 1 dans lequel ledit deuxième dispositif de détermination (105) est adapté pour appliquer ledit modèle audit ou auxdits deuxièmes cumulants. 3. Système selon la revendication 1 ou 2 dans lequel ledit premier dispositif de détermination est en outre adapté pour déterminer lesdits deuxièmes cumulants d'ordre 1 à 5, et dans lequel ledit deuxième dispositif de détermination est en outre adapté pour déterminer lesdits premiers cumulants d'ordre 1 à 5. 4. Système selon l'une des revendications 1 à 3 dans lequel ledit troisième dispositif est adapté en outre pour utiliser un développement d'Edgeworth. 5. Procédé de détermination d'une distribution d'une erreur de position d'un récepteur (102) de signaux de localisation, lesdits signaux étant envoyés par au moins un satellite, ledit procédé comportant ; - une étape de réception, par un récepteur (102), desdits signaux de localisation satellitaire, une position dudit récepteur étant dite première position et étant affectée par une erreur, dite première erreur, présentant une distribution, dite première distribution, - une première étape de détermination, par un premier dispositif de détermination (104), d'au moins une position dudit ou d'au moins un dit ou de chaque dit satellite, dite deuxième position, - une étape de transmission, par un dispositif de transmission (103), de ladite ou d'au moins une dite ou de chaque dite deuxième position dudit premier dispositif de détermination vers ledit récepteur, ledit procédé étant caractérisé en ce que - ladite première distribution est définie par au moins un premier cumulant, d'ordre supérieur à 2, - ladite première étape de détermination est en outre adaptée pour la détermination d'au moins un deuxième cumulant, d'ordre supérieur à 2, représentatif d'une deuxième distribution d'une deuxième erreur sur ladite deuxième position,- ladite étape de transmission est en outre adaptée pour la transmission dudit ou d'au moins un dit ou de chaque dit deuxième cumulant, associé à ladite deuxième position dudit premier dispositif de détermination vers ledit récepteur, - de plus ledit procédé comporte : o une deuxième étape de détermination, par un deuxième dispositif de détermination dudit récepteur, dudit ou d'au moins un dit ou de chaque dit premier cumulant, à partir de ladite deuxième position, dudit deuxième cumulant et d'un modèle de détermination de ladite première position du récepteur à partir de distances entre ledit récepteur et ledit ou lesdits satellites et o une troisième étape de détermination par un troisième dispositif de détermination dudit récepteur, de ladite première distribution, à partir dudit ou d'au moins un dit ou de chaque dit premier cumulant. 6. Procédé selon la revendication 5 dans lequel ladite deuxième étape de détermination est adaptée pour appliquer ledit modèle audit ou auxdits deuxièmes cumulants. 7. Procédé selon la revendication 5 ou 6 dans lequel ladite première étape de détermination de détermination est en outre adaptée pour déterminer lesdits deuxièmes cumulants d'ordre 1 à 5, et dans lequel ladite deuxième étape de détermination est en outre adaptée pour déterminer lesdits premiers cumulants d'ordre 1 à 5. 8. Procédé selon l'une des revendications 5 à 7 dans lequel ladite troisième étape de détermination est adaptée en outre pour utiliser un développement d'Edgeworth.30