FR 3038067 A1 20161230 FR 1555805 A 20150624 1 Procédé de localisation d'un récepteur au sein d'un système de positionnement La présente invention a pour objet un procédé de localisation d'au moins un récepteur au sein d'un système de positionnement ainsi qu'un tel système de positionnement. L'invention s'applique plus particulièrement à des systèmes utilisant le multiplexage par le code, encore appelés en anglais « code division multiple access » (CDMA), ce qui est par exemple le cas du GPS (« global positioning system » en anglais) et du GNSS (« global navigation satellite system » en anglais). Ces systèmes mettent en oeuvre des émetteurs émettant un signal comprenant une porteuse modulée par un code. Les systèmes de positionnement connus reposent, pour la très grande majorité d'entre eux, sur une approche de triangulation ou de multilatération afin de définir, dans un repère global ou local, les coordonnées d'un terminal mobile, appelé récepteur. D'autres approches reposent quant à elle sur un ensemble de mesures en coordonnées polaires afin d'obtenir les coordonnées du terminal cible. Ces approches sont mises en oeuvre au niveau du terminal, ou à un niveau plus général, tel qu'un serveur distant avec lequel communique le terminal, ou à la fois sur le terminal et sur le serveur. Les mesures conduisent à des distances ou à des angles. Les systèmes de positionnement connus sont efficaces en milieu extérieur, à ciel ouvert, milieu dit « outdoor » en anglais, mais leurs performances se dégradent en milieu intérieur, dit « indoor » en anglais ou lorsque le ciel est obstrué par des immeubles ou obstacles hauts. De surcroît, le niveau de précision requis en intérieur est supérieur à celui requis en extérieur, l'environnement étant plus petit et plus dense en informations en intérieur. Il est en outre important d'assurer la continuité de positionnement entre les milieux extérieur et intérieur. Les réseaux de capteurs offrent des performances excellentes mais leur infrastructure présente une très forte contrainte. Les capteurs employés sont les capteurs à ultrasons, infrarouge, de pression ou encore des étiquettes électroniques. Ces dernières, notamment du type « Near Field Communication » (NFC), sont de plus en plus utilisées, mais ne permettent pas un positionnement continu, ni dans le temps ni dans l'espace. Les techniques utilisant les réseaux de télécommunication mobiles, tels que GSM, UMTS, ou de quatrième génération dits « 4G », ou locaux, tels que Bluetooth, 3038067 2 Bluetooth Low Energy (BLE), Wi-Fi ou Ultra Large Band (ULB), nécessitent une infrastructure pouvant être déployée par ailleurs, pour des besoins autres que la localisation. Il est cependant nécessaire de mettre en oeuvre également des moyens spécifiques à la localisation. Une approche fondée sur l'identification de cellule dans le cas 5 de réseaux mobiles, appelée « Cell-Id » en anglais, a été développée pour les balises locales, évitant le recours aux diverses phases de calibration souvent requises. Cette technique est peu précise. Des approches dites de «positionnement symbolique » ont été développées, comme décrit dans l'article de Magda Chelly et Nel Samama, « WiFi Indoor Symbolic 10 Positioning System Implemented On a PDA », GNSS 2007, fondées sur une solution hybride utilisant l'identification de cellule dans les réseaux mobiles et les signaux GNSS, et proposant un positionnement continu, peu précis et fiable à ce jour. Des systèmes inertiels sont également utilisés, mettant en oeuvre une infrastructure réduite mais un terminal mobile complexe et une modélisation non triviale. 15 Ces capteurs sont intégrés dans certains téléphones intelligents actuels, « smartphones » en anglais. Les approches utilisant des capteurs détournent ces derniers de leur utilisation naturelle pour réaliser une fonction de localisation. Les systèmes utilisant l'analyse des variations locales du champ magnétique 20 permettent de bonnes performances, mais ne sont pas encore prêts techniquement pour être intégrés dans des smartphones. Des approches similaires reposent sur l'utilisation de balises magnétiques, et se rapprochent des techniques de type « empreintes » ou « fingerprinting» en anglais, couramment utilisées par les réseaux locaux de télécommunications. 25 Des systèmes optiques sont également connus, utilisant des caméras intégrées aux smartphones, mais conduisant à une complexité de calcul importante pour les processeurs de ces derniers. Les systèmes de reconnaissance d'images sont performants mais n'assurent pas une totale continuité : des contraintes de mise en oeuvre persistent, ces systèmes requérant une action consciente de l'utilisateur pour effectuer la prise de vue, et 30 nécessitant une image non floue, ce qui limite le mouvement du terminal mobile durant la prise de vue. 3038067 3 La technique dite «Li-Fi», pour «Light Fidelity» en anglais, émerge également. Cette approche s'apparente à la technique « Cell-Id » mentionnée plus haut, mais utilise les détecteurs de luminosité des smartphones, codant chaque luminaire en fonction de sa position par l'intermédiaire d'une modulation de la lumière, le smartphone 5 décodant cette modulation pour déterminer sa position. L'ergonomie d'usage d'une telle approche n'est pas optimale. Les systèmes de type radar reposent sur des mesures de temps de parcours aller-retour d'une onde, permettant d'obtenir une distance, et sur des mesures de direction d'arrivée en coordonnées sphériques. Ces systèmes représentent une simplification plutôt 10 élevée des contraintes habituelles d'un système de positionnement, à savoir la synchronisation entre les parties en présence et la diversité spatiale nécessaire. Un unique émetteur est suffisant pour déterminer une position. Néanmoins des terminaux mobiles dédiés sont nécessaires. Certaines approches utilisent des antennes doubles espacées d'une longueur d'onde du signal, et mesurent la différence de phases des signaux reçus des deux 15 antennes pour définir l'angle d'arrivée du signal émis. Des solutions utilisant les signaux GNSS ont été développées, notamment AGPS pour « Assisted GPS» en anglais, et HSGNSS pour « High Sensibility GNSS» en anglais, offrant une amélioration des performances en positionnement dans des environnements auparavant non couverts par des satellites, et pouvant être mises en oeuvre 20 sans infrastructure dédiée. Ces solutions sont néanmoins limitées en termes de détection, de précision et de temps nécessaire à l'obtention d'une position. Il est connu du brevet US 5 899 957 de corriger des signaux GPS porteurs de phase à une position arbitraire en utilisant des mesures de phase reçues de différentes stations GPS de base dont la position est connue, pour améliorer l'efficacité du 25 positionnement. Le brevet US 4 646 096 décrit un procédé utilisant le principe du mouvement de la constellation de satellites autour d'un récepteur en mesurant et stockant les variations de phase des porteuses correspondant aux signaux provenant desdits satellites, afin de déterminer l'ambiguïté entière des mesures du code des signaux correspondants pour des 30 récepteurs monocanal. Il est connu, notamment du brevet US 7 212 155, d'utiliser des mesures de variation de la phase de porteuses de signaux provenant de satellites, sans utiliser de station 3038067 4 de base, pour corriger progressivement une position hypothétique de départ d'un récepteur en mouvement. Il est connu des brevets US 7 671 794 et US 6 549 165 des systèmes de positionnement utilisant des mesures de différence de phase réalisées soit sur deux 5 antennes distinctes en réception, soit entre deux signaux émis à des fréquences différentes. La demande WO 2014/060777 décrit le positionnement d'un terminal mobile, utilisant un signal envoyé par ledit terminal et reçu par des antennes multiples de points d'accès multiples sans fil, et se fondant sur des calculs de différences de phase entre les différentes antennes pour obtenir un angle d'arrivée du signal sur les antennes. 10 La plupart des méthodes précédemment décrites nécessitent de coupler ou fusionner plusieurs techniques et/ou technologies afin de parvenir à la continuité de positionnement entre les milieux extérieur et intérieur. Les approches connues à base de répélites mettent à profit les récepteurs et les signaux GNSS transmis, afin de permettre à un unique récepteur GNSS, présent dans le 15 terminal mobile, d'assurer une fonction de localisation. Ces approches utilisent le concept de pseudo-satellites, ou « pseudolites », c'est-à-dire un réseau d'émetteurs terrestres transmettant des signaux de même structure que les signaux envoyés par un satellite, comme décrit dans l'article de Kee C. et al. "Centimeter-Accuracy Indoor Navigation Using GPS-Like Pseudolites" , GPS World 2001, 20 pages 14-20. Les mesures les plus précises des signaux GNSS, qui permettent en particulier un positionnement centimétrique dans de bonnes conditions de réception en extérieur, sont utilisées dans ces approches, offrant des performances expérimentales de grande qualité, et étant fiables et simples à utiliser. Elles permettront à court terme de réduire le nombre d'émetteurs nécessaires et de simplifier la gestion du système, 25 notamment en supprimant la synchronisation ou en utilisant des émetteurs bas coût. Les signaux élaborés dans le cadre du projet européen de système de positionnement par satellites « Galileo » permettent une réduction importante des niveaux d'interférence, laissant à penser qu'il sera possible de conserver les bandes utilisées par les signaux GNSS pour assurer la continuité complète de la fonction de localisation. 30 La géolocalisation utilisant des pseudolites calcule des intersections géométriques de sphères centrées en des points connus, les positions des pseudolites, en utilisant le temps de propagation des signaux pour mesurer les rayons des sphères, comme 3038067 5 décrit dans le chapitre 2 du livre de Kaplan E. et Hegarty C. "Understanding GPS Principles and Applications", Artech House, 2006, 2nde édition. Le temps de propagation d'un signal peut être mesuré grâce au code ou à la porteuse du signal. Cependant, le phénomène de chemins multiples, dit « multipath » en anglais, affecte fortement les temps 5 de mesure en extérieur, et encore plus en intérieur, à cause de la présence de nombreux obstacles, comme décrit dans le chapitre 6 du livre susmentionné et dans l'article de Fluerasu A. et al. "Multipath modelisation of typical indoor environments optimisation of GNSS based indoor positioning", ENC-GNSS 2008. Il a été démontré qu'il est difficile d'obtenir des mesures précises à plus de 2 ou 3 mètres en utilisant uniquement des mesures 10 de code, dans les articles de Jardak N. et Samama N. "Short Multipath Insensitive Code Loop Discriminator", IEEE Transactions on Aerospace and Electronic Systems 2010, vol 46, pages 278-295, de Vervisch-Picois A. et al. "2D Indoor Dynamic Positioning Using GNSS Based Repeaters", ION GNSS 2006, et de Jee G. I. et al. " Indoor Positioning Using Time Synchronised Switching GPS Repeater", ION GNSS 2005. 15 Les mesures de la phase de la porteuse ont l'avantage d'être moins sensibles au problème de chemins multiples, et permettent une précision de quelques dizaines de centimètres, comme présenté dans les articles de Rizos C. et al. "LocataNet: Intelligent Time-Synchronised Pseudolite Transceivers for cm-Level Stand-Alone Positioning", TAIN World Congress 2003, et de I. Selmi et al. "Experimental Positioning Results of the 20 Repealite Based Indoor Positioning System", IEEE International Conference on Indoor Positioning and Indoor Navigation 2012. Cependant, il est difficile de réaliser un positionnement absolu avant de connaitre l'ambiguïté de la phase de la porteuse issue de chaque pseudolite. Cette ambiguïté Nk correspond au résultat de la division euclidienne de la distance de 25 propagation par la longueur d'onde de la porteuse : DPLk = Ok ).Nk, avec DPLK la distance entre le récepteur et le pseudolite PLk, Ok la valeur de la phase de la porteuse, et X, la longueur d'onde, par exemple égale à 0,19 m pour le signal Ll GPS. Pour résoudre ce problème en intérieur, des méthodes utilisant les principes liés au milieu extérieur sont connues, décrites dans les articles de Kee C. et al. "Development of 30 indoor Navigation System using asynchronous pseudolites" , ION GPS 2000, de Xiaoguang Wan et al. "The Pseudolite-based Indoor Navigation System Using Ambiguity Resolution On The Fly" , IEEE 3rd International Symposium on Systems and Control in Aeronautic and 3038067 6 Astronautics 2010, pages 212-217, et de Puengnim A. et al. "Precise Positioning for Virtually Synchronized Pseudolite System", IEEE International Conference on Indoor Positioning and Indoor Navigation 2013, pages 1-8. Néanmoins, l'hypothèse de linéarité du milieu extérieur n'étant pas valable en intérieur, la méthode décrite par Kee et al. repose 5 sur une version non linéaire de l'algorithme de positionnement décrit dans le chapitre 2 du livre susmentionné, conduisant à une convergence de l'algorithme non garantie dans certaines situations. Les méthodes décrites dans les articles de Xiaoguang Wan et al. et Puengnim A. et al. reposent sur l'utilisation de filtres de Kalman étendus, la première mettant en oeuvre des mesures du code des signaux et l'approche « Least square 10 AMBiguity Adjustement » (LAMBDA) avec une position initiale statique, et la deuxième utilisant le mouvement du récepteur pour linéariser le problème. Les approches décrites dans les articles de Petrovski I. et al. "Precise Navigation Indoor", IEEE SICE Annual Conference 2004, vol 2, pages 1739-1744, et de Kao Wei-Wen et Tsai Chin-Lang, "Carrier Phase Indoor Positioning using Pseudolites 15 And INS" ION GNSS 2003, se fondent sur le même principe mais en utilisant le mouvement, de la station de base dans le premier cas, ou du récepteur, à l'aide d'un capteur inertiel dans le deuxième cas. D'autres approches utilisant uniquement des mesures Doppler sont connues, notamment décrites dans l'article de Schelkshorn S. et Detlefsen J. "Indoor Navigation 20 Based On Doppler Measurements", IEEE Workshop on Positioning, Navigation and Communication 2007, pages 37-40. Cet article présente une approche effectuant simultanément des mesures Doppler d'un objet en mouvement à localiser par quatre capteurs situés à différentes positions. L'article de Sakamoto Y. et al. "Doppler Positioning with a Movable Receiver Antenna and a Single Pseudolite for Indoor 25 Localization", IEEE/ASME International Conference on Advanced Intelligent Mechatronics 2011, pages 19-24, présente une méthode de positionnement par Doppler uniquement appliquée à un robot ayant une antenne rotative. Il existe un besoin pour obtenir une localisation précise au sein d'un système de positionnement, permettant la continuité entre les milieux intérieur et extérieur, et 30 mettant en oeuvre une infrastructure peu complexe et compatible avec l'ensemble des terminaux mobiles actuels. 3038067 7 L'invention a pour but de répondre à ce besoin et elle y parvient, selon l'un de ses aspects, grâce à un procédé de localisation d'au moins un récepteur au sein d'un système de positionnement, le système comprenant : - au moins deux émetteurs, chaque émetteur émettant un signal comprenant 5 une porteuse modulée par un code, et un récepteur mobile au sein du système, configuré pour détecter les signaux émis par les émetteurs, procédé dans lequel : - durant le déplacement du récepteur, à partir d'une position initiale estimée 10 prédéfinie de ce dernier, on effectue des mesures successives de la phase de la porteuse du signal émis par chaque émetteur, pour différentes positions du récepteur, on calcule les variations de la phase de la porteuse des signaux entre chaque position ultérieure du récepteur pour laquelle la phase a été mesurée et la position initiale estimée du récepteur, pour chaque émetteur, et 15 - on utilise ces variations de phase pour calculer la variation de distance entre le récepteur et les émetteurs, afin de déterminer la position initiale réelle du récepteur au sein du système de positionnement. La méthode selon l'invention s'affranchit, lors des mesures, de la connaissance de la valeur absolue de la distance séparant l'émetteur du récepteur. On obtient ainsi un 20 positionnement absolu, grâce au mouvement du récepteur et à partir de mesures relatives des variations de phase de la porteuse, fournissant des résultats proches de ceux d'un Doppler intégré. Les mesures de variation de distances entre les émetteurs et le récepteur sont possibles, avec une précision élevée, par exemple décimétrique ou centimétrique, par des 25 mesures de la phase de la porteuse de signaux radio, notamment de type GNSS. Ces mesures traduisent la vitesse de déplacement radiale relative entre les émetteurs, qui sont fixes selon l'invention, et le récepteur. La méthode de localisation selon l'invention a une sensibilité moins grande aux perturbations dues à la propagation des signaux. L'invention ne nécessite qu'une infrastructure réduite, notamment un nombre 30 réduit d'émetteurs, indépendants les uns des autres, et ne nécessitant pas d'être synchronisés entre eux. 3038067 8 L'invention permet la continuité du positionnement entre les milieux extérieurs et intérieurs. Le procédé selon l'invention prend en compte le fait que, en fonction de l'endroit où l'on se déplace, les distances entre les émetteurs et le récepteur varient 5 différemment. Une fois que la position initiale du récepteur est déterminée, un positionnement selon les méthodes connues peut être réalisé. L'invention peut être mise en oeuvre avec tous les récepteurs de terminaux mobiles actuels, notamment de téléphones portables. 10 Le procédé peut être mis en oeuvre en intérieur (en indoor), servant par exemple pour le positionnement d'objets dans des zones non couvertes par le GPS ou le GNSS, par exemple des bâtiments. Le procédé selon l'invention se fonde avantageusement sur la génération de signaux de type GNSS transmis par voie radio, par exemple sur des fréquences RNSS, 15 «Radio Navigation Satellite Services» en anglais, susceptibles d'être traités par des récepteurs GNSS, au moins en ce qui concerne leur traitement numérique. La fréquence des signaux émis par les émetteurs peut être quelconque, mais est de préférence égale à 1,575 GHz, les signaux étant par exemple de type GNSS Galileo. Le terme « émetteur » peut se substituer au terme « générateur », et vice-versa. 20 Pseudo-distance Comme évoqué précédemment, et comme représenté à la figure 1, un pseudolite PL1, PL2, PL3, PL4 est un émetteur synchronisé sur un système global 31 de temps de référence, commun à l'ensemble de la constellation. On peut mesurer la pseudodistance PRk(t) d'un pseudolite PLk à un instant t, ou «pseudorange » en anglais, 25 correspondant à une mesure indirecte de distance par le repérage de l'instant de réception d'un signal provenant du pseudolite et daté à l'émission, lorsque les horloges du pseudolite et du récepteur 30 ne sont pas synchronisées. Cette mesure peut s'exprimer par : V(x(t) - x pik)2 (Y(t)- pik)2 (z(t) z pik)2 + c.b(t) = PR k(t) , avec (x(t), y(t), z(t)) les coordonnées du récepteur à l'instant t, (xpik, ypik, zpik) les 30 coordonnées du pseudolite PLk, supposées connues, c la vitesse de la lumière et b(t) le biais d'horloge entre le temps de référence de la constellation et l'horloge du récepteur à 3038067 9 l'instant t. Il est à noter que la hauteur z est considérée comme connue dans le cas d'un positionnement en deux dimensions. La pseudo-distance peut être mesurée en utilisant le délai du code ou la phase de la porteuse du signal entrant. Dans le cas où la phase est utilisée, l'expression devient : 5 .N1(X(t)- Xpik )2 (Y(t) Yp/k )2 ± (Z(t) Zpik )2 c.b(t) = k (t) k avec Ok la phase de la porteuse dépendant du temps t, la longueur d'onde et Nk l'ambiguïté de la phase introduite précédemment, ne dépendant pas du temps. Variations de la distance récepteur-émetteurs Comme l'ambiguïté de la phase ne dépend pas du temps, il est avantageux, 10 pour l'éliminer, de mesurer la pseudo-distance d'un émetteur selon l'invention à deux moments successifs t1 et t2 et de soustraire les équations obtenues : ( X (t 2) X plk )2 + (Y (t2 - Y pik )2 (Z (t2) Z plk )2 - x t1)- x pik )2 + (Y (ti ) y pik ) 2 + Z t ) Z plk ) 2 + c .( b ( t 2 ) - b ( 1'1 )) = k (t2 ) k (t1) Dans un souci de simplicité, dans la suite, un positionnement en deux 15 dimensions sera considéré, et les inconnues x(t1), y(t1), x(t2), y(t2) et b(t2)-b(t1) seront notées respectivement x1, yt, x2, yz et Abu, et la phase it^k(-tt) sera notée (Plik. Au moins un troisième instant t3 peut être considéré, correspondant à une position du récepteur x3, yz et à un biais d'horloge b(t3), pour obtenir une nouvelle mesure et augmenter le nombre d'équations. 20 Le nombre d'émetteurs no est avantageusement fonction de la dimension m du positionnement et du nombre de mesures kpt réalisées pour des positions différentes du récepteur : m.(1 + -1) n k PI - pt Cette relation permet de connaître le nombre minimum d'émetteurs nécessaires 25 pour la bonne mise en oeuvre du procédé selon l'invention. 3038067 10 La dimension m du positionnement peut être égale à deux ou trois, correspondant respectivement à un positionnement en deux ou en trois dimensions. Dans le cas d'un nombre limité d'émetteurs, notamment deux émetteurs, des indications supplémentaires peuvent faciliter la détermination de la position initiale du 5 récepteur, par exemple une observation de la zone de l'environnement dans lequel évolue le récepteur et de sa trajectoire, permettant de discriminer certaines symétries, et/ou de déterminer si la trajectoire du récepteur coupe ou non la droite passant par les deux émetteurs. Le procédé selon l'invention permet de s'adapter aux conditions requises et 10 disponibles : un nombre plus élevé d'émetteurs peut être utilisé si peu de mesures de la phase doivent être réalisées, par exemple pour obtenir rapidement la position du récepteur, ou, au contraire, un plus faible nombre d'émetteurs peut être utilisé si une infrastructure réduite est souhaitée, un nombre plus élevé de mesures de la phase étant alors réalisé. La longueur du déplacement du récepteur entre deux positions est également à déterminer, 15 selon la géométrie de la région couverte par le système de positionnement. Le système d'équations obtenues selon le nombre de mesures effectuées peut être linéarisé et résolu itérativement par l'algorithme de Newton ou par l'algorithme des moindres carrés, appelé « Least-Square » en anglais. On peut calculer la variation de la phase de la porteuse entre une position 20 ultérieure (xi, yi, zi) et la position initiale estimée prédéfinie (xi, yi, zi) du récepteur, correspondant à l'instant ti et à l'instant ultérieur ti, par la relation : Xj - Xpik )2 + (Yi Yp/k )2 + (1 - Ip/k )2 - V(X1 - Xplk\) 2 - 2 + (z - zplk)2 (Y1 Y plk) =0 - j 01 k k On peut utiliser une série de Taylor du premier ordre pour calculer la variation de la phase de la porteuse : 25 - aykdyi + ai dx. + ai dy + c. dAb - ,+,i 61,1 xk Yk ' `11( 3038067 11 \ 2 f avec pi, = - xplk) 2-F^yu Yplk)2 + - Zpik)2 + C. b(tu), au = (xpik-szu) , et Xk 2 2 r \i(Ssu-xplk) +07u-Yplk) +(z-zpik) \ \ \ 2 (37p1k-37u) aYk \ 2 r,, \ 2 \ 2 -\/(Ss lk) +0u-Yplk) +(z-zpik) u p 5 Les valeurs dxi, dyi, dxi, dy et dAbij, correspondant à la variation de distance entre la position initiale estimée et la position initiale réelle du récepteur, sont avantageusement déterminées à chaque itération de l'algorithme. Un produit de matrices est avantageusement utilisé pour calculer la variation de distance AX entre le récepteur et les émetteurs : H. AX = d0(1), --Al A2 0 0 0 - -A1 0 A3 0 0 10 avec H = ai xi = a1 xnpl -A1 0 0 Ai 0 0 -A1 0 0 - - - Akpt- aixl pour j = 1, = Ynpl-npl ai 1 yi ai 1 pour j > 1, et Ynpi dxi dx2 dy2 c. dAbi2 dxj dyj c. dAbii dxkpt dykpt c. dAbikpt_ 'hl (-A2 (g_ - FU) (g pi (po 2n p1no (g_ - (1)1- - - Pl) ^n1 (g pl (rnpl rnpl - ipt - 1 _ (-1(pt - Pi) AX = ,d041 = A,k t ((pt i- Ynpi epin nn pi)_ 3038067 12 La matrice H peut être inversée pour obtenir la matrice de variations de distance AX = . Le procédé selon l'invention est de préférence réitéré tant que la variation de distance entre le récepteur et les émetteurs est supérieure à un premier seuil prédéfini. 5 Le choix de la position initiale estimée du récepteur est fondamental pour la convergence de l'algorithme. Le procédé selon l'invention comporte avantageusement une étape dans laquelle, si la variation de distance entre le récepteur et les émetteurs est supérieure à un deuxième seuil prédéfini, on modifie la position initiale estimée prédéfinie du récepteur. Les mesures de phase peuvent être réitérées à partir de cette nouvelle position 10 initiale estimée. Lorsque la variation de distance entre le récepteur et les émetteurs est inférieure au premier seuil prédéfini, on vérifie si la position initiale du récepteur ainsi déterminée appartient à la région couverte par le système de positionnement, délimitée par les positions des émetteurs. 15 Si la position initiale du récepteur déterminée appartient à la région couverte par le système de positionnement, la position initiale estimée du récepteur est avantageusement retenue comme position réelle, sinon, on peut réitérer le procédé selon l'invention. Système de positionnement 20 L'invention a encore pour objet, selon un autre de ses aspects, un système de positionnement, comportant : - au moins deux émetteurs, chaque émetteur émettant un signal comprenant une porteuse modulée par un code, et - un récepteur mobile au sein du système, configuré pour détecter les 25 signaux émis par les émetteurs, le récepteur étant configuré pour : - durant son déplacement au sein du système, à partir d'une position initiale estimée prédéfinie, effectuer des mesures successives de la phase de la porteuse du signal émis par chaque émetteur, pour différentes positions du récepteur, 30 - calculer les variations de la phase de la porteuse des signaux entre chaque position ultérieure du récepteur pour laquelle la phase a été mesurée et la position initiale estimée du récepteur, pour chaque émetteur, et 3038067 13 - utiliser ces variations de phase pour calculer la variation de distance entre le récepteur et les émetteurs, afin de déterminer la position initiale réelle du récepteur au sein du système de positionnement. Chaque émetteur émet avantageusement sur une fréquence différente. Cela 5 permet d'éliminer le problème du phénomène d'éblouissement ou d'interférence intrinsèque entre deux signaux émettant sur une même fréquence. Récepteur L'invention a encore pour objet, selon un autre de ses aspects, un récepteur destiné à être utilisé au sein d'un système de positionnement comprenant au moins deux 10 émetteurs, chaque émetteur émettant un signal comprenant une porteuse modulée par un code, le récepteur étant mobile au sein du système et configuré pour détecter les signaux émis par les émetteurs, le récepteur étant configuré pour : - durant son déplacement au sein du système, à partir d'une position initiale 15 estimée prédéfinie, effectuer des mesures successives de la phase de la porteuse du signal émis par chaque émetteur, pour différentes positions du récepteur, - calculer les variations de la phase de la porteuse des signaux entre chaque position ultérieure du récepteur pour laquelle la phase a été mesurée et la position initiale estimée du récepteur, pour chaque émetteur, et 20 - utiliser ces variations de phase pour calculer la variation de distance entre le récepteur et les émetteurs, afin de déterminer la position initiale réelle du récepteur au sein du système de positionnement. Le récepteur comporte avantageusement une boucle à verrouillage de phase configurée pour mesurer la phase 0, de la porteuse des signaux émis par les émetteurs. 25 Le récepteur est avantageusement multicanal. On peut utiliser comme récepteur un récepteur standard, GNSS par exemple, aucune modification d'un tel récepteur n'étant nécessaire pour mettre en oeuvre le procédé selon l'invention. Grandeur géométrique représentative de la position du récepteur 30 En variante ou en combinaison, l'invention a encore pour objet un procédé de localisation d'au moins un récepteur au sein d'un système de positionnement, le système comprenant : 3038067 14 - au moins deux générateurs, chaque générateur émettant sur une même porteuse au moins deux signaux ayant chacun un code différent, et - un récepteur configuré pour détecter les signaux émis par les générateurs, procédé dans lequel : 5 - le récepteur mesure, pour chacun des générateurs, la différence de phase entre les deux signaux émis par le générateur, et - en fonction de ces mesures de différence de phase, on calcule au moins une grandeur géométrique représentative de la position du récepteur par rapport aux générateurs, afin de localiser le récepteur au sein du système de positionnement. 10 Les mesures sont faites pour une seule position du récepteur, qui n'a ainsi pas besoin d'être mobile dans l'environnement. La méthode selon l'invention s'affranchit, lors des mesures, de la connaissance de la valeur absolue de la distance séparant le générateur du récepteur. L'invention produit des mesures instantanées non ambiguës permettant d'obtenir un positionnement absolu. 15 Les mesures des grandeurs géométriques représentatives de la position du récepteur par rapport aux générateurs sont possibles, avec une précision élevée, par exemple décimétrique ou centimétrique, par des mesures de différence de phase de signaux radio, notamment de type GNSS. En fonction du signe des différences de phase mesurées, les symétries 20 potentielles sont éliminées et un unique point est obtenu. L'utilisation de générateurs à deux codes permet d'utiliser un récepteur n'ayant qu'une antenne, correspondant à un « radar inversé ». Les deux signaux émis par chaque générateur sont de préférence émis à partir de deux zones d'émission du générateur distantes l'une de l'autre d'une distance prédéfinie 25 d12. Les deux signaux émis par chaque générateur sont avantageusement émis à la même fréquence. Dans un mode de réalisation, la grandeur géométrique représentative de la position du récepteur par rapport à un générateur est l'angle d'arrivée a des signaux émis 30 par le générateur j sur le récepteur, fonction de la distance prédéfinie d12 entre les deux 3038067 15 zones d'émission du générateur et des mesures de différence de phase 39' , et défini par : a1 = Arc cos( 391 L'angle est avantageusement pris dans le plan défini par les deux zones d'émission du générateur et le récepteur. On considère que l'angle au milieu des deux 5 zones d'émission est identique aux angles issus directement de chaque zone vers le récepteur, ce qui est le cas lorsque le récepteur est infiniment éloigné desdites zones d'émission. Il est avantageux que la distance entre les zones d'émission et le récepteur soit supérieure à au moins dix fois la distance prédéfinie entre les zones, étant par exemple supérieure à 2 mètres. 10 Au moins un deuxième calcul d'angle est avantageusement réalisé avec au moins un deuxième générateur, afin de réaliser une intersection géométrique. L'utilisation de l'angle comme grandeur géométrique rend le procédé simple, et fonctionne même à bonne distance, en supposant que les lignes d'émission des signaux sont parallèles. 15 Dans une variante ou en combinaison, la grandeur géométrique représentative de la position (xi-, yr) du récepteur par rapport à un générateur est la distance les séparant, fonction des distances di, d2 entre la position (xai, Yai, xa2, ya2) de chaque zone d'émission du générateur et la position du récepteur, définies par exemple par : d2 - d1= .^/(xa2 xr )2 ± (Ya2 - Yr )2 - V (xai/2±(Yai Yr)2 = &p pour un 20 positionnement en deux dimensions (2D). Après développement, on peut obtenir: 4392042 = (d22 d12 394 - 239 2 (0/22 042 ), d22 d12 = (X,2 Xr )2 (Ya2 Yr )2 -k(ai xr )2 +(J'ai - Yr )2 .1 2 2 a2 X al )± 2()Cal Xa2 ))Cr + ya22 Y al2 )± 2(Yal Ya2 )Yr - Et, par suite, on peut obtenir les expressions suivantes : d22 d12 2AX12Xr 2AY12Yr + g X21 + A21721 (0/22 - 042 = 43X122 Xr 2 ± 43,17122Y,2 +(2X21 +A21721)2 +84X1241712x,y, + 43,X12(A2X-21 gY21 )Xr +43X12(g -X-21 ± g1721)Yr' 0112 ) avec: 25 3038067 16 dans lesquelles: 3.X12 = Xal - Xa2 et AY12 - v .z al Ya2' 3.2X21 Xa 22 - xal2 et i A2 1 v 21 Ya 22 Yal2 - 5 On obtient la relation générale suivante, fonction des distances entre la position de chaque zone d'émission du générateur et la position du récepteur : Y21 )1 A2y21 -.-/z-21 )2 694 4 =0 10 dont il est possible de simplifier la forme globale comme suit : xr 2 169 yr 2 2 - AX12 2 .1+ 2 ] - X12 (A2 x21 ± ± 2AX12 [39 2 yu yr Xr [6(P 2 AX12 269 2 Xa1 (A2 xll - ± 39 A1712 xr 2 (cal 2 ± yai + yr [39 2 Ayi2 239 2yai (LI 2 'A,2 2 2 ± 4 2 )± "P2 (A2 x2i 3.2 y2i ) A(x r - ef )2 +B(yr-yref)2 +C)cryr +D =0, avec A, B C et D d'une part, et Xref et yref d'autre part, des coefficients déterminés par la 15 géométrie du problème et la mesure de différence de phase 69 . Cela revient à déterminer une hyperbole pour un positionnement en deux dimensions, ou un hyperboloïde en trois dimensions (3D) sur lequel se trouve le récepteur, les zones d'émission des générateurs étant leurs foyers. Au moins un deuxième générateur est avantageusement utilisé pour obtenir un 20 système d'équations et une intersection des hyperboles, pour obtenir la position (x,-, yr) du récepteur. L'utilisation comme grandeur géométrique de la distance séparant le récepteur d'un générateur, et l'expression sous la forme d'hyperboles, ne nécessite aucune approximation. 25 Il faut au moins deux générateurs pour un positionnement en deux dimensions et trois générateurs pour un positionnement en trois dimensions pour réaliser les intersections géométriques et obtenir un point unique. La position des générateurs au sein du système de positionnement est supposée connue. Les générateurs sont de préférence statiques. 3038067 17 Les deux signaux émis par chaque générateur peuvent être issus d'au moins deux antennes appartenant au générateur et distantes l'une de l'autre d'une distance prédéfinie. La distance prédéfinie dont les deux antennes d'un même générateur sont 5 distantes l'une de l'autre est avantageusement égale à la longueur d'onde de la porteuse des signaux émis par les antennes. Cela permet d'obtenir des mesures non ambiguës au niveau du récepteur. Les signaux émis par un même générateur sont de préférence synchronisés entre eux, par exemple en interne par le générateur. 10 Chaque générateur émet avantageusement sur une fréquence différente. Cela permet d'éliminer le problème du phénomène d'éblouissement ou d'interférence intrinsèque entre deux signaux émettant sur une même fréquence. Le récepteur est avantageusement multicanal, au moins un premier canal étant configuré pour détecter le premier signal d'un générateur, et au moins un deuxième canal 15 étant configuré pour détecter le deuxième signal du même générateur. La combinaison des deux procédés définis ci-dessus présente un certain intérêt. Un système complet de positionnement mettant en oeuvre ces deux procédés permet, à partir des mêmes mesures de phases, d'avoir au moins trois modes de fonctionnement : un positionnement absolu en effectuant des différences, à un instant donné, entre les deux 20 antennes d'un générateur et cela pour plusieurs générateurs, puis un positionnement relatif à partir de la position obtenue, d'une précision décimétrique, par mesures de différence de phases à des instants successifs, puis un positionnement absolu en mouvement par analyses de l'évolution de la phase des porteuses en provenance des diverses antennes à plusieurs instants successifs. 25 L'invention pourra être mieux comprise à la lecture de la description qui va suivre, d'exemples non limitatifs de mise en oeuvre de celle-ci, et à l'examen du dessin annexé, sur lequel : - la figure 1 représente de façon schématique un système de positionnement, selon l'art antérieur, utilisant des pseudolites, 30 - la figure 2 illustre des étapes de mise en oeuvre du procédé selon l'invention, 3038067 18 - la figure 3 représente de façon schématique un système de positionnement dans lequel peut être mis en oeuvre le procédé selon l'invention, - la figure 4 illustre un résultat de trajectoire d'un récepteur obtenue par la mise en oeuvre du procédé selon l'invention, 5 - la figure 5 représente de façon schématique un système de positionnement dans lequel peut être mis en oeuvre une variante de procédé selon l'invention, - la figure 6 représente la géométrie mise en oeuvre dans un système de positionnement selon une variante de l'invention, et - les figures 7 à 9 sont des courbes montrant les performances d'une variante 10 de procédé selon l'invention. On a représenté à la figure 3 un exemple de système 1 dans lequel peut être mis en oeuvre le procédé selon l'invention, dont des étapes sont illustrées à la figure 2. Le système 1 comprend un récepteur 2 mobile au sein du système 1 et une pluralité d'émetteurs PL1, PL2, PL3, PL4 formant une constellation locale. Comme 15 représenté sur la figure 3, le système 1 peut être mis en oeuvre en intérieur (en indoor), par exemple à l'intérieur d'un bâtiment, ou à l'extérieur dans une zone très urbaine, entre des murs très hauts par exemple. Chaque émetteur PL1, PL2, PL3, PL4 émet de préférence un signal comprenant une porteuse modulée par un code. 20 Les signaux émis par les émetteurs PL1, PL2, PL3, PL4 sont par exemple des signaux GNSS de fréquence égale à 1,575 GHz. Le récepteur 2 est configuré pour détecter les signaux émis par les émetteurs. Ces derniers sont par exemple reçus par l'antenne du récepteur 2 puis amplifiés et convertis en fréquence intermédiaire FI plus basse que leur fréquence initiale. 25 Ces signaux sont dans l'exemple considéré échantillonnés, puis numérisés avant d'être traités par les canaux de réception du récepteur 2. Ces canaux de réception peuvent mettre en oeuvre des boucles de poursuite. Le récepteur 2 comporte avantageusement une boucle à verrouillage de phase configurée pour mesurer la phase c de la porteuse des signaux émis par les émetteurs. 30 Comme expliqué précédemment, et comme représenté à la figure 2, lors d'une étape 11, une position initiale estimée prédéfinie X est choisie pour le récepteur 2. Durant son déplacement au sein du système de positionnement 1, à partir de cette position initiale 3038067 19 estimée, lors d'une étape 12, le récepteur 2 effectue des mesures successives de la phase Ok de la porteuse du signal émis par chaque émetteur PL1, PL2, PL3, PL4, pour différentes positions ultérieures du récepteur (xi, n, zi). Lors d'une étape 13, on calcule les variations 31:1 de la phase de la porteuse des 5 signaux entre chaque position ultérieure du récepteur 2 pour laquelle la phase Ok a été mesurée et la position initiale estimée du récepteur, pour chaque émetteur PL1, PL2, PL3, PL4. La matrice d31:1 définie précédemment est créée lors d'une étape 13bis. Afin de calculer la variation de distance AX entre le récepteur 2 et les émetteurs et de déterminer la position initiale réelle du récepteur 2 au sein du système de 10 positionnement 1, la matrice H définie précédemment est créée lors d'une étape 14 et est inversée lors d'une étape 15. La variation de distance AX est calculée lors d'une étape 16. La variation de distance AX ainsi calculée et la position initiale estimée prédéfinie X du récepteur 2 sont additionnées lors d'une étape 17 afin de former une nouvelle position initiale X. 15 Lors d'une étape 18, la variation de distance AX est comparée à un deuxième seuil prédéfini £2. Si AX est supérieure à ce deuxième seuil prédéfini, on modifie la position initiale estimée prédéfinie X du récepteur 2, les mesures de phase et de variations 301:1 de phase étant réitérées à partir de cette nouvelle position initiale estimée. Si AX est inférieure à ce deuxième seuil prédéfini £2, on la compare à un 20 premier seuil prédéfini ci, le procédé selon l'invention étant avantageusement réitéré, à partir de l'étape 14, tant que la variation de distance AX entre le récepteur et les émetteurs est supérieure à ce premier seuil prédéfini. Lorsque la variation de distance AX entre le récepteur 2 et les émetteurs est inférieure au premier seuil prédéfini ci, on vérifie, lors d'une étape 19, si la position initiale 25 du récepteur ainsi déterminée appartient à la région couverte par le système de positionnement 1, délimitée par les positions des émetteurs. Si la position initiale du récepteur déterminée appartient à la région couverte par le système de positionnement 1, la position initiale estimée X du récepteur 2 est retenue comme position réelle lors d'une étape 20, sinon, on réitère le procédé. 3038067 20 Le premier seuil prédéfini est compris entre 10-5 m et 10-1 m, étant par exemple égal à 10-2 m. Le deuxième seuil prédéfini est compris entre 102 m et 105 m, étant par exemple égal à 103 m. La figure 4 représente des résultats obtenus en mettant en oeuvre l'invention 5 dans un canyon urbain en forme de U de 20 m par 30 m, formés par des immeubles hauts. Dans cet exemple, quatre émetteurs sont déployés au dernier étage des immeubles, à une altitude de 18 m environ. Un récepteur 2 est en mouvement dans cet environnement. Le récepteur 2 a comme position initiale estimée P0, puis effectue un aller-retour entre les positions PO et 10 P3 en passant par les positions ultérieures P1 et P2, puis va de la position PO à la position P4. Le procédé de localisation du récepteur 2 selon l'invention est mis en oeuvre et les positions obtenues sont comparées avec la trajectoire réelle, comme visible à la figure 4. La trajectoire reconstituée obtenue est très proche de la trajectoire réelle, connue au préalable pour les besoins de l'expérience, l'écart restant en dessous de 50 cm, comme 15 représenté dans le tableau ci-dessous montrant l'erreur et l'erreur moyenne pour chaque position. Position X(m) Y(m) Erreur(m) PR(m) EHDOP Erreur moyenne PO 4.04 14.82 0.44 0.05 26.7 P1 4.39 17.56 0.43 0.05 P2 5.61 20.78 0.40 0.06 P3 8.8 24.15 0.13 0.28 P4 4.09 3.37 1.08 0.23 P4c 3.72 4.00 0.73 0.23 4.2 Les positions P3 et P4 donnent les erreurs les plus grandes, car elles sont plus difficiles à estimer de par la configuration de la trajectoire, correspondant à un chemin 20 multiple et des glissements de cycle cumulant des sources d'erreur. Le procédé selon l'invention a ainsi été utilisé pour corriger la position P4, en utilisant une mesure supplémentaire à une position intermédiaire entre P3 et P4, conduisant au résultat P4c dans le tableau ci-dessus, donnant une erreur plus faible. La dernière colonne du tableau montre les résultats obtenus pour la dilution 25 horizontale étendue de la précision, ou « extended horizontal dilution of precision » en 3038067 21 anglais, valeur spécifiant l'influence de la géométrie de l'environnement sur la précision du système de positionnement, et adaptée au procédé selon l'invention, représentant ainsi l'influence de l'erreur de mesures sur la position initiale du récepteur. En considérant que les erreurs sur chaque mesure de phase sont des gaussiennes centrées en zéro, et que ces 5 erreurs sont distribuées de manière identique pour chaque émetteur et indépendantes les unes des autres, la dilution horizontale étendue de la précision peut être calculée à partir de la formule : cov(Ex) = (11'11)- 1 'ARE/ avec cov(Ex) la matrice de covariance de l'erreur sur la position estimée du récepteur, et 0-1JERE la variance de l'erreur sur la mesure, ou « User Equivalent Range Error » en anglais. En notant o les éléments de cov(cx) et Dij VOE11+OE22 OEUERE 10 les éléments de (Eltil)-1, on obtient : EHDOP =+ D v - il + - 22 = Une valeur de dilution horizontale étendue de la précision comprise entre 10 et 40 est considérée comme acceptable, les résultats obtenus pour les positions PO à P4 sont donc corrects. Néanmoins, la valeur de dilution horizontale étendue de la précision obtenue pour la position P4c définie ci-dessus est excellente. Cela montre qu'il est intéressant et 15 nécessaire d'utiliser le procédé selon l'invention au cours de l'évolution du récepteur au sein du système pour diminuer les erreurs s'accumulant le long de la trajectoire. Comme expliqué précédemment, le nombre d'émetteurs npi du système 1 est fonction de la dimension m du positionnement et du nombre de mesures kpt réalisées pour 1 des positions différentes du récepteur 2 : m.(1 + n pl . Par exemple, pour m = 3 11' , pt 20 correspondant à un positionnement en 3D, le nombre minimum d'émetteurs npi est égal à 5 si trois positions différentes sont utilisées pour faire trois mesures kpt : 3(1+1/3) = 4. Si seulement une mesure kpt est utilisée, pour avoir la position du récepteur rapidement, il faut : 3(1+1) = 6 émetteurs. On a représenté à la figure 5 un exemple de système 1 dans lequel peut être 25 mise en oeuvre une variante de procédé selon l'invention. Le système 1 comprend un récepteur 2 et trois générateurs PL1, PL2, PL3 formant une constellation locale. Chaque générateur PL1, PL2, PL3 émet, dans l'exemple décrit, sur une même porteuse, deux signaux ayant chacun un code différent. Les émetteurs peuvent être les 30 générateurs décrits précédemment. Ces deux signaux sont émis à partir de deux zones 3038067 22 d'émission PL11, PL12 du générateur PL1 distantes l'une de l'autre d'une distance prédéfinie d12. Dans l'exemple décrit et de préférence, les deux signaux émis par chaque générateur PL1, PL2, PL3 sont issus de deux antennes appartenant au générateur et distantes l'une de l'autre d'une distance prédéfinie égale par exemple à la longueur d'onde 5 de la porteuse des signaux émis par les antennes. Les signaux émis par un même générateur sont de préférence synchronisés entre eux. Chaque générateur PL1, PL2, PL3 émet de préférence sur une fréquence différente. Les signaux émis par les générateurs PL1, PL2, PL3 sont par exemple des signaux GNSS de fréquence égale à 1,575 GHz. 10 Le récepteur 2 est configuré pour détecter les signaux émis par les générateurs. Le récepteur 2 est avantageusement configuré pour mesurer, pour chacun des générateurs PL1, PL2, PL3, la différence de phase &p' entre les deux signaux émis par le générateur PLj. Comme défini précédemment, en fonction de ces mesures de différences de 15 phase 39' , le récepteur 2 est avantageusement configuré pour calculer au moins une grandeur géométrique représentative de la position du récepteur par rapport aux générateurs PL1, PL2, PL3, afin d'obtenir sa position (xr, yr) au sein du système 1. Dans l'exemple de la figure 6, et comme défini précédemment, la grandeur géométrique représentative de la position (xr, yr) du récepteur 2 par rapport à un générateur 20 PLj est l'angle d'arrivée ai des signaux émis par le générateur PLj sur le récepteur, fonction de la distance prédéfinie d1 entre les deux zones d'émission du générateur et des mesures de différences de phase &p' : = Arc cos91 \d12 ) Dans une variante, la grandeur géométrique représentative de la position (xr, yr) du récepteur 2 par rapport à un générateur PLj est la distance les séparant, fonction des 25 distances (dl, d2) entre la position (xai, yai, xa2, ya2) de chaque zone d'émission PLji, PLj2 du générateur PLj et la position (xr, yr) du récepteur 2: d2 dl =./(xa2 Xr )2 + (Ya2 Yr )2 - V(Xal Xr )2 + (Yal Yr )2 Les figures 7 à 9 représentent des résultats obtenus en mettant en oeuvre l'invention dans un hall de 20 m par 20 m, dans lequel deux générateurs PL1 et PL2 sont 30 positionnés aux positions (0, 0) et (20, 0), étant ainsi positionnés du même côté du hall. 3038067 23 Le procédé selon l'invention est mis en oeuvre, comme défini précédemment, pour trouver un récepteur 2, se trouvant aux coordonnées (5;15) du hall. Environ 660 mesures de différence de phase sont réalisées afin d'évaluer les performances du procédé selon l'invention. 5 La figure 7 montre la distribution des erreurs de position par mesure, correspondant à la distance entre la position du récepteur 2 calculée et sa position réelle. Selon la figure 7, la précision du positionnement semble être centrée sur une trentaine de centimètres. La figure 8 représente la distribution des erreurs de position en fonction des 10 erreurs de mesure de la phase des signaux. La figure 9 montre les positions du récepteur 2 calculées selon l'invention dans le plan Oxy, le point recherché étant représenté par un carré aux coordonnées (5, 15) et les divers points calculés étant représentés par des losanges foncés. Des moyennes glissantes sur 10 points ont été réalisées, représentées par des triangles, permettant de confirmer que 15 la précision du positionnement est de l'ordre d'une trentaine de centimètres. Le losange clair représente la moyenne générale du point obtenu sur les 650 mesures, qui se trouve être à 7,5 cm de la position réelle du récepteur 2. L'invention n'est pas limitée aux exemples qui viennent d'être décrits. D'autres signaux que les signaux GPS et GNSS peuvent être utilisés, par 20 exemple des signaux Wifi, 4G, GSM ou radio. Les procédés selon l'invention peuvent être couplés à d'autres procédés de géolocalisation, tels que des procédés utilisant un lecteur mobile, des balises et une application de géolocalisation et de guidage exécutable sur ledit lecteur. L'invention peut être mise en oeuvre dans des lieux où sont pratiqués des sports 25 collectifs, notamment en intérieur, par exemple en fixant ou en incorporant un système dans des cages de football ou handball afin de calculer en temps réel la position des joueurs, ou sur des chantiers de construction, par exemple pour la pose de faux plafonds, en utilisant une ou plusieurs réglettes constituées d'un réseau d'une dizaine d'antennes, ou dans des grands centres commerciaux. 30 L'invention peut être utilisée en extérieur, permettant de s'affranchir de station de base. 3038067 24 L'expression « comportant un » doit être comprise comme synonyme de l'expression « comprenant au moins un », sauf lorsque le contraire est précisé. REVENDICATIONS1. Procédé de localisation d'au moins un récepteur (2) au sein d'un système de positionnement (1), le système comprenant : - au moins deux émetteurs (PL1, PL2), chaque émetteur émettant un signal comprenant une porteuse modulée par un code, et - un récepteur (2) mobile au sein du système (1), configure pour détecter les signaux émis par les émetteurs (PLI, PL2), procédé dans lequel : - durant le déplacement du récepteur (2), à partir d'une position initiale estimée prédéfinie (xi, yi, zi) de ce dernier, on effectue des mesures successives de la phase (0,) de la porteuse du signal émis par chaque émetteur (PLk), pour différentes positions ultérieures (xi, yi, zi) du récepteur (2), - on calcule les variations de la phase (d34) de la porteuse des signaux entre chaque position ultérieure (xj, yi, zi) du récepteur (2) pour laquelle la phase (0,) a été mesurée et la position initiale estimée (xi, yi, zi) du récepteur (2), pour chaque émetteur (PLk), et - on utilise ces variations de phase (d34) pour calculer la variation de distance (AX) entre le récepteur (2) et les émetteurs (PL1, PL2), afin de déterminer la position initiale réelle du récepteur (2) au sein du système de positionnement (1). 2. Procédé selon la revendication 1, dans lequel on calcule la variation de la phase (0k -0,1) de la porteuse entre une position ultérieure (xj, yi, zi) et la position initiale estimée prédéfinie (xi, yi, zi) du récepteur (2) par la relation : Xj- Xplk )2 ± (Yi Yp/k )2 +(1 - Ip/k )2 - V(X1 - Xplk )2 + (Y1 - Yp/k )2 + (11 - Zplk )2 =Aj avec (xpik, ypik, zpik) les coordonnées de l'émetteur (PLk), (c) la vitesse de la lumière, et Ab, le biais d'horloge entre l'émetteur (PLk) et le récepteur (2). 3038067 26 3. Procédé selon l'une quelconque des revendications précédentes, dans lequel on utilise une série de Taylor du premier ordre pour calculer la variation de la phase de la porteuse : - aykdyi + aix dx. + ai dy + c. dAb (2.,i 2,1 k Yk ' - F'k), \ 2 f avec pi, = - xplo 2-F^yu Yplk) + - Zpik)2 + C. b(tu), 5 au = (xpik-szu) , et xk .\/(Ss ) -x lk -Fl ) Yu-Yplk u p (Z-zpik) \ 2 r,, \ 2 \ 2 (37p1k-37u) aYk \ 2 f, \ 2 .\i(5s(u-xplk) +Uu-Yplk) +(z-zplk)2. 10 4. Procédé selon l'une quelconque des revendications précédentes, dans lequel un produit de matrices est utilisé pour calculer la variation de distance (OA) entre le récepteur (2) et les émetteurs (PLk) : H. AX = -A1 A2 0 0 0- -A1 0 A3 0 0 avec H = -A1 0 --- 0 Ai 0 . 0 -A1 0 0 - - - O Akpt a a axe a Yi xi Y 1 1 i i x1 1 Al = [ al xnpl al Ynpli pour j = 1, Ai = [ . : aj xn i pour j > 1, et i pi aYnpl 1 15 3038067 27 dx1 dy1 dx2 dy2 c. dAbi2 dxi dyi c. dAbli dxkpt dykpt c. dAbikpt_ (g_7,2 ;;U)1) - 61,1 ( F 1:12no (pno (pno pno) - (1)1 - - Pl) (g pl (rnpl rnpl) - (pipt - _ (-1(pt -11 - Pi) AX = , d3,41 = pt (,k t Pn13131 Pn1)- - pl la matrice (H) étant de préférence inversée pour obtenir la matrice de variations de distance : AX = d31). 1/-1. 5. Procédé selon l'une quelconque des revendications précédentes, étant réitéré 5 tant que la variation de distance (0A) entre le récepteur (2) et les émetteurs (PLk) est supérieure à un premier seuil prédéfini (Ci). 6. Procédé selon l'une quelconque des revendications précédentes, comportant en outre une étape dans laquelle, si la variation de distance (0A) entre le récepteur (2) et les émetteurs (PLk) est supérieure à un deuxième seuil prédéfini (c2), on modifie la position 10 initiale estimée prédéfinie du récepteur (2), les mesures de phase étant notamment réitérées à partir de cette nouvelle position initiale estimée. 7. Procédé selon l'une quelconque des revendications précédentes, dans lequel, lorsque la variation de distance (0A) entre le récepteur (2) et les émetteurs (PLk) est inférieure au premier seuil prédéfini (Ci), on vérifie si la position initiale du récepteur (2) 15 ainsi déterminée appartient à la région couverte par le système de positionnement (1), délimitée par les positions des émetteurs (PLk). 8. Procédé selon la revendication précédente, dans lequel, si la position initiale du récepteur (2) déterminée appartient à la région couverte par le système de positionnement (1), la position initiale estimée du récepteur (2) est retenue comme position 20 réelle, sinon, on réitère le procédé. 9. Procédé selon l'une quelconque des revendications précédentes, dans lequel le nombre d'émetteurs (no) est fonction de la dimension (m) du positionnement et du nombre de mesures (4) réalisées pour des positions différentes du récepteur (2) : 3038067 28 m.(1 + -1) n P1 - pt 10. Procédé selon l'une quelconque des revendications précédentes, dans lequel la fréquence des signaux émis par les émetteurs (PLk) est égale à 1,575 GHz, les signaux étant notamment de type GNSS Galileo. 5 11. Procédé selon l'une quelconque des revendications précédentes, étant mis en oeuvre en intérieur (indoor). 12. Système (1) de positionnement, comportant : au moins deux émetteurs (PL1, PL2), chaque émetteur émettant un signal comprenant une porteuse modulée par un code, et 10 - un récepteur (2) mobile au sein du système (1), configuré pour détecter les signaux émis par les émetteurs (PLI, PL2), le récepteur (2) étant configuré pour : - durant son déplacement au sein du système (1), à partir d'une position initiale estimée prédéfinie (xi, yi, zi), effectuer des mesures successives de la phase (¢k) 15 de la porteuse du signal émis par chaque émetteur (PLk), pour différentes positions ultérieures (xi, yi, zi) du récepteur (2), - calculer les variations de la phase (d34) de la porteuse des signaux entre chaque position ultérieure (xi, yi, zi) du récepteur (2) pour laquelle la phase (¢k) a été mesurée et la position initiale estimée (xi, yi, zi) du récepteur (2), pour chaque émetteur 20 (PLk), et utiliser ces variations de phase (d34) pour calculer la variation de distance (OA) entre le récepteur (2) et les émetteurs (PL1, PL2), afin de déterminer la position initiale réelle du récepteur (2) au sein du système de positionnement (1). 13. Récepteur (2), destiné à être utilisé au sein d'un système (1) de 25 positionnement comprenant au moins deux émetteurs (PL1, PL2), chaque émetteur émettant un signal comprenant une porteuse modulée par un code, le récepteur (2) étant mobile au sein du système (1) et configuré pour détecter les signaux émis par les émetteurs (PLI, PL2), le récepteur (2) étant configuré pour : 3038067 29 - durant son déplacement au sein du système (1), à partir d'une position initiale estimée prédéfinie (xi, yi, zi), effectuer des mesures successives de la phase (¢k) de la porteuse du signal émis par chaque émetteur (PLk), pour différentes positions ultérieures (xj, yj, zj) du récepteur (2), 5 calculer les variations de la phase (d34) de la porteuse des signaux entre chaque position ultérieure (xj, yj, zj) du récepteur (2) pour laquelle la phase (0,) a été mesurée et la position initiale estimée (xi, yi, zi) du récepteur (2), pour chaque émetteur (PLk), et utiliser ces variations de phase (d34) pour calculer la variation de distance 10 (OA) entre le récepteur (2) et les émetteurs (PL1, PL2), afin de déterminer la position initiale réelle du récepteur (2) au sein du système de positionnement (1). 14. Récepteur selon la revendication précédente, comportant une boucle à verrouillage de phase configurée pour mesurer la phase (0,) de la porteuse des signaux émis par les émetteurs (PLI, PL2).