FR 2734900 A1 19961206 FR 9506528 A 19950601 La présente invention se rapporte à la détermination d'une orientation liée à un système mobile, par exemple à un viseur de casque et a plus particulièrement pour objet un procédé, de détermination d'une orientation liée à un système mobile permettant d'obtenir une grande précision. Le problème consiste à repérer avec précision la direction d'un axe orienté sans avoir de contact matériel avec lui, en particulier dans les systèmes viseurs de casque utilisés dans les avions de combat ou les hélicoptères. Le fonctionnement d'un tel système est rappelé brièvement ci- après: par un dispositif collimateur annexe, le pilote voit, à travers une glace semi-réfléchissante solidaire du casque et interposée sur un axe de vision, l'image d'un réticule projeté à l'infini en superposition avec la scène extérieure. Lorsqu'il veut désigner une cible, le pilote amène ce réticule à coïncider avec la cible et signale que la coïncidence est réalisée au moyen par exemple d'une commande par bouton poussoir prévue à cet effet. A condition de repérer, au moment de la signalisation, la position angulaire exacte du casque, il est possible de déterminer par rapport à I'avion, la direction de visée et de désigner l'objectif à un système d'arme, ou de pointer dans cette direction un système optique ou autre. On connaît notamment par le brevet français n 79 14441 un dispositif de mesure d'orientation, notamment pour système viseur de casque, dans lequel, pour ramener la ligne de visée déterminée par la position du casque du pilote dans un repère lié à l'aéronef, un capteur de position constitué de trois bobines électromagnétiques orthogonales est placé sur le casque supposé rigide; par ailleurs, en un point fixe de la cabine, sont placées trois autres bobines électromagnétiques, formant un émetteur. Le procédé consiste alors à faire passer successivement, dans chaque bobine de l'émetteur, formant le trièdre fixe, un courant électrique. Ces courants font apparaître successivement, trois champs magnétiques associés qui sont captés par les bobines du capteur, formant le trièdre mobile lié au casque. Pour chaque champ magnétique émis les trois composantes orthogonales correspondantes sont captées dans les bobines du capteur. Les neuf composantes résultantes permettent de passer du trièdre mobile au trièdre fixe. En effet, ces neuf composantes dépendent: - de la position du capteur dans l'espace; - de l'orientation du capteur. En fait, la transformation résultant d'une telle méthode peut être assez différente de la transformation des coordonnées permettant de passer effectivement du trièdre de référence mobile au trièdre de référence fixe par rapport à l'aéronef, notamment du fait que l'on utilise des champs magnétiques. En effet, ces champs peuvent être notablement perturbés du fait de courants électriques circulant dans tous les éléments métalliques de l'environnement de la cabine, par exemple des courants de Foucault. De plus, des défauts inhérents à la réalisation des bobines d'émission ou de réception magnétiques créent également des modifications par rapport au champ dipolaire habituellement pris comme modèle de source rayonnée. La détermination d'une direction en utilisant les lois de l'électromagnétisme en espace libre, sans perturbations n'est donc pas possible quand l'espace considéré est la cabine de pilotage d'un aéronef, ou un environnement quelconque formé par des pièces conductrices. C'est pourquoi certains dispositifs connus prévoient, avant chaque mesure, la mise en mémoire de mesures préalables et la comparaison de mesures effectuées à un instant donné avec les mesures mises en mémoire afin de déterminer la rotation qui fait passer du trièdre fixe au trièdre mobile. Une telle méthode nécessite la mise en mémoire d'une quantité très importante de mesures qui doivent être reprises intégralement dès que des paramètres tels que, par exemple, la position exacte de l'émetteur mobile par rapport à la ligne de visée, changent. D'autres procédés permettant d'améliorer les performances ont été réalisés comme dans le brevet français 90 08230. Toutefois, comme il est mentionné dans ce brevet, la méthode de modélisation prend globalement en compte trois sources de génération des champs magnétiques a savoir: - les antennes d'émission et de réception, y compris leurs défauts - les perturbations dues aux courants de Foucault circulant dans les parties conductrices du cockpit, parfois décrites comme des réflexions des champs émis par l'émetteur - les perturbations dues aux parties conductrices et magnétiques fixées sur le casque du pilote, et donc liées au repère capteur magnétique. Le principal inconvénient des méthodes de l'état antérieur est le manque de précision et de robustesse, c'est-à-dire une sensibilité trop importante, vis-à-vis des variations des champs émis par les trois sources précédentes dépendant elles-mêmes: - des variations de position et orientation de l'émetteur dans un repère de référence - des variations géométriques des cockpits se traduisant par des variations de position et orientation des perturbations, de plaques métalliques par exemple, par rapport à un repère mécanique de référence - des variations des positions ou des orientations des pièces perturbantes fixées sur le casque (casques de différentes tailles, personnalisation des réglages inter-pupillaires pour chaque pilote...). En ce qui concerne la cabine (cockpit) des avions ou hélicoptères, il faut distinguer deux types de variation: - des variations dans le temps d'une même cabine (vibrations, efforts au cours des vols,...) - les variations par rapport à une cabine de référence. La seconde cause est la plus importante et la plus difficile à résoudre. La méthode indiquée dans le brevet 90 08230 ne résout pas ces problèmes et de plus suppose que la cartographie magnétique faite initialement est applicable sur toutes les cabines de même type. Or, il est très difficile d'obtenir une précision très importante en restant robuste par rapport à des changements de casques ou seulement leur démontage/remontage pour un même ou plusieurs cockpits, et à des changements de cockpits pour un ou plusieurs casques différents, c'est-à- dire un ou plusieurs pilotes notamment. Les précédentes méthodes connues ne résolvent pas les problèmes d'interchangeabilité, ou ne présentent pas toutes les garanties de maintien des performances par changement de cabines ou/et changement de casque. La présente invention a pour but de pallier les inconvénients précités. A cet effet, I'invention a pour objet un procédé de détermination d'une position et d'une orientation d'un système mobile par rapport à un émetteur de champ électromagnétique, le système mobile placé dans un porteur étant lié à un capteur de champ magnétique, caractérisé en ce qu'il comporte: - une première étape de modélisation analytique des champs électromagnétiques en fonction des coordonnées du capteur, un premier champ (Bo) étant crée par l'émetteur, un deuxième champ (B1) étant crée par les courants électriques induits dans le porteur par le premier champ (Bo) et un troisième champ (B2) étant crée par les courants électriques induits dans le système mobile par les deux premiers champs (Bo +Bl), l'effet électromagnétique de chaque champ (Bo, B1, B2) étant caractérisé indépendamment des effets des autres champs par les coefficients de son propre modèle; - une deuxième étape de calcul, en temps réel, de la position et de l'orientation du capteur en utilisant la mesure courante du champ électromagnétique au niveau du capteur et en utilisant les modèles (Bo, Bl, B2) de la première étape, la position et l'orientation du capteur étant définie à partir du champ mesuré (B,,,) auquel le troisième champ (B2) a été retranché. L'invention a pour principaux avantages qu'elle permet: - d'améliorer les performances des méthodes précédentes et aussi celle décrite dans le brevet 90 08230, tout en étant plus robuste par rapport aux déformations temporelles du cockpit; - de rendre le système insensible aux variations d'environnement (cockpit, casque). La méthode du brevet 90 08230 nécessite la mise en oeuvre de mesures cartographiques pouvant durer 24 à 30 heures, plus de 20 000 points de mesure étant nécessaires, pour utiliser des ordres de polynômes élevés afin d'assurer la précision requise, I'aéronef étant alors immobilisé plus d'une semaine; - d'identifier séparément les 3 sources de perturbation de manière simple et surtout quasi "transparentes" pour l'utilisateur, c'est-à-dire sans que cela augmente le temps d'immobilisation de l'aéronef nécessaire à une harmonisation. L'harmonisation étant définie comme l'opération permettant soit de régler un référentiel de mesure parallèle au référentiel de l'aéronef, soit de mesurer les coordonnées angulaires de ses vecteurs unitaires dans ce même référentiel de l'aéronef. Une cartographie d'une demi-heure à deux heures suffisant pour la mise en oeuvre du procédé objet de l'invention. Pour des cabines différentes d'une même série, I'harmonisation de l'émetteur par des moyens optiques par rapport au repère de référence avion est remplacée par une mesure magnétique (mesures permettant de déterminer la position et l'orientation de l'émetteur). Une caractéristique de l'invention est notamment d'associer un modèle polynominal exprimant globalement et de façon approchée la position/orientation du système mobile en fonction de mesures par exemple du type du brevet français n 90.08230 et un modèle des champs non linéaire par rapport à la position/orientation du capteur. Le premier modèle permet une résolution approchée non itérative pour donner une position/orientation approchée. Le second modèle utilise les valeurs approchées de la résolution précitée pour identifier très précisément, ou du moins plus précisément la position/orientation en un temps très court par une méthode itérative d'un critère non linéaire. D'autres caractéristiques et avantages de l'invention apparaîtront à l'aide de la description qui suit, faite en regard de dessins annexés qui représentent: - La figure 1, un premier système de cartographie pour la mise en oeuvre du procédé selon l'invention, - La figure 2, un deuxième système de cartographie utilisé pour plusieurs opérations de modélisation ou de caractérisation de divers sous- ensembles, - La figure 3, une illustration d'une mise en oeuvre possible de la première étape du procédé selon l'invention, - La figure 4, un casque de pilote équipé d'un capteur de champ magnétique et d'éléments produisant des perturbations électromagnétiques, - La figure 5, une modélisation possible d'une source de perturbation électromagnétique précitee, - La figure 6, une représentation d'une perturbation électromagnétique causée par la source précitée, - La figure 7, une illustration d'une mise en oeuvre possible de la deuxième étape du procédé selon l'invention. - La figure 8, une illustration d'une sous-étape possible du procédé selon l'invention. La modélisation du champ électromagnétique dans un cockpit est fondée sur les hypothèses suivantes: Le champ global est la somme de trois champs le champ Bó produit par les bobines d'émission; - le champ BI, induit par Bo qui provoque des courants de Foucault dans les conducteurs fixes du cockpit; - le champ B2, induit par Bo + B, qui provoque des courants de Foucault dans les conducteurs mobiles dans le champ et liés à l'équipement de tête (casque). En général I B2 I La somme Bo + B, est le champ existant dans le cockpit, en l'absence de conducteurs dans le casque. Les champs sont exprimés par des fonctions analytiques, dont les paramètres sont les coordonnées du capteur; ces fonctions sont solutions des équations des champs selon les équations de Maxwell connues. Les solutions des équations des champs ne dépendent théoriquement que des conditions aux limites du volume, c'est-à-dire à sa surface. Par la suite, on recherchera une modélisation séparée des champs (Bo + BI) et B2. Pour la modélisation de (Bo + BI), on utilise la connaissance du modèle de Bo (espace libre), pour en déduire le champ B1. La détermination de l'orientation du capteur fixé sur le casque s'opère en effectuant une cartographie magnétique à partir de laquelle un modèle du champ magnétique est déterminé pour le volume de débattement du capteur. Sur la figure 1, est présenté le synoptique général d'un premier système de cartographie pour la mise en oeuvre du procédé selon l'invention comportant: une antenne d'émission pour l'émetteur 1 liée à un repère RE de centre OE; 30. une antenne de réception pour le capteur 2 liée a un repère Rc de centre Oc; a un équipement de tête 5 (casque équipé) composé d'une coque sur laquelle sont fixés notamment le capteur 2, son préamplificateur 3, le (ou les) tube (s) CRT 6, le module optique 19. Le module optique est constitué d'un système de génération d'un symbole et d'un système optique de collimation à l'infini, et enfin d'un système semi réfléchissant permettent de superposer ce symbole au paysage extérieur, définissant la ligne de visée Lv. L'équipement de tête comporte en outre une mémoire 7 associée au capteur 2, et une mémoire 8 associée au casque. Ces mémoires sont par exemple des PROM, EPROM ou EEPROM. une unité centrale 10 de l'équipement; un circuit d'interface 30 et un calculateur de dialogue 60 connectés pour la gestion des mémoires associées aux-sous-ensembles interchangeables, capteur, casque par exemple, ou pour le téléchargement de programmes. Un module de connexion 12 relie l'équipement de tête à l'unité centrale 10. Un bus d'alimentation est relié aux circuits. Une mémoire 9, placée par exemple sur un amplificateur 11 câblé entre l'antenne d'émission 1 et l'unité centrale 10, est associée à l'antenne d'émission 1. Une mémoire 14 associée au porteur, c'est-à-dire par exemple au cockpit ou cabine de l'aéronef est reliée par exemple à l'unité centrale 10. Ces mémoires sont par exemple des PROM, EPROM ou EEPROM. La mémoire 9 associée à l'antenne d'émission 1 contient des coefficients définissant un modèle du champ Bo, précité. La mémoire 14 associée au porteur contient des coefficients définissant un modèle du champ Bl précité. La mémoire 8 associée au casque contient des coefficients définissant un modèle du champ B2 précité. Les mémoires associées aux éléments ou sousensembles 7, 8, 9, 14 chargées de leurs coefficients font, par exemple, parties intégrantes de ces éléments ou sous-ensembles. Elles sont par exemple chargées en usine avant d'être associées à leurs éléments respectifs. Chaque élément ou sous-ensemble 7, 8, 9, 14 est monté ou changé avec sa mémoire associée, le sous ensemble 14 est associé avec le cockpit. Sur la figure 2, est présenté pour la mise en oeuvre du procédé selon l'invention le synoptique d'un second système de cartographie, pour la mise en oeuvre du procédé selon l'invention, utilisé pour plusieurs des opérations de modélisation ou caractérisation des divers sous-ensembles ou cockpit. En effet, il peut être utilisé de façon indifférente: - pour caractériser les antennes en usine, en utilisant un modèle espace libre B0, c'est à dire sans perturbations le capteur étant alors monté sur un support spécifique, ou une coque de casque permettant de positionner le capteur conformément à l'utilisation réelle sur le casque; - pour caractériser les casques équipés en utilisant un modèle B2 défini par la suite; - pour caractériser les perturbations du porteur, la configuration étant identique à la première mentionnée ci-dessus utilisée en usine pour caractériser les antennes. L'ensemble est harmonisé dans le repère de l'avion de façon connue. Une cartographie magnétique est réalisée en translations et aussi par rotations si nécessaire, du capteur dans le volume de la cabine de l'aéronef. Selon l'invention, toutes les matrices de mesures sont mémorisées. A titre indicatif, mille matrices de 9 termes suffisent pour obtenir un bon modèle de l'environnement du cockpit. Le capteur 2 est monté sur un support mécanique 21 qui permet de le faire tourner selon trois rotations orthogonales R, S, G. Le principe de l'invention consiste en une première étape de modélisation initiale comprenant: - une modélisation analytique des champs électromagnétiques en fonction des coordonnées du capteur, et satisfaisant les lois de l'électromagnétisme; - une modélisation séparée des effets, en particulier ceux * liés au matériel: en créant un modèle B0 des antennes d'émission et réception en espace libre avec leurs corrections; et un modèle B2 des effets perturbateurs du casque et équipement de tête * liés à l'environnement dans lequel est utilisé le matériel: en créant un modèle B1 des effets perturbateurs dus aux pièces conductrices de la cabine. -une caractérisation de chaque source de perturbation indépendamment l'une de l'autre, par les coefficients de son propre modèle, coefficients qui sont ensuite mémorisés dans des mémoires, par exemple PROM, EPROM, ou EEPROM, associées aux éléments modélisés considérés comme "interchangeables", ces éléments sont par exemple l'émetteur, le capteur, le casque, le porteur (cabine ou cockpit par exemple). Dans une deuxième étape le procédé selon l'invention consiste en un calcul, en temps réel, de la position et l'orientation notamment du capteur ainsi que du casque et de la ligne de visée, en utilisant la mesure courante et le modèle mémorisé. Cette seconde étape fait appel, par exemple, à: -un algorithme non itératif permettant de fournir une valeur approchée de la position et l'orientation du capteur en un temps très court; -un algorithme itératif permettant d'optimiser le calcul de la position et l'orientation du capteur qui interviennent dans le modèle analytique dont les paramètres sont mémorisés dans les mémoires associées 7, 8, 11, 14. L'utilisation du modèle du casque B2 doit notamment satisfaire aux conditions suivantes: - les modèles d'antenne Bo et du cockpit BI étant déterminés, le modèle de B2 s'applique indifféremment à tous les champs (Bo + BI) pourvu qu'ils soient connus même, les coefficients de B2 étant déterminés sans connaître B1 à l'origine. Notamment, les coefficients de B2, s'identifient indépendamment de BI c'est à dire en usine par exemple, Bo étant généré et connu. Par contre, la compensation de B2 dans le cockpit, qui dépend de (Bo + B1) est valable pour tout (Bo + Bi) pourvu que Bo et BI soient connus et utilisés dans le modèle de copensation de B2; - L'identification ou caractérisation des antennes et du casque, c'est-à-dire l'ensemble de l'équipement est faite en usine; celle du cockpit est faite indépendamment sur le site avec un outillage de cartographie sans casque avec des antennes préalablement caractérisées. Tout cela est rendu réalisable par la complète interchangeabilité des sous-ensembles (cartes, amplificateurs, antennes après corrections). Un outillage de cartographie est par exemple un capteur étalonné sur son support mécanique. Ce support peut être harmonisé dans le repère avion de référence, par des mesures magnétiques, par exemple, ou bien son orientation peut être déterminée optiquement dans le repère avion. Avant de décrire plus en détail les étapes du procédé selon l'invention, les notations suivantes doivent être prises en considération: (j)kX X, Y Z étant les coordonnées Il l v = avec X =Y du capteur estimées Rot Z G G,S,R étant les angles définissant Rot = S la rotation du capteur R Le vecteur v est l'estimée du vecteur v. Les modèles B0, B1, B2 sont définis par: Bo = Bo Ao, X) B1 = B1 AI,X) B2 = B2 A2, X, RotI (Ao) = est le vecteur de paramètres caractérisant les antennes: capteur i et émetteur j. On peut donc définir Ao par Ao0 J ij pour relier les jeux de paramètres aux antennes correspon- dantes. (A I) = est le vecteur de paramètres du champ Bl de la cabine t numéro e. A2 = est le vecteur de paramètres du champ B2 relatif au casque k numéro k. AMK1 = est le vecteur dont les composantes sont les paramètres du modèle (Bo + Bl) identifiés par la méthode décrite dans le brevet 90.08230 pour une configuration standard. La figure 3 illustre la première étape du procédé selon l'invention. Cette première étape, de modélisation, comprend une première sous-étape 31 de détermination des coefficients du modèle des effets du capteur n 1 et de l'émetteur n j en espace libre. L'obtention du modèle en espace libre Bo est obtenue de la manière suivante: * A chaque émission, on dispose des trois composantes du champ émis, soient les trois vecteurs: - - -_. Ml, M2, M3 que l'on regroupe en une matrice de mesure [M]=[MI, M2, M3] * L'objet du modèle B0 en espace libre est d'identifier l'ensemble des paramètres dont dépend ce modèle, à savoir: + les coordonnées angulaires RE/R et de position XE de l'émetteur dans un repère de référence; + les coordonnées angulaires du capteur dans un repère de référence RCO/R lorsque l'outil c'est à dire le support capteur est en position angulaire de référence. + les matrices de correction de facteur d'échelle et de couplage de l'émetteur [HR KR]et du capteur [HCii Ko] + le gain global émetteur capteur: kCR. * L'écriture du modèle de la mesure dans le repère du capteur est le suivant: P t M = kcR [HcKc] (Rco/1t) (RClCO)A (RE/1r)f p3 [H3] [HR KR] Cela revient à exprimer: BE P[H2]P' p3 expression connue du champ magnétique d'un dipôle dans le repère orthonormé de l'émetteur avec classiquement: -2 0 O [H3]= 0 -1 0 matrice de couplage magnétique dipolaire. 0 0 -1 _] [P]= UXVE.: U -E,UxEII -'-' '-4 --'9 UXE est le vecteur unitaire de XE = Xoc - XOE, c'est-à-dire du vecteur émetteur-capteur (Xoc, XOE: vecteurs dont les composantes sont respectivement les coordonnées du centre du capteur et celles du centre de l'émetteur) dans le repère émetteur. UXE' et UXE" sont les vecteurs unitaires perpendiculaires formant le repère orthonormé dit repère radial. La mesure corrigée s'écrit: Bl = [Hc Kc]-' M[HRKR] kcR dans laquelle: Hc et HR sont les matrices définies par les vecteurs unitaires des axes respectivement capteur et émetteur exprimés dans un repère --- orthonormé lié à des deux éléments (l'axe X de ces repères est colinéaire avec un axe pris comme référence, les deux autres axes perpendiculaires forment un repère orthonormé); kCR est un coefficient de normalisation des facteurs d'échelles émetteur et capteur globalement. On a donc la relation dans le repère de référence mécanique: B,, = (Rc01, jj)t (RClCO)"1 (RE1,À)kî BE = (RCIA) (RE/AI )AI BE BAI = (&/,f)Y (BE)AL BE etBAi sont les composantes du vecteur champ exprimé respectivement dans le repère de l'émetteur et le repère mécanique. L'identification du modèle espace libre dipolaire consiste à calculer à partir d'un certain nombre de mesures M et les rotations (Rc/cO)A , les paramètres suivant du modèle kCR, [HcKC], [HR KR], [RE/II], ['ZOE]A,, [RCOIM], faisant respectivement 1 terme, 5 termes, 5 termes, 3 termes, 3 termes, et 3 termes ce qui fait 20 termes au total. En utilisation en espace libre, il s'agit de trouver à partir de la mesure M, la position et l'orientation du capteur (de laquelle on déduit la ligne de visée) de la façon suivante: Calcul de: a) BM (connaissant kcR,[Hc KCJ,[HR,KR]) b) XE soit pet P dont on déduit BE = P[H3 P Connaissant RE/,I, on en déduit (BE) t champ BE exprimé dans le repère de référence (mécanique M). On en déduit [Pcl/ ] Af = [BE],Af [BB1] En présence d'émission non dipolaire, les termes de la matrice de couplage magnétique ne sont plus constants mais dépendent du vecteur XE selon les fonctions associées de Legendre connues de l'homme du métier. Les coefficients, qui sont les composantes du vecteur Ao sont mémorisés dans la mémoire 7 associée au capteur n i et dans la mémoire 11 associée à l'émetteur n j. La première sous-étape 31 est suivie d'autres sous-étapes 32, 33 qui utilisent par exemple les coefficients A o ij Une sous-étape suivante 32 consiste en la modélisation d'un casque n0 k équipé avec un capteur n i' par des mesures en Espace Libre. Les figures 4, 5 et 6 illustrent un exemple d'obtention du modèle B2 du casque n k. Sur la figure 4, un capteur est représenté par un repère orthonormé C d'origine Oc, le capteur est par exemple fixé sur le casque 41. Des sources de perturbation électromagnétiques 42, 43, 44 sont par ailleurs elles aussi fixées sur le casque 41. A titre d'exemple, elles sont limitées à trois. Elles sont représentées par des repères orthonormés Si,S2,S3,...Sk d'origines respectives 01, 02, 03,... 0*. et d'axes respectifs X1, YI, Z1, X2, Y2, Z2 X3, Y3,Z3 et Xk,Xk, Zk. Les directions de (Xl, YI, Zl),(X2,Y2,Z2),(X3, Y3,Z3),(Xk,Xk,Xk.) sont quelconques. Les origines 01,02, 03...0k des repères sont par exemple les centres des sources de perturbations 42, 43, 44. Par la suite, celles-ci seront confondues avec leurs repères. Le repère orthonormé C du capteur est mobile dans un repère de référence Réf d'origine OR. Les repères SI, S2, S3, et plus généralement les repères Sk lorsque k > 3, des sources de perturbation sont fixes par rapport au repère C du capteur. Un repère orthonormé E d'origine OE, fixe dans le repère de référence Réf est par exemple lié à un émetteur d'induction magnétique, servant notamment à l'aide du capteur fixé sur le casque 1 à déterminer la position et l'orientation de ce dernier comme expliqué précédemment. Dans une première phase, les sources de perturbation S, S2,S3,...Sk étant plongées dans un champ magnétique donné ambiant, par exemple par l'antenne d'émission [Bo], pour chaque source, le modèle de la perturbation d'une source au niveau du capteur C est exprimé comme une fonction explicite du champ moyen existant, par exemple à l'origine Ok du repère définissant cette source. Cela permet notamment de définir un modèle de la source dépendant explicitement du champ dans lequel est plongé chaque source Sk, et par conséquent de pouvoir déterminer les paramètres du modèle de la source en dehors de son environnement électromagnétique d'utilisation, ceux-ci devenant de cette façon indépendants du champ ambiant. Dans une deuxième phase, le capteur C étant plongé dans le champ magnétique à mesurer, la perturbation causée par chaque source Sk est déterminée à l'aide de son modèle et du champ magnétique moyen existant à son niveau. Dans une troisième phase, la somme des perturbations dues aux sources est effectuée. Enfin, dans une quatrième phase, cette somme est retranchée à la mesure réalisée par le capteur C. De cette manière, I'ensemble des paramètres du modèle de source représentatif du phénomène de perturbation produit par cette dernière, est indépendant du champ dans lequel est plongé le capteur et l'ensemble des sources, donc dans l'exemple d'application, cet ensemble de paramètres du modèle est indépendant du champ magnétique dans lequel est plongé le casque. C'est àdire que la détermination du modèle peut se faire notamment dans un champ différent du champ utilisé de manière opérationnelle par le pilote. La première phase précitée peut ainsi être par exemple réalisée en usine en dehors de l'environnement opérationnel, c'est à dire en dehors de la cabine d'aéronef par exemple. Dans le cas d'application au viseur de casque notamment, les modèles de source peuvent être identifiés en usine par des moyens spécifiques, tels qu'un émetteur quasi-dipolaire rayonnant en espace libre par exemple, et être utilisés dans une cabine d'aéronef quels que soient les champs électromagnétiques existant dans une bande de fréquence contenant la fréquence d'oscillation du champ produit par l'émetteur E fixé dans la cabine, pourvu que leur modèle soit connu par ailleurs et identifié dans une étape antécédente. La création d'un modèle représentatif de la perturbation d'une source peut être réalisée notamment en considérant de manière simplifiée qu'une source, et notamment une pièce conductrice, se comporte comme une ou plusieurs boucles élémentaires de courant. En espace libre, par exemple, le champ magnétique mesuré par le capteur C est la somme du champ émis par l'émetteur E et des champs produits par chaque source Sk. du casque. Ce champ magnétique mesuré peut s'exprimer par la relation suivante: NS Bc (Oc)= (RCIE) BE (OC)+ E R k (Bk (0k)) Sk (1)' k=l c/Sk fkO - Bc (Oc) étant le champ total mesuré exprimé dans le repère C du capteur et au centre Oc du capteur; - BE (Oc) étant le champ ambiant, produit par exemple par l'émetteur au centre Oc du capteur et exprimé dans le repère E de l'émetteur, ce champ BE (Oc) est le champ à mesurer; - (Rc/E)t étant la matrice transposée de la matrice de passage du repère E de l'émetteur au repère C du capteur; - Rc/E s'écrit (RcIRef)E (RE/Ref)tE = (RE/Ref)'Rtf (Rc/Ref)Rdf O Rc/Ref est la matrice de passage du repère de référence Ref au repère capteur (C), exprimée dans le repère E, elle s'écrit (Rc/Ref)E ou dans le repère Ref, elle s'écrit (RCRc f) Re fr - Rt/E BE (Oc) = R' c/Re f [BE (Oc) Ref [BE (Oc)] Ref est le champ ambiant exprimé dans le repère de Référence Ref. - Bk(Ok) est le champ magnétique moyen au centre Ok et exprimé par exemple dans le repère Sk; la fonction fk représente la relation existant entre le champ ambiant au point Ok et le champ perturbant au point 0c produit par la source k et exprimé dans le repère (Sk). - RcI/Sk est la matrice de passage du repère St. de la source Sk au repère C du capteur; - NS est le nombre de sources perturbatrices. La figure 5 illustre un modèle représentatif de la perturbation d'une source Sk.. Ce modèle est constitué d'une boucle ou spire élémentaire SPk de centre 0k., origine du repère orthonormé Sk de la source, et de rayon rk. Il est supposé que le champ magnétique moyen Bk(Ok) au centre de la spire SPk est sensiblement constant dans une sphère de rayon rk. Le champ moyen est donc supposé réparti de façon sensiblement constante dans la spire. Ce champ moyen est variable dans le temps et provoque dans la spire un courant ik proportionnel à sa variation selon la loi connue de Lentz. Ce courant provoque à son tour un champ magnétique selon les lois connues de l'électromagnétisme. La source peut aussi être une aimantation induite par le champ ambiant dans des matériaux perméables (effets ferromagnétiques). En considérant dans un premier temps une spire élémentaire de faibles dimensions, I'expression du courant ik circulant dans la spire élémentaire plongée dans le champ moyen Bk(Ok) est donnée par la relation suivante: Rk ik + Lk d ik = dk (2)' dt dt Rk et Lk étant respectivement la résistance ohmique et l'auto- inductance de la spire SPk, Ak étant le flux magnétique traversant la spire. Pour un champ magnétique Bk(Ok.) au niveau de la spire variant sinusoïdalement dans le temps à la pulsation co, il vient: ik [Rk + j Lk] = A- i' Ak (Bk(Ok). k) (3)' - Ak étant l'aire de la spire SPk; - /ik étant la normale à la spire; - Bk(Ok.). it. représentant le produit scalaire du vecteur champ magnétique ambiant au centre de la spire par la normale /7k à cette spire. De la relation (3), il vient: ik =+ Ak Bk (Ok). k (4)' Rk + j Lk cW En notant I le module du courant alternatif ik, il vient: I =kkI IaBk(Ok)."k =a IIBk (Ok). Cos( (Ok), k) (5)' avec a=- j Ak o (6)' Rk + j Lk co En considérant un repère radial p,uO,U- centré au point Oc, ok Oc centre du capteur, comme l'illustre la figure 5, avec,,= i, ioet Ok Oc étant deux autres vecteurs unitaires tels que up, uO, ui forment un repère (uk)orthonormé direct avec o dans le plan (un, ,), le vecteur champ (BPk)uk crée par la spire SPk. s'exprime dans le repère radial up, u0, u par trois composantes, une composante Bp selon l'axe ỉp, une composante Bo selon l'axe ỏ et une composante Bo selon l'axe iu. Ce champ (BPk),k modélise la perturbation de la source k. Les lois connues de l'électromagnétisme enseignent que pour [Ok 0óc > rk, ., les composantes Bp, Bo et Bo sont données par les relations suivantes Bp= I [Uao 2 zrk2 (-) (2e + 1)!! 2e2 P2e+(C k)(7)' 4gr P e=o 2()' 0 2e(e!)pe--5Pel (ó') B I|o 2 (-0e (2e(+1)!! 1 + l(rk)2(C O ' (8)' B,0, =-')e=o o= (9)' avec Po.o| I = lik | - I'k =(k,,),nk =Ok Zk sur la figure 5; - Pt"' (x) = polynôme de Legendre d'ordre e,m -(k)!! = k (k-2)(k-4)x... 5 x 3 x 1, k étant impair. Les expressions pour p modifient pas le modèle final et de plus dans la réalité p> rk. L'expression du champ (Bpk) de la source Skou de la spire SPk, exprimée dans le repère C du capteur, est donnée par la relation suivante: (BPk)c = (IRk/c) (BPk) 11k (10)' o Ru,/C est la matrice de passage du repère C du capteur, précédemment défini, au repère u21/k. = ip, io, ait précédemment défini. De façon plus générale, si l'on considère un système de trois spires d'axes de symétrie de révolution orthonormés x,,niy, iz, définissant par exemple le repère (Sk), le champ résultant au niveau du capteur Oc est la somme des trois champs produits par les spires ci-dessus. Dans le repère radial (tak), le champ dû à la source k est: (k (Oc)),. = + BPk + BPk] k (11)' B",2,3 obéissent aux équations (7)', (8)', (9)'. Pk Dans un premier temps, on supposera pour simplifier l'exposé que Pk = Oc Ok "> rXrVrz rayons des spires. Avec ces hypothèses, les formules (7)', (8)', (9)' donnant les B1'2'3 ci-dessus se simplifient et s'écrivent Pk de façon connue de l'homme de l'art par: (BPk(Oc)),k 2 p 0 ++P +stij (12) (Bp (j)uk=--- o - 0Pl o,,. [xrir,+SV Yi+ Szlzhz]i (1 2)' Pk. 0 0 -I avec P=[ [p,/i, ii0,]: matrice de passage du repère (Sk), au repère radial (11k)- Sx = rr2, Sy = rr., Sz = rr2. IX, Iy,]z sont les courants I des bobines d'axes six, ly, nz. Sx Ix, Sy Iv, Sz Iz sont respectivement les moments magnétiques d'axes zx nix, iz. Cette matrice P est constante pour la source k puisqu'elle ne dépend que de Oc Ok, fixe dans le repère (C) du capteur ou du repère (Sk) de la source. O10 Dans la relation (12)', les courants IxIv, Iz s'expriment selon les relations (4), (5), (6), soit: Ix =IaxBk(Ok)I COs(0x) Sk Iy =|ayl Bk(0k)ICos(0y) Sk (13)- z = laZ.l BPk(Ok)l Cos (Oz) Sk Comme l'illustre la figure 6, le vecteur B,. (Ok))= Vk est repéré dans le repère (Sk) par deux angles Ok et 0k, Ok repérant la projection du vecteur Vk représentant le champ magnétique dans le plan Ok., Xk, Yk et Ok l'angle du vecteur V7., avec Z. =iz. On prend les conventions suivantes: Cos ( O)Sk = Cos Ok Cos(Oy)sk =Sin 0k Sin qk (14)' Cos (0x)Sk = Sin 0k Cos Ok En substituant les relations (13)' et (14)' dans la relation (12)', l'expression du champ perturbateur (Bpk (Oc))c au point Oc, origine du repère C du capteur, exprimé dans ce repère C, est obtenue par la relation suivante en utilisant la relation (10)': -->o -Sisl 0k (ax Cos Ok + bx Si, k k)+ c Cos ok Bpk(OC) = IIBPkP(Ok)|| IRUkICI Si 6k (ay Cos çk + by Si" Ok)+ cy COS Ok,( L Si, Ok (a, Cos Ok + b: S ink)+ * Cos Ok j )' IBPk(Ok)il étant le module du vecteur Bpk(Ok). Bk étant le champ ambiant au point Ok eta, a,, y, baz, b, bz, ex, cv, c, représentent les coefficients de la matrice 2 0 0 laxisx 0 ax bx cx j 0 -I 0 [P1 0 o - YSy ay by cy o o -1 o o |:J.:z %z bz: Si l'hypothèse Pk >>{rX,rv,,r} n'est pas respectée, la matrice -2 0 0 0 -1 0 dépend de Pk et de 0'k et ó'k définissant les coordonnées de 0 0 -1 oc Ok selon les formules (7)', (8)', (9)' que l'homme du métier saurait appliquer. Quoi qu'il en soit, cette matrice reste constante pour la source Sk. Ainsi l'hypothèse simplificatrice ne change rien à la formule (15)' du champ perturbateur. L'expression du champ BpPk (0ó) dans le repère capteur est du même type que celui exprimé dans le repère (11k), seuls les coefficients ax, ay, az, bx, by, bz, cx, Cy, cz seraient modifiés dans le modèle. Les coefficients sont identifiés dans le repère capteur (C) par exemple. La formule (15)' montre en fait que chaque composante du modèle de la perturbation, le vecteur Bpk, s'écrit comme le développement limité à l'ordre 1 des parties réelles et imaginaires de la fonction complexe de Legendre associée définie par: L e 2e+1 (e-m)!, m Zj'( 4, O) (= m) P (Cos)e jm (16)' e=l m=o 4 (e m)!e avece Jm = Cosmo+jSinm M et avec Pm?(x=(-) (l 2)m2 de+m (x e (17)' e (X) 2e e! adxt+m (x2) A I'ordre 1,L= 1 etm=oetm= 1. De façon plus générale, les calculs et expériences réalisés par la Déposante ont montré qu'il est possible d'adopter pour chaque composante du vecteur champ Bpk, modèle de la source Sk, exprimée dans le repère C du capteur, la définition suivante: I k e k (BPk ())1 =V | k Pe (CoS Ok)A[ i Cos(m k) + C k i Sin (mók)] (18) (Bpk(c))i=Pky Z Pl'(Osk[rCsó) e,m, e=l t=o Dans cette formule, Ok est l'angle de Vk avec l'axe 0z du capteur, et,k. repère dans le plan 0c Xc Zc de ce même capteur la projection de Vk par rapport à Oc selon la figure 6. Lk. est l'ordre du développement du polynôme de Legendre relatif à la source k. Cette nouvelle définition ne change rien aux hypothèses précédentes puisque le repère (Sk;) a été supposé quelconque. On peut dans ce cas le considérer confondu avec le repère capteur (C). - (Bpk, (OC)), étant la ième composante du vecteur (Bpk(Oc)), i= 1, 2,3, cette composante est produite par trois spires de modélisation de la source Sk.; _ k k - Ami et Ce,,i étant des constantes à identifier pour la source k. Ce modèle (18)' tient compte du fait que le modèle par 3 spires filiformes (15)' est une vue simplifiée des phénomènes physiques réels. En fait, il faut parler plutôt de pseudo- spires qui correspondent aux lignes de courants circulant dans les éléments conducteurs. On vient d'établir le modèle (18)' de la perturbation et on constate qu'il dépend notamment, et de façon explicite, de Bk (0k), c'est à dire du champ extérieur au point Ok., centre de la source k. Ce champ est supposé connu: - lors de l'identification du modèle, ce champ est le champ BE émis en espace libre par l'émetteur; - en utilisation, dans un champ ambiant quelconque, par exemple dans une cabine d'aéronef, c'est le champ émis par l'émetteur BE auquel s'ajoute éventuellement le champ perturbateur de la cabine. Dans la première étape de l'invention qui consiste à déterminer les coefficients inconnus du modèle (18), le champ ambiant BE émis par l'émetteur est connu selon par exemple une méthode enseignée dans le brevet français n* 90.08230: la mesure du champ B, (Oc) résultant de l'émission de chacune des trois bobines d'émission permet selon cette méthode de fournir une valeur très proche de la position du centre du capteur (Xc, YZc) dans le repère (E) de l'émetteur ou le repère de référence Ref lorsqu'on fait subir au casque des rotations et des translations. Connaissant une estimation de Xcc, Zc, on en déduit la valeur du champ ambiant BE au centre du capteur. On exprime ensuite le champ Vk =Bk (OÈ,) en fonction des coordonnées de O. dans le repère de référence: OR Ok = OR Oc + Oc Ok o oc/Ok)= &/Ref oc Ok - = ROT. (P) c (o)c est constant dans le repère capteur ou dans le repère (Sk). 5. ROT est la rotation du casque égale à Rc/Ref matrice de passage du repère de référence Ref au repère capteur (C). Dans la phase d'identification du modèle, ROT est connue, on en déduit les coordonnées de la source k dans le repère de référence (Ref) en fonction des coordonnées X'k de la source dans le repère capteur, inconnues qu'il s'agit d'identifier au cours de cette première étape de l'invention. Les coordonnées de la source k dans le repère de référence s'expriment par Xk = X+ ROT.C) ou encore Xl. = Xc + ROT. (?'k) si X'k = (P)C En phase d'identification, ROT est la rotation connue du casque, (Xc), coordonnées du centre du capteur dans le repère de référence, sont connues par exemple de façon approchée par le procédé du brevet français 90.08230. On veut identifier 3 composantes de X'k = (P)c On en déduit (Bk(Ok))c = ROT' '[BE(Xk)]Re f que l'on substitue dans la formule (18)'. Après plusieurs itérations sur X'k et sur les Ak et Ck mni elú,m,i l'estimation de Xk devient meilleure et les valeurs des paramètres du modèle c'est à dire les composantes (XC)1,2,3 du vecteur Xc et les k k (Ak; C;,,,;Ce,,) convergent vers leur vraie valeur. La convergence est acquise lorsque le champ total mesuré par le capteur et son modèle (1)' sont sensiblement égaux en tout point de l'espace, c'est à dire lorsque (X'k)i, Ak etCk minimisent l'écart quadratique moyen entre mesure et champ modélisé selon les relations (1)' et (18)'. Le modèle doit être réidentifié pour chaque champ BE émis par chaque bobine de l'émetteur, les bobines d'émission étant au nombre de trois par exemple. Ces champs sont émis séquentiellement ou non. Ils sont par exemple linéairement indépendants, c'est le cas notamment si les bobines d'émission sont quasi-orthogonales. Dans le cas d'un émetteur à trois bobines, il y a donc 9 termes du type de la relation (18)' à définir, trois composantes d'espace du vecteur Bk étant définies pour les trois bobines d'émission. Cela est à réaliser pour toutes les sources de perturbation, k variant de là NS,NS étant le nombre total de ces sources. Ainsi, pour l'identification du modèle de chaque source Sk, il s'agit d'identifier en fait les coefficients A ki et Ck de la relation (18)', pour les trois bobines d'émission ainsi que les coordonnées (X'k). du centre Ok de k k la source SA. Les coefficients AkeX et Cki et (X'k. sont par exemple des paramètr úmjes rmepré k isentpatif e xml Ces, paramètres représentatifs du casque et de la source Sk pour un casque donné. Il est possible de placer le support liant rigidement le capteur C et la source de perturbation Sk, c'est à dire le casque, dans le champ des bobines d'émission et d'effectuer des translations et des rotations du casque. Connaissant au préalable le champ, ou son modèle, émis par les bobines, il est possible par exemple de mesurer par différence la perturbation et d'identifier par une méthode classique itérative les paramètres (X'k)i, A k" et Ck qui minimisent l'écart quadratique moyen entre la mesure et le champ préalablement connu ou son modèle. Cette opération, peut par exemple être réalisée en usine, de façon définitive. Les paramètres (X'k), Ak,,j et Ceki sont par exemple mémorisés dans un [' C,ml,i e., circuit 8 résident sur le casque, ce circuit pouvant être une mémoire de type PROM, EPROM, EEPROM notamment. L'identification de ces paramètres Ck A,,,,, et Ce,,,i et (X'k) et leur mémorisation dans les casques permet de rendre ceux-ci facilement interchangeables. Une autre sous-étape suivante 33 consiste en la modélisation d'une cabine n l, avec des mesures dans cette cabine équipée d'un outillage de cartographie, d'un capteur n i", et d'un émetteur n j". Avant de décrire, par un exemple de réalisation, cette sous-étape 33, on peut rappeler que l'on cherche à connaître avec précision le champ perturbateur induit par l'ensemble des conducteurs de la cabine sous l'effet d'un champ inducteur appelé B0 (x y z) connu dans le repère de l'émetteur. Le volume balayé par le capteur est intérieur à la cabine et ne contient pas de sources de courants. Pour identifier le modèle du champ perturbateur induit ci-dessus, appelé Bl (x,y, z), on connaît un certain nombre de mesures du champ total [Bo (x y z) + B1 (x y z)] mesuré en des points (x,,,my,,, z,,,) du volume, m = 1, M, correspondant au centre du capteur d'orientation noté par sa matrice de rotation RPc/A par rapport à un repère mécanique de référence. Ces points de mesure s'obtiennent par des translations connues dans un repère mécanique d'un outillage portant le capteur et par des rotations de ce même outillage connues dans le repère de référence (RC, ")= (Rcico) ij (Rco/if),X\ (C i) La notation CO correspond à l'orientation du capteur lorsque la rotation de l'outil est nulle (C'est à dire Rco/o=II (I étant la matrice identité). A partir de ces mesures, on identifie les paramètres inconnus du modèle connaissant à une translation et une rotation près les coordonnées et les orientations du capteur en chaque point de mesure. Puisque le problème est du type magnétostatique, les deux équations de base sont les suivantes: V.B=O0 V x H = J (C2) ou encore div B = 0 rot H = J Puisque dans le volume d'intérêt, c'est à dire le volume débattu par le capteur, la densité de courant J(xyz) est nulle, les équations se réduisent à V x H = O avec la condition V. B = 0. Puisque V x H = O, H dérive d'un potentiel scalaire 01 (x y z): H = - V i. Il vient également la relation B = po H. La relation sur le potentiel s'écrit: V uo V 01 = o > v2 01 = o (C3) Cela est l'équation de Laplace bien connue. Or, les solutions de l'équation de Laplace sont aussi bien connues, ce sont des fonctions harmoniques - pour un volume parallèlépipédique, on choisit des coordonnées rectangulaires et les fonctions sont du type: 9)(xyz)= Z An,, Sin (a x+Ox Sin (p. yv+ o) Sinh(yr.mz+0.m,.) n, im=l Les paramètres An,,, Ox.n, Oyn, Oznn sont identifiés grâce à la connaissance du champ mesuré en des points (x y z)capteur pour lesquels on a: [B(x yZ)capteurlnesure =-,P. V (x.vz) (C5) O et sachant que Jr m rL, m (C6) a -, p.1 - /;,r.,, - V+ b (c6) a b 'b,1= - pour un volume sphérique pour lequel on choisirait un système de coordonnées sphérique (p, 0, ç), la solution générale est fondée sur la série de fonctions associées de Legendre: O(, o. ço) = [Aep +Be, p-e+, p- Yl),M(O, ç) e=o 0II=--e o les fonctions harmoniques sphériques Ye,,, (S, 0) sont définies par x =,opcosOcos ço;y=,opcosOsin ço;z= + p sin 0(C7) Ye, ",I ( O ) = e2- n)! (cos) eim (C8) et o Pé" (x) (-l)"' ( -x2)2 d"P e(x) (9) d x"' * avec Pe (x): Polynôme de Legendre d'ordre e. Souvent e, m sont des entiers positifs. Mais, on peut très bien, sans changer la méthode, utiliser des solutions plus générales: fonction de Legendre de première espèce et d'ordre y non entier déterminé par les conditions aux limites. La fonction Ye, devient YY, v o y et v sont fonctionnaires. L'identification des coefficients se fait comme décrit par la suite. s De la même façon que précédemment, les coefficients Ae, m, Be, m qui sont des nombres complexes sont déterminés par la minimisation d'un critère de distance (moindres carrés) du type: g = MIN de 1 [ |V[E A e." loi + Be",? Pi.]] Ye" (0,P)i] iM = N =O,i=-eú | B(P 0S A)i]2 (CIO) -[B(, c, capteur mnesuré) avec M: nombre de mesures L: ordre du polynôme de Legendre. v = -1 ó ip +e n 1 e (cil) Vó=el 0+ e2 --±e3 sn ----lI ep p 90 p sin 0 9ço "', e2, e3 sont les vecteurs unitaires de OR Oc et des 2 vecteurs formant un trièdre direct. L'ordre du modèle est déterminé par tâtonnement jusqu'à ce que l'erreur s2 soit suffisamment faible. En fait, I'ordre dépend de la valeur spécifiée de l'erreur de la ligne de visée ST qui doit être inférieure à une certaine valeur dans 95 % des cas. L'erreur de la ligne de visée est calculée à partir de l'erreur de rotation du capteur comme cela sera vu par la suite. Deux possibilités existent pour identifier les coefficients du modèle dans la cabine: - la première consiste à identifier par le modèle (C10) I'ensemble des champs B(x.vz)=Bo (xyz)+Bl(xyz). Cette méthode donne des résultats corrects - la seconde qui sera précisée ci-après consiste à identifier le seul champ BI (x y z). On opère en deux phases: a Le champ Bo (x yZ) est identifié comme en espace libre en utilisant la mesure du capteur contenant Bo +BI. En fait, on identifie les paramètres de Bo et Bl. Cette solution fonctionne dans des zones plutôt proches de l'émetteur o Bl Les paramètres RCO/.f, RE/Ig XOE/M donnant respectivement l'orientation du capteur dans le repère mécanique et la position et I'orientation de l'émetteur dans le même repère sont correctement identifiés dans les cockpits d'avions ou d'hélicoptères, cela bien que non évident, a été vérifié par les résultats obtenus par la Demanderesse. a De même de la mesure corrigée [B,,,], on retire le modèle Bo précédemment identifié pour obtenir l'estimée de B,: [1 XYZ] = [B,,,] - [ (X Y)] (CI2) Les coefficients A.,,,, Be,,, du modèle sont estimés en remplaçant [B(PO, )i]capteur) de (CI 0) par (B1 (x, y, z)) estimé de (C12) x, y, z. Cette solution est préférée à la première parce que les positions et orientations de l'émetteur sont ainsi connues dans le repère du porteur. Il est ainsi possible d'en tenir compte lors de changements d'émetteurs puisque ces paramètres sont mémorisés dans les mémoires 9 liés à l'émetteur. La première étape du procédé selon l'invention comprend une dernière sous-étape 34 de génération d'un modèle global B, somme des deux modèles Bo, B1 identifiés indépendamment et à partir duquel des coefficients globaux du type de cout identifiés selon le procédé du brevet français 90 08230 sont obtenus, c'est à dire un modèle polynominal linéaire par rapport aux variables d'entrée représentatives des coefficients du champ global B, obtenu par une phase d'identification 341, selon par exemple le procédé du brevet français 90 08230. Les coefficients sont placés dans une mémoire 35 appelée mémoire équipement. La figure 7 illustre la deuxième étape du procédé selon l'invention, étape de la position et de l'orientation du capteur en fonctionnement opérationnel effectués dans une première phase ou sous-étape 71. A partir de la mesure des 3 champs émis par les 3 axes du capteur, c'est à dire des 9 mesures effectuées dans une première sous-étape 71, à l'instant courant m, on doit déterminer la position et l'orientation du capteur et en déduire l'orientation de la ligne de visée dans le repère de l'aéronef. Pour cela, on dispose, comme le montre la figure 7, outre les 9 mesures de [M],, les paramètres des modèles: - espace libre: (Ào)ij relatifs au capteur i et à l'émetteur j mémorisés dans la mémoire capteur 7 n i et dans la mémoire émetteur 9 n j - cabine: (j,)e relatifs à la cabine e, mémorisés dans la mémoire cabine 14 n e. - casque: (À2),k relatifs au casque k, mémorisés dans la mémoire casque n k. Ces paramètres ont été obtenus dans la première étape illustrée par la figure 3. On dispose par ailleurs, grâce à la mémoire équipement 35 du modèle du champ global approché de (Bo + B1). Le problème revient à trouver le vecteur v, défini précédemment, dont les composantes sont la position et l'orientation du capteur et qui permettent de trouver une distance enclidienne minimum entre modèle fonction de v et mesure courante corrigée B,,,. La mesure courante corrigée B,,, est obtenue dans une deuxième sous- étape 72 de correction des antennes du capteur et de l'émetteur à partir des coefficients stockés dans leurs mémoires associées 7, 9 et des mesures de la première sous-étape 71. Le principe consiste à ensuite, dans une troisième sous-étape 73, déterminer le vecteur v optimal par une méthode d'optimisation non linéaire. Il s'agit de déterminer les composantes de V qui minimisent un critère C de distance euclidienne entre la mesure corrigée BM et le modèle de cette mesure C = [Bt [Ro] [Bo (c)] + B2 (Rot, Xx) i--! 9 étant le nombre de composantes. La principale contrainte à cette méthode itérative est le temps de convergence. L'optimum doit être trouvé dans un temps fixé à l'avance garantissant le temps réel, par exemple à 100 Hz, on doit fournir l'information v de Position/Orientation. Pour garantir un temps faible de convergence il faut aussi s'assurer qu'à chaque nouveau cycle de calcul, toutes les 10 ms par exemple, la valeur initiale de v soit très proche de l'optimum. Cela suppose qu'aucune perte de la position/orientation ne puisse arriver suite à des coupures dans les mesures. Pour garantir d'avoir à chaque cycle la bonne position/orientation initiale, on utilise par exemple la méthode dite MK1 du brevet 90 08230 pour initialiser v (position/orientation) ou tout autre modèle polynominal donnant la position et l'orientation du capteur par une combinaison linéaire des valeurs définies par la position captée. On a montré dans la première étape comment les coefficients AMK1 du modèle MK1 ont été identifiés. Ces coefficients sont par exemple stockés dans la mémoire équipement 35. La figure 8 illustre la troisième sous-étape 73 en montrant les différentes phases de la détermination de vm à partir de [BM],,,. La sous-étape 73 se décompose en une phase d'initialisation et une phase d'optimisation non linéaire 74 permettant de déterminer en un minimum de temps la position et l'orientation du capteur. Une première approximation v1 est par exemple fournie à partir d'un modèle MK1 représentatif d'un ensemble moyen ou standard (antennes et cabines), sachant que cette méthode de détermination de v n'est pas itérative en vertu de la linéarité du modèle. Un autre modèle de type polynominal pourrait être utilisé à base d'éléments finis par exemple. Cette première identification permet de calculer un modèle approché du champ B2:[B2], que l'on retranche au champ mesuré corrigé [BM],,,. Le modèle du champ B2 calculé est par exemple celui défini par les coefficients obtenus dans la sous-étape 32 de la première étape de modélisation relative au modèle du casque. La méthode itérative qui permet de déterminer V, à partir de V"n fondée sur la minimisation d'un critère C peut répondre par exemple à l'aquation suivante: Min C = AMin v [Bç] (J) - ROT-' [BR (X) + Bi (XC] + B2 (ROT, Xc)} V V M.tYI IWI > j=l ROT-I= ROT', 9 étant le nombre de composantes mesurées à chaque instant m, j est l'indice de chaque terme des matrices B,,Bo,BI,B2 (j=1à9) A partir de l'orientation du capteur ROT, on détermine l'orientation de la ligne de visée LV qui est un vecteur fixe lié au capteur. Le modèle B2 peut recouvrir également, outre les effets dus à des courants de Foucault, des effets purement magnétiques de type fer doux, relatifs par exemple à une aimantation ferromagnétique. L'invention a été décrite pour une application à la détermination de la position et de l'orientation de la ligne de visée d'un casque dans une cabine d'aéronef. Elle peut cependant être appliquée à tous systèmes mobiles dans un environnement électromagnétique. Dans un cas d'application à un viseur de casque, les antennes et le casque sont modélisés indépendamment les uns des autres et de la cabine de pilotage, cette modélisation ayant lieu par exemple en usine. Tous les sous-ensembles de l'équipement sont interchangeables, de même que les cabines, tous ces éléments possédant leur propre modèle. REVENDICATIONS 1. Procédé de détermination d'une position et d'une orientation d'un système mobile (Lv) par rapport à un émetteur (1) de champ électromagnétique, le système mobile placé dans un porteur étant lié à un capteur (2) de champ magnétique, caractérisé en ce qu'il comporte: - une première étape de modélisation analytique des champs électromagnétiques en fonction des coordonnées du capteur (2), un premier champ (Bo) étant créé par l'émetteur (1), un deuxième champ (BI) étant créé par les courants électriques induits dans le porteur par le premier champ (Bo) et un troisième champ (B2) étant créé par les courants électriques ou des aimantations induits dans le système mobile par les deux premiers champs (Bo + BI), I'effet électromagnétique de chaque champ (Bo, BI, B2) étant caractérisé indépendamment des effets des autres champs par les coefficients de son propre modèle; - une deuxième étape de calcul, en temps réel, de la position et de l'orientation du capteur en utilisant la mesure courante du champ électromagnétique au niveau du capteur (2) et en utilisant les modèles (Bo, B1, B2) de la première étape, la position et l'orientation du capteur étant définie à partir du champ mesuré (B,,,) auquel le troisième champ (B2) a été retranché. 2. Procédé selon la revendication 1, caractérisé en ce que l'effet électromagnétique de chaque champ (Bo, B1, B2) est caractérisé séparément des effets des autres champs. 3. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il associe un modèle polynominal (B) exprimant globalement et de façon approchée la position/orientation du système et un modèle des champs (Bo, B1, B2) non linéaire par rapport à la position/orientation du capteur, le premier modèle permettant une résolution approchée non itérative et le second modèle utilisant les valeurs approchées de cette résolution pour identifier plus précisément la position/orientation en un temps très court par une méthode itérative de résolution d'un critère non linéaire. 4. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que le modèle des effets du premier champ (Bo) est stocké dans une mémoire associée à l'émetteur (1) et du capteur (2), le modèle des effets du deuxième champ (BI) dans une mémoire liée au porteur et le modèle des effets du troisième champ (B2) dans une mémoire associée au casque (2). 5. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que les coefficients d'un modèle global générique sont stockés dans une mémoire (35). 6. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la première étape comporte une première sous-étape (31) de détermination des coefficients du modèle des effets d'un capteur et d'un émetteur, en espace libre, sans perturbation. 7. Procédé selon la revendication 6, caractérisé en ce que le modèle en espace libre est obtenu de la façon suivante: * A chaque émission, on dispose des trois composantes du champ émis, soient les trois vecteurs: MI, M2, M3 que l'on regroupe en une matrice de mesure [M] = MI M2, M3]. * L'objet du modèle B0 en espace libre est d'identifier l'ensemble des paramètres dont dépend ce modèle, à savoir: + les coordonnées angulaires REIR et de position XE de l'émetteur dans un repère de référence; + les coordonnées angulaires du capteur dans un repère de référence RCOIR lorsque l'outil c'est à dire le support est en position angulaire de référence; + les matrices de correction de facteur d'échelle et de couplage de l'émetteur [HR KR] et du capteur [HC KC]; + le gain global émetteur capteur: KCR 8. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la première étape comporte une sous-étape (32) de modélisation en espace libre du système mobile équipé avec un capteur (2) et comprenant des matériaux perturbateurs du champ. 9. Procédé selon la revendication 8, caractérisé en ce que des sources de perturbation (S1,S2, S3, Sk) étant liées au capteur (C), la sous- étape (32) consiste dans une première phase, les sources (S1,S2, S3.... Sk) étant plongées dans un champ magnétique ambiant quelconque donné, à exprimer la perturbation causée par chaque source comme une fonction explicite du champ magnétique moyen existant au niveau de la source pour créer un modèle représentatif de sa perturbation au centre du capteur (C), dans une deuxième phase, le capteur étant plongé dans le champ magnétique à mesurer, à déterminer la perturbation causée par chaque source à l'aide de son modèle et du champ magnétique moyen existant à son niveau, dans une troisième phase à sommer les perturbations dues aux sources (SI,S2,S3,....Si.) et dans une quatrième phase à retrancher cette somme à la mesure effectuée par le capteur (C), chaque source (Sk) étant modélisée par au moins une spire (Spk), les variations du champ ambiant (Bk) au niveau de la spire créant un courant (i.) dans la spire, lequel produit un champ de perturbation. 10. Procédé selon la revendication 9, caractérisé en ce que la perturbation causée par une source (Sk.) est modélisée par un vecteur champ magnétique (Bpk. (Oc)) au niveau du capteur (C), ce champ étant exprimé de la façon suivante pour chaque composante: Lk e kC (Bpk( C)) = |Bk (0k)| 2 ptm (Cos Ok)[At,;Cos(m 0k)+Cti Sin(m ók)] m=o -(BPk (Oc))" étant la ième composante du vecteur (Bpk (Oc)) au niveau du capteur (C) est produite par au moins une spire de modélisation de la source (Sk), i = 1, 2, 3; - (BPk, (Oc)) étant la valeur du champ magnétique moyen ambiant au centre O,. de la source (S,. ) et (B,. (Oc)) son module; - Ok étant l'angle entre un axe Oc c du repère capteur (C) et le vecteur champ magnétique Bk (Ok) au centre Ok de la spire; - s,. étant l'angle entre un axe OCc du repère capteur et la projection du vecteur champ magnétique précité Bk(Ok) sur le plan passant par le centre Oc du repère capteur et contenant les axes Oc-- et Oc Yc, (Oc, Xc, Yc, Zc) étant le repère du capteur. - P,"(x) étant un polynôme de Legendre d'ordre e, m défini par 0 pm"(x) =(-1) (I-x) de+)e k! dXe+'" (2x -1) - Ak, et C,,, i étant des constantes. 11. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la première étape comporte une sous-étape (33) de modélisation du porteur, par des mesures dans ce dernier le porteur étant équipé d'un outillage de cartographie, d'un capteur et d'un émetteur. 12. Procédé selon la revendication 11, caractérisé en ce qu'à partir des mesures, on identifie les coefficients inconnus du modèle connaissant à une translation et une rotation près les coordonnées et les orientations du capteur en chaque point de mesure. 13. Procédé selon la revendication 11, caractérisé en ce que l'identification des coefficients (A,,,, Be,,,) du modèle se fait par la minimisation d'un critère de distance du type des moindres carrés: ú 2=MIN de V E [Atn,, + BR,,,-P(e+i)1 Ye,n, (Oep)J Atm, Btm i= p=o m=-t avec M: nombre de mesures L: ordre du polynôme de Legendre. V)= ey+ e2 - - - (ci1) ep p' +e psin 8 {op L'ordre du modèle est déterminé par tâtonnement jusqu'à ce que l'erreur s2 soit suffisamment faible. 14. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la première étape comprend une dernière sous-étape 34) de génération d'un modèle global (B), somme des deux modèles Bo, BI) identifiés indépendamment. 15. Procédé selon la revendication 14, caractérisé en ce que le modèle global (B) est un modèle polynominal. 16. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la deuxième étape comprend une première sous- étape (71) de mesures des champs émis, une deuxième sous-étape (72) de correction des antennes du capteur et de l'émetteur à partir des coefficients stockés dans leurs mémoires associées et des mesures de la prmeière sous-étape (71) et une troisième sous-étape (73) de détermination de la position/orientation du capteur. 17. Procédé selon la revendication 16, caractérisé en ce que la troisième sous-étape (73) se décompose en une phase d'initialisation et une phase d'optimisation non linéaire (74) permettant de déterminer la position/orientation du capteur. 18. Procédé selon la revendication 17, caractérisé en ce que la phase d'optimisation non linéaire (74) est itérative, une détermination V," de la position/orientation étant définie à partir d'une détermination V2,, par la minimilisation d'un critère C: vn C =i v E [B] - ROT Bi (,) + Bi () + B2= (ROT,,XC)} j=l ROT - Bo, B1, B2 étant les champs émis - B étant le champ global - 9 étant le nombre de composantes mesurées à chaque instant m. 19. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que le système mobile est la ligne de visée (Lv) d'un viseur de casque. 20. Procédé selon la revendication 19, caractérisé en ce que les antennes et le casque sont modélisés en usine indépendamment les uns des autres et de la cabine de pilotage.