Brevet FR3015049A1 — publié le 20150619

PROCEDE DE LOCALISATION D'OBJETS MAGNETIQUES MOBILES PRESENTES DEVANT UN RESEAU DE MAGNETOMETRES [001] L'invention concerne un procédé et un dispositif de localisation d'objets magnétiques mobiles présentés devant un réseau de magnétomètres. L'invention a également pour objet un support d'enregistrement d'informations pour la mise en oeuvre de ce procédé. [002] Les demandes de brevet FR2 988 862 et FR2 988 874 divulguent des procédés de localisation d'objets magnétiques mobiles présentés devant un réseau de magnétomètres comportant N magnétomètres tri-axes reliés mécaniquement les uns aux autres sans aucun degré de liberté pour conserver une distance connue entre chacun de ces magnétomètres, où N est un nombre entier supérieur ou égal à cinq. Ces procédés connus comportent : a) la mesure par chacun des magnétomètres de l'amplitude du champ magnétique le 15 long de chacun de ses axes de mesure, b) l'estimation à partir des mesures des magnétomètres, des valeurs de plusieurs variables en résolvant un système actuel d'équations reliant ces variables à chaque mesure d'un magnétomètre tri-axe du réseau, chaque variable correspondant à la position ou à l'orientation ou à l'amplitude d'un moment magnétique d'un aimant 20 permanent de l'objet magnétique, - l'itération des étapes a) et b) à des instants successifs dans le temps pour obtenir les valeurs estimées de chaque variable à ces différents instants successifs. [003] Ces procédés fonctionnent particulièrement bien. En particulier, ces procédés sont capables de localiser de nombreux objets magnétiques mobiles différents. De 25 plus, ils sont aussi capables de localiser simultanément plusieurs objets magnétiques mobiles mécaniquement indépendants les uns des autres. A cause de ces capacités, on dit que ces procédés sont « flexibles ». [004] Toutefois, lorsque le nombre d'objets magnétiques mobiles à localiser simultanément augmente, la précision de la localisation de chacun de ces objets 30 magnétiques mobiles diminue. [005] L'invention vise à améliorer la précision de ces procédés connus de localisation d'objets magnétiques mobiles, notamment lorsque le nombre d'objets magnétiques mobiles à localiser simultanément augmente, tout en conservant leur flexibilité. 35 [006] Elle a donc pour objet un tel procédé comportant : c) l'identification, en fonction des objets mobiles présentés devant le réseau de magnétomètres, d'une équation supplémentaire reliant, par une relation d'égalité, une première variable du système actuel d'équations à un terme, ce terme étant une relation prédéfinie entre une ou plusieurs autres variables du système actuel d'équations ou une valeur numérique constante, d) le remplacement, dans le système actuel d'équations, de cette première variable par le terme auquel elle est égale, identifié lors de l'étape c), pour obtenir un nouveau 5 système d'équations dans lequel le nombre de variables dont la valeur doit être estimée est plus petit que dans le système actuel d'équations, e) l'utilisation du nouveau système d'équations à la place du système actuel d'équations lors des itérations suivantes des étapes a) et b), le nouveau système d'équations devenant ainsi, pour les itérations suivantes des étapes a) et b), le 10 système actuel d'équations. [7] Le déposant a constaté que réduire le nombre de variables dans le système d'équations utilisé lors de l'étape b) permet de réduire le bruit d'estimation sur les positions et orientations estimées des objets magnétiques mobiles. Cela se traduit donc par une amélioration de la précision de localisation des objets magnétiques 15 mobiles. Toutefois, cela réduit également la flexibilité du procédé. [8] Pour réduire ce nombre de variables sans nuire à la flexibilité du procédé, le procédé ci-dessus commence par utiliser un système d'équations ayant de nombreuses variables indépendantes. Ensuite, en fonction des objets mobiles présentés devant le réseau de magnétomètres, il identifie des équations 20 supplémentaires reliant entre elles certaines variables du système actuel d'équations ou reliant une variable du système actuel d'équations à une valeur constante. Ces équations supplémentaires sont alors utilisées pour réduire le nombre de variables du système actuel d'équations. Ainsi, ce procédé adapte automatiquement et dynamiquement le nombre de variables du système d'équations en fonction des 25 objets magnétiques mobiles à localiser. Cela permet de conserver la flexibilité du procédé tout en améliorant la précision de localisation. [9] De plus, lorsque le nombre de variables pour lesquelles une valeur doit être estimée lors de l'étape b) est réduit, le procédé de localisation est plus rapide. [0010] Le procédé ci-dessus ne nécessite pas de stocker un grand nombre de 30 systèmes d'équations prédéfinis mais simplement un système d'équations initial et des équations supplémentaires. [0011] Enfin, diminuer le nombre de variables à estimer permet aussi d'accroître la robustesse du procédé. [0012] Les modes de réalisation de ce procédé peuvent comporter une ou plusieurs 35 des caractéristiques suivantes : ^ l'identification de l'équation supplémentaire est réalisée en fonction des valeurs estimées, lors d'une précédente itération des étapes a) et b), des variables ; ^ l'identification de l'équation supplémentaire comporte : - le calcul d'une grandeur représentative de l'amplitude des variations de la première variable au cours de précédentes itérations des étapes a) et b) à partir des valeurs de cette première variable estimée au cours de ces précédentes itérations des étapes a) et b), - la comparaison de cette grandeur calculée à un seuil prédéterminé pour établir si cette première variable a varié au cours de ces précédentes itérations des étapes a) et b), et - si ce seuil prédéterminé est franchi, l'établissement d'une nouvelle équation supplémentaire reliant, par une relation d'égalité, la valeur de cette première variable à une valeur constante, la valeur constante étant fonction d'au moins une valeur de cette première variable estimée lors de ces précédentes itérations des étapes a) et b) au cours desquelles il a été établi que sa valeur ne varie pas, puis l'utilisation de cette nouvelle équation supplémentaire lors de l'étape d), et - si ce seuil prédéterminé n'est pas franchi, l'absence d'établissement de cette nouvelle équation supplémentaire ; ^ l'étape c) comporte : - l'acquisition d'une caractéristique distinctive de l'objet magnétique mobile actuellement présenté devant le réseau de magnétomètres, et - la comparaison de la caractéristique distinctive acquise à des conditions de sélection préenregistrées dans une base de données, cette base de données associant à chaque condition de sélection au moins une équation supplémentaire reliant, par une relation d'égalité, une variable du système actuel d'équations à un terme, ce terme étant une relation prédéfinie entre une ou plusieurs autres variables du système actuel d'équations ou une valeur numérique constante, et - uniquement si la caractéristique distinctive acquise correspond à l'une des conditions de sélection préenregistrées, alors l'utilisation lors de l'étape d) de l'équation supplémentaire associée à cette condition de sélection ; ^ l'acquisition d'une caractéristique distinctive comporte l'obtention de cette caractéristique distinctive à partir des valeurs des variables estimées à l'aide du système actuel d'équations ; ^ le procédé comporte : - le calcul d'une erreur d'estimation avec le système actuel d'équations, l'erreur d'estimation étant représentative de l'écart entre : des valeurs estimées des mesures des magnétomètres lorsque la position, l'orientation et l'amplitude du moment magnétique de chaque objet magnétique sont égales à celles estimées en utilisant le système actuel d'équations lors d'une itération des étapes a) et b), et les mesures des magnétomètres relevées à l'étape a) lors de la même itération des étapes a) et b), et - si l'erreur d'estimation calculée franchit un seuil prédéterminé : la sélection d'au moins une des équations supplémentaires précédemment utilisées lors de l'étape d) pour obtenir le système actuel d'équations, le remplacement, dans le système actuel d'équations, du terme de l'équation supplémentaire sélectionnée, par la variable que ce terme a remplacé pour obtenir un nouveau système d'équations dans lequel le nombre de variables dont la valeur doit être estimée est plus grand que dans le système actuel d'équations et l'utilisation du nouveau système d'équations à la place du système actuel d'équations lors des itérations suivantes des étapes a) et b), et - si l'erreur calculée ne franchit pas ce seuil prédéterminé, le maintien de l'utilisation du système actuel d'équations lors des itérations suivantes des étapes a) et b) ; ^ la sélection d'au moins une équation supplémentaire comporte : - la sélection d'au moins deux équations supplémentaires précédemment utilisées lors de l'étape d) pour obtenir le système actuel d'équations, - pour chacune des équations supplémentaires sélectionnées, le remplacement, dans le système actuel d'équations, du terme de l'équation supplémentaire sélectionnée, par la variable que ce terme a remplacée pour obtenir des troisième est quatrième systèmes d'équations différents, chacun de ces quatrième et troisième systèmes d'équations comportant un nombre de variables dont la valeur doit être estimée plus grand que dans le système actuel d'équations, - le calcul de l'erreur d'estimation avec ces troisième et quatrième systèmes d'équations, et - la sélection de l'équation supplémentaire, parmi celles initialement sélectionnées, qui minimise l'erreur d'estimation calculée. 30 [0013] Ces modes de réalisation du procédé de localisation présentent en outre les avantages suivants : - utiliser les valeurs précédemment estimées pour identifier l'équation supplémentaire permet d'utiliser les magnétomètres à la fois pour localiser les objets magnétiques et pour sélectionner l'équation supplémentaire, 35 - en cas d'absence de variation d'une des variables, établir l'équation supplémentaire selon laquelle cette variable est égale à une valeur constante permet de diminuer le nombre de variables à estimer sans pour cela utiliser des connaissances préenregistrées sur le ou les comportements possibles de cette variable, - utiliser un seuil pour déterminer si une variable varie, ce seuil étant déterminé à partir de la variation maximale observée pour cette variable, permet d'augmenter la précision du procédé, - l'identification d'une équation supplémentaire en comparant les caractéristiques distinctives d'un objet mobile à des conditions de sélection préenregistrées permet de limiter le nombre de variables à estimer en prenant en compte une base de connaissances prédéfinie sur les déplacements possibles des objets magnétiques mobiles, - obtenir la caractéristique distinctive à partir des mesures réalisées par le réseau de 10 magnétomètres permet d'identifier l'équation supplémentaire en utilisant les mêmes magnétomètres que ceux utilisés pour localiser l'objet magnétique mobile ; - ajouter à nouveau une variable au système actuel d'équations lorsque l'erreur d'estimation dépasse un certain seuil permet d'adapter automatiquement le procédé à un changement d'objet magnétique mobile ou à un changement dans son utilisation ; 15 - sélectionner la variable à réintroduire dans le système actuel d'équations à partir de l'erreur d'estimation calculée dans le cas où cette variable est absente du système d'équations et dans le cas où cette variable est présente dans le système d'équations permet de systématiquement diminuer l'erreur à chaque fois qu'une variable est réintroduite dans le système d'équations. 20 [0014] L'invention a également pour objet un support d'enregistrement d'informations comportant des instructions pour l'exécution du procédé ci-dessus, lorsque ces instructions sont exécutées par un calculateur électronique. [0015] Enfin, l'invention a également pour objet un dispositif de localisation d'objets magnétiques mobiles, ce dispositif comportant : 25 - un réseau de magnétomètres comportant N magnétomètres tri-axes reliés mécaniquement les uns aux autres sans aucun degré de liberté pour conserver une distance connue entre chacun de ces magnétomètres, où N est un nombre entier supérieur ou égal à cinq, - une unité de traitement apte: 30 a) à acquérir les mesures de chacun des magnétomètres du réseau de magnétomètres, b) à estimer, à partir des mesures des magnétomètres, des valeurs de plusieurs variables en résolvant un système actuel d'équations reliant ces variables à chaque mesure d'un magnétomètre tri-axe du réseau, chaque variable 35 correspondant à la position ou à l'orientation ou à l'amplitude d'un moment magnétique d'un objet magnétique, - à réitérer les étapes a) et b) à des instants successifs dans le temps pour obtenir les valeurs estimées de chaque variable à ces différentes instants successifs, dans lequel l'unité de traitement est également apte : c) à identifier, en fonction des objets mobiles présentés devant le réseau de magnétomètres, une équation supplémentaire reliant, par une relation d'égalité, une première variable du système actuel d'équations à un terme, ce terme étant une relation prédéfinie entre une ou plusieurs autres variables du système actuel d'équations ou une valeur numérique constante, d) à remplacer, dans le système actuel d'équations, cette première variable par le terme auquel elle est égale, identifié lors de l'étape c), pour obtenir un nouveau système d'équations dans lequel le nombre de variables dont la valeur doit être estimée est plus petit que dans le système actuel d'équations, e) à utiliser le nouveau système d'équations à la place du système actuel d'équations 15 lors des itérations suivantes des étapes a) et b), le nouveau système d'équations devenant ainsi, pour les itérations suivantes des étapes a) et b), le système actuel d'équations. [0016] L'invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple non limitatif et faite en se référant aux dessins, 20 sur lesquels : - la figure 1 est une illustration schématique d'une interface homme-machine permettant de commander un appareil électrique. - la figure 2 est une illustration partielle, en vue de dessus, d'un réseau de magnétomètres mis en oeuvre dans l'interface homme-machine de la figure 1 ; 25 - la figure 3 est une illustration schématique, en coupe verticale, d'un curseur utilisable dans l'interface de la figure 1 ; - la figure 4 est une illustration schématique, en coupe verticale, d'un bouton utilisable dans l'interface de la figure 1 ; - la figure 5 est une illustration schématique d'une base de données utilisée dans 30 l'interface de la figure 1 ; - la figure 6 est un organigramme d'un procédé de commande d'un appareil électrique à l'aide de l'interface homme-machine de la figure 1. [0017] Dans ces figures, les mêmes références sont utilisées pour désigner les mêmes éléments. 35 [0018] Dans la suite de cette description, les caractéristiques et fonctions bien connues de l'homme du métier ne sont pas décrites en détails. [0019] La figure 1 représente une interface homme-machine 2 permettant de commander un appareil électrique 4. Ici, l'appareil électrique comporte un écran et une unité 5 de commande capable de commander l'affichage d'une image sur cet écran. [0020] Ici, le fonctionnement de l'interface 2 est illustré dans le cas où l'unité 5 est une console de jeux vidéo. Par exemple, l'unité 5 commande le déplacement d'un 5 personnage 6 sur l'écran. Toutefois, l'interface 2 est utilisable dans de nombreuses autres applications comme cela est décrit à la fin de cette description. [0021] L'interface 2 comporte plusieurs ustensiles actionnables directement à la main par un être humain, dénommé « utilisateur » par la suite. Chacun de ces ustensiles comporte au moins un objet magnétique. Ici chaque objet magnétique est 10 un aimant permanent. Pour simplifier la figure 1, seul un ustensile 10 est représenté sur cette figure. D'autres ustensiles de l'interface 2 sont décrits en référence aux figures 3 et 4. L'interface 2 comporte également un dispositif 12 de localisation de chaque aimant permanent. [0022] Dans ce mode de réalisation, pour modifier la configuration de l'interface 2, 15 chaque ustensile est librement déplaçable, directement par la main de l'utilisateur, dans un repère orthogonal Rm fixé sans aucun degré de liberté au dispositif 12. Ici, les directions X et Y du repère Rm sont horizontales et la direction Z est verticale. A cet effet, chaque ustensile pèse moins d'un kilo et, de préférence, moins de 200g. Les dimensions de chaque ustensile sont suffisamment réduites pour qu'il puisse être 20 saisi et déplacé par une seule main de l'utilisateur. [0023] Dans ce mode de réalisation, l'ustensile 10 comprend une manette 14, un aimant permanent 16 et un support 18. La manette 14 a une forme oblongue pour former une poignée facilement saisissable par l'utilisateur. Elle s'étend le long d'un axe longitudinal 17. La manette 14 est destinée à être utilisée comme un joystick. A 25 cet effet, une extrémité inférieure de la manette 14 est fixée au support 18 par l'intermédiaire d'une rotule 19. Le support 18 est pourvu d'une face plane inférieure destinée à être fixée ou à reposer sur un plan de manière à immobiliser, dans le repère Rm, la position de la rotule 19 lorsque l'ustensile 10 est utilisé. Ainsi, la manette 14 peut être inclinée librement par l'utilisateur autour du centre de rotation de 30 la rotule 19. Lors de son utilisation, la manette 14 présente donc trois degrés de liberté en rotation. [0024] Par exemple, la manette 14, le support 18 et la rotule 19 sont entièrement réalisés dans un matériau non-magnétique, c'est-à-dire un matériau ne présentant pas de propriété magnétique mesurable par le dispositif 12. Ce matériau est par 35 exemple du plastique. La position de la manette 14 autour des axes X, Y et Z est localisée à partir de la position de l'aimant 16. L'aimant 16 présente un moment magnétique non nul même en absence de champ magnétique extérieur [0025] Typiquement, dans cette description, le champ magnétique coercitif de chaque aimant est supérieur à 100 A.m-1 ou 500 A.m-1. Par exemple, il est réalisé en matériau ferro- ou ferrimagnétique. La puissance de chaque aimant permanent est typiquement supérieure à 0,01 A.m2 ou 0,1 A.m2. [0026] L'aimant 16 est fixé sans aucun degré de liberté sur la manette 14. La direction du moment magnétique de l'aimant 16 est ici confondue avec l'axe longitudinal 17 de la manette 14. Sur la figure 1 et les suivantes, la direction du moment magnétique d'un aimant est représentée par une flèche. La plus grande longueur de cet aimant est notée L par la suite. [0027] Lors de l'utilisation de la manette 14, les coordonnées de l'aimant 16 sont reliées les unes aux autres par la relation vectorielle suivante : P + Rd*(M/ M II) = Po, 10 où - Po est le vecteur codant la position du centre de rotation de la rotule 19 de l'ustensile 10 dans le repère Rm, - P est le vecteur codant la position du centre géométrique de l'aimant 16 dans le repère Rm, 15 - M est le moment magnétique de l'aimant 16, - II M II est l'amplitude du moment magnétique M, - Rd est la distance qui sépare le centre de rotation de la rotule 19 du centre géométrique de l'aimant 16. [0028] A partir de cette relation vectorielle, il est possible d'établir trois équations 20 supplémentaires reliant entre elles les variables codant la position, l'orientation et l'amplitude du moment magnétique de l'aimant 16. Par exemple, ces équations supplémentaires sont les suivantes : ,\/m, x2+m- y2+m- z2 Rd *Mx = Px- Pox (1) 25 -Rd *My = Py- Poy (2) ~IMx2+M y2+m z2 Rd *Mz = Pz- Poz (3) ,\Liv y2+m- z2 où : 30 - Mx, My et Mz sont les projections orthogonales du moment magnétique M sur les axes X, Y et Z du repère Rm. [0029] Ces équations expriment le fait : - que le centre géométrique de l'aimant 16 est toujours situé sur un sphère centrée sur le centre de rotation de la rotule 19 et de rayon Rd, et - que le moment magnétique de l'aimant 16 pointe en permanence sur le centre de rotation de la rotule 19. [0030] Le centre de géométrie d'un objet est le barycentre de l'ensemble des points de cet objet en affectant le même poids à chacun de ces points. [0031] Pour simplifier la description, dans la suite de cette description, les équations supplémentaires sont exprimées dans un repère local lié sans aucun degré de liberté à l'ustensile. Il est à la portée de l'homme du métier d'exprimer ces mêmes équations supplémentaires dans le repère Rm en tenant compte des coordonnées et de l'orientation du repère local par rapport au repère Rm. Il s'agit simplement d'un changement de repère. [0032] Le dispositif 12 permet de localiser un ou plusieurs aimants dans le repère Rm. Par localisation d'un aimant, on entend ici la détermination de la position et de l'orientation de l'aimant 16 dans le repère Rm. La position est définie sans ambiguïté par les valeurs de trois variables, par exemple les coordonnées x, y et z dans le repère Rm. Plus précisément, dans le cas d'un aimant, les variables x, y et z sont les coordonnées du centre géométrique de l'aimant. L'orientation du moment magnétique d'un aimant est définie dans le repère Rm par les valeurs de deux variables ey et 0,. Ici, les variables ey et 0, sont les angles du moment magnétique de l'aimant, respectivement, par rapport aux axes Y et Z du repère Rm. Le dispositif 12 détermine également un sixième variable A. La variable A est l'amplitude du moment magnétique de l'aimant. [0033] Le dispositif 12 comporte à cet effet un réseau de N magnétomètres tri-axes Mi. Sur la figure 1, les traits ondulés verticaux indiquent qu'une partie du dispositif 12 n'a pas été représentée. [0034] Typiquement, N est supérieur à cinq et, de préférence, supérieur à seize ou trente-deux. Ici, N est supérieur ou égal à soixante-quatre. [0035] Dans ce mode de réalisation, les magnétomètres Mi; sont alignés en lignes et en colonnes pour former une matrice. Ici, cette matrice comporte huit lignes et huit colonnes. Les indices i et j identifient, respectivement, la ligne et la colonne de cette 30 matrice à l'intersection de laquelle se trouve le magnétomètre ni. Sur la figure 1, seuls les magnétomètres Mil, -M2, -Mi3, -NI Mi8 d'une ligne i sont visibles. La position - des magnétomètres Mi; les uns par rapport aux autres est décrite plus en détail en référence à la figure 2. [0036] Chaque magnétomètre Mi; est fixé sans aucun degré de liberté aux autres 35 magnétomètres. A cet effet, les magnétomètres Mi; sont fixés sans aucun degré de liberté sur une face arrière 22 d'une plaque rigide 20. Cette plaque rigide présente une face avant 24 tournée vers l'aimant 16. La plaque 20 est réalisée dans un matériau non magnétique rigide. Par exemple, la plaque 20 est en verre. [0037] Chaque magnétomètre Mu mesure la direction et l'intensité du champ magnétique généré par les aimants présents devant la face 24. Pour cela, chaque magnétomètre Mu mesure la norme de la projection orthogonale du champ magnétique au niveau de ce magnétomètre Mu sur trois axes de mesure de ce magnétomètre. Ici, ces trois axes de mesure sont orthogonaux entre eux. Par exemple, les axes de mesure de chacun des magnétomètres Mu sont, respectivement, parallèles aux axes X, Y et Z du repère. La sensibilité du magnétomètre Mu est par exemple de 4*10-7T. [0038] Chaque magnétomètre Mu est raccordé par l'intermédiaire d'un bus 28 de 10 transmission d'informations à une unité de traitement 30. [0039] L'unité de traitement 30 est capable, pour chaque aimant, de le localiser dans le repère Rm et d'estimer l'amplitude de son moment magnétique à partir des mesures des magnétomètres Mu. A cet effet, l'unité 30 comporte un calculateur électronique 32 programmable apte à exécuter des instructions enregistrées sur un 15 support d'enregistrement d'informations. L'unité 30 comporte donc aussi une mémoire 34 contenant les instructions nécessaires pour l'exécution par le calculateur 32 du procédé de la figure 6. Dans ce mode de réalisation, pour chaque nombre P d'objets magnétiques susceptibles d'être simultanément utilisés dans l'interface 2, l'unité 30 implémente un modèle mathématique Mp associant chaque mesure d'un 20 magnétomètre Mu aux positions, orientations et amplitudes des moments magnétiques de P objets magnétiques dans le repère Rm. Chaque modèle M p se présente sous la forme d'un système d'équations dans lequel un premier jeu de variables représente les positions et orientations des P objets magnétiques ainsi que les amplitudes des moments magnétiques de ces objets. Un second jeu de variables 25 représente les mesures des magnétomètres Mu. Pour obtenir les positions, orientations et amplitudes des moments magnétiques des P objets magnétiques, les variables du premier jeu sont les inconnues et les valeurs des variables du second jeu sont connues. Ce modèle est typiquement construit à partir des équations physiques du magnétisme. Ce modèle est paramétré par les distances connues 30 entres les magnétomètres Mu. Ici, les objets magnétiques sont les aimants permanents. Pour construire ce modèle, chaque aimant permanent est approximé par un dipôle magnétique. Cette approximation n'introduit que très peu d'erreurs si la distance entre l'aimant permanent et le magnétomètre Mu est supérieure à 2L et, de préférence, supérieure à 3L, où L est la plus grande dimension de l'aimant 35 permanent. Typiquement, L est inférieur à 20 cm et, de préférence, inférieur à 10 ou 5 cm. [0040] Ici, le modèle Mp est non linéaire. L'unité 30 le résout en mettant en oeuvre un algorithme d'estimation de sa solution. Par exemple, l'algorithme utilisé est un filtre de Kalman d'ensemble plus connu sous le terme anglais de « Unscented Kalman Filter ». [0041] Étant donné que chaque objet magnétique est caractérisé par trois variables pour connaître sa position, deux variables pour connaître son orientation et une 5 variable pour connaître l'amplitude de son moment magnétique, le nombre maximum d'objets magnétiques simultanément localisables par le réseau de N magnétomètres est inférieur à N/2. Dès lors, la valeur du nombre P est inférieure ou égale à N/2 et, de préférence, inférieure à N/5 ou N/10 ou N/20 pour avoir des mesures redondantes. La redondance des mesures permet d'améliorer la précision de la localisation des 10 objets magnétiques. [0042] Dans ce mode de réalisation, seul cinq modèles mathématiques M1, M2, M3, M4 et M5 sont implémentés dans l'unité 30 pour, respectivement, 1, 2, 3, 4 et 5 aimants permanents simultanément présents devant la face 24. [0043] L'unité 30 est également capable de transmettre une commande à l'appareil 15 4 par l'intermédiaire d'une interface 36 reliée à cet appareil 4. [0044] La mémoire 34 comporte aussi une base 38 de données dans laquelle sont enregistrées plusieurs lois de commande de l'appareil 4. Chaque loi de commande permet de générer la commande de l'appareil 4 correspondant à l'état actuel de l'ustensile auquel elle est associée. A cet effet, chaque loi de commande associe : 20 - à au moins un état possible de l'ustensile, une commande de l'appareil 4, et - à un autre état possible du même ustensile, une autre commande de l'appareil 4 ou une absence de commande de l'appareil 4. [0045] La loi de commande détermine donc comment fonctionne l'appareil 4 en réponse à l'actionnement, par l'utilisateur, de l'ustensile qui est associé à cette loi de 25 commande. Cette base 38 est décrite plus en détail en référence à la Figure 5. [0046] La figure 2 représente une partie des magnétomètres Mu du dispositif 12. Ces magnétomètres Mu sont alignés en lignes i parallèles à la direction X. Ces magnétomètres sont également alignés en colonnes j parallèles à la direction Y pour former une matrice. Les lignes i et les colonnes j sont disposées dans l'ordre des 30 indices croissants. [0047] Le centre du magnétomètre Mu se trouve à l'intersection de la ligne i et de la colonne j. Le centre du magnétomètre correspond au point où le champ magnétique est mesuré par ce magnétomètre. Ici, les indices i et j appartiennent à l'intervalle [1 ; 8]. 35 [0048] Les centres de deux magnétomètres Mu et Mu., immédiatement consécutifs le long d'une ligne i sont séparés par une distance connue d,,,,,.1. De façon similaire, le centre de deux magnétomètres Mu et M,+1,,, immédiatement consécutifs le long d'une même colonne j sont séparés par une distance connue dJ,',.1. [0049] Dans le cas particulier décrit ici, quelle que soit la ligne i, la distance d,,,,J+1 est la même. Cette distance est donc notée d,. De façon similaire, quelle que soit la colonne j, la distance dJ,',.1 entre deux magnétomètres est la même. Cette distance est donc notée d,. [0050] Ici, les distances d, et dj sont toutes les deux égales à d. [0051] Typiquement, la distance d est inférieure, et de préférence deux fois inférieure, à la plus petite distance qui peut exister entre deux objets magnétiques simultanément présent devant la face 24 lors d'une utilisation normale de l'interface 2. Ici, la distance d est comprise entre 1 et 4 cm lorsque : - la puissance de l'aimant permanent est de 0,5 A.m2, - la sensibilité des magnétomètres est de 4*10-7T, et - le nombre de magnétomètres Mu est de soixante-quatre. [0052] Les figures 3 et 4 représentent d'autres ustensiles utilisables à la place de l'ustensile 10 ou simultanément avec l'ustensile 10 dans l'interface 2. [0053] La figure 3 représente un curseur 90. Ce curseur permet, par exemple, de régler la valeur d'un paramètre de l'appareil 4. Ce curseur 90 comporte à cet effet un coulisseau 92 monté déplaçable en translation dans une fente rectiligne 94 ménagée dans une face supérieure d'un boîtier 96. L'utilisateur peut faire coulisser à la main le coulisseau 92 pour modifier l'état du curseur 90. Dans le cas du curseur 90, chaque état correspond à une position particulière du coulisseau 92 le long de la fente 94. [0054] Pour mesurer la position du coulisseau 92 le long de la fente 94, le curseur 90 est équipé : - d'un aimant permanent 98 fixé sans aucun degré de liberté au coulisseau 92, et - d'un aimant permanent 100 fixé sans aucun degré de liberté au boîtier 96. 25 [0055] On définit un repère orthogonal local Rgo lié sans aucun degré de liberté au boîtier 96 est dont : - l'axe Xgo est parallèle à la fente 92, - l'axe Ygo est parallèle au fond du boîtier 96, - l'axe Zgo est perpendiculaire au fond du boîtier 96, et 30 - l'origine est confondue avec le centre géométrique de l'aimant 100. [0056] Pour des raisons de clarté, le repère Rgo est représenté à côté du curseur 90 sur la figure 3. [0057] Dans ce mode de réalisation, la direction du moment magnétique de l'aimant 98 est parallèle à la direction x9o. La direction du moment magnétique de l'aimant 100 35 est parallèle à la direction z90. [0058] Le nombre de degré de liberté du curseur 92 par rapport au boîtier 96 est limité à un degré de liberté en translation le long de l'axe Xgo. Dans ces conditions, les équations supplémentaires suivantes existent entre les coordonnées des aimants 98 et 100 dans le repère Rgo Z98 = Z100 + h98 et y98 = y100, où : - x98, y98, z98 sont les coordonnées dans le repère Rgo du centre géométrique de l'aimant 98, - xioo, yloo, zioo l'aimant 100, - h98, est une valeur numérique constante connue fonction des dimensions de l'ustensile 90. [0059] On a également les équations supplémentaires suivantes : e -y98 = ey100 et ez100 = ez98 + 90°, où : - ey98 et ez98 sont les angles entre la direction du moment magnétique de l'aimant 98 est les axes y90 et z90, et - ey100 et ez100 sont les angles entre la direction du moment magnétique de l'aimant 100 est les axes y90 et z90. [0060] Les amplitudes des moments magnétiques des aimants 98 et 100 sont uniques dans l'interface 2 et différentes l'une de l'autre. [0061] La figure 4 représente un bouton 110. Ce bouton 110 est déplaçable à la main par un utilisateur entre un état de repos (représenté sur la figure 4) et à un état enfoncé. Typiquement, le bouton 110 est destiné à déclencher une action de l'appareil 4 uniquement lorsqu'il atteint son état enfoncé. [0062] Ce bouton 110 comporte une touche 112 uniquement déplaçable en translation, le long d'un axe 114 solidaire d'un boîtier 116, entre une position de repos (représenté sur la figure 4) et une position enfoncée. Dans la figure 4, l'axe 114 est vertical. Les états de repos et enfoncé correspondent, respectivement, aux positions de repos et enfoncée de la touche 112. [0063] Le bouton 110 comporte un ressort 118 interposé entre une partie inférieure du bouton 112 et une face supérieure du boîtier 116. Une tige rectiligne 120 s'étend le long de l'axe 114, cette tige 120 est fixée, d'un côté, sans aucun degré de liberté, à la touche 112 et, de l'autre côté, à une butée 124. Cette tige 120 est montée à coulissement à l'intérieur d'un orifice 122 ménagé dans la face supérieure du boîtier 116. La butée 124 permet de retenir l'extrémité inférieure de cette tige 120 à l'intérieur du boîtier 116. [0064] Pour déterminer l'état dans lequel se trouve le bouton 110, la tige 120 comporte un aimant permanent 126 fixé sans aucun degré de liberté sur cette tige. Dans ce mode de réalisation, la direction du moment magnétique de l'aimant 126 est confondue avec l'axe 114. Un autre aimant permanent 128 est également fixé sans aucun degré de liberté sur le fond du boîtier 116. [0065] Dès lors, l'état dans lequel se trouve le bouton 110 peut être établi à partir de la position relative de l'aimant 126 par rapport à l'aimant 128. Par exemple, l'état est établi à partir de la valeur d'une distance d110 et d'un seuil So. La distance d110 est la distance la plus courte qui sépare les centres géométriques des aimants 126 et 128. sont les coordonnées dans le repère Rgo du centre géométrique de Le seuil So est utilisé pour discriminer l'état enfoncé de l'état de repos. Si la distance duo est inférieure au seuil So, cela signifie que le bouton 110 est dans son état enfoncé. A l'inverse, si la distance d110 est supérieure à ce seuil So, le bouton 110 est dans son état de repos. [0066] On définit un repère orthogonal local R110 fixé sans aucun degré de liberté au bouton 110. Pour des raisons de clarté, il est représenté à côté du bouton 110 sur la figure 4. L'origine du repère R110 est confondue avec le centre géométrique de l'aimant 128. Les axes X110 et Z110 de ce repère R110 sont confondus avec les directions des moments magnétiques, respectivement, des aimants 128 et 126. L'axe Z110 de ce repère R110 est donc confondu avec l'axe 114. [0067] L'aimant 126 ne présente qu'un seul degré de liberté en translation par rapport à l'aimant 128. Dans ces conditions, les équations supplémentaires suivantes relient entre elles les coordonnées des aimants 126 et 128 dans le repère R110 : x126 = X128, Y126 = Y128, eY126 = eY128, et ez126 = OZ128 + 90°, où : - x126, Y126, 2126 sont les coordonnées dans le repère R110 du centre géométrique de l'aimant 126, - X128, Y128, 2128 sont les coordonnées dans le repère R110 du centre géométrique de l'aimant 128. [0068] Ici, les aimants 126 et 128 sont identiques. [0069] La figure 5 représente la base de données 38. Cette base 38 comporte pour chaque ustensile U,: - une condition C, de sélection permettant de sélectionner, à partir d'une caractéristique distinctive de l'ustensile U' seulement cet ustensile U, et non pas les autres ustensiles ayant un agencement d'aimants différents et susceptibles d'être 25 utilisés dans l'interface 2, - une loi de commande L, préenregistrée, - un ensemble E, d'équations supplémentaires définissant la ou les relations prédéterminées qui existent entre les coordonnées et les orientations des moments magnétiques des aimants de cet ustensile. 30 [0070] Par la suite, on suppose que la caractéristique distinctive de chaque ustensile est l'amplitude du ou des moments magnétiques de ses aimants permanents. [0071] De préférence, la condition de sélection s'exprime sous la forme d'une plage de valeurs à laquelle doit appartenir la caractéristique distinctive de l'ustensile pour être vérifiée. 35 [0072] Si la condition C, est vérifiée par un aimant permanent ou une paire d'aimants permanents, alors la loi L, et l'ensemble E, sont associés à l'ustensile équipé de cet aimant ou paire d'aimants. L'indice « i » identifie la loi L, et l'ensemble E, associés à cette condition C, par la base 38. Ici, chaque condition C, porte sur les valeurs des variables estimées pour les aimants permanents qui équipent un ustensile et qui, généralement, ne sont pas utilisées pour établir l'état dans lequel se trouve actuellement cet ustensile. [0073] A titre d'illustration, un exemple de contenu de la base 38 va maintenant être décrit. [0074] La loi Li est une loi de commande qui utilise uniquement les valeurs des variables ey et Oz. Plus précisément, elle transforme chaque valeur des angles Oy et Oz en une commande respective de l'appareil 4. Par exemple, cette commande peut provoquer un déplacement du personnage 6 dans une direction spécifiée par les valeurs des angles ey et Oz. [0075] Cette loi Li est associée à une condition Ci. La condition Ci est la suivante : la valeur de la variable A doit être comprise entre Si et 52, où Si et 52 sont des limites prédéfinies telles que seule l'amplitude du moment magnétique de l'aimant 16 soit comprise entre ces limites. Cette loi Li est donc destinée à être associée à l'ustensile 10. Si le dispositif 12 détermine que la valeur de la variable A d'un aimant permanent est comprise entre Si et S2, alors il utilise la loi Li pour convertir les valeurs des angles ey et Oz de ce même aimant en une commande l'appareil 4. Dès lors, en modifiant l'inclinaison de l'ustensile 10, il est possible de commander l'appareil 4. Ceci reste vrai quelle que soit la position de l'ustensile 10 par rapport au dispositif 12 à partir du moment où la présence de l'aimant 16 peut être détectée. L'emplacement de l'ustensile 10 peut donc librement être choisi par l'utilisateur, ce qui correspond à une infinité de configuration possible pour l'interface 2. L'ensemble E1 comporte l'équation supplémentaire décrite en référence à la figure 1. [0076] La loi L2 est une loi de commande qui utilise la valeur d'un distance entre deux aimants permanents pour générer une commande de réglage d'un paramètre 25 de l'appareil 4 en fonction de cette distance. Cette loi L2 précise que la distance est mesurée entre : - un aimant permanent dont l'amplitude du moment magnétique est comprise entre des limites S71 et S81 prédéfinies telle que seul l'aimant 98 puisse satisfaire cette condition, et 30 - un aimant permanent dont l'amplitude du moment magnétique est comprise entre des limites S72 et S82 prédéfinies telle que seul l'aimant 100 puisse satisfaire cette condition. [0077] Cette loi de commande est typiquement conçue pour être associée au curseur 90 de sorte que la distance correspond à la distance d90. A cet effet, la loi L2 35 est associée à une condition C2. Cette condition C2 est la suivante : la valeur de la variable A d'un aimant permanent est comprise entre les limites S72 et S82. [0078] Avec une telle condition C2, la loi L2 est associée au curseur 90. L'ensemble E2 comporte les équations supplémentaires décrites en référence à la figure 3. [0079] La loi L3 est une loi de commande qui utilise uniquement la distance la plus courte entre deux aimants permanents identiques pour générer une commande de l'appareil 4. Ici, la loi L3 est destinée à être associée au bouton 110 de sorte que la distance la plus courte correspond à la distance d110. Par exemple, la loi de commande L3 est la suivante : - si la distance d110 est supérieure au seuil So, alors l'unité 30 ne génère aucune commande de l'appareil 4, et - si la distance d110 est inférieure ou égale au seuil So, alors l'unité 30 génère et transmet une commande à l'appareil 4. [0080] La loi L3 est associée à une condition C3. Ici, cette condition C3 est la suivante : - les valeurs des variables A de deux d'aimants permanents distincts sont comprises toutes les deux entre des limites S9 et Slo, et - la distance la plus courte entre ces deux aimants est inférieure à un seuil D max110. [0081] Les limites S9 et Slo sont des constantes prédéfinies telles que seule l'amplitude du moment magnétique des aimants 126 et 128 soit comprise entre ces limites. La valeur du seuil Dmax110 est choisi égale à la plus grande valeur possible de la distance d110. L'ensemble E3 comporte les équations supplémentaires décrites en référence à la figure 4. [0082] Le fonctionnement de l'interface 2 va maintenant être décrit plus en détail en référence au procédé de la figure 6. [0083] Ce procédé débute par une phase 140 d'initialisation lors de laquelle les lois L' les ensembles E, et les conditions C, sont enregistrées dans la base 38. [0084] Ensuite, on procède à une phase 142 de configuration de l'interface 2. Lors de cette phase 142, l'utilisateur choisit un ou plusieurs ustensiles parmi les ustensiles 10, 90 et 110. Ici, on suppose qu'il ne dispose que d'un seul exemplaire de l'ustensile 10 et de deux exemplaires des ustensiles 90 et 110. Ensuite, il les dispose librement sur la face 24 du dispositif 12. De plus, il s'assure que la configuration réalisée ne conduit pas à avoir simultanément plus de cinq aimants permanents devant la face 24. L'interface homme-machine 2 est alors configurée. La suite de cette description est faite dans le cas particulier où l'utilisateur a disposé sur la face 24 un exemplaire de chacun des ustensiles 10, 90 et 110. [0085] Il est alors possible de procéder à l'utilisation de cette interface 2 pour commander l'appareil 4. [0086] Cela débute par une étape 146 lors de laquelle les magnétomètres Mu mesurent simultanément le champ magnétique du ou des aimants simultanément présents devant la face 24. [0087] Ensuite, lors d'une étape 148, l'unité 30 estime à partir des mesures des magnétomètres Mu la position, l'orientation et l'amplitude du moment magnétique de chacun des aimants présents. [0088] Pour cela, lors d'une opération 150, l'unité 30 résout le système d'équations 5 du modèle M1 à un dipôle magnétique. Elle obtient un jeu de coordonnées x,, yi, z1, 01 et cp, et une amplitude Al correspondant à l'estimation de la position, l'orientation et l'amplitude d'un seul aimant permanent. [0089] Ensuite, lors d'une opération 152, l'unité 30 calcule une erreur d'estimation E1 représentative de l'écart entre : 10 - des valeurs estimées des mesures des magnétomètres, calculées à partir du système d'équations M1 et à partir des position, orientation et amplitude obtenues à l'issue de l'opération 150, et - les mesures des magnétomètres relevées lors de l'étape 146. [0090] L'erreur est estimée pour chaque magnétomètre. Ici, l'erreur d'estimation 15 globale résultant de l'utilisation du modèle M1 est par exemple obtenue en moyennant les erreurs calculées pour chacun des magnétomètres. [0091] Dans le cas où l'algorithme utilisé pour résoudre le modèle est un filtre de Kalman d'ensemble ou étendu, la résolution de ce système d'équations lors de l'opération 150 fournit également une estimation de cette erreur El. 20 [0092] Les étapes 150 et 152 sont réalisées pour P = 1 jusqu'à P = 5. De préférence, les étapes 150 et 152 pour chaque valeur de P sont réalisées en parallèle. [0093] Ensuite, lors d'une opération 154, l'unité 30 sélectionne le modèle Mp qui donne l'erreur Ep la plus petite. Ainsi, s'il n'y a qu'un seul aimant permanent devant la face 24, l'unité 30 sélectionne automatiquement le modèle M1. Si, au contraire, il y a 25 deux aimants permanents, l'unité 30 sélectionne automatiquement le modèle M2 et ainsi de suite. Ici, c'est le modèle M5 qui est automatiquement sélectionné car l'interface 2 comporte cinq aimants permanents. Ce modèle M5 comporte trente variables dont les valeurs sont inconnues. [0094] A l'issue de cette étape 148, le calculateur dispose d'une première estimation 30 des valeurs des variables x, y, z, ey, 02 et A pour chacun de cinq aimants. Cette première estimation est celle obtenue suite à l'exécution de l'opération 150 en utilisant le modèle M5. [0095] L'étape 148 se poursuite alors par une première phase 166 d'adaptation automatique du modèle M5 à la configuration réalisée par l'utilisateur lors de l'étape 35 142. [0096] Cette phase 166 débute par une étape 168 d'identification des ensembles E, qui peuvent être utilisés pour réduire le nombre de variables à estimer dans le modèle M5. A cet effet, lors de l'étape 168, le calculateur 32 identifie les ustensiles présents devant la face 24. [0097] Par exemple, lors d'une opération 170, le calculateur acquière la caractéristique distinctive de chaque ustensile présent devant la face 24. Ici, cette caractéristique distinctive est obtenue à partir de la précédente estimation de la position, orientation et amplitude de chaque aimant permanent. Dans ce cas particulier, cette caractéristique distinctive est soit l'amplitude du moment magnétique estimé d'un aimant permanent, soit les amplitudes estimées des moments magnétiques d'une paire d'aimants permanents. Les amplitudes estimées sont par exemple celles estimées lors de l'exécution de l'opération 150 avec le modèle M5. [0098] Ensuite, lors d'une opération 172, le calculateur 32 compare chaque 10 caractéristique distinctive acquise à chaque condition C, de la base 38. [0099] Si la caractéristique distinctive acquise correspond à l'une des conditions C, alors, lors d'une opération 176, un ustensile est identifié. A chaque fois qu'un ustensile est identifié, un identifiant de cet ustensile est ajouté à une liste d'ustensiles identifiés. Cette liste d'ustensiles identifiés contient également pour chaque ustensile 15 identifié l'ensemble E, et la loi L, associés à cet ustensile par la base 38. [00100] Si la caractéristique distinctive acquise ne correspond à aucune des conditions C' alors l'ustensile n'est pas identifié. [00101] A l'issue de l'étape 168, les équations supplémentaires utilisables pour diminuer le nombre de variables du modèle M5 sont donc identifiées. Toutefois, 20 chacune de ces équations supplémentaires, exprimée dans le repère Rm, dépend des coordonnées de l'origine et de l'orientation des axes du repère locale par rapport au repère Rm. [00102] Par conséquent, le procédé se poursuit par une étape 178 lors de laquelle, pour le repère local de chaque ustensile identifié, le calculateur 32 détermine : 25 - la position de son origine dans le repère Rm, et si nécessaire, - l'orientation de ses axes dans le repère Rm. [00103] Pour cela, le calculateur 32 utilise la position et l'orientation estimées du et des aimants permanents de l'ustensile dans le repère Rm et des informations préenregistrées sur la position du repère local par rapport aux aimants permanents 30 de cet ustensile. [00104] Par exemple, dans le cas de l'ustensile 10, la position de l'origine du repère R10 est déterminée en translatant la position estimée de l'aimant 16 dans le repère Rm d'une distance Rd dans la direction de son moment magnétique et en se rapprochant de la face 24. 35 [00105] Dans le cas de l'ustensile 90, la position de l'origine du repère Rgo dans le repère Rm est égale à la position estimée du centre géométrique de l'aimant permanent 100. Les orientations des axes Zgo et Xgo sont égales aux orientations estimées des moments magnétiques, respectivement, des aimants 100 et 98. [00106] De façon similaire, le calculateur 32 détermine la position et l'orientation du repère local R11o. [00107] A l'issue de l'étape 178, chacune des équations supplémentaires identifiées exprime une égalité entre une première variable du modèle M5 et un terme. La valeur de ce terme est exprimée à l'aide d'une ou plusieurs autres variables du modèle M5 différentes de la première variable. Typiquement, ces autres variables du modèle M5 sont celles qui codent la position et l'orientation des aimants permanents du même ustensile. [00108] Ensuite, pour chaque équation supplémentaire, lors d'une étape 180, le calculateur 32 remplace, dans le système actuel d'équations du modèle M5, la première variable par le terme auquel elle est égale. Ainsi, à l'issue de l'étape 180 on obtient un nouveau système d'équations et donc un nouveau modèle courant, noté M5c, comportant moins de variables à estimer que le système actuel. Pour les étapes qui suivent le système d'équations du modèle Mc devient le système actuel d'équations. [00109] Plus précisément le modèle M5 initial est un modèle générique qui permet de localiser simultanément cinq aimants permanents sans tenir compte des relations mécaniques qui peuvent exister entre ces aimants permanents. Ce modèle M5 est donc un système d'équations qui comporte trente variables à estimer. Après l'exécution de l'étape 180, le nouveau modèle Ms, est un modèle spécifique qui dépend des ustensiles actuellement disposés sur la face 24. Ce modèle M5, prend notamment en compte les relations mécaniques qui limitent le nombre de degrés de liberté en déplacement des aimants permanents les uns par rapport aux autres. Il comporte donc moins de variables à estimer que le modèle M5. Par exemple, ici, le modèle Ms, comporte seulement 19 variables à estimer au lieu de trente. En effet, avant même d'avoir déterminé si les ustensiles 10, 90 et 110 sont immobiles où non dans le repère Rm et d'avoir remplacé leur moment magnétique par des constantes, il est possibles de simplifier le système actuel d'équations en utilisant les trois équations supplémentaires de l'ustensile 10 et les quatre équations supplémentaires de chacun des ustensiles 90 et 110. [00110] Par la suite, les variables remplacées sont directement calculées à partir des équations supplémentaires identifiées et des valeurs estimées pour les variables du modèle M5c. [00111] Le procédé se poursuit par une étape 190 lors de laquelle les magnétomètres 35 Mu mesurent le champ magnétique présent. [00112] Lors d'une étape 192, le calculateur 32 acquiert ces mesures. [00113] Lors d'une étape 194, le calculateur 32 estime la position et l'orientation de chaque aimant permanent en résolvant le système d'équations du modèle M5c. Ce système d'équations est par exemple résolu comme décrit précédemment lors de l'opération 150. [00114] Lors d'une étape 196, le calculateur 32 calcule également l'erreur d'estimation actuelle Eerr liée à l'utilisation du modèle M5c. Cette erreur Eerr est calculée 5 comme décrit lors de l'opération 150 sauf que c'est le modèle Msc qui est utilisé. [00115] Ensuite, lors d'une étape 198, l'erreur Eerr est comparée à un seuil prédéterminé Serr. Si l'erreur Eerr est inférieure à ce seuil Serr, alors le calculateur 32 procède à une étape 200. [00116] Lors de l'étape 200, les lois L, sélectionnées et les positions et orientations 10 estimées lors de l'étape 194 sont utilisées pour générer une ou plusieurs commandes transmises à l'appareil 4. Ensuite, à l'issue de l'étape 200, le procédé retourne à l'étape 190. [00117] Si l'erreur Eerr est supérieure au seuil Serr, alors on procède à une étape 202 de réintroduction d'une ou plusieurs variables dans le modèle courant Msc pour faire 15 diminuer cette erreur Eerr. [00118] L'étape 202 peut être réalisée de nombreuses façons différentes. En effet, il existe de nombreuses méthodes pour réintroduire dans le modèle courant une ou plusieurs variables de manière à limiter l'erreur Eerr. Par exemple, une première méthode consiste à retourner systématiquement à l'étape 146. Dans ce cas, 20 l'ensemble des étapes à partir de l'étape 146 est réitéré à chaque fois que l'erreur Eerr dépasse le seuil Serr. Cette méthode présente l'avantage de permettre de réintroduire l'ensemble des variables des modèles M1 à M5. De plus, elle est simple à mettre en oeuvre. [00119] Une deuxième méthode consiste à réintroduire les variables les unes après 25 les autres dans l'ordre inverse où elles ont été éliminées du modèle courant. Ainsi, le calculateur 32 réintroduit dans le système actuel d'équations la dernière variable éliminée de manière à obtenir un nouveau système d'équations contenant une variable de plus. Ensuite, ce nouveau système d'équations est utilisé lors des mises en oeuvre suivantes des étapes 192 à 198. Cette deuxième méthode présente 30 l'avantage qu'il n'est pas nécessaire de réitérer systématiquement l'étape 148 et la phase 166 d'adaptation automatique du modèle. Un autre avantage de cette méthode est que la dernière variable fixée est la plus susceptible d'être modifiée. [00120] Une troisième méthode consiste à tester la réintroduction dans le modèle courant de chaque variable et de ne réintroduire que la variable ou les variables qui 35 permettent de diminuer l'erreur Eerr. Par exemple, pour cela, le calculateur 32 réintroduit d'abord une première variable dans le système actuel d'équations, puis vérifie si l'erreur Eerr calculée avec le nouveau système d'équations ainsi obtenu est inférieure à celle obtenue lors de la dernière itération de l'étape 196. Si l'erreur diminue, alors cette variable est sélectionnée. A l'inverse, si l'erreur n'est pas modifiée ou au contraire augmente, alors cette variable n'est pas sélectionnée. Le calculateur 32 teste ainsi chacune des variables du modèle courant Ms, les unes après les autres. Ensuite, il réintroduit dans le modèle courant Ms, toutes les variables sélectionnées. [00121] Si la mise en oeuvre de la deuxième ou de la troisième méthode conduit à réintroduire dans le modèle courant toutes les variables et que malgré cela, l'erreur Eerr reste supérieure au seuil Serr, alors, le procédé retourne à l'étape 148. En effet, cela peut signifier que l'augmentation de l'erreur Eerr est alors peut-être causée par le retrait ou l'ajout d'un nouvel ustensile sur la face 24 qui modifie le nombre d'aimants permanents présents. [00122] En parallèle des étapes 196 à 200, lors d'une étape 206, l'utilisateur manipule les ustensiles 10, 90 et 110 pour commander l'appareil 4. [00123] Egalement en parallèle des étapes 196 à 200, le calculateur 32 exécute une seconde phase 210 d'adaptation automatique du modèle courant. [00124] Cette phase 210 débute par une étape 211 d'identification d'une équation supplémentaire. Dans ce cas, lors de l'étape 211, l'équation supplémentaire est construite et non pas obtenue à partir d'informations préenregistrées sur les objets mobiles. Pour cela, lors d'une opération 212, le calculateur 32 calcule l'amplitude des variations de chaque variable du modèle Ms, au cours des Q précédentes itérations des étapes 190 à 194. Typiquement Q est un nombre entier supérieur à dix et de préférence supérieur à cent ou mille. Par exemple, lors de l'opération 212, le calculateur 32 calcule la variance de chacune des variables du modèle M5c. [00125] Ensuite, lors d'une opération 214, le calculateur 32 compare chaque variance calculée à un seuil prédéterminé respectif Svar. Par exemple, le seuil Svar est égal à une fraction de la variance maximale mesurée dans le passé pour cette variable. Par exemple, ce seuil Svar est pris égal à 10 % ou 5 % ou 1 % de cette variance maximale mesurée. Le seuil Svar peut aussi être une constante numérique prédéterminée. [00126] Si la variance calculée lors de l'opération 212 est supérieure à ce seuil Svar, alors on considère que la valeur de cette variable varie et qu'aucune nouvelle 30 équation supplémentaire n'est identifiée. [00127] Dans le cas contraire, le calculateur 32 considère que cette variable ne varie plus, ce qui permet d'établir une nouvelle équation supplémentaire. Cette équation supplémentaire est une relation d'égalité entre cette variable et une constante numérique. La valeur de cette constante numérique est par exemple prise égale à la 35 moyenne de plusieurs des précédentes valeurs estimées pour cette variable ou simplement prise égale à la dernière valeur estimée pour cette variable. [00128] Ensuite, lors d'une étape 216, les variables qui sont considérées comme étant des constantes sont remplacées dans le modèle Ms, par les valeurs de ces constantes. On obtient alors un nouveau système d'équations comportant encore moins de variables que le système actuel d'équations. Ce nouveau système d'équations est alors utilisé lors des exécutions suivantes des étapes 190 à 198. Il remplace donc le précédent système actuel d'équations. [00129] Ainsi, la phase 210 permet d'adapter le modèle courant en exploitant le fait qu'en pratique, un degré de liberté d'un ustensile n'est pas nécessairement utilisé par l'utilisateur. Les équations supplémentaires ainsi construites sont ajoutées à dans la liste des ustensiles et des équations supplémentaires déjà identifiés. Ces variables éliminées lors de la phase 210 peuvent être réintroduite dans le système d'équations lors de l'étape 202. [00130] Lorsque l'utilisateur souhaite reconfigurer l'interface 2. Il retourne à la phase 142. Lors de cette nouvelle exécution de la phase 142, il peut modifier la configuration de l'interface 2, par exemple : - en enlevant ou en ajoutant un ustensile, et/ou - en déplaçant un ustensile. [00131] Ainsi, l'interface 2 est facilement configurable tout en restant simple à fabriquer. [00132] De nombreux autres modes de réalisation sont possibles. Par exemple, la face avant 24 n'est pas nécessairement plane. Par exemple, en variante, il s'agit d'une forme en relief. Typiquement, cette forme en relief peut comporter des logements en creux destinés à recevoir et à immobiliser les ustensiles dans ces logements. [00133] D'autres ustensiles que ceux précédemment décrits sont réalisables et utilisables dans l'interface 2. Par exemple, on pourra se référer à ce sujet à la demande de brevet FR1354160 déposé le 7 mai 2013 par le Commissariat à l'énergie 25 atomique et aux énergies alternatives. [00134] De nombreuses autres caractéristiques distinctives et conditions de sélection associées peuvent être conçues. Par exemple, lorsque l'ustensile comporte K aimants permanents, avec K strictement supérieur à un, la ou les caractéristiques distinctives de l'ustensile sont choisies dans le groupe composé de : 30 - la position relative des K aimants permanents les uns par rapport aux autres, - l'orientation des moments magnétiques des K aimants permanents les uns par rapport aux autres, - les amplitudes des moments magnétiques des K aimants permanents ou les amplitudes relatives des moments magnétiques des K aimants permanents les unes 35 par rapport aux autres. [00135] Lorsque la caractéristique distinctive est la position d'un aimant permanent d'un ustensile par rapport à la position d'un autre aimant permanent du même ustensile, il n'est pas nécessaire d'utiliser l'amplitude des moments magnétiques pour sélectionner l'ensemble Ei à associer à cet ustensile. [00136] Lorsque l'ustensile comporte seulement un aimant permanent, la caractéristique distinctive de l'ustensile est l'amplitude du moment magnétique de son aimant permanent. [00137] La caractéristique distinctive peut aussi être la position d'un aimant 5 permanent dans le repère Rm. Dans ce cas, la condition de sélection vérifie si cette aimant permanent est situé à l'intérieur d'une zone prédéfinie dans le repère Rm. Dans l'affirmative une loi de commande et un ensemble d'équations supplémentaires sont associés à cet ustensile. Il est également possible de définir des conditions vérifiant si la position d'un aimant est non pas située à l'intérieur d'un volume 10 prédéfini mais dans une position particulière par rapport à un point fixe du repère Rm ou à un point dont la position est définie par un autre aimant permanent distinct de l'ustensile. [00138] La caractéristique distinctive peut aussi être l'inclinaison du moment magnétique d'un aimant permanent par rapport à un axe de référence du repère Rm, 15 par exemple la verticale. Dans ce cas, la condition de sélection teste si l'inclinaison estimée d'un moment magnétique d'un aimant permanent de l'ustensile est comprise dans une plage de valeurs prédéterminée. Dans l'affirmative, le calculateur identifie l'ensemble E, d'équations supplémentaires associé à cette condition de sélection comme permettant de réduire le nombre de variables du système actuel d'équations. 20 [00139] L'acquisition d'une caractéristique distinctive d'un objet peut être réalisée différemment de ce qui a été décrit. Par exemple, chaque objet est équipé d'une étiquette RFID, (Radio Frequency Identification). Dans ce cas, le dispositif 12 comporte en plus un lecteur d'étiquettes RFID. Ainsi, lorsqu'un ustensile est présenté devant le dispositif 12, celui-ci est identifié à partir de la lecture de l'identifiant de cet 25 ustensile contenu dans son étiquette RFID. Ensuite, cet identifiant est utilisé pour sélectionner dans une base de données préenregistrée l'ensemble E, associé à cet ustensile. Ce même identifiant peut aussi être utilisé pour sélectionner la ou les lois de commande à utiliser avec cet ustensile. Ainsi, la caractéristique distinctive n'est pas nécessairement une caractéristique construite à partir des mesures des 30 magnétomètres Mu. [00140] L'acquisition de la caractéristique distinctive n'est pas non plus nécessairement automatique. Par exemple, en variante, l'utilisateur indique, par l'intermédiaire d'une interface homme-machine l'identifiant du ou des objets ajoutés ou enlevés. 35 [00141] Dans un mode de réalisation où l'interface comporte uniquement un ou plusieurs exemplaires identiques du même ustensile, les conditions de sélection peuvent être omises. En effet, dans ce cas, l'ensemble E, d'équations supplémentaires associé à chaque ustensile est toujours le même. La base de données 38 peut alors contenir qu'un seul enregistrement. Il suffit simplement de dénombrer le nombre d'aimants permanents. [00142] La base de données 38 ne contient pas nécessairement, en plus des ensembles E' les lois de commandes L. Par exemple, la base 38 contient uniquement les conditions C, et les ensembles E, et une autre base de données contient des conditions C', et les lois L. Dans ce cas les conditions C, et C', ne sont pas nécessairement les mêmes. Ainsi, un ensemble E, n'est pas nécessairement et systématiquement associé à la même loi de commande. [00143] Si le nombre P d'aimants permanents simultanément présents devant la face 24 est connu à l'avance, le procédé de commande peut être simplifié en utilisant seulement le modèle Mp correspondant pour déterminer la position, l'orientation et l'amplitude du moment magnétique de ces P aimants permanents. Par exemple, le nombre P d'aimants permanents est saisi par l'utilisateur lors de la phase de configuration. L'étape 148 peut donc être simplifiée. [00144] Le changement de configuration de l'interface 2 n'est pas nécessairement réalisé directement à la main par un utilisateur. Dans une autre variante, l'interface comporte des actionneurs électriques permettant de modifier la configuration de l'interface. Par exemple, à cet effet, le déplacement des ustensiles ou le remplacement de ces ustensiles par d'autres ustensiles est motorisé, de sorte que l'utilisateur n'a pas lui-même à manipuler chacun de ces ustensiles. [00145] Il est possible de limiter le nombre d'emplacements où un ustensile peut être disposé par rapport à la face avant 24. Pour cela, l'interface 2 comporte plusieurs emplacements, équipé chacun de détrompeur propre à empêcher la fixation à cet emplacement d'un certain nombre d'ustensiles et, au contraire, à autoriser la fixation à cet emplacement d'autres ustensiles différents. [00146] L'approximation utilisée pour construire le modèle Mp peut aussi être une approximation quaternaire ou supérieure, c'est-à-dire que les équations du magnétisme sont approximées à un ordre supérieur à celui correspondant à l'approximation dipolaire. [00147] De nombreuses méthodes différentes sont utilisables pour déterminer la position et l'orientation de l'objet magnétique. Par exemple, la méthode décrite dans US6269324 est utilisable. Ces méthodes n'utilisent pas nécessairement un filtre de Kalman. Par exemple, les méthodes décrites dans US2002/171427A1 ou US6263230B1 sont possibles. [00148] Les magnétomètres du réseau de magnétomètres ne sont pas nécessairement rangés en colonnes et en lignes. Ils peuvent être agencés aussi selon d'autres motifs. Par exemple, les magnétomètres sont disposés sur chaque sommet de chaque maille triangulaire ou hexagonale d'un maillage d'un plan. [00149] La disposition des magnétomètres les uns par rapport aux autres peut aussi être aléatoire ou non régulière. Ainsi, la distance entre deux magnétomètres immédiatement consécutifs dans le réseau n'est pas nécessairement la même pour toutes les paires de deux magnétomètres immédiatement consécutifs. Par exemple, la densité de magnétomètres dans une zone donnée du réseau peut être plus importante qu'ailleurs. Augmenter la densité dans une zone donnée peut permettre d'augmenter la précision de la mesure dans cette zone. [00150] Le réseau de magnétomètres peut aussi s'étendre dans trois directions non colinéaires de l'espace. Dans ces conditions, les magnétomètres sont répartis à 10 l'intérieur d'un volume à trois dimensions. [00151] Le nombre N de magnétomètres peut également être supérieur ou égal à soixante-quatre ou quatre-vingt dix. [00152] Tous les magnétomètres du réseau de magnétomètres ne sont pas nécessairement identiques les uns aux autres. En variante, les magnétomètres n'ont 15 pas tous la même sensibilité. [00153] L'appareil 4 peut être remplacé par tout type d'appareil électrique qui doit être commandé en réponse à une action d'un être humain. Par exemple, l'appareil commandé peut être un robot, une machine-outil ou autre. [00154] Dans tous les modes de réalisation décrits ici, l'aimant permanent peut être 20 remplacé par un objet magnétique non alimenté en permanence et qui acquiert un moment magnétique en présent d'un champ magnétique continu extérieur, tel que le champ magnétique terrestre. Par exemple, l'aimant permanent est remplacé par une pièce en matériau magnétique doux. Un matériau magnétique est considéré comme doux si son champ magnétique coercitif est inférieur à 10 ou 1 A.m-1. Une telle pièce 25 présente un moment magnétique créé par l'interaction entre le champ magnétique terrestre et la pièce en matériau magnétique doux. [00155] De nombreux autres modes de réalisation du procédé de la figure 6 sont possibles. Par exemple, l'une ou l'autre des phases 166 et 210 est omise. [00156] Le nombre de mesures utilisées pour calculer la variance lors de l'opération 30 212 peut être fonction de la durée d'une fenêtre temporelle glissante. Typiquement, seules les mesures réalisées à l'intérieur de cette fenêtre temporelle glissante sont prises en compte pour calculer la variance. REVENDICATIONS1. Procédé de localisation d'objets magnétiques mobiles présentés devant un réseau de magnétomètres comportant N magnétomètres tri-axes (Mu) reliés mécaniquement les uns aux autres sans aucun degré de liberté pour conserver une distance connue entre chacun de ces magnétomètres, où N est un nombre entier supérieur ou égal à cinq, ce procédé comportant : a) la mesure (190) par chacun des magnétomètres de l'amplitude du champ magnétique le long de chacun de ses axes de mesure, b) l'estimation (194) à partir des mesures des magnétomètres, des valeurs de plusieurs variables en résolvant un système actuel d'équations reliant ces variables à chaque mesure d'un magnétomètre tri-axe du réseau, chaque variable correspondant à la position ou à l'orientation ou à l'amplitude d'un moment magnétique d'un objet magnétique, - l'itération des étapes a) et b) à des instants successifs dans le temps pour obtenir les valeurs estimées de chaque variable à ces différents instants successifs, caractérisé en ce que le procédé comporte : c) l'identification (168 ; 211), en fonction des objets mobiles présentés devant le réseau de magnétomètres, d'une équation supplémentaire reliant, par une relation d'égalité, une première variable du système actuel d'équations à un terme, ce terme étant une relation prédéfinie entre une ou plusieurs autres variables du système actuel d'équations ou une valeur numérique constante, d) le remplacement (180 ; 216), dans le système actuel d'équations, de cette première variable par le terme auquel elle est égale, identifié lors de l'étape c), pour 25 obtenir un nouveau système d'équations dans lequel le nombre de variables dont la valeur doit être estimée est plus petit que dans le système actuel d'équations, e) l'utilisation du nouveau système d'équations à la place du système actuel d'équations lors des itérations suivantes des étapes a) et b), le nouveau système d'équations devenant ainsi, pour les itérations suivantes des étapes a) et b), le 30 système actuel d'équations. 2. Procédé selon la revendication 1, dans lequel l'identification (168 ; 211) de l'équation supplémentaire est réalisée en fonction des valeurs estimées, lors d'une précédente itération des étapes a) et b), des variables. 35 3. Procédé selon la revendication 2, dans lequel l'identification de l'équation supplémentaire comporte : - le calcul (212) d'une grandeur représentative de l'amplitude des variations de la première variable au cours de précédentes itérations des étapes a) et b) à partir des 5 valeurs de cette première variable estimée au cours de ces précédentes itérations des étapes a) et b), - la comparaison (214) de cette grandeur calculée à un seuil prédéterminé pour établir si cette première variable a varié au cours de ces précédentes itérations des étapes a) et b), et 10 - si ce seuil prédéterminé est franchi, l'établissement (214) d'une nouvelle équation supplémentaire reliant, par une relation d'égalité, la valeur de cette première variable à une valeur constante, la valeur constante étant fonction d'au moins une valeur de cette première variable estimée lors de ces précédentes itérations des étapes a) et b) au cours desquelles il a été établi que sa valeur ne varie pas, puis l'utilisation de cette 15 nouvelle équation supplémentaire lors de l'étape d), et - si ce seuil prédéterminé n'est pas franchi, l'absence d'établissement de cette nouvelle équation supplémentaire. 4. Procédé selon l'une quelconque des revendications précédentes, dans lequel 20 l'étape c) comporte : - l'acquisition (170) d'une caractéristique distinctive de l'objet magnétique mobile actuellement présenté devant le réseau de magnétomètres, et - la comparaison (172) de la caractéristique distinctive acquise à des conditions de sélection préenregistrées dans une base de données, cette base de données 25 associant à chaque condition de sélection au moins une équation supplémentaire reliant, par une relation d'égalité, une variable du système actuel d'équations à un terme, ce terme étant une relation prédéfinie entre une ou plusieurs autres variables du système actuel d'équations ou une valeur numérique constante, et - uniquement si la caractéristique distinctive acquise correspond à l'une des 30 conditions de sélection préenregistrées, alors l'utilisation lors de l'étape d) de l'équation supplémentaire associée à cette condition de sélection. 5. Procédé selon la revendication 4, dans lequel l'acquisition d'une caractéristique distinctive comporte l'obtention (170) de cette caractéristique distinctive à partir des 35 valeurs des variables estimées à l'aide du système actuel d'équations. 6. Procédé selon l'une quelconque des revendications précédentes, dans lequel le procédé comporte : - le calcul (198) d'une erreur d'estimation avec le système actuel d'équations, l'erreur d'estimation étant représentative de l'écart entre : - des valeurs estimées des mesures des magnétomètres lorsque la position, l'orientation et l'amplitude du moment magnétique de chaque objet magnétique sont égales à celles estimées en utilisant le système actuel d'équations lors d'une itération des étapes a) et b), et - les mesures des magnétomètres relevées à l'étape a) lors de la même itération des étapes a) et b), et - si l'erreur d'estimation calculée franchit un seuil prédéterminé : - la sélection (202) d'au moins une des équations supplémentaires précédemment utilisées lors de l'étape d) pour obtenir le système actuel d'équations, - le remplacement (202), dans le système actuel d'équations, du terme de l'équation supplémentaire sélectionnée, par la variable que ce terme a remplacé pour obtenir un nouveau système d'équations dans lequel le nombre de variables dont la valeur doit être estimée est plus grand que dans le système actuel d'équations et l'utilisation du nouveau système d'équations à la place du système actuel d'équations lors des itérations suivantes des étapes a) et b), et - si l'erreur calculée ne franchit pas ce seuil prédéterminé, le maintien de l'utilisation du système actuel d'équations lors des itérations suivantes des étapes a) et b). 7. Procédé selon la revendication 6, dans lequel la sélection d'au moins une équation supplémentaire comporte : - la sélection d'au moins deux équations supplémentaires précédemment utilisées lors de l'étape d) pour obtenir le système actuel d'équations, - pour chacune des équations supplémentaires sélectionnées, le remplacement, dans le système actuel d'équations, du terme de l'équation supplémentaire sélectionnée, par la variable que ce terme a remplacée pour obtenir des troisième est quatrième systèmes d'équations différents, chacun de ces quatrième et troisième systèmes d'équations comportant un nombre de variables dont la valeur doit être estimée plus grand que dans le système actuel d'équations, - le calcul de l'erreur d'estimation avec ces troisième et quatrième systèmes d'équations, et- la sélection de l'équation supplémentaire, parmi celles initialement sélectionnées, qui minimise l'erreur d'estimation calculée. 8. Support (34) d'enregistrement d'informations, caractérisé en ce qu'il comporte 5 des instructions pour l'exécution d'un procédé conforme à l'une quelconque des revendications précédentes, lorsque ces instructions sont exécutées par un calculateur électronique. 9. Dispositif de localisation d'objets magnétiques mobiles, ce dispositif comportant : 10 - un réseau de magnétomètres comportant N magnétomètres tri-axes (M,J) reliés mécaniquement les uns aux autres sans aucun degré de liberté pour conserver une distance connue entre chacun de ces magnétomètres, où N est un nombre entier supérieur ou égal à cinq, - une unité (30) de traitement apte: 15 a) à acquérir les mesures de chacun des magnétomètres du réseau de magnétomètres, b) à estimer, à partir des mesures des magnétomètres, des valeurs de plusieurs variables en résolvant un système actuel d'équations reliant ces variables à chaque mesure d'un magnétomètre tri-axe du réseau, chaque variable correspondant à la 20 position ou à l'orientation ou à l'amplitude d'un moment magnétique d'un objet magnétique, - à réitérer les étapes a) et b) à des instants successifs dans le temps pour obtenir les valeurs estimées de chaque variable à ces différentes instants successifs, caractérisé en ce que l'unité (30) de traitement est également apte : 25 c) à identifier, en fonction des objets mobiles présentés devant le réseau de magnétomètres, une équation supplémentaire reliant, par une relation d'égalité, une première variable du système actuel d'équations à un terme, ce terme étant une relation prédéfinie entre une ou plusieurs autres variables du système actuel d'équations ou une valeur numérique constante, 30 d) à remplacer, dans le système actuel d'équations, cette première variable par le terme auquel elle est égale, identifié lors de l'étape c), pour obtenir un nouveau système d'équations dans lequel le nombre de variables dont la valeur doit être estimée est plus petit que dans le système actuel d'équations, e) à utiliser le nouveau système d'équations à la place du système actuel d'équations 35 lors des itérations suivantes des étapes a) et b), le nouveau système d'équations devenant ainsi, pour les itérations suivantes des étapes a) et b), le système actuel d'équations.