FR 3038064 A1 20161230 FR 1555810 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, 3038064 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. 3038064 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 3038064 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érence 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 3038064 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 k, 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 3038064 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. 3038064 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 générateurs, chaque générateur émettant sur une même 5 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 : - 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 10 - 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. 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 15 d12. Les deux signaux émis par chaque générateur sont avantageusement émis à la même fréquence. 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 20 des mesures instantanées non ambiguës permettant d'obtenir un positionnement absolu. 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. La méthode de localisation selon l'invention a une 25 sensibilité moins grande aux perturbations dues à la propagation des signaux. Les mesures sont faites pour une seule position du récepteur, qui n'a ainsi pas besoin d'être mobile dans l'environnement. En fonction du signe des différences de phase mesurées, les symétries potentielles sont éliminées et un unique point est obtenu. 30 L'invention ne nécessite qu'une infrastructure réduite, notamment un nombre réduit de générateurs, indépendants les uns des autres, et ne nécessitant pas d'être 3038064 8 synchronisés entre eux. L'utilisation de générateurs à deux codes permet d'utiliser un récepteur n'ayant qu'une antenne, correspondant à un « radar inversé ». L'invention permet la continuité du positionnement entre les milieux extérieurs et intérieurs. 5 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. Le procédé peut être mis en oeuvre en intérieur (en indoor), servant par 10 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, « Radio Navigation Satellite Services » en anglais, susceptibles d'être traités par des 15 récepteurs GNSS, au moins en ce qui concerne leur traitement numérique. La fréquence des signaux émis par les générateurs 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. Angle d'arrivée des signaux sur le récepteur 20 Dans un premier mode de réalisation de l'invention, 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 par le générateur j sur le récepteur, fonction de la distance prédéfinie d12 entre les deux zones d'émission du générateur et des mesures de différence de phase &p', et défini par : a = Arc cos( 25 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 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 30 supérieure à au moins dix fois la distance prédéfinie entre les zones, étant par exemple supérieure à 2 mètres. 0112 ) 3038064 9 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 5 sont parallèles. Distance entre un générateur et le récepteur 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 d1, d2 entre la position (xai, Yai, xa2, ya2) de chaque zone d'émission 10 du générateur et la position du récepteur, définies par exemple par : d2 di= ^1(xa2- xr )2 (Ya2 Yr )2 - V(Xal Xr )2 (Yal Yr )2 pour un positionnement en deux dimensions (2D). Après développement, on peut obtenir: 43 2 di 2 = (d22 - d12 +894 - 239 2 (0/2 2 - 0/12 15 avec: \ 2 d22 0/12 = (X,2 Xr )2 + (Ya2 Yr )2 -krai-xr)2 +(J'ai- Yr) .1 2 2 a2 Xal )± 2(Xal Xa2 )Xr a22 Y al2 )± 2(Yal Ya2 )Yr Et, par suite, on peut obtenir les expressions suivantes : d22 d12 = 232(12 xr 2AY12 Yr + g X21 + A21721 (0/22 - 042 = 43X122 Xr 2 ± 4A17122Y,2 (2X21 g1721)2 + 8AXI2AYux,..y,. + 44X12 -X21 ± A21721 r 4 Ay12 (g X21 ± gY21 r' dans lesquelles: AX12 = Xal - xa2 et AY v 12 - al Ya2 42 X21 Xa22 - xal et i A2 1 v 21 Ya22 Ya12 - 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 : 20 25 3038064 10 Xr [69 2 - AX12 2 .1+ Xr [6(P 2 AX12 - 269 2 Xal 4X12 (g X21 ± g Y21 )1 ± yr 2 [39 2 yu 2 ]+ yr [39 2 Ayi2 239 2yai (A2 xll +42y21)1 - 24/V12417u xr yr ± 39 2 (cal 2 ± yai 2 )± 39 2 (g x2i 3,2 y2i 1 (g x2 ± gym )2 (59 4 2 \ 1 j 4 =0 dont il est possible de simplifier la forme globale comme suit : 5 A(xi. - ef )2 +B(yr-yref)2 +Cx,yr +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 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 10 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 système d'équations et une intersection des hyperboles, pour obtenir la position (x,-, yr) du récepteur. 15 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. Générateurs L'invention a encore pour objet, selon un autre de ses aspects, un générateur 20 destiné à être utilisé au sein d'un système de positionnement, ledit système comprenant au moins un autre générateur et un récepteur configuré pour détecter les signaux émis par ledit générateur et l'autre générateur, ledit générateur étant configuré pour émettre sur une même porteuse au moins deux signaux ayant chacun un code différent, le récepteur étant configuré pour mesurer la différence de phase entre lesdits deux signaux émis par le 25 générateur. 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. 3038064 11 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. 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 5 prédéfinie. La distance prédéfinie dont les deux antennes d'un même générateur sont 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. 10 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. 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. 15 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 générateurs, chaque générateur émettant sur une même porteuse au moins deux signaux ayant chacun un code différent, le récepteur détectant les signaux émis par les générateurs, 20 le récepteur étant configuré pour : - mesurer, 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, que soit calculée au moins une grandeur géométrique représentative de la position du récepteur par rapport aux 25 générateurs, afin de localiser le récepteur au sein du système de positionnement. Le récepteur comporte avantageusement une boucle à verrouillage de phase configurée pour mesurer les phases des signaux émis par les générateurs. 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 30 étant configuré pour détecter le deuxième signal du même générateur. Le récepteur peut être statique ou mobile. 3038064 12 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. Système de positionnement 5 L'invention a encore pour objet, selon un autre de ses aspects, un système de positionnement, comportant : 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 détectant les signaux émis par les générateurs, 10 le système étant configuré pour que : - 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, soit calculée au moins une grandeur géométrique représentative de la position du récepteur par rapport aux 15 générateurs, afin de localiser le récepteur au sein du système de positionnement. Variations de la distance récepteur-émetteurs 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 : 20 - 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 signaux émis par les émetteurs, procédé dans lequel : 25 durant le déplacement du récepteur, à partir d'une position initiale estimée 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 ultérieures 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 30 initiale estimée du récepteur, pour chaque émetteur, et 3038064 13 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 combinaison des deux procédés définis ci-dessus présente un certain intérêt. 5 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 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 10 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. Comme évoqué précédemment, et comme représenté à la figure 6, un pseudolite PL1, PL2, PL3, PL4 est un émetteur synchronisé sur un système global 31 de 15 temps de référence, commun à l'ensemble de la constellation. On peut mesurer la pseudo- distance PRk(t) d'un pseudolite plk à un instant t, ou «pseudorange » en anglais, 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 : 20 V(x(t) - xpll, )2 ± (Y(0 Y plk )2 ± (Z(t)- Zpik )2 ± c.b(t) = PRk(t) avec (x(t), y(t), z(t)) les coordonnées du récepteur à l'instant t, (xpik, ypik, zpik) les coordonnées du pseudoliteplk, 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 à l'instant t. Il est à noter que la hauteur z est considérée comme connue dans le cas d'un 25 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 : .N1(X(t)- Xpik )2 (y(t) Yp/k )2 ± (Z(t)- Zpik )2 c.b(t) = k (t) k, avec ¢k la phase de la porteuse dépendant du temps t, la longueur d'onde et Nk 30 l'ambiguïté de la phase introduite précédemment, ne dépendant pas du temps. 3038064 14 Comme l'ambiguïté de la phase ne dépend pas du temps, il est avantageux, pour l'éliminer, de mesurer la pseudo-distance d'un émetteur selon l'invention à deux moments successifs tt et tz et de soustraire les équations obtenues : V(x(t 2) - x pik ) 2 + ( Y ( t2 - Y pll, ) 2 + Z t2 ) Z p/k ) 2 - VGX(ti)- X plk 2 + (Y (t1) plk 2 + Z(t1) Z p/k ) 2 C b (t2 b(tl)) = k t2 ) k tl ) 5 Dans un souci de simplicité, dans la suite, un positionnement en deux dimensions sera considéré, et les inconnues x(t1), y(t1), x(tz), y(tz) et b(t2)-b(t1) seront notées respectivement x1, yt, )(2, 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 10 position du récepteur x3, yz et à un biais d'horloge b(t3), pour obtenir une nouvelle mesure et augmenter le nombre d'équations. 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 : 15 m.(1 + 1 n pl - pt Cette relation permet de connaître le nombre minimum d'émetteurs nécessaires pour la bonne mise en oeuvre du procédé selon l'invention. La dimension m du positionnement peut être égale à deux ou trois, correspondant respectivement à un positionnement en deux ou en trois dimensions. 20 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 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 25 émetteurs. Le procédé selon l'invention permet de s'adapter aux conditions requises et disponibles : un nombre plus élevé d'émetteurs peut être utilisé si peu de mesures de la 3038064 15 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, 5 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 10 ultérieure (xj, yj, zj) et la position initiale estimée prédéfinie (xi, yi, zi) du récepteur, correspondant à l'instant ti et à l'instant ultérieur tj, par la relation : xj- Xplk)2 + (y - plk)2 (1 - Ip/k )2 )2 + (Y1 -Yplk)2 _L - plk )2 On peut utiliser une série de Taylor du premier ordre pour calculer la variation de la phase de la porteuse : 15 - aykdyi + ai dx. + ai dy + c. dAb - ,+,i ,+,1 (7,i 2,1- xk Yk ' F'k), \ ' avec pi, = - xplo 2-F^yu Yplk)2 f + - Zpik)2 + C. b(tu), au - (xpik-jcu) et xk \ 2, \ f,, \ 2 Ni(SU-X131k) -Flu-Yplk) -E(Z-Zpik) 20 aYk = (Ypik- u) N2 f, \ 2 .\i(5s(u-xplk) -FlYu-Yplk) +(z-zplk)2. Les valeurs dxi, dyi, dxj, 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. 3038064 16 Un produit de matrices est avantageusement utilisé pour calculer la variation de distance AX entre le récepteur et les émetteurs : H. AX = d041, -A1 A2 0 0 0 - -A1 0 A3 0 0 avec H = -A1 0 - - - 0 Ai 0 0 -A1 0 0 Akpt_ e, 1 axi ai J3Ti 1 Al= pour j = 1, = pour j > 1, et xnpi aynpl xnpl -Ynpl 1 5 10 2n (pino (p2no pino) - (1)1 - - Pl) AX = , d3,41 = 11:111Pt rhpt `rnpl 61,1 (2 ;;1_ (g - FU) dx1 dy1 dx2 dy2 c. dAbi2 dxi dyi c. dAbu dxkpt dykpt c. dAbikpt_ ^n1 (g pl (rnpl rnpl - - - (-,/(Pt - P1) (hl Ynpl (Pnpi Pflpl La matrice H peut être inversée pour obtenir la matrice de variations de distance AX = d31). H-1. 10 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. 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 15 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 initiale estimée. 3038064 17 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. 5 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. Les émetteurs peuvent être tels que les générateurs définis précédemment. 10 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 dans lequel peut être mis en oeuvre un procédé selon l'invention, 15 - la figure 2 représente la géométrie mise en oeuvre dans un système de positionnement selon l'invention, - les figures 3 à 5 sont des courbes montrant les performances du procédé selon l'invention, - la figure 6 représente de façon schématique un système de positionnement, 20 selon l'art antérieur, utilisant des pseudolites, - la figure 7 illustre des étapes de mise en oeuvre d'une variante de procédé selon l'invention, - la figure 8 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, et 25 - la figure 9 illustre un résultat de trajectoire d'un récepteur obtenue par la mise en oeuvre d'une variante de procédé selon l'invention. On a représenté à la figure 1 un exemple de système de positionnement 1 dans lequel peut être mise en oeuvre l'invention. Le système comprend un récepteur 2 et trois générateurs PL1, PL2, PL3 30 formant une constellation locale. 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. 3038064 18 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. Ces deux signaux sont émis à partir de deux zones 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 5 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 à la longueur d'onde 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 PU, PL2, PL3 émet de préférence sur une fréquence 10 différente. Les signaux émis par les générateurs PU, PL2, PL3 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 générateurs. 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. 15 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, étant multicanal. 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 Ok de la porteuse des signaux émis par les générateurs. 20 Le récepteur 2 est avantageusement configuré pour mesurer, pour chacun des générateurs PU, 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érence de phase 39' , le récepteur 2 est avantageusement configuré pour calculer au moins une 25 grandeur géométrique représentative de la position du récepteur par rapport aux générateurs PU, PL2, PL3, afin d'obtenir sa position (xr, yr) au sein du système 1. Dans l'exemple de la figure 2, 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 PLj est l'angle d'arrivée ai des signaux émis par le générateur PLj sur le récepteur, 3038064 19 fonction de la distance prédéfinie d1 entre les deux zones d'émission du générateur et des mesures de différence de phase 39 : a = Arc cos 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 5 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 3 à 5 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 10 positionnés aux positions (0, 0) et (20, 0), étant ainsi positionnés du même côté du hall. 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. 15 La figure 3 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 3, la précision du positionnement semble être centrée sur une trentaine de centimètres. La figure 4 représente la distribution des erreurs de position en fonction des 20 erreurs de mesure de la phase des signaux. La figure 5 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 25 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. On a représenté à la figure 8 un exemple de système 1 dans lequel peut être mise en oeuvre une variante de procédé selon l'invention, dont des étapes sont illustrées à 30 la figure 7. d12j ) 3038064 20 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 représenté sur la figure 6, 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 5 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. Les émetteurs peuvent être les générateurs décrits précédemment. Dans ce cas, une seule antenne par générateur est avantageusement considérée. 10 Comme expliqué précédemment, et comme représenté à la figure 7, 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 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 15 positions ultérieures du récepteur (xj, yj, zj). Lors d'une étape 13, on calcule les variations 31:1 de la phase de la porteuse des 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. 20 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 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 25 prédéfinie X du récepteur 2 sont additionnées lors d'une étape 17 afin de former une nouvelle position initiale X. 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 30 phase étant réitérées à partir de cette nouvelle position initiale estimée. 3038064 21 Si AX est inférieure à ce deuxième seuil prédéfini £2, on la compare à un 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. 5 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 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 10 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é. 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. 15 La figure 9 représente des résultats obtenus en mettant en oeuvre l'invention 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 20 comme position initiale estimée P0, puis effectue un aller-retour entre les positions PO et 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 9. La trajectoire reconstituée obtenue est très proche de la trajectoire réelle, connue au 25 préalable pour les besoins de l'expérience, l'écart restant en dessous de 50 cm, comme représenté dans le tableau ci-dessous montrant l'erreur et l'erreur moyenne pour chaque position. 30 3038064 22 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 multiple et des glissements de cycle cumulant des sources d'erreur. Le procédé selon 5 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 horizontale étendue de la précision, ou « extended horizontal dilution of precision » en 10 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 erreurs sont distribuées de manière identique pour chaque émetteur et indépendantes les 15 unes des autres, la dilution horizontale étendue de la précision peut être calculée à partir de la formule : cov(Ex) = (HtH)-1 aT '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 aii les éléments de cov(cx) et Dij les éléments de «11)-1, on obtient : EHDOP = - il + + - D 22 = VOEU. +OE22 - ,UERE 20 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 3038064 23 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 5 des positions différentes du récepteur 2 : m.(1 + n pl . Par exemple, pour m = 3 11' , pt 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. 10 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 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 15 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 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, 20 en utilisant une ou plusieurs réglettes constituées d'un réseau d'une dizaine d'antennes, ou dans des grands centres commerciaux. L'invention peut être utilisée en extérieur, permettant de s'affranchir de station de base. L'expression « comportant un » doit être comprise comme synonyme de 25 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 générateurs (PLI, PL2), chaque générateur émettant sur une même porteuse au moins deux signaux ayant chacun un code différent, et - un récepteur (2) configuré pour détecter les signaux émis par les générateurs (PLI, PL2), procédé dans lequel : - le récepteur (2) mesure, pour chacun des générateurs (PL1, PL2), la différence de phase (&9' ) entre les deux signaux émis par le générateur (PLj), et - en fonction de ces mesures de différence de phase (&9' ), on calcule au moins une grandeur géométrique représentative de la position (xr, yr) du récepteur (2) par rapport aux générateurs (PL1, PL2), afin de localiser le récepteur (2) au sein du système de positionnement (1). 2. Procédé selon la revendication 1, dans lequel les deux signaux émis par chaque générateur (PLj) sont émis à partir de deux zones d'émission (PLji, PLj2) du générateur distantes l'une de l'autre d'une distance prédéfinie (d12). 3. Procédé selon la revendication précédente, dans lequel la grandeur géométrique représentative de la position (xr, yr) du récepteur (2) par rapport à un générateur (PLj) est l'angle d'arrivée (agi) des signaux émis par le générateur (PLj) sur le récepteur (2), fonction de la distance prédéfinie (di)) entre les deux zones d'émission (PLji, PLj2) du générateur (PLj) et des mesures de différence de phase (39' ), et défini par : a = Arc cos( 391 4. Procédé selon la revendication 1 ou 2, dans lequel 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 distances (d1, 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), définies par : d2 - = .\/()C,2 X, )2 ± (Ya2 Yr )2 (x ai - xr )2 (Yal Yr )2 = 89 - 0112 ) 3038064 25 5. Procédé selon l'une quelconque des revendications précédentes, dans lequel les deux signaux émis par chaque générateur (PLj) sont issus d'au moins deux antennes (PLji, PLj2) appartenant au générateur et distantes l'une de l'autre d'une distance prédéfinie (d12). 5 6. Procédé selon la revendication précédente, dans lequel la distance prédéfinie (d12) selon laquelle les deux antennes (PLji, PLj2) d'un même générateur (PLj) sont distantes l'une de l'autre est égale à la longueur d'onde (X) de la porteuse des signaux émis par les antennes (PLji, PLj2). 7. Procédé selon l'une quelconque des revendications précédentes, dans lequel 10 les signaux émis par un même générateur (PLj) sont synchronisés entre eux. 8. Procédé selon l'une quelconque des revendications précédentes, dans lequel chaque générateur (PLj) émet sur une fréquence différente. 9. Procédé selon l'une quelconque des revendications précédentes, dans lequel la fréquence des signaux émis par les générateurs (PL1, PL2) est égale à 1,575 GHz, les 15 signaux étant notamment de type GNSS Galileo. 10. Procédé selon l'une quelconque des revendications précédentes, étant mis en oeuvre en intérieur (indoor). 11. Système de positionnement (1), comportant : - au moins deux générateurs (PLI, PL2), chaque générateur émettant sur une 20 même porteuse au moins deux signaux ayant chacun un code différent, et - un récepteur (2) détectant les signaux émis par les générateurs (PL1, PL2), le système (1) étant configuré pour que : - le récepteur (2) mesure, pour chacun des générateurs (PL1, PL2), la différence de phase (39' ) entre les deux signaux émis par le générateur (PLj), et 25 - en fonction de ces mesures de différence de phase (39' ), soit calculée au moins une grandeur géométrique représentative de la position (xi-, yr) du récepteur (2) par rapport aux générateurs (PL1, PL2), afin de localiser le récepteur (2) au sein du système de positionnement (1). 12. Récepteur (2) destiné à être utilisé au sein d'un système de positionnement 30 (1) comprenant au moins deux générateurs (PL1, PL2), chaque générateur émettant sur une même porteuse au moins deux signaux ayant chacun un code différent, le récepteur (2) détectant les signaux émis par les générateurs (PLI, PL2), 3038064 26 le récepteur (2) étant configuré pour : - mesurer, pour chacun des générateurs (PL1, PL2), la différence de phase ( &p') entre les deux signaux émis par le générateur (PLj), et - en fonction de ces mesures de différence de phase (39' ), que soit calculée 5 au moins une grandeur géométrique représentative de la position (x,, yr) du récepteur (2) par rapport aux générateurs (PL1, PL2), afin de localiser le récepteur (2) au sein du système de positionnement (1). 13. Récepteur selon la revendication précédente, comportant une boucle à verrouillage de phase configurée pour mesurer la phase (0) des signaux émis par les 10 générateurs (PL 1, PL2). 14. Générateur (PL1) destiné à être utilisé au sein d'un système de positionnement (1), ledit système comprenant au moins un autre générateur (PL2) et un récepteur (2) configuré pour détecter les signaux émis par ledit générateur et l'autre générateur (PL1, PL2), ledit générateur (PLI) étant configuré pour émettre sur une même 15 porteuse au moins deux signaux ayant chacun un code différent, le récepteur (2) étant configuré pour mesurer la différence de phase (60) entre lesdits deux signaux émis par le générateur (PLI).