FR 2966001 A1 20120413 FR 1003997 A 20101008 Procédé d'identification d'émetteurs par un terminal dans un réseau iso-fréquence L'invention concerne un procédé d'identification d'émetteurs par un terminal dans un réseau iso-fréquence. Elle s'applique notamment au domaine de la géo-localisation. Les terminaux mobiles de télécommunications, comme par exemples les téléphones portables ou les ordinateurs portables, sont utilisés 1 o de fréquemment pour la mise en oeuvre d'applications requérant des informations précises de localisation. Ainsi, il est habituel que ces terminaux comportent des moyens leur permettant d'estimer le plus précisément possible leur position. Des systèmes de géo-localisation par satellites sont particulièrement utilisés. On peut citer à titre d'exemple le système GPS, 15 acronyme venant de l'expression anglo-saxonne « Global Positionning System ». Pour qu'un système de géo-localisation par satellites fonctionne correctement, il est requis qu'un nombre suffisant de satellites soient visibles par le terminal. Le choix du nombre et de la position des satellites dans la constellation d'un tel système permet de satisfaire à cette contrainte. II existe 20 néanmoins des situations dans lesquelles le terminal ne peut recevoir des signaux en provenance d'un nombre suffisant de satellites. C'est par exemple le cas dans des centres urbains comprenant des bâtiments de grande taille empêchant la visibilité des satellites par le terminal ou quand le terminal se trouve dans des zones partiellement ou totalement couvertes. 25 D'autres techniques de géo-localisation peuvent être utilisées et s'appuient notamment sur les caractéristiques des formes d'ondes et des réseaux d'accès de systèmes radio-mobiles tels que le GSM et l'UMTS ou bien s'appuyant sur des systèmes de télécommunications sans fil tels que le Wifi ou le WiMax. Cependant, la précision d'estimation est moins bonne que 30 celle obtenue par les techniques existantes de géo-localisation par satellites. Les réseaux iso-fréquence, appelés aussi réseaux SFN, acronyme venant de l'expression anglo-saxonne « Single Frequency Network » sont plus rarement utilisés dans le cadre de la géo-localisation. Un réseau SFN est un réseau comprenant une pluralité d'émetteurs et de récepteurs. Dans 35 ce type de réseaux, chaque émetteur transmet le même signal de manière synchronisée et sur une même fréquence porteuse. En l'absence de trajets multiples, un récepteur reçoit plusieurs répliques du même signal provenant de plusieurs émetteurs. Les différentes répliques sont affectées d'une atténuation et d'un retard différents, ladite atténuation et ledit retard dépendant notamment de la distance entre le terminal récepteur et chacun des émetteurs. Ces caractéristiques des réseaux SFN rendent difficiles pour un terminal récepteur d'identifier la provenance d'une réplique du signal, c'est-à-dire d'identifier l'émetteur. La localisation du terminal en utilisant des signaux transmis par les émetteurs d'un réseau est en principe possible sous réserve de pouvoir localiser et identifier lesdits émetteurs. Le tatouage numérique est une technique connue permettant notamment aux terminaux d'identifier les émetteurs d'un réseau SFN. Cette technique est souvent désignée par le mot anglais « watermarking ». Le principe de cette technique est qu'un signal d'identification de l'émetteur est superposé au signal comportant les données applicatives à destination de l'utilisateur du terminal récepteur. Le signal comportant ces données utilisateur est appelé dans la suite de la description signal utile. Lesdites données utilisateur correspondent par exemple à des flux vidéo associés à des chaînes de télévision. Le signal d'identification est habituellement superposé au signal utile avec une puissance relative très faible, ce qui ne perturbe pas la réception des données utilisateur. Cependant, le signal d'identification est brouillé lorsque l'on s'approche d'un autre émetteur. Cela correspond au phénomène bien connu de l'homme du métier et désigné par l'expression anglo-saxonne « near-far effect ». Dans ce cas, seuls les émetteurs les plus proches peuvent être identifiés. Comme explicité précédemment, pour estimer la position d'un terminal, il est requis d'identifier autant d'émetteurs que possible. Un but de l'invention est notamment de pallier les inconvénients 30 précités. A cet effet l'invention a pour objet un procédé d'identification d'émetteurs par un terminal dans un réseau iso-fréquence comprenant une pluralité d'émetteurs, lesdits émetteurs étant synchronisés et émettant avec un retard artificiel ri propre à chaque émetteur. Le procédé comporte au 35 moins : une étape d'acquisition de la position approximative du terminal p , de la position pi d'une liste d'émetteurs {Tx} au voisinage du terminal et des retards des retard Ti leurs étant associés ; une étape de mesures de pseudo-distances pi entre les émetteurs et le terminal ; une étape d'association des mesures p; aux émetteurs de positions connues pi en minimisant une fonction de coût, ladite fonction de coût v(pi, p,a) correspondant à la norme de l'erreur entre les mesures pi et un modèle de mesures des pseudo-distances appliqué à une permutation de la position des émetteurs a . Selon un aspect de l'invention, le procédé comporte une étape 15 d'estimation de la position affinée du terminal en utilisant les positions des émetteurs et les pseudo-mesures leurs étant associées. Les positions p; et les retards artificiels Ti des émetteurs peuvent être mémorisés dans une base de données localisée dans le terminal. Alternativement, les positions pi et les retards artificiels Ti des 20 émetteurs peuvent être transmises au terminal en utilisant des moyens de télécommunications. L'acquisition de la position approximative p du terminal est, par exemple, le résultat d'estimations réalisées sur des systèmes radio-mobiles tels que le GSM ou l'UMTS. 25 II est aussi possible que la position approximative du terminal récepteur corresponde à une position acquise par GPS. La fonction de coût correspond par exemple à l'expression suivante : CüY v(pi, p,a)= PI =1 I 11 (x-xi )2 +(ÿ-y~)2 -cXzi ) . II 30 dans laquelle : pi représente la i-ème mesure de pseudo-distances ; p = [x, ÿ] représente la position approximative du terminal ; 20 25 30 a représente l'ensemble des combinaisons possibles des émetteurs de la liste {Tx} ; card(a) représente le cardinal de a ; c est une constante représentant la vitesse de la lumière ; Ti correspond au délai artificiel associé au i-ème émetteur. Selon un aspect de l'invention ; la meilleure combinaison émetteurs-mesures 6 est déterminée en utilisant l'expression : 6 = arg min[v(p, p,a)] 6 dans laquelle : arg min( ) est une fonction indiquant que la recherche du minimum ne se fait k que sur la variable k. L'invention a aussi pour objets un terminal de télécommunication comprenant des moyens pour recevoir des signaux émis pas les émetteurs d'un réseau iso-fréquence. Ce terminal comporte des moyens pour mettre en oeuvre le procédé d'identification d'émetteurs tel que décrit précédemment. D'autres caractéristiques et avantages de l'invention apparaîtront à l'aide de la description qui suit donnée à titre illustratif et non limitatif, faite en regard des dessins annexés parmi lesquels : la figure 1 illustre le principe d'identification des émetteurs et de géo-localisation d'un terminal dans un réseau iso-fréquence selon l'invention ; la figure 2 donne un exemple simple d'application du procédé selon l'invention lorsqu'il n'y a pas d'erreurs d'estimations ; la figure 3 donne un exemple simple d'application du procédé selon l'invention en présence d'erreurs d'estimations. La figure 1 le principe d'identification des émetteurs et de géolocalisation d'un terminal dans un réseau iso-fréquence selon l'invention. Afin d'éviter que deux répliques de signal en provenance d'un même émetteur arrivent en même temps au niveau du terminal récepteur, le procédé selon l'invention peut être mis en ceuvre pour un réseau SFN dans lequel les émetteurs sont synchronisés, mais où chaque émetteur transmet le signal avec un retard différent, une valeur prédéfinie de retard étant propre à chacun desdits émetteurs. Ces retards sont appelés retards ou délais artificiels dans la suite de la description. Le concept de retards artificiels a été introduit dans la demande de brevet FR1001868. Afin d'estimer sa position, un terminal réalise des mesures de pseudo-distances pi à partir de signaux reçus en provenance des différents émetteurs du réseau iso-fréquence. Ces pseudo-distances mesurées 1 o peuvent être décrites à l'aide de l'expression suivante : pi = h(p, pi , z~ ) + n, avec i E {Tx} (1) Dans laquelle : 15 p = [x, y] correspond à la position exacte du terminal récepteur ; p = [x; , yi ] correspond à la position du i-ème émetteur ; Ti correspond au délai artificiel associé au i-ème émetteur ; n, correspond au bruit affectant la mesure de pseudo-distance ; h(p, pi ,Ti ) = -\j(x - x, )z +(y - y, )z + c x z, est la fonction de modélisation de la 20 pseudo-distance ; {Tx} correspond à la liste des émetteurs pour lesquels une pseudo distance est mesurée. Sur la base des mesures de pseudo-distances pi seules, il n'est pas possibles d'associer les signaux reçus aux différents émetteurs. En 25 d'autres termes, il n'est pas possible d'identifier quel émetteur a émis une réplique de signal reçue par le terminal récepteur. Afin d'identifier les émetteurs, le procédé utilise d'une part les mesures de pseudo-distances pi 100 et d'autre des données d'assistance 101 permettant d'identifier la configuration d'émetteurs correspondant à la 30 collection de pseudo-distances mesurées. En effet, de par l'introduction de retards artificiels, une seule combinaison émetteurs/pseudo-distances est possible pour une position donnée du terminal. Chaque mesure de pseudodistance est ainsi associée sans ambiguïté à un émetteur. II est alors possible de localiser avec précision le terminal par un simple calcul bien connu de l'homme du métier, les positions des émetteurs et les pseudodistances leurs étant associées étant connues du terminal. Plusieurs types de données d'assistance 101 sont utilisés. Ainsi une position approximative du récepteur et la position des émetteurs au voisinage du terminal sont utilisées par le procédé. La fourniture de ces données d'assistance peut être réalisée de différentes manières. Par exemples, les positions et les retards artificiels des émetteurs du réseau isofréquence peuvent être mémorisés dans une base de données accessible par le terminal ou bien lui être transmis par voie radio. Pour ce qui est de la position approximative du terminal récepteur, elle peut être le résultat d'estimations réalisées en GSM ou en UMTS. Celle-ci peut aussi être obtenue en utilisant la dernière position estimée par un module GPS avant que celui-ci ne soit plus capable de fonctionner correctement, par exemple lorsqu'il ne dispose pas d'un nombre suffisant de signaux satellites pour effectuer ses estimations. En résumé, un récepteur dispose alors d'un ensemble de mesures de pseudo-distances 100. Il connaît sa position approximative ainsi que les caractéristiques des émetteurs qui sont dans son voisinage 101, c'est-à-dire leur position ainsi que les valeurs de retards artificiels leurs étant associés. Sur la base de ces données, pour chaque combinaison d'associations possibles entre les émetteurs et les mesures de pseudodistances, le récepteur détermine un coût, ce coût étant obtenu à l'aide d'une fonction de coût v( ). Le récepteur sélectionne ensuite l'une des combinaisons 25 d'émetteurs possibles. Pour cela, la combinaison d'émetteurs associée au coût minimum est sélectionnée. La fonction de coût est égale, par exemple, à la norme de l'erreur entre les mesures de pseudo-distances et un modèle de mesures appliqué à une permutation de la position des d'émetteurs. 30 Une fonction de coût pouvant être utilisée dans le cadre de l'invention est donnée par l'expression suivante : ca v(pi , p,a)= =1~ - ~(x - xi )2 + (ÿ - y )2 - c x Ti l I (2) dans laquelle : p = [x, ÿ] représente la position approximative du terminal ; a représente une des combinaisons possibles des émetteurs de la liste {Tx} ; card(a) représente le cardinal de a ; c est une constante représentant la vitesse de la lumière. La position approximative du terminal p peut être obtenue par d'autres moyens de positionnement. Elle correspond par exemples à une estimation GPS antérieure réalisée avant d'entrer dans un bâtiment ou à une estimation sur la base de signaux GSM ou UMTS. A titre de remarque, si la fonction de coût v(pt, p,a) est calculée pour la bonne combinaison 6,,, lorsque la position approximative du terminal est égale à sa position réelle, c'est-à-dire lorsque p = p , et en prenant l'hypothèse que les pseudo-distances sont estimées sans erreur, le coût vaut 15 v(pi, p,a)=0. L'estimation de la meilleure combinaison émetteurs-mesures 6 revient à un problème de minimisation de la fonction de coût v(pi, p,a) tel que : 20 6 = argmin[v(pi, p,a)] (3) 6 expression dans laquelle : arg min( ) est une fonction indiquant que la recherche du minimum ne se fait k que sur la variable k. 25 La figure 2 donne un exemple simple d'application du procédé selon l'invention lorsqu'il n'y a pas d'erreur d'estimation. Dans cet exemple, deux émetteurs Tx1 et Tx2 appartenant à un réseau iso-fréquence sont visibles d'un terminal récepteur 201. Le premier émetteur Tx1 est à une 30 distance d, du terminal et le second émetteur est à une distance dz du terminal. Dans cet exemple, il est supposé que la position estimée du terminal est égale à la position vrai, c'est-à-dire p = p. Il est également pris comme hypothèse que les mesures de pseudo-distances p, et p2 sont obtenues sans erreur d'estimation. Dans cet exemple, il y a deux combinaisons d'émetteurs appelées 6, et u2 à tester : 6, correspond à la combinaison : émetteur Tx1 associé à la pseudo-distance mesurée p, et émetteur Tx2 associé à la pseudo-distance mesurée p2. 62 correspond à la combinaison : émetteur Tx2 associé à la pseudo-distance mesurée p, et émetteur Tx1 associé à la pseudo-distance 10 mesurée p2. En prenant les valeurs suivantes : - d, = 500 m, d2 = 200, - positions des émetteurs : x, _ -400, x2 = 200, y, = y2 = 0 , - position réelle du terminal : x= y = 0 m, 15 - position estimée du terminal x = 0 m et = 0 m car p = p , - c x z, = 100 et c x a2 = 0, - p, = d2 et p2 = d,, on obtient alors les deux coûts suivants déterminés en utilisant la fonction de coût telle que définie par l'expression (2) : 20 v(pi , p, 61 ) = 180000 et v(pi , p, 62) = 0 (4) La fonction de coût est minimisée pour la combinaison 62 , les mesures de pseudo-distances sont alors associées sans ambiguïté aux différents 25 émetteurs du réseau iso-fréquence. Il est alors possible de déterminer avec précision la position du terminal. On peut noter qu'il est possible d'utiliser l'identification des émetteurs dans un réseau iso-fréquence avec un but autre que celui de positionner des terminaux. A titre d'exemple, cette identification peut être 30 utilisée dans le cadre de la surveillance d'interférences intra-système dans un réseau iso-fréquence à retards artificiels. Cela permet ainsi de déterminer quel émetteur brouille la réception TV dans une zone donnée. La figure 3 donne un exemple simple d'application du procédé 35 selon l'invention en présence d'erreurs d'estimations. Dans cet exemple, deux émetteurs Tx1 et Tx2 appartenant à un réseau iso-fréquence sont visibles d'un terminal récepteur 301. Le premier émetteur Tx1 est à une distance d', du terminal et le second émetteur est à une distance d'2 du terminal. La position estimée du terminal 302 n'est pas exacte et est différente de la position vrai, c'est-à-dire p p. II est également pris comme hypothèse que les mesures de pseudo-distances p, et p2 sont obtenues avec des erreurs d'estimation. Comme dans l'exemple de la figure 2, il y a deux combinaisons 10 d'émetteurs 6, et 62 à tester : ul correspond à la combinaison : émetteur Tx1 associé à la pseudo-distance mesurée p, et émetteur Tx2 associé à la pseudo-distance mesurée p2. 62 correspond à la combinaison : émetteur Tx2 associé à la 15 pseudo-distance mesurée p, et émetteur Tx1 associé à la pseudo-distance mesurée p2. En prenant les valeurs suivantes : 20 on obtient alors les deux coûts déterminés en utilisant la fonction de coût de 25 l'expression (2) : v(pi, p, ul ) = 212500 et v(pi, p, 62 ) = 1700 (5) Malgré les erreurs d'estimation, la fonction de coût est minimisée pour la 3o combinaison 62. A titre d'exemple, il est possible d'estimer la performance de ce procédé en calculant la différence entre la seconde valeur de coût la plus faible et la valeur de coût la plus faible, c'est-à-dire celle correspondant à la 35 combinaison retenue. Lorsque cette différence est maximisée, cela signifie - d, = 500 m, d2 = 200 m, - positions des émetteurs : x, = -400, x2 = 200, y, = y2 = 0 , - position réelle du terminal : x= y= 0 m, - position estimée du terminal x = 20 m, y' = 0 m, - c x zl = 100 et c x z2 = 0, - pl = 220 m et p2 = 530 m. que les combinaisons émetteurs-mesures sont identifiées avec un faible taux d'erreur. Cette différence est maximisée lorsque les retards artificiels sont choisis de manière à rejeter les zones d'ambigüité en dehors de la couverture des émetteurs. En d'autres termes, l'introduction de retards artificiels rend la signature des mesures du récepteur unique par rapport à une position possible. Avantageusement, ce procédé est robuste même lorsque les données d'assistance sont de qualité médiocre, comme par exemple lorsque la position approximative des récepteurs comprend une erreur de l'ordre io d'une centaine de mètres. REVENDICATIONS1- Procédé d'identification d'émetteurs par un terminal dans un réseau iso-fréquence comprenant une pluralité d'émetteurs, lesdits émetteurs étant synchronisés et émettant avec un retard artificiel Ti propre à chaque émetteur, le procédé étant caractérisé en ce qu'il comporte au moins : une étape (100) d'acquisition de la position approximative du terminal p , de la position pi d'une liste d'émetteurs {Tx} au voisinage du terminal et des retards des retard z, leurs étant associés ; une étape (101) de mesures de pseudo-distances p, entre les émetteurs et le terminal ; une étape (102) d'association des mesures p; aux émetteurs de positions connues p; en minimisant une fonction de coût, ladite fonction de coût v(p,, p,6) correspondant à la norme de l'erreur entre les mesures p; et un modèle de mesures des pseudo- distances appliqué à une permutation de la position des émetteurs 6. 2- Procédé selon la revendication 1 caractérisé en ce qu'il comporte une étape d'estimation de la position affinée du terminal en utilisant les 25 positions des émetteurs et les pseudo-mesures leurs étant associées. 3- Procédé selon l'une quelconque des revendications précédentes caractérisé en ce que les positions pi et les retards artificiels z, des émetteurs sont mémorisées dans une base de données localisée 30 dans le terminal. 4- Procédé selon l'une quelconque des revendications 1 ou 2 caractérisé en ce que les positions pi et les retards artificiels ri des émetteurs sont transmises au terminal en utilisant des moyens de 35 télécommunications. 15 20 5- Procédé selon l'une quelconque des revendications précédentes caractérisé en ce que l'acquisition de la position approximative p du terminal est le résultat d'estimations réalisées sur des systèmes radio- mobiles tels que le GSM ou l'UMTS. 6- Procédé selon l'une quelconque des revendications 1 à 4 caractérisé en ce que la position approximative du terminal récepteur correspond à une position acquise par GPS. 7- Procédé selon l'une quelconque des revendications précédentes caractérisé en ce que la fonction de coût correspond à l'expression suivante : carda 2 v(pi , p, 6) = ( pi - ~(x - xi )2 + (ÿ - yi )2 - c x Ti I . i=1 dans laquelle : p, représente la i-ème mesure de pseudo-distances ; p = [.z, ÿ] représente la position approximative du terminal ; 20 a représente l'ensemble des combinaisons possibles des émetteurs de la liste {Tx} ; card(6) représente le cardinal de 6 ; c est une constante représentant la vitesse de la lumière ; z, correspond au délai artificiel associé au i-ème émetteur. 25 8- Procédé selon l'une quelconque des revendications précédentes caractérisé en ce que la meilleure combinaison émetteurs-mesures 6 est déterminée en utilisant l'expression : 30 6 = arg min[v(p, p, a)] 6 dans laquelle : arg min( ) est une fonction indiquant que la recherche du minimum ne k se fait que sur la variable k. 15 9- Terminal de télécommunication comprenant des moyens pour recevoir des signaux émis pas les émetteurs d'un réseau iso-fréquence, caractérisé en ce qu'il comporte des moyens pour mettre en ceuvre le procédé d'identification d'émetteurs selon l'une quelconque des revendications précédentes.