FR 3031194 A1 20160701 FR 1463374 A 20141229 3031194 PROCEDE ET DISPOSITIF DE DETECTION ET DE LOCALISATION EN INTERIEUR DE PERSONNES OU D'OBJETS Domaine de l'invention L'innovation se situe dans le domaine de la détection et de la localisation de personnes à l'intérieur d'un bâtiment, en particulier à l'aide d'ondes radio. Etat de la Technique Il existe deux grandes approches pour localiser une personne, selon que la personne dispose ou non d'un moyen électronique de communication, tel qu'un Smartphone par exemple. Si c'est le cas l'approche est dite à «système actif », sinon elle est dite à « système passif ». De nombreuses solutions de localisation pour système actif sont connues. Elles utilisent un réseau d'antennes disposées dans la zone de localisation, telles les antennes des points d'accès Wifi. Certaines techniques utilisent les signaux reçus à partir de l'antenne de l'usager pour déterminer la distance séparant l'usager des antennes réparties et pour déterminer par triangulation la position de l'usager. D'autres techniques vont utiliser une base de données constituée au préalable des empreintes (« fingertyping » en anglais) des intensités des ondes reçues sur les différentes antennes distribuées dans la zone, en provenance d'une antenne du même type que celle de l'usager. La position de l'usager sera déterminée ensuite en comparant les signaux émis par celui -ci et ceux figurant sur la base de données. Il existe par contre peu de solutions de localisation à système passif. Ces techniques vont utiliser l'impact de la présence de la personne ou de l'objet sur les signaux qu'un ensemble d'antennes émettent en 2 303 1 1 9 4 direction d'un autre groupe d'antennes. La difficulté de la localisation précise d'une personne (ou d'un objet) non équipée d'un moyen électronique de communication, située dans un bâtiment ou dans un espace clos ou semi-clos, résulte de la faible variation de ces signaux 5 surtout lorsque des murs séparent la personne des antennes. Par ailleurs, l'emplacement et le nombre d'objets dans une pièce peuvent évoluer. Les objets peuvent être déplacés, enlevés, d'autres objets peuvent être ajoutés, le tout modifiant la configuration et la structure initiale d'un environnement. Ces changements dans le temps 10 de l'environnement vont poser problème car ils induisent sur les signaux collectés par les antennes des variations à l'instar de celles provoquées par la personne à localiser. Certaines techniques utilisent des antennes disposées à la périphérie de la zone à surveiller et mesurent les variations d'intensité sur 15 des liaisons point à point en présence de la personne à détecter. Il est possible de trouver dans l'art antérieur des solutions pour localiser des personnes ou des objets. Ainsi, dans la demande de brevet WO 2011/129983 A2 de R. W. Schlub, un dispositif électronique permet de déterminer si une personne 20 ou un objet est à une distance donnée d'une antenne, en faisant des mesures d'impédance d'antenne. Cette solution est adaptée quand la personne ou l'objet est à proximité de l'antenne. L'article de Chang et al. intitulé People « Tracking with UWB Radar Using a Multiple-Hypothesis Tracking of Clusters (MHTC) Method » publié 25 en 2010 divulgue des techniques d'ondes radio à large spectre pour la détection d'humains en mouvement, dans une configuration avec une seule antenne en émission et en réception. Il est aussi connu un projet du Massachussetts Institue of Technology intitulé WiTrack qui a donné lieu à plusieurs publications de 30 F.Adib et al. tel l'article intitulé « 3D Tracking via Body Radio Reflections » ou l'article « See Through Walls with Wi-Fi ! ». Le procédé qui est 3 3031194 présenté mesure à des temps réguliers des « Time-of-Flight » (TOF) en utilisant une technique dite `FMCW' pour « Frequency Modulated Carrier Wave » basée sur le calcul des temps mis par les ondes à parcourir la distance entre les antennes et les obstacles. Puis par comparaisons 5 itératives des temps, un modèle géométrique de l'environnement est établi. Cependant, le procédé requiert que la personne soit en mouvement pour permettre sa localisation, et est limité à la détection derrière un seul mur. Il existe alors le besoin d'une solution qui pallie aux inconvénients 10 des approches connues. La présente invention répond à ce besoin. Résumé de l'invention Un objet de la présente invention est de proposer un dispositif et un procédé de détection et de localisation en intérieur de personnes ou 15 d'objets non équipés de moyen électronique de communication. La solution proposée utilise deux réseaux linéaires d'antennes, un en émission et un en réception, et une méthode de calcul pour déterminer la position de la/des personne(s) ou objet(s), qui utilise les intensités et les phases des éléments des matrices de transfert à plusieurs entrées, 20 plusieurs sorties, dites matrices MIMO « Multiple Input Multiple Output ». Chaque élément de matrice est obtenu à fréquence fixe en faisant le rapport entre le voltage reçu sur une des antennes de réception et le voltage transmis par une des antennes d'émission. Avantageusement, la solution proposée utilise les effets de la 25 présence d'un corps ou d'une masse (humain ou non) sur les canaux des matrices MIMO, pour déterminer la position de celui-ci, même dans les cas où les canaux peuvent être partiellement obstrués par des murs. Avantageusement, la solution proposée permet de déterminer la position d'une personne ou d'un objet qui peut être relativement éloignés 30 du dispositif de localisation, de 1 mètre à plusieurs dizaines de mètres, 4 3031194 dans un espace intérieur clos ou semi-clos pouvant comporter un ensemble de pièces. Avantageusement, le dispositif de l'invention utilise des réseaux d'antennes à une seule fréquence. Toujours avantageusement, les deux 5 réseaux d'antennes peuvent se trouver dans des pièces différentes. La solution proposée s'appuie sur un résultat d'étude menée par les inventeurs qui a permis de démontrer qu'il existe une relation entre une matrice de transfert MIMO et les propriétés géométriques d'objets (leur position, leur taille) constituant un environnement. De manière 10 avantageuse, cette relation est appliquée à la détermination des coordonnées spatiales d'une personne ou d'un objet situé dans un bâtiment, en effectuant une seule mesure à fréquence fixe de la matrice de transfert MIMO. De manière générale, la méthode de détection et de localisation 15 proposée consiste premièrement dans un espace intérieur clos ou semiclos vide c'est-à-dire en l'absence d'une personne ou d'un objet à localiser, à déterminer une première matrice de transfert MIMO pour un ensemble de NT antennes en émission et NR antennes en réception, les antennes étant disposées en deux réseaux d'antennes linéaires distants 20 l'un de l'autre et où à l'intérieur de chaque réseau, les antennes sont distantes d'une longueur d'onde divisée par deux. Ensuite, une pluralité de matrices de transfert sont calculées à intervalles réguliers `kT'. A partir de ces relevés, une fonction de localisation h(uR, uT ;k) périodique par rapport à des variables `LIR' et 'UT' est définie, les coefficients de Fourier 25 de cette fonction étant les éléments des matrices de transfert. Le maximum de la fonction est ensuite déterminé. Les coordonnés spatiales x(kT), y(kT) de la position de la personne ou de l'objet à un temps `kT' donné sont calculées à partir de ce maximum. La présente invention trouvera une application avantageuse dans les domaines de la téléassistance, pour les personnes âgées vivant 5 3031194 seules, pour la surveillance des enfants ou encore dans des applications dites de « vision à travers les murs ». D'autres applications sont envisageables pour la détection d'intrusion dans les bâtiments ou encore pour l'éclairage des pièces ou de bureaux, en fonction de la présence de 5 personnes. Pour obtenir les résultats recherchés, un procédé et un dispositif sont proposés. En particulier, un procédé pour déterminer l'ajout de personnes ou d'objets dans un espace intérieur est proposé, l'espace 10 étant clos ou semi-clos et équipé d'un premier réseau d'antennes en émission distant d'un second réseau d'antennes en réception, le procédé comprenant les étapes suivantes : - mesurer une valeur initiale de champ électromagnétique créé par les réseaux d'antennes pour l'espace intérieur ; 15 - mesurer à intervalles successifs des valeurs du champ électromagnétique au voisinage du réseau d'antennes en réception ; et - comparer les valeurs mesurées à la valeur initiale pour déterminer une différence, ladite différence étant représentative d'un ajout dans l'espace intérieur. 20 Dans un mode de réalisation préférentiel, les étapes de mesures du champ électromagnétique consistent à déterminer et enregistrer une matrice de transfert MIMO initiale constituée par les deux réseaux d'antennes pour l'espace intérieur vide, puis à mesurer des matrices de 25 transfert MIMO successives, et l'étape de comparaison consiste à déterminer s'il existe des différences entre les éléments de la matrice de transfert initiale et une matrice ultérieure. Avantageusement, à l'intérieur de chaque réseau, les antennes 30 sont distantes d'une longueur d'onde divisée par deux. 6 3031194 Dans une variante, le procédé comprend de plus, une étape de détermination de la position de la/les personne(s) ou objet(s) s'il existe une différence entre les éléments de la matrice de transfert initiale et une matrice ultérieure. 5 Dans un mode de réalisation, l'étape de détermination de la position de la/les personne(s) ou objet(s) consiste à définir une fonction de localisation dont les coefficients de Fourier sont les éléments des matrices de transfert, et à calculer le/les maxima de ladite fonction. 10 Avantageusement, le procédé comprend de plus une étape de calcul des coordonnés spatiales de ladite personne ou dudit objet à un temps `kT' à partir du maximum de la fonction de localisation. 15 Avantageusement, les intervalles successifs de mesure sont dans une plage de 1 à 10 secondes. L'invention couvre aussi un système pour déterminer l'ajout de personnes ou d'objets dans un espace intérieur clos ou semi-clos, le 20 système comprenant des moyens pour mettre en oeuvre les étapes du procédé. Avantageusement, le système comprend un circuit d'émission distant d'un circuit de réception, le circuit d'émission comprenant au 25 moins un réseau linéaire d'antennes en émission, le circuit de réception comprenant au moins un réseau linéaire d'antennes en réception, les antennes à l'intérieur de chaque réseau étant distantes d'une longueur d'onde divisée par deux. 7 3031194 Dans une implémentation, le circuit de réception comprend de plus une unité de traitement pour enregistrer les mesures et effectuer les calculs des matrices de transfert et des coordonnées de position. 5 Avantageusement, une partie des éléments de l'invention peut opérer sous la forme d'un produit programme d'ordinateur qui comprend des instructions de code permettant d'effectuer les étapes du procédé revendiqué. En particulier, l'unité de traitement peut comprendre un processeur apte à effectuer les calculs après les étapes de mesures. 10 Description des figures Différents aspects et avantages de l'invention vont apparaitre en appui de la description d'un mode préféré d'implémentation de l'invention 15 mais non limitatif, avec référence aux figures ci-dessous : La figure 1 illustre un environnement pour implémenter le dispositif de l'invention dans un mode de réalisation; La figure 2 illustre un réseau d'antennes en émission et réception selon un mode de réalisation de l'invention ; 20 La figure 3 illustre un enchaînement des étapes du procédé de localisation de l'invention selon un mode de réalisation ; La figure 4 illustre les blocs fonctionnels des circuits d'émission et de réception dans un mode de réalisation de l'invention. 25 Description détaillée de l'invention Référence est faite à la figure 1 qui illustre de manière schématisée un environnement (100) pour implémenter et opérer la présente invention. 8 3031194 Un ensemble de pièces ou bureaux (102-1 à 102-10) sont distribués de part et d'autre d'un couloir central (103). Chaque pièce peut être vide ou contenir du mobilier (104) de forme et taille différente, et positionné à des endroits identiques ou différents dans chaque pièce. Il est à noter que la 5 configuration présentée est simplifiée pour faciliter la description de l'invention, mais n'est en rien limitative et que toute variation sur le nombre, la forme, l'emplacement, tant des pièces, des parois ou du mobilier, peut être considérée pour opérer l'invention. Ainsi, l'espace intérieur peut être clos ou semi-clos, tel qu'un espace n'ayant pas de 10 plafond ou ne présentant qu'un nombre limité de parois latérales. L'invention s'applique avantageusement dans un espace intérieur considéré semi-clos qui peut ne comporter que quelques parois verticales ou des piliers, pouvant avoir des parois vitrées sur l'extérieur. L'invention s'applique à tout espace qui présente des éléments de paroi ayant une 15 influence sur la propagation des ondes radioélectriques, tels les murs, plafond, plancher ou mobilier. Un réseau (106) d'antennes de transmission TX et un réseau (108) d'antennes de réception RX sont respectivement disposés aux extrémités 20 du couloir, un premier groupe (106) constituant `NT' antennes de transmission TX, un second groupe (108) constituant `NR' antennes de réception RX. Le nombre d'antennes de transmission et de réception est de manière préférentielle le même, mais l'homme de l'art pourra adapter la méthode de calcul décrite à un nombre d'antennes en émission et en 25 réception qui soit différent. Ainsi pour des raisons de clarté et de simplification de la description, le nombre d'antennes en émission et en réception est pris identique et égal à 'N'. De même ce nombre peut être pair ou impair sans que cela impacte le principe de la méthode. Dans une implémentation standard pour des appartements ou un ensemble de 30 bureaux de diamètre inférieur à 30 mètres, N peut être égal à 5 ou 7 antennes. 9 3 03 1 194 Dans un mode préférentiel, les antennes sont de type omnidirectionnel, par exemple des antennes dipôle, à fréquence fixe de l'ordre du Gigahertz. Dans l'exemple décrit, les réseaux d'antennes sont linéaires et les directions des réseaux d'émission et de réception sont 5 dans un plan parallèle au sol. Pour simplifier les calculs décrits ci-après, les deux réseaux sont positionnés à la même hauteur et leur direction est horizontale. Les antennes se situent de manière préférentielle à une hauteur du sol d'au moins 1 mètre et à une distance du plafond supérieure à 1 m. 10 La figure 1 illustre le cas de deux personnes (P1, P2) situées dans deux bureaux différents, respectivement 102-3 et 102-9, et qui ne sont pas en vue directe des antennes. Tel qu'il va être décrit plus loin, l'invention permet de détecter la présence des deux personnes et de déterminer leur position au sein du bâtiment. 15 Référence est faite à la figure 2 qui illustre une vue théorique des réseaux d'antennes en émission TX et en réception R. Dans la suite de la description, les éléments suivants désignent : `0' pour l'origine des coordonnées et le milieu du réseau des antennes d'émission TX ; 20 'I' pour le point correspondant au milieu du réseau des antennes de réception RX ; `aT' pour le vecteur correspondant au pas du réseau TX, entre deux antennes de transmission ; `aR' pour le vecteur correspondant au pas du réseau RX, entre 25 deux antennes de réception ; 'e-r' pour l'angle entre la direction de aT et une ligne arbitraire issue de O et qui est aussi définie par une direction azimutale 1:1-r' ; et '6)F1' pour l'angle entre la direction de aR et une ligne arbitraire issue de I et qui est aussi définie par une direction azimutale '4)R'. 10 3 0 3 1 1 9 4 Le procédé de l'invention se base sur la mesure des éléments d'une matrice de transfert `H°' d'un système MIMO constitué du réseau des 'NT' antennes de transmission TX et du réseau des `KIR' antennes de réception RX. Un élément Hmn de la matrice de transfert, est obtenu en 5 mesurant la tension complexe vm, obtenue sur la m-ième antenne RX lorsqu'une tension complexe vn est appliquée sur la n-ième antenne TX, et en effectuant leur rapport, tel que : Hmn = vm/vn. Il a été établi empiriquement que les éléments d'une matrice de transfert 'H' sont donnés par l'expression suivante : 10 Hmn = JS2xS2 C(OR)t A(QR, SET) C(QT) dQR dOT où : S2 représente la sphère de rayon unité ; QT= (4)T, 0T) est une direction issue du milieu du réseau TX (4)T, eT étant respectivement les angles d'azimut et d'élévation) ; 15 OR= (4)R, OR) est une direction issue du milieu du réseau RX ; c(4T) and c(4R) sont respectivement des matrices 2 x NT and 2 x NR dont les éléments dépendent seulement des diagrammes d'antennes correspondants et de la géométrie des réseaux ; A(OR,OT) est une matrice 2 x 2 , dite fonction d'étalement, qui 20 dépend seulement de l' environnement. Tel que décrit dans les articles de F. Bentosela et al. « On the transfer matrix of a MIMO system » dans Mathematical Methods in the Applied Sciences, vol. 34, no. 8, 963-976, 2011, et « Influence of 25 environment richness on the increase of MIMO capacity with number of antennas » dans IEEE Transactions on Antennas and Propagation, Vol. 62, No. 7, 2014, p. 3786-3796, il existe un lien entre la fonction A(CIR,QT) et les propriétés physiques, comme la conductivité, la permittivité, et les propriétés géométriques de l'environnement. 30 11 3031194 Selon les principes de base de la propagation des ondes électromagnétiques dans un environnement physique, le mouvement des électrons pulsé par un oscillateur dans une antenne d'émission TX, va induire à l'intérieur des objets environnants et de l'antenne de réception 5 RX un mouvement de charges électriques. Le champ total vu par l'antenne RX est la somme des champs crées par tous les courants: ceux dans l'antenne TX qui créent les champs dits incidents (Ei, Hi) et ceux induits dans les objets de l'environnement qui vont créer à leur tour les champs dits diffractés (Es, Ha). Pour déterminer ces derniers, tel que cela 10 est connu de l'homme de l'art, il est montré qu'il existe un opérateur intégral 'T' qui transforme les champs incidents sur l'objet en courants dans l'objet et qui dépend à la fois de la géométrie de l'objet et de ses caractéristiques électriques (permittivité, conductivité). Une fois ces courants déterminés, l'objet est considéré comme se comportant comme 15 une antenne et les champs diffractés (Es, Hs) se calculent en n'importe quel point de l'espace, de manière usuelle, par des intégrales dont l'intégrant contient la densité de courants. Dans les articles précédemment cités, les auteurs ont justifié 20 mathématiquement la procédure de calcul ci-dessus et montré le lien entre la fonction A et le noyau de l'opérateur T. Il est expliqué que si on considère les angles solides issus des centres (O, I) des réseaux d'antennes TX et RX et dont l'enveloppe est tangente à l'objet, comme cela est illustré sur la figure 2, alors l'ensemble des directions CIR,CIT à 25 l'intérieur des angles solides forme le support de la fonction A(SIR,C2T). Le support d'une fonction est l'ensemble des valeurs des variables pour lesquelles la fonction ne prend pas la valeur zéro. Il est à noter que si "l'objet " entoure complètement les antennes comme c'est le cas si celles-ci sont dans un environnement intérieur clos, alors les angles solides sont 30 les sphères entières. 12 3031194 Si un objet est ajouté à l'environnement, et s'il est vu des réseaux d'antennes TX et RX sous des angles solides relativement petits, alors les champs électromagnétiques au voisinage des antennes RX sont en première approximation la somme des champs en l'absence de l'objet 5 rajouté et des champs dus aux courants créés dans l'objet rajouté. Ce dernier peut être vu comme une antenne qui émet dans l'environnement initial un champ électromagnétique qui est la somme d'un champ direct dans l'espace libre et de champs images en général sensiblement plus faibles résultant des diverses réflexions, diffractions dans 10 l'environnement. Sur cette base, la méthode proposée par la présente invention permet de déterminer la présence en intérieur de personnes ou d'objets dans un espace clos ou semi-clos en mesurant à des intervalles 15 successifs les valeurs électromagnétiques au voisinage des antennes, puis en les comparant afin de déterminer une variation, qui est alors représentative d'un changement dans l'environnement, cette variation peut être l'absence, la présence, l'arrivée, le départ d'une ou plusieurs personnes ou objets de l'environnement intérieur. 20 La méthode permet ensuite par le calcul des matrices de transfert de l'environnement, sans et avec la/les personne(s) ou objet(s) ajouté(s), de déterminer une/des position(s). En effet, la fonction d'étalement Ai+2(QR,u -T, ) correspondant à un environnement avec une personnes ou 25 un objet ajouté, et la fonction d'étalement Al(CIR,QT) correspondant à cet environnement initial seul, diffèrent approximativement par une fonction dont le support est dans les angles solides sous lesquels la personne ou l'objet ajouté est « vu » depuis les réseaux d'antennes TX et RX. 13 303 1 1 9 4 Dans le cas où les antennes sont identiques et forment des réseaux linéaires, les éléments de la matrice de transfert 'H' peuvent s'écrire sous la forme : -2 Tt i Hmn = fs2xS2 dORdOT e (naRnR+narnT) f(QR, OT) 5 ou encore en procédant aux intégrations par rapport aux variables cp-r et (PR et en introduisant les variables UT = COS(e-r) et UR = COS(eR) dans le cas où la distance entre antennes est À/2 : ff 1 Hmn = duRduT e'i(nuR+ nuT ) KuR, uT La mesure des matrices de transfert Hml+n2 et Hml n correspondant 10 aux situations avec ou sans personne ou objet ajouté, permet de déterminer les fonctions correspondantes h1+2(uR, u-r) et hl (uR, u-r) puisque les Hmn sont les coefficients de Fourier de h(uR, UT). La différence « h1+2(uR, LIT) - hi (uR ) >> permet de localiser la 15 position de l'objet par une recherche du maximum de cette fonction de localisation. Avantageusement, plus les réseaux comportent d'antennes, plus la détermination des fonctions à partir des coefficients de Fourier est précise. De manière préférentielle, les réseaux peuvent être formés de 20 NT = NR = 5 ou NT = NR = 7 antennes. La figure 3 illustre un enchainement des étapes du procédé (300) de localisation de l'invention selon un mode de réalisation. Dans une première étape (302), le procédé permet de générer la 25 matrice de transfert initiale H°. Pour cela, chaque antenne d'émission TX émet un signal sinusoïdal de fréquence fixe de l'ordre du GHz, choisie dans des plages règlementairement autorisées (2,4-2,5 GHz , 5,725-5,875GHz, 24-24,25 GHz,...). 14 303 1 1 9 4 Après un délai qui permet la stabilisation du signal au niveau des antennes de réception RX, l'amplitude et la phase du signal reçu sont mesurées. Le délai est dans une implémentation préférentielle, inférieur à la micro-seconde. Il est à considérer que la phase du signal 5 devant être mesurée avec précision, le procédé doit déterminer et prendre en compte les déphasages introduits par les différents composants électroniques du réseau de réception : switches, amplificateur, détecteur de phase. Le procédé permet d'obtenir un ensemble de « NRxNT » 10 nombres complexes qui constituent la matrice de transfert initiale du système MIMO. De manière avantageuse, le processus total d'acquisition doit être assez rapide afin que la personne à localiser puisse être considérée comme statique pendant la durée de l'acquisition. 15 Plusieurs modes d'acquisition sont possibles. Dans un mode préférentiel, la matrice de transfert initiale est d'abord enregistrée en l'absence de toute personne au sein de l'espace intérieur. Puis des matrices de transfert sont successivement enregistrées (304) à intervalles réguliers - kT - qui vont permettre de détecter la présence 20 d'une ou plusieurs personnes qui seront ensuite localisées. La période d'échantillonnage T peut-être de l'ordre de la seconde ou de la dizaine de secondes si l'on veut suivre le déplacement des personnes. La matrice initiale est désignée par H°, les matrices successives par Hk, le procédé permet dans une étape suivante (306) 25 d'établir la différence Hk - H° = ,8,Hk. A l'aide des éléments OH ,n, une fonction de localisation 'ha' est définie pour des variables uT et uR prises dans un intervalle [-1,1] par : 15 3 03 1 194 N+1 N+1 2 2 11A (uR, uT; := ejirmuR e jirnuT en N-1 N-1 m=- n= 2 2 Puis le procédé permet de déterminer (308) le maximum de la fonction . Les valeurs de uR, uT pour lesquelles le maximum est atteint sont désignées par : 5 uiV (k),ur. (k) Avec les notations suivantes : - le centre du réseau d'antennes TX de coordonnées (xi-,.117-) ; - le centre du réseau d'antennes RX de coordonnées (xR, yR) ; 10 - la direction du réseau d'antennes TX est donnée par l'angle ar; - la direction du réseau d'antennes RX est donnée par l'angle aR; et en notant : PI ° arccos iVl PR ---= arecos(uIti k le procédé permet de calculer (310) les coordonnées (x,y) de la 15 personne et déterminer sa position au temps kT selon les équations suivantes : 16 303 1 1 9 4 Dans la situation où il y a plusieurs personnes suffisamment 5 séparées, la fonction de localisation présentera plusieurs maxima et le même algorithme sera appliqué à chaque maximum de la fonction pour déterminer la position de chaque personne. Référence est maintenant faite à la figure 4 qui illustre les différents 10 blocs fonctionnels des circuits d'émission et de réception dans un mode de réalisation de l'invention. Le circuit d'émission comprend un oscillateur (402) qui génère une tension sinusoïdale à la fréquence fixe requise. Il est couplé à un switch 404, dont le rôle est de diriger successivement la tension à la sortie de l'oscillateur vers les différentes antennes de 15 transmission du réseau TX (406). Le circuit de réception comprend un réseau (408) d'antennes de réception RX, couplé à un switch (410). Un amplificateur (412) permet d'amplifier le signal reçu, puis un bloc de détection de phase et d'amplitude (414) permet de mesurer l'amplitude et la phase du signal 20 reçu. Une unité de traitement (416) permet d'enregistrer les mesures faites, d'effectuer les calculs des matrices de transfert telles que décrites précédemment, et de calculer les coordonnées de position d'une ou plusieurs personnes situées dans une ou plusieurs pièces, selon les équations décrites précédemment. 25 17 3031194 Dans l'exemple illustré, l'oscillateur (la source du signal RF) est placé du côté du circuit d'émission. Alternativement, il peut être placé du côté du circuit de réception RX auquel cas un câble coaxial est ajouté pour permettre d'envoyer le signal sur le circuit interrupteur (404) du 5 réseau d'émission. Avantageusement, pour chaque appareil, les déphasages sont calibrés, en particulier ceux introduits par les interrupteurs. Ils sont ensuite mémorisés dans un bloc mémoire de l'unité de calcul pour être pris en compte lors de la détermination des coordonnées d'une personne. 10 Ainsi, il a été décrit une méthode permettant de déterminer la présence en intérieur de personnes ou d'objets dans un espace clos ou semi-clos, en effectuant une détection différentielle d'une masse additionnelle dans un champ électromagnétique au moyen d'antennes multiples et de mesures de matrices de transfert. La localisation et les 15 coordonnées de la position de la ou des personne(s) ou objet(s) sont calculées à partir des mesures différentielles des matrices de transfert, selon les équations décrites. Les coordonnées de la personne localisée peuvent ensuite être utilisées selon la nature de l'application, soit en étant affichées sur un 20 écran local ou à distance, soit pour déclencher des actions comme lancer une alerte par exemple. L'homme du métier peut par des moyens courants coupler le résultat obtenu par le procédé décrit à des opérations postérieures selon les usages requis. 25 Ainsi, la présente description illustre une implémentation préférentielle de l'invention, mais n'est pas limitative. Un exemple a été choisi pour permettre une bonne compréhension des principes de l'invention, et une application concrète, mais il n'est en rien exhaustif et doit permettre à l'homme du métier d'apporter des modifications et des 30 variantes d'implémentation en conservant les mêmes principes. 18 REVENDICATIONS1. Un procédé pour déterminer l'ajout de personnes ou d'objets dans un espace intérieur, l'espace étant clos ou semi-clos et équipé d'un premier réseau d'antennes en émission distant d'un second réseau d'antennes en réception, le procédé comprenant les étapes suivantes : - mesurer une valeur initiale de champ électromagnétique créé par les réseaux d'antennes pour l'espace intérieur ; - mesurer à intervalles successifs des valeurs du champ électromagnétique au voisinage du réseau d'antennes en réception ; et - comparer les valeurs mesurées à la valeur initiale pour déterminer une différence, ladite différence étant représentative d'un ajout dans l'espace intérieur. 2. Le procédé selon la revendication 1, dans lequel : - les étapes de mesures du champ électromagnétique consistent à déterminer et enregistrer une matrice de transfert MIMO initiale constituée par les deux réseaux d'antennes pour l'espace intérieur vide, puis à mesurer des matrices de transfert MIMO successives ; et - l'étape de comparaison consiste à déterminer s'il existe des différences entre les éléments de la matrice de transfert initiale et une matrice ultérieure. 3. Le procédé selon la revendication 1 ou 2 dans lequel à l'intérieur de chaque réseau, les antennes sont distantes d'une longueur d'onde divisée par deux. 19 3031194 4. Le procédé selon l'une quelconque des revendications 1 à 3 comprenant de plus, s'il existe une différence entre les éléments de la matrice de transfert initiale et une matrice ultérieure, une étape de détermination de la position desdites 5 personnes ou objets ajoutés. 5. Le procédé selon la revendication 4 dans lequel l'étape de détermination de la position desdites personnes ou objets consiste à définir une fonction de localisation dont les 10 coefficients de Fourier sont les éléments des matrices de transfert, et à calculer le ou les maxima de ladite fonction. 6. Le procédé selon la revendication 5 comprenant de plus une étape de calcul des coordonnés spatiales desdites personnes 15 ou objets à un temps `IcT' à partir du maximum de la fonction de localisation. 7. Le procédé selon l'une quelconque des revendications 1 à 6 dans lequel les intervalles successifs de mesure sont dans une 20 plage de 1 à 10 secondes. 8. Un système pour déterminer l'ajout de personnes ou d'objets dans un espace intérieur, le système comprenant des moyens pour mettre en oeuvre les étapes du procédé selon l'une 25 quelconque des revendications 1 à 7. 9. Le système selon la revendication 8, dans lequel les moyens comprennent un circuit d'émission distant d'un circuit de réception, le circuit d'émission comprenant au moins un réseau 30 linéaire d'antennes en émission, le circuit de réception 20 3031194 comprenant au moins un réseau linéaire d'antennes en réception, les antennes à l'intérieur de chaque réseau étant distantes d'une longueur d'onde divisée par deux. 5 10.Le système selon la revendication 9 dans lequel le circuit de réception comprend de plus une unité de traitement pour enregistrer des mesures et effectuer des calculs de matrices de transfert et de coordonnées de position. 10 11. Le système selon la revendication 10 dans lequel l'unité de traitement comprend un processeur apte à effectuer lesdits calculs.