FR 3038729 A1 20170113 FR 1556501 A 20150709 1 Procédé de localisation d'un objet magnétique L'invention concerne un procédé permettant de déterminer la localisation et/ou l'orientation d'une source de champ magnétique dans l'espace. Les appareils de détection des mouvements dans l'espace par des 5 capteurs de champ magnétiques, ou magnétomètres, se composent généralement d'une série de capteurs de champ magnétiques, d'un ou plusieurs aimants et d'un ou plusieurs objets dont on veut détecter les mouvements. Une première méthode consiste à détecter le champ magnétique émis par 10 un aimant fixe. Le mouvement d'un l'objet induit des variations du champ magnétique qui sont détectées par des capteurs fixes. Ces variations sont ensuite ajustées à un modèle pour en déduire le mouvement de l'objet. Dans US 4,785,242, un mouvement à un degré de liberté (rotation) est détecté de cette manière. 15 Une seconde méthode consiste à lier de manière solidaire les sources de champ magnétique et les objets à détecter. Yan, L. et al. présentent dans An Orientation Measurement Method Based on Hall-effect Sensors for Permanent Magnet Spherical Actuators with 3D Magnet Array (Scientific reports, 4., 2014) une méthode pour détecter la rotation à trois degrés de liberté d'un rotor lié à des 20 aimants permanents. Le champ magnétique est détecté par des capteurs à effet Hall, et le mouvement est calculé par un ajustement des données basé sur une approximation exponentielle. Une limitation de ces méthodes de détection provient de la difficulté de 25 corréler les signaux des capteurs avec un modèle analytique ou numérique du champ magnétique en trois dimensions. W02013144342 divulgue un procédé et un dispositif de détermination de la localisation et de l'orientation d'un objet magnétique en utilisant un réseau de 30 capteurs de champ magnétique tri-axes, permettant au capteur de déterminer le sens et l'amplitude de l'induction magnétique au point de mesure. Lorsque le champ magnétique a été mesuré par tous ou une partie des capteurs, le dispositif 3038729 2 implémente un modèle mathématique associant chaque mesure de capteur à la position et à l'orientation de la source de champ magnétique dans l'espace. Le modèle utilisé effectue une approximation de la source de champ magnétique par un dipôle magnétique et est implémenté sous la forme d'un filtre de Kalman. L'utilisation d'un modèle dipolaire présente l'avantage d'être compatible avec la charge de calcul réalisable par un processeur usuel. En revanche, ce modèle possède une précision limitée. Il introduit en particulier des erreurs importantes si la source magnétique, de dimension maximum L, est située à une distance inférieure à 10L d'un capteur de champ magnétique. L'erreur entre le champ magnétique mesuré par le capteur et le champ magnétique selon un modèle dipolaire à l'endroit du capteur augmente quand la source de champ magnétique se rapproche du capteur, pour diverger quand la source de champ magnétique et le capteur magnétique tendent vers le même point. L'implémentation d'un modèle prenant en compte la géométrie exacte de la source de champ magnétique serait avantageuse pour calculer avec précision la localisation et l'orientation d'une source proche d'un capteur, située par exemple à une distance inférieure à 10L d'un capteur, mais peut représenter une charge de calcul incompatible avec la charge de calcul imposable au processeur pour un traitement des données en temps réel. L'invention a pour but de surmonter les inconvénients précités de l'art antérieur. Plus particulièrement, elle vise à surmonter les limitations des modèles dipolaires tout en maintenant un niveau de complexité computationnelle acceptable, notamment pour des applications en temps réel utilisant un processeur embarqué dans un appareil portable. Un objet de l'invention permettant d'atteindre ce but est un procédé de détermination de la position relative dans l'espace d'au moins une source de champ magnétique par rapport à au moins un capteur de champ magnétique comprenant 30 les étapes consistant à : a) acquérir des mesures du champ magnétique généré par au moins une dite source de champ magnétique au moyen d'au moins un capteur de champ magnétique ; 3038729 3 b) calculer au moyen d'au moins un processeur une solution de l'expression du champ magnétique généré par au moins une dite source de champ magnétique en modélisant chaque dite source de champ magnétique par un élément choisi parmi 5 une superposition de solénoïdes et une superposition de plans chargés, puis en estimant la valeur d'une intégrale elliptique complète liée au dit modèle par un algorithme utilisant une transformation de Landen ; c) calculer au moins un élément choisi parmi la position et l'orientation de 10 chaque dite source de champ magnétique en utilisant au moins une méthode choisie parmi : une minimisation du coût entre lesdites mesures desdits capteurs de champ magnétique et lesdites solutions calculées à l'étape b) ; un filtrage bayésien, à des instants successifs, de chaque dite mesure de 15 chaque dit capteur de champ magnétique, ledit filtrage étant fonction de chaque dite solution calculée à l'étape b). Avantageusement, on choisit d'utiliser au moins un dit filtrage bayésien dans l'étape c) , ledit filtre bayésien étant un filtre de Kalman. 20 Avantageusement, au moins une dite source de champ magnétique est en mouvement par rapport à au moins un dit capteur de champ magnétique et au moins un élément choisi parmi la position et l'orientation d'au moins une dite source de champ magnétique et la position et l'orientation d'au moins un dit capteur de 25 champ magnétique est calculé à partir d'un filtrage de Kalman desdites mesures d'au moins un dit capteur de champ magnétique, ledit filtrage de Kalman étant fonction de chaque dite solution calculée, à chacun desdits instants successifs. Avantageusement, au moins une dite source de champ magnétique peut 30 être placée à moins de dix centimètres d'au moins un dit capteur de champ magnétique. L'invention concerne aussi une interface homme/machine comportant au moins une dite source de champ magnétique, au moins un dit capteur de champ 3038729 4 magnétique, et au moins une unité de traitement, au moins une unité de traitement étant programmée de manière opportune pour la mise en oeuvre de chaque étape d'un dit procédé. 5 Avantageusement, l'interface homme/machine comporte une pluralité de dits capteurs de champ magnétique, dans laquelle l'ensemble des dits capteurs de champ magnétique est un réseau non plan. Avantageusement, au moins une dite source de champ magnétique de 10 ladite interface homme/machine est choisie parmi un solénoïde monocouche et un aimant cylindrique plein. Avantageusement, l'interface homme/machine comporte au moins deux dites sources de champ magnétique, au moins une dite source de champ 15 magnétique étant un solénoïde, et une dite unité de traitement de ladite interface est programmée pour commander une émission périodique de l'intensité du champ magnétique dudit solénoïde à une fréquence différente d'au moins une autre dite source de champ magnétique. 20 Avantageusement, au moins une dite source de champ magnétique de ladite interface est un aimant permanent composé d'un alliage d'éléments choisis parmi le fer, le vanadium, le néodyme, le bore, l'aluminium, le nickel, le titane, le cuivre, le samarium et le cobalt. 25 On rappelle que, dans les matériaux paramagnétiques ou diamagnétiques, le vecteur H, mesuré en ampère par mètre (A/m) et qui est aussi appelé « champ magnétique » dans de nombreux ouvrages, est lié au vecteur B, appelé induction magnétique, mesurée en Tesla (T), par la relation B=p(H).H + Br, où Br est l'induction rémanente. La plupart des développements qui suivent demeurent valables si on s'intéresse au champ H au lieu du champ B dans la zone de linéarité entre B et H. La description suivante présente plusieurs exemples de réalisation du dispositif de l'invention : ces exemples sont non limitatifs de la portée de l'invention. 3038729 5 Ces exemples de réalisation présentent à la fois les caractéristiques essentielles de l'invention ainsi que des caractéristiques additionnelles liées aux modes de réalisation considérés. Par souci de clarté, les mêmes éléments porteront les mêmes repères numériques dans les différentes figures. 5 L'invention sera mieux comprise et d'autres avantages, détails et caractéristiques de celle-ci apparaîtront au cours de la description explicative qui suit, faite à titre d'exemple en référence aux dessins annexés dans lesquels : 10 - la figure 1 présente une vue schématique de interface homme/machine ; la figure 2 présente la description d'un modèle ampérien d'une source de champ magnétique ; 15 - la figure 3 illustre l'utilisation du théorème de superposition dans une modélisation utilisée dans l'invention ; la figure 4 présente la description d'un modèle coulombien d'une source de champ magnétique ; 20 la figure 5 illustre l'utilisation du théorème de superposition dans une modélisation utilisée dans l'invention ; la figure 6 présente la mise en oeuvre du problème du calcul de champ 25 magnétique d'un solénoïde dans un référentiel cartésien ; la figure 7 présente une vue schématique du fonctionnement d'un filtre de Kalman. 30 Dans la suite du texte, le terme capteur de champ magnétique est définit par un capteur de champ magnétique tri-axe : la mesure du champ magnétique est effectuée sur trois axes différents de l'espace à l'endroit du capteur de champ magnétique. 3038729 6 La figure 1 présente une vue schématique de l'interface homme/machine 1. L'interface homme/machine 1 peut comporter une ou plusieurs sources de champs magnétique 3, un ou plusieurs capteurs de champ magnétique s 4, chacun reliés par un bus 8 à une unité de traitement 7. Le bus 8 est représenté sur la figure 1 par une double flèche : il peut être réalisé par une liaison de données filaire ou sans fil. Chacun des capteurs de champ magnétique 4 et des sources de champ magnétique 3 peut être mécaniquement indépendant, ou liés de manière solidaire. Dans la réalisation particulière illustrée par la figure 1, l'élément de repère 10 2 comporte un capteur de champ magnétique 4 et une source de champ magnétique 3, solidairement liés et fixes dans un repère (i,j,k) orthogonal. Dans la réalisation illustrée par la figure 1, une des deux sources de champ magnétique 3 n'appartient pas à l'élément de repère 2 et peut être déplacée librement par rapport au repère (i,j,k). De même, un des deux capteurs de champ magnétique 4 15 n'appartient pas à l'élément de repère 2 et peut être déplacé librement par rapport au repère (i,j,k). Pour réaliser une interface homme/machine 1, on peut par exemple lier de manière solidaire l'élément de repère 2 à une partie du corps humain, et lier le capteur de champ magnétique 4 et la source de champ magnétique 3 n'appartenant pas à l'élément de repère 2, à d'autres parties du corps humain, afin 20 d'en mesurer les déplacements respectifs par rapport au repère (i,j,k). Dans une réalisation particulière de l'invention, la source de champ magnétique 3 peut être un aimant permanent c'est-à-dire un aimant dont le moment magnétique est non nul en l'absence de champ magnétique extérieur. La forme de 25 la source de champ magnétique 3 peut être un cylindre plein ou un cylindre troué. La source de champ magnétique 3 peut être fabriquée dans un alliage d'éléments choisis parmi le fer, le vanadium, le néodyme, le bore, l'aluminium, le nickel, le titane, le cuivre, le samarium et le cobalt. 30 Dans une réalisation particulière de l'invention, les capteurs de champ magnétique 4 d'un même élément de repère 2, lorsque le dispositif en utilise un nombre supérieur ou égal à quatre, ne sont pas coplanaires. Cette disposition permet d'obtenir des informations sur le gradient de champ magnétique pour chacun des axes définis en figure 1. 3038729 7 La source de champ magnétique 3 peut aussi être un solénoïde parcouru par un courant électrique alternatif. 5 L'interface homme/machine 1 peut comporter plusieurs sources de champ magnétique 3. Dans le cas où plusieurs de ces sources de champ magnétique 3 sont des solénoïdes parcourus par des courants alternatifs, on peut imposer une fréquence de variation du courant électrique particulière à chacun des solénoïdes, dans le but de différentier les sources de champ magnétiques 3 lors de 10 leur détection. Plus généralement, dans des réalisations de l'invention, l'interface homme/machine peut comporter au moins une source de champ 3 réalisée par un solénoïde; l'unité de traitement 7 peut être programmée pour commander une 15 émission périodique de l'intensité du champ magnétique dudit solénoïde, à une fréquence différente d'au moins une autre fréquence d'émission du champ magnétique d'une autre source de champ magnétique 3. Cette caractéristique permet de distinguer l'origine du champ magnétique reçu par le ou les capteurs de champ magnétique (4) en analysant les différentes fréquences composant un ou 20 plusieurs signaux. Elle permet la localisation simultanée de plusieurs sources de champ magnétique 3. Dans une autre réalisation particulière de l'invention, la source de champ magnétique 3 est réalisée par le passage de courant électrique dans un 25 solénoïde monocouche. L'élément de mesure 2 peut, dans une réalisation particulière de l'invention, comporter un réseau de capteurs de champ magnétique 4. Les capteurs de champ magnétique 4 peuvent être placés en ligne et en colonne pour former une 30 matrice. Chaque capteur de champ magnétique 4 peut être fixé sans aucun degré de liberté aux autres capteurs de champ magnétique 4. A cet effet, les capteurs de champ magnétique 4 sont fixés sur la face supérieure d'un substrat 3038729 8 rigide. Chaque capteur de champ magnétique 4 mesure le sens et l'intensité du champ magnétique généré par une ou plusieurs sources de champ magnétique 3 à l'emplacement du capteur. Pour cela, chaque capteur de champ magnétique 4 mesure la norme de la projection orthogonale du champ magnétique généré par la 5 ou les sources de champ magnétique 3, à l'emplacement de ce capteur de champ magnétique 4, sur les trois axes de ce capteur de champ magnétique 4. Dans un exemple de réalisation, ces trois axes sont orthogonaux entre eux. Par exemple, un capteur de champ magnétique 4 va mesurer les composantes du champ magnétique sur des axes colinéaires à i,j et k. 10 Chaque capteur de champ magnétique 4 est raccordé à l'unité de traitement 7 par l'intermédiaire d'un bus 8 de transmission d'information. L'unité de traitement 7 est capable de déterminer la position et/ou 15 l'orientation de la ou des sources de champ dans le repère (i,j,k) de l'élément de mesure 2 à partir des mesures du ou des capteurs de champ magnétique 4. L'unité de traitement 7 comporte, entre autres, un ou plusieurs processeurs et mémoires. Le ou les processeurs sont programmés pour exécuter des instructions enregistrées dans la ou les mémoires. 20 En particulier, l'unité de traitement 7 associe un modèle mathématique avec l'ensemble ou à une partie des mesures des capteurs de champ magnétique 4 pour en déduire la position et/ou l'orientation de la ou des sources de champ magnétique 3. 25 L'unité de traitement 7 est aussi capable de restituer la localisation et/ou l'orientation d'une ou plusieurs sources de champ magnétique 3 à des plateformes numériques par l'intermédiaire de la sortie 9. 30 La figure 2 présente la description d'un modèle ampérien d'une source de champ magnétique 3. Ce modèle est utilisé par l'unité de traitement 7 de l'invention pour calculer la position et/ou l'orientation d'une ou plusieurs sources de champ magnétique 3 en fonction des mesures effectuées par les capteurs de champ magnétique 4. 3038729 9 Trémolet de Lacheisserie, É. (Magnétisme: Fondements, tome I de Collection Grenoble Sciences. EDP Sciences, 2000) montre que l'induction B générée par la matière aimantée et/ou par des courants électriques imposés de 5 densité jo , en tout point de l'espace se ramène à un problème d'électromagnétisme du vide dans lequel deux types de courants sont à considérer : les courants réels (ou libres) de densité volumique jo les courants associés à la matière aimantée (ou courants liés), de densité volumique j, et surfacique jrns données par : 10 jm = rot M (1) Ems = Mxn (2) 15 où n est le vecteur unitaire normal à la surface de la matière aimantée illustrée en figure 2 et orienté vers l'extérieur de la source de champ magnétique 3 et M est l'aimantation en tout point de la matière aimantée. Cette équivalence est appelée modèle ampérien de l'aimantation. 20 La relation (1) indique que dans le cas d'une aimantation uniforme, jm=0. Ainsi, un cylindre uniformément aimanté est équivalent à un solénoïde parcouru par une densité de courant surfacique Ijms1=IMI- La figure 3 illustre l'utilisation du théorème de superposition dans une 25 modélisation utilisée dans l'invention. Dans un exemple de réalisation de l'invention, la source de champ magnétique 3 est un cylindre troué. On peut utiliser le théorème de superposition en ramenant le problème au calcul du champ magnétique généré par deux solénoïdes emboités (modèle ampérien) et parcourus par des courants surfaciques opposés. Le modèle ampérien peut être utilisé par l'unité de traitement 30 7 de l'invention pour calculer la position et/ou l'orientation d'une ou plusieurs sources de champ magnétique 3 en fonction des mesures effectuées par les capteurs de champ magnétique 4. Le calcul de l'induction générée par l'aimant cylindrique troué peut être fait en sommant les inductions générées par la 3038729 10 superposition des inductions magnétiques des deux solénoïdes présentés dans la figure 3. La figure 4 présente une description du modèle coulombien. Ce modèle 5 peut être utilisé par l'unité de traitement 7 de l'invention, alternativement au modèle ampérien, pour calculer la position et l'orientation d'une ou plusieurs sources de champ magnétique 3 en fonction des mesures effectuées par les capteurs de champ magnétique 4. Trémolet de Lacheisserie, É. (Magnétisme: Fondements, tome de Collection Grenoble Sciences. EDP Sciences, 2000) montre que le champ 10 magnétique Hm produit par la matière aimantée est identique en tout point de l'espace à celui créé par des masses magnétiques fictives de densités volumique p, et surfacique am données par: 15 = -div M (3) = M.n (4) Le champ magnétique considéré est ici le champ magnétique dû à la matière (appelé parfois champ démagnétisant dans la matière aimantée, et champ de la matière hors de la matière aimantée). L'induction magnétique B s'obtient en 20 prenant en compte le champ magnétique ambiant Ho et l'aimantation M par la relation B=po(Ho+H,+M). Dans le cas d'un aimant de type FeNdB, le champ ambiant, supposé de l'ordre du champ magnétique terrestre dans cette réalisation de l'invention, sera négligé vis-à-vis de l'aimantation rémanente. Quand l'aimantation est uniforme, d'après (3), p,=0 ; on peut caractériser l'aimantation par 25 des densités de charge superficielle am.Par exemple, le champ Hm créé (dans la matière et en dehors) par un cylindre uniformément aimanté le long de son axe est le même que celui engendré par deux disques de densités surfaciques uniformes am=± M et situés sur chacune des faces planes du cylindre. 30 La figure 5 illustre l'utilisation du théorème de superposition dans une modélisation utilisée dans l'invention. Dans un exemple de réalisation de l'invention, la source de champ magnétique 3 est un cylindre troué. On peut utiliser le théorème de superposition en ramenant le problème au calcul du champ généré par deux paires de disques chargés uniformément (am=± , modèle coulombien). 3038729 11 Durand E. (Electrostatique et Magnétostatique, Tome I Les distributions, Edition Masson et Cie, p. 246 - 250) a présenté les calculs des inductions magnétiques générées pour des distributions de courants ou de masses 5 magnétiques fictives, en particulier pour un solénoïde parcouru par un courant surfacique kis. Comme décrit précédemment, un solénoïde d'axe principal z parcouru par un courant surfacique js est l'équivalent ampérien d'un aimant permanent cylindrique d'axe z, et d'aimantation uniforme M=Ijslz. Ainsi, la connaissance de l'induction magnétique générée en tout point de l'espace par un 10 solénoïde parcouru par un courant surfacique js permet de calculer l'induction magnétique générée en tout point de l'espace par un aimant cylindrique d'aimantation uniforme M=Ijslz. En coordonnées cylindriques, le centre du repère étant le centre 15 géométrique du solénoïde O, un point M étant repéré par ses coordonnées (r,0,z), l'induction produite en tout point de l'espace par le solénoïde est indépendante de la variable 0 en raison de la symétrie du problème. Les composantes Br et B, sont données par : goJs 2n-k [(k2 - 2)«k2) + 2E(k 2 Br = 20 Bz = e+ i 1 s v7r Zk [ K(k2) a+r n(h2, k2)1 (6) e- avec : 25 k2 = 4 (7) (a+ra)2r+z2 r1- = z ± L/2 (8) (5) 30 et K, E et H les intégrales elliptiques complètes de la première, seconde et troisième espèce: 3038729 12 5 i; , 1 K(m) - f v -0 t-msin29d0 Tt E(m) = fo \/1 - msin20 clO ir 1 II(n, m) = f 2 dB 0 (1-nsin20b11-msin2 (9 Une des caractéristiques de l'invention est de combiner l'utilisation des modèles ampérien et/ou coulombien présentés ci-dessus avec une résolution des 10 intégrales elliptiques (9), (10) et (11) par des algorithmes dit de Bulirsch (Bulirsch R., Numerical Calculation of Elliptic Integrals and Elliptic Functions, Numerische Mathematik 7, 78-90, 1965). Il s'agit d'algorithmes itératifs permettant l'évaluation de fonctions mathématiques communes, en l'occurrence des intégrales elliptiques, en utilisant le moins d'itérations possible. 15 Une implémentation avantageuse de cet algorithme, en termes d'efficacité, est présentée dans Derby, N., & Olbert, S. (2010), « Cylindrical magnets and ideal solenoids », American Journal of Physics, 78(3), 229-235. Cette implémentation permet de calculer la fonction C (pour Complete Elliptic Integral, 20 c'est-à-dire intégrale elliptique complète) définie par : C(kc,p,c,$) = fj/2 c cos2 yo+ s sin2 yo chp (cos2yo+psin2yo)lcos2cp+ksin2cp On peut alors réécrire les équations (5) et (6) selon les notations de 25 Derby et al.: Br(M) = Bo[a+C(k+,1,1, -1) - a_C(k_,1,1, -1)] Bz(M) = Bo a+a p[ig+C(k+, y2,1, y) - fi_ C(k_,Y2,1,Y)] (14) (12) (13) 30 3038729 13 11BrIl Ill/on/II Bo = = W 71- 7 a ,\Iz_2E+(p+a)2 z+ ,\Iz_2E+(p+a)2 a-p Y = a+p k± jz_+(a-p)2 z_+(a+p)2 10 La fonction C peut être calculée efficacement en utilisant les algorithmes de Bulirsch, traduisant numériquement une transformation de Landen. La publication de Bulirsch considère trois cas distincts. Dans le cas du 15 calcul d'une intégrale complète de première espèce, qui peut s'écrire sous la forme : cell(kc) = f ( g o (20) ,\1(1-g2)(1-Ficc22) l'algorithme traduit un calcul intégral numérique utilisant une transformation de Landen. Dans le cas d'une intégrale complète de seconde espèce (incluant le cas 20 des intégrales de première espèce), qui peut s'écrire sous la forme : ce12(kc,a,b) = fo g (1-F2),\1(1-g-2)(1+kc22) a+g 2 (21) l'algorithme traduit aussi un calcul intégral numérique utilisant une transformation de 25 Landen. On peut par ailleurs noter que pour un argument complexe de l'intégrale elliptique, la résolution des intégrales elliptiques de première et seconde espèce utilise une transformation de Gauss car la transformation de Landen n'est pas stable numériquement dans ce cas, mais la plupart des réalisations de l'invention utilisent des transformations de Landen. a+ = 5 fi+ = 3038729 14 Enfin, un troisième algorithme est proposé dans le cas d'une intégrale elliptique complète de troisième espèce, c'est-à-dire de la forme : 5 ce12(kc,m,p) = foc° . (22) (1-F7g2).\1(1-g.2)(1+kc22) Ce troisième algorithme traduit lui aussi un calcul intégral numérique utilisant une transformation de Landen (alors que la forme plus générale, dans le cas d'une intégrale elliptique non complète est calculée en utilisant une 10 transformation de Gauss). Dans les réalisations de l'invention, le calcul des intégrales K(m), E(n) et I-I(n,m) ou C(kc, p, c, s) correspond au calcul d'intégrales complètes : les réalisations de l'invention utilisent des transformations de Landen pour calculer avantageusement ces intégrales. 15 Le calcul de la fonction C demande en moyenne cinq itérations pour converger à une précision de 10-12 ce qui peut permettre d'utiliser cette méthode pour modéliser, par exemple, le champ magnétique généré par des solénoïdes de tailles différentes sans utiliser l'approximation dipolaire. Les inventeurs ont découvert que l'utilisation de cette méthode de calcul du champ magnétique pouvait 20 permettre de rendre la charge de calcul acceptable pour une unité de traitement usuelle, ou un processeur usuel, tout en modélisant le champ magnétique à des distances très proches de la ou des sources de champ magnétique. Par exemple, dans le cas d'une source magnétique d'une dimension caractéristique d'un centimètre, une réalisation de l'invention permet de modéliser le champ magnétique, 25 grâce au calcul des intégrales elliptiques complètes par une transformation de Landen, à une distance inférieure à dix centimètres, ce qui est impossible en utilisant un modèle dipolaire des sources de champ magnétique. Lors de la mesure de la position et/ou de l'orientation relative d'une source de champ magnétique 3 par rapport à un capteur 4, il est possible de rapprocher la source 3 et le capteur 4 à 30 une distance inférieure à 10 cm. Dans d'autres réalisations de l'invention, la modélisation du champ magnétique 3 est faite par un modèle coulombien. On peut calculer des intégrales elliptiques complètes équivalentes à K, E et H par les algorithmes de Bulirsch. 3038729 15 La figure 6 présente la mise en oeuvre du problème de calcul d'une source champ magnétique 3 dans un référentiel cartésien. Le panneau A de la figure 6 présente la configuration générique du problème de calcul du champ 5 magnétique. Les résultats des expressions du champ magnétique obtenus par un algorithme de Bulirsch peuvent être mis en oeuvre pour calculer les coordonnées de l'induction magnétique B en un point M quelconque de l'espace à partir de la connaissance de paramètres de la source de champ magnétique 3 (moment magnétique, dimensions et position). 10 Les formules (13) à (19) sont exprimées dans le repère cylindrique local de la source de champ magnétique 3. Dans la figure 3, on traite un cas général où le repère absolu est quelconque. 15 Pour passer du repère global w(O,i,j,k) à un repère local relié au solénoïde S(S,i',j',k'), il faut effectuer une translation et deux rotations. La translation est effectuée avec le vecteur ps puis une rotation autour de l'axe k permet de placer le solénoïde dans le plan (1k). Enfin une rotation autour de j permet de dresser le solénoïde de telle façon que son axe soit confondu avec l'axe 20 k. La position du point M dans le repère S est donnée par Ms : MS = Ry(-6)R,(-(p)(0M - OS) (23) Le point M est décrit dans le repère cartésien lié à la source magnétique. 25 Le panneau B de la figure 6 illustre l'alignement du moment magnétique avec l'axe k. On peut placer M sur le plan (1k) en appliquant une rotation à M. On obtient alors le point M' de coordonnées (p,0,z) données par : (p, 0, z) = Rz(-(pOms (24) L'effet de cette rotation est présenté par le panneau C de la figure 6. 30 Une fois que p et z sont calculés, il est possible d'utiliser (16) et (17) pour calculer le 3038729 16 champ magnétique au point projeté sur (ik). Une fois le champ magnétique calculé, on peut effectuer les rotations inverses afin de déterminer l'induction magnétique Bw dans le repère initial : 13(') = Rz((p)Ry((p)R,((p0(Bp, 0, Bz) (25) Dans des réalisations de l'invention, les matrices de rotation sont 5 générées par des opérations géométriques ne faisant pas intervenir des fonctions trigonométriques de type arc cosinus ou arc tangente pour des raisons de stabilité numérique et de rapidité d'exécution matérielle. On peut avantageusement utiliser une seconde méthode de passage du repère global w au repère local S relié à la source de champ magnétique 3. On définit par Pmis = OM - OS la position relative du point de mesure M par rapport au centre du solénoïde S. On définit le vecteur M comme le moment de la source de champ magnétique 3 (colinéaire à l'axe de la source de champ magnétique 3). On définit Mn = M/IIMII le moment normalisé de la source de champ magnétique 3. On 15 définit Pm/s,' = Pm/s/ II PII le vecteur unitaire de S vers M. Dans le repère cylindrique lié à la source de champ magnétique 3, les coordonnées du point M sont (x,y,z) = (0,0,P. Mn). Le vecteur d'induction magnétique B, dans le repère global w est alors défini après le calcul du vecteur d'induction magnétique Bs dans le repère local S par B, = Bs. Mn + Bs.M x (-n - P m/sn X Mn). 20 L'invention utilise des méthodes d'estimation qui utilisent d'une part les modèles analytiques d'induction introduits précédemment, et d'autre part, des mesures issues de capteurs. L'utilisation conjointe de ces méthodes permet alors de calculer les positions et/ou les orientations relatives d'une ou plusieurs source(s) de 25 champ magnétique 3 par rapport à un ou plusieurs capteurs de champ magnétique 4. Elle peut aussi être utilisée pour estimer la position et/ou l'orientation d'un ou plusieurs capteurs magnétiques 4 par rapport à une ou plusieurs sources de champ magnétique 3. 30 Dans les réalisation de l'invention, une étape intervient après l'acquisition des mesures du ou des capteurs de champ magnétique 4 et après le 3038729 17 calcul par l'unité de traitement 7 de l'expression du champ magnétique 3 généré par au moins une source de champ magnétique 3 en modélisant les sources 3 par une superposition de solénoïdes et/ou une superposition de plans uniformément chargés, comme décrit précédemment. Cette étape consiste à estimer la 5 localisation et/ou de l'orientation d'une source de champ magnétique 3 et/ou d'un capteur de champ magnétique 4 suivant une méthode choisie entre une optimisation, un filtrage bayésien ou une méthode hybride entre une optimisation et un filtrage bayésiens. 10 Dans le cas de l'optimisation, on cherche à trouver un vecteur de paramètres x (par exemple la position et/ou l'orientation de la source de champ magnétique 3, le moment magnétique de la source de champ magnétique 3, ou la position et/ou l'orientation du capteur de champ magnétique 4) qui minimisent un coût donné. En l'occurrence, l'optimisation consiste à minimiser les coûts entre les 15 mesures desdits capteurs et les solutions calculées lors de la modélisation du champ magnétique. Par exemple, on peut utiliser la méthode des moindres carrées: x = argmin,((z - h(x))T (z - h(x))) (26) Ici z est la mesure réalisée par le/les capteurs de champ magnétique 4. 20 h(x) = 13,0(x) est le modèle proposé dans la section précédente. Pour résoudre ce problème on peut par exemple utiliser la méthode de Gauss-Newton ou la méthode de Levenberg-Marquardt. Dans le cas du filtre bayésien, on peut utiliser l'approximation 25 Markovienne d'ordre 1 pour estimer la position d'une ou plusieurs sources de champ magnétique 3 par rapport à ou aux capteur(s) de champ magnétique ou inversement. Cette opération peut être effectuée en deux étapes : - une prédiction utilisant l'état à l'instant précédent utilisant un 30 modèle dynamique d'évolution d'état pour prédire l'état actuel, 3038729 18 une mise à jour basée sur l'innovation apportée par les mesures données par les capteurs de champ magnétique 4 par rapport aux mesures prédites, calculées sur la base du modèle ampérien exposé précédemment. 5 De façon générale, pour estimer l'état a posteriori à partir d'un a priori et d'une mesure, le théorème de Bayes peut être utilisé: P(zt+i ixt+i, zi1)13(xt+i I zi.t) vraisemblance x a priori (27) P (zt+ilzil) 10 Le calcul de la vraisemblance fait intervenir une fonction h qui lie l'état à la mesure. Ici l'état est la position et/ou l'orientation de la source magnétique et la mesure est le champ magnétique mesuré par le ou les capteurs de champ magnétique 4. Ce principe s'applique à tous les estimateurs Bayésiens. En 15 particulier, un filtre de Kalman peut être utilisé pour résoudre ce problème. Lors de ce calcul, on peut estimer le vecteur moment magnétique au lieu de l'orientation pour éliminer certains calculs trigonométriques. La figure 7 présente une vue schématique du fonctionnement d'un filtre 20 de Kalman. Le filtre de Kalman, permet de trouver l'état à l'instant t+1 en fonction de l'état à l'instant précédant en deux opérations (prédiction et mise à jour). En particulier dans le cas où on utilise une fonction de mesure non linéaire (tel est le cas du modèle proposé), on peut recourir à la version étendue du filtre. 25 Dans un filtre de Kalman, l'état est décrit par les deux premiers moments de sa distribution de probabilité. Comme cette dernière est supposée Gaussienne, les deux moments sont suffisants pour la décrire complètement. On appelle xt+i l'espérance de l'état et Pt+1 sa matrice de covariance. évidence 3038729 19 Un filtre de Kalman est dit étendu s'il réutilise les mêmes équations utilisées dans un filtre de Kalman classique en linéarisant localement la fonction de mesure : 11(x+1) = Ht+1Xt+1 (28) 5 La fonction est calculée à partir de (25). Ht.0 est la matrice jacobienne évaluée en lt+iit. La phase de prédiction pour un modèle d'évolution statique est : xt+iit Fxtit (29) Pt+lit = FPtitFT + Qt 10 où F est une matrice identité et Q est la matrice de covariance du bruit d'état. Elle décrit l'incertitude sur l'évolution de l'état. La phase de mise à jour permet d'utiliser le modèle proposé pour corriger la prédiction à partir de l'information apportée par la mesure actuelle, appelée innovation 1+1. On a : Xt+iit+i Kt-E1S1t+1 (30) Pt+11t+1 Pt+ii t Pt+iit 15 Ici, Hun = Vx[h(xt+i)] est la matrice jacobienne de la fonction de mesure h évaluée en xt+iit. Kt+1 est le gain de Kalman. L'innovation et le gain de Kalman sont calculés comme suit : 20 Slt+i = zt+1 h(xt+i Kt+1 = Pt+iitHt+1TIjitl+131t+i Ii1+137-t+i Ht+iPt+iitHt+iT t Rt+1 3038729 20 Dans des réalisations de l'invention, la position et/ou l'orientation relative d'une source de champ 3 par rapport à la position et/ou l'orientation d'un capteur 4 sont calculées à partir d'un filtrage de Kalman tel que décrit précédemment. Le même principe s'applique pour le filtre de Kalman de type UKF, ainsi que les filtres à 5 particules et tous les autres filtres Bayésiens. On peut avantageusement utiliser un algorithme appelé LMA (décrit dans Aloui, S., Villien, C., & Lesecq, S. (2014), « A new approach for motion capture using magnetic field: models, algorithms and first results. » International Journal of 10 Adaptive Control and Signal Processing) à la place du filtre de Kalman. Cette approche est basée sur l'utilisation d'un algorithme hybride Bayésien/Optimisation de type Levenberg-Marquard. Lors de l'utilisation de l'algorithme LMA, on effectue dans un premier temps une étape de prédiction, similaire à l'étape de prédiction décrite précédemment lors de l'utilisation d'un filtre de Kalman. Dans un deuxième 15 temps, on effectue une étape de mise à jour, similaire à la méthode de Kalman décrite précédemment mais en remplaçant les équations d'optimisation du coût par : Chaque itération i est calculée par : 11 t+1 xt+0+1 + At-F-0+1 (35) où : (P +111,-1) A'diag((P t+1 (36) 25 et : 20 Z . -- h( Xf,_ ) - lit Rï+11 Xt+:1[D+1 = 11/i Htt- 1)-rRt-11 (zt - h (xi-1 )) ( Xt +11t+1) (37) _ H2tir-Ri+11H2t-+11 , (r-t+11t) (38) où À est le facteur d'amortissement, non constant. Le facteur d'amortissement est initialisé à .1° = 1. Pour chaque itération i, l'erreur ei e2112 est calculée. Si 3038729 21 xti lit est laissé tel quel et À est diminué (par exemple `' A). -01-1 _ -1 Sinon, est rejeté (Xt+1It = xt+11t) et À est augmenté (par exemple ,V+1 = 10)x\ i- 5 La matrice de covariance de la distribution estimée a posteriori est Pt+lit+1 = Elle est calculée comme dans l'exemple précédent du filtre de Kalman. Cette méthode fait converger le calcul vers la borne de Cramér-Rao (borne minimale de l'erreur d'estimation). REVENDICATIONS1. Procédé de détermination de la position relative dans l'espace d'au moins une source de champ magnétique (3) par rapport à au moins un capteur de champ magnétique (4) comprenant les étapes consistant à : a) acquérir des mesures du champ magnétique généré par au moins une dite source de champ magnétique (3) au moyen d'au moins un capteur de champ magnétique (4); b) calculer au moyen d'au moins un processeur une solution de l'expression du champ magnétique généré par au moins une dite source de champ magnétique (3) en modélisant chaque dite source de champ magnétique (3) par un élément choisi parmi une superposition de solénoïdes et une superposition de plans chargés, puis en estimant la valeur d'une intégrale elliptique complète liée au dit modèle par un algorithme utilisant une transformation de Landen ; c) calculer au moins un élément choisi parmi la position et l'orientation de chaque dite source de champ magnétique (3) en utilisant au moins une méthode choisie parmi : une minimisation du coût entre lesdites mesures desdits capteurs de champ magnétique (4) et lesdites solutions calculées à l'étape b) ; un filtrage bayésien, à des instants successifs, de chaque dite mesure de chaque dit capteur de champ magnétique (4), ledit filtrage étant fonction de chaque dite solution calculée à l'étape b). 2. Procédé selon la revendication précédente dans lequel au moins le filtrage bayésien est choisi comme méthode pour l'étape c), ledit filtre bayésien étant un filtre de Kalman. 3038729 23 3. Procédé selon l'une des revendications précédentes dans lequel au moins une dite source de champ magnétique (3) est en mouvement par rapport à au moins un dit capteur de champ magnétique (4) et dans lequel au moins un élément choisi parmi la position et l'orientation d'au moins une dite source de champ magnétique (3) et la position et l'orientation d'au moins un dit capteur de champ magnétique (4) est calculé à partir d'un filtrage de Kalman desdites mesures d'au moins un dit capteur de champ magnétique (4), ledit filtrage de Kalman étant fonction de chaque dite solution calculée, à chacun desdits instants successifs. 4. Procédé selon l'une des revendications précédentes dans lequel au moins une dite source de champ magnétique (3) peut être placée à moins de dix centimètres d'au moins un dit capteur de champ magnétique (4). 5. Interface homme/machine comportant au moins une source de champ magnétique (3), au moins un capteur de champ magnétique (4), et au moins une unité de traitement (7) programmée de manière opportune pour la mise en oeuvre de chaque étape de l'une des revendications précédentes. 6. Interface homme/machine selon la revendication précédente comportant une pluralité de capteurs de champ magnétique (4) formant dans leur ensemble un réseau non plan. 7. Interface homme/machine selon l'une des revendications 5 à 6 dans laquelle au moins une dite source de champ magnétique (3) est choisie parmi un solénoïde monocouche et un aimant cylindrique plein. 8. Interface homme/machine selon l'une des revendications 5 à 7 comportant au moins deux sources de champ magnétique (3), dans laquelle au moins une desdites sources de champ magnétique est un solénoïde, et dans laquelle l'unité de traitement (7) est programmée pour commander une émission périodique de l'intensité du champ magnétique dudit solénoïde à une fréquence différente d'au moins une autre desdites sources de champ magnétique (3). 3038729 24 9. Interface homme/machine selon l'une des revendications 5 à 8 dans laquelle au moins une dite source de champ magnétique (3) est un aimant permanent composé d'un alliage d'éléments choisis parmi le fer, le vanadium, le néodyme, le bore, l'aluminium, le nickel, le titane, le cuivre, le samarium et le cobalt.