FR 2732773 A1 19961011 FR 9504241 A 19950410 La présente invention concerne un procédé d'identification simultanée et de correction d'erreurs dans les mesures d'un magnétomètre monté à bord d'un véhicule, ainsi qu'un dispo- sitif pour la mise en oeuvre de ce procédé et divers systè-5 mes utilisant ledit dispositif. On sait qu'un magnétomètre monté à bord d'un véhicule, en particulier un aéronef, est généralement soumis à divers types de perturbations, par exemple des perturbations magnétiques et/ou des perturbations mécaniques, qui faussent10 les mesures de champ magnétique effectuées par ce magnéto- mètre. Afin d'obtenir des mesures exactes, il est donc nécessaire, le cas échéant, d'identifier les erreurs engendrées par ces perturbations et de les corriger. Différentes méthodes sont15 actuellement connues pour effectuer de telles corrections. Toutefois, pour l'essentiel, ces méthodes connues sont uniquement destinées à corriger des erreurs de mesure dues à des perturbations magnétiques et ne permettent pas de corriger des erreurs présentant d'autres causes et/ou dues à20 des perturbations non magnétiques. Ainsi, on connaît par exemple une première méthode usuelle utilisée pour corriger des erreurs de mesure dues à des perturbations magnétiques, qui nécessite pour sa mise en oeuvre une installation précise du magnétomètre sur le25 véhicule. Selon cette première méthode connue, on définit un modèle de perturbation magnétique de la forme Hm = [B].H + HB, dans lequel: - Hm représente le champ perturbé mesuré par le magnétomè- tre; - H représente la valeur effective du champ magnétique; - [B] est une matrice symétrique tenant compte des perturba- tions engendrées par les "fers doux", l'absence d'une partie antisymétrique étant due à l'installation précise et appropriée du magnétomètre sur le véhicule; et - HB représente un champ perturbateur dû aux perturbations engendrées par les "fers durs" On notera qu'en l'absence de perturbations magnétiques, la matrice [B] représentant alors la matrice identité et le champ perturbateur HB étant nul, le module du champ magnéti-10 que mesuré est constant et l'extrémité décrit une sphère. Sous l'effet de perturbations, cette sphère est transformée en un ellipsoïde décentré. Ladite première méthode de correction connue préconise d'effectuer les opérations suivantes:15 - effectuer un grand nombre de mesures du champ magnétique au cours de variations du cap et de l'assiette du véhi- cule; - déterminer l'ellipsoïde correspondant au mieux aux mesures effectuées; et - calculer à partir de l'ellipsoïde ainsi déterminée les composantes de la matrice [B] et du champ perturbateur HB. Le modèle obtenu alors permet de calculer le champ effec- tif H à partir de mesures ultérieures effectuées par le magnétomètre et soumises à des perturbations magnétiques. Cette méthode exige toutefois, comme indiqué précédemment, un agencement particulier et très précis du magnétomètre sur le véhicule, ce qui nécessite une fabrication et/ou un réglage spécifiques et coûteux. En outre, la précision obtenue par cette méthode est peu satisfaisante. Ces in-30 convénients sont partiellement résolus par une seconde méthode connue qui est soumise à des contraintes beaucoup moins importantes que la méthode précédente concernant l'agencement du magnétomètre sur le véhicule. Cette seconde méthode connue, qui permet également de corriger uniquement des erreurs de mesures dues à des perturbations magnétiques, est mise en oeuvre en deux étapes successives: - dans une première étape, on détermine, à partir d'un procédé semblable à celui de la méthode précédente, une mesure conservant la constance de la norme du champ, ladite mesure n'étant plus affectée que par un mésaligne- ment; et - dans une seconde étape, on détermine la composante non symétrique de la matrice de transformation semblable à la matrice [B] précédente, à partir des variations de la composante verticale du champ magnétique, cette étape étant mise en oeuvre en utilisant des mesures d'assiettes. Cette seconde méthode connue présente de nombreux inconvé- nients. En particulier, la mise en oeuvre des deux étapes précitées suppose une procédure de manoeuvres complexe et longue du véhicule, au cours de laquelle sont effectuées les mesures utilisées. Afin de remédier aux inconvénients des deux méthodes connues précitées, la demanderesse a décrit, dans la demande de brevet 94 15231 déposée le 19 décembre 1994, un procédé de correction d'erreurs de mesure d'un magnétomètre. Toutefois, ce procédé avantageux permet uniquement la correction des erreurs dues à des perturbations magnétiques, ainsi que des erreurs dues à un mauvais alignement du trièdre de référence du magnétomètre par rapport au trièdre de référence du véhicule. Par conséquent, il n'existe actuellement aucune méthode connue permettant de corriger des erreurs dues, au moins en30 partie, à d'autres types de perturbations que celles décri- tes précédemment. La présente invention a pour objet de pallier cette lacune de l'état de la technique. Elle concerne un procédé permet- tant de corriger, de façon générale et très précise, des erreurs dues à des perturbations déterminées et variées dans les mesures d'un magnétomètre monté à bord d'un véhicule, en particulier un aéronef. A cet effet, selon l'invention, ledit procédé est remarqua- ble en ce que: - on définit un modèle théorique correcteur des erreurs de mesure du magnétomètre de la forme: i=n Hc = [A].Hm + E ([C]i. [E]i) + Hp, i=l n étant un nombre entier supérieur ou égal à 1, dans lequel: Hc est le champ corrigé, 15. [A] est une matrice à déterminer, Hm est la valeur mesurée par le magnétomètre du champ magnétique dans un référentiel lié audit véhicule, les éléments [C]i sont des matrices à déterminer de dimensions respectives (3xmi), mi étant un nombre entier supérieur ou égal à 1, les éléments [E]i sont des matrices de dimensions respectives (mixl) dont les coefficients sont des mesures de paramètres perturbateurs ou des combinaisons de telles mesures, et 25. Hp est un champ perturbateur à déterminer; - on considère que le champ corrigé Hc est le champ terres- tre défini dans le référentiel du véhicule de sorte que: i=n [A].Hm + E ([C]i.[E]i) + Hp = [M].H i=l avec: H la valeur effective du champ magnétique dans un référentiel de base, et [M] une matrice de passage dudit référentiel de base audit référentiel du véhicule; - on effectue des mesures Hm successives du champ magnétique au moyen du magnétomètre lors d'évolutions du véhicule; - on détermine les coefficients des différentes matrices [E]i, à partir des mesures des paramètres perturbateurs, effectuées simultanément aux mesures Hm; - on détermine les coefficients de la matrice de passage [M]; - on définit un vecteur d'erreur E de composantes Ex, Ey et Ez, défini par la relation: i=n E = [M].H - ([A].Hm + E ([C]i.[E]i) + Hp) i=l - on définit une erreur globale E2 = Ex2 + Ey2 + Ez2; - on définit un système d'équations formé de m équations du type (DEE2)/3cl = 0, 1 = 1 à m, (âEE2)/cl correspondant à la dérivée partielle de la somme des erreurs globales EE2, pour toutes les mesures, par rapport à un coefficient cl, les différents coefficients cl représentant les coeffi- cients à déterminer dudit modèle théorique, à savoir les coefficients de la matrice [A], des matrices [C]i et du vecteur Hp, ainsi qu'au moins une composante du champ magnétique effectif H; - on résout ledit système d'équations de manière à obtenir lesdits coefficients cl, dont ladite composante du champ magnétique effectif H; et - on utilise le modèle correcteur obtenu à partir desdits coefficients pour corriger les erreurs dues à des pertur- bations déterminées dans les mesures dudit magnétomètre. Ainsi, grâce à l'invention, on obtient une méthode générale permettant d'identifier et de corriger de façon précise des erreurs de mesure d'un magnétomètre, dues à des perturba-10 tions variées et multiples, comme on le verra plus en détail ci-après à partir de différents modes de réalisation parti- culièrement avantageux. De plus, grâce à l'invention, la détermination du modèle correcteur est effectuée en une seule étape, ce qui entraîne15 une mise en oeuvre simplifiée, en particulier en ce qui concerne la procédure de manoeuvres du véhicule, au cours de laquelle sont effectuées les mesures. En outre, le procédé conforme à l'invention permet de déter- miner, simultanément au calcul dudit modèle correcteur, au20 moins une composante du champ magnétique effectif dans le référentiel de base. On notera que les calculs effectués lors de la mise en oeuvre du procédé sont simplifiés dans le cas o on ne détermine qu'une seule composante dudit champ magnétique. En effet, dans les modes de réalisation pour lesquels on ne détermine que la composante verticale du champ magnétique, on utilise avantageusement comme erreur globale E2 = Ez2, de sorte que seuls les coefficients de la troisième ligne de la -1. -1 matrice [M] sont utilisés, ladite matrice [M] étant la matrice de passage dudit référentiel du véhicule audit référentiel de base, lesdits coefficients étant déterminés à partir des assiettes longitudinale et transversale dudit véhicule. En revanche, dans les modes de réalisation pour lesquels on détermine les trois composantes du champ magnétique, les coefficients de la matrice [M] sont déterminés de façon avantageuse à partir des assiettes longitudinale et trans- versale, ainsi qu'à partir du cap dudit véhicule. Avantageusement, la valeur de cap utilisée dans ce cas peut correspondre soit à une valeur de cap géographique, soit à une valeur de cap affectée d'une erreur constante quel- conque, telle que celle délivrée par une centrale d'attitude et de cap en mode gyroscopique libre. Pour bien mettre en évidence la variété et la multiplicité de corrections possibles par le procédé conforme à l'inven- tion, on présente ci-dessous plusieurs modes de réalisation particulièrement avantageux, ces différents modes de réali- sation pouvant être réalisés simultanément lors d'une mise en oeuvre unique dudit procédé. A titre d'exemple, la présente invention permet de corriger des erreurs de mesure du magnétomètre dues à des perturba- tions mécaniques, en particulier des déformations de l'élé- ment porteur du magnétomètre qui sont susceptibles de modifier l'alignement de ce dernier, ce qui provoque alors25 des erreurs de mesure. En premier lieu, dans le cas o le magnétomètre est monté sur un élément porteur, du type poutre de queue d'un héli- coptère, susceptible d'être déformé lorsqu'il est soumis à des contraintes mécaniques variables, selon l'invention,30 pour identifier et corriger les erreurs dues à de telles déformations dans les mesures du magnétomètre, on introduit de façon avantageuse dans le modèle théorique le produit d'une matrice [C] de dimension (3,5) et d'une matrice [E] de dimension (5,1), telles que: 0 -Cez 0 0 Cey [C]= Cez 0 -Cex 0 0 0 Cex 0 -Cey 0 et Fv.Hmx Fv.Hmy [E] = Fv.Hmz Fh.Hmx Fh. Hmz avec: - Fv: l'effort latéral engendré sur l'élément porteur, - Fh: l'effort vertical engendré sur l'élément porteur, - Cex: le coefficient d'élasticité de torsion de l'élément porteur, - Cey: le coefficient d'élasticité de flexion en tangage de l'élément porteur, Cez: le coefficient d'élasticité de flexion en lacet de l'élément porteur, et - Hmx, Hmy et Hmz: les composantes respectives du champ mesuré Hm dans le référentiel du véhicule. En second lieu, dans le cas o le magnétomètre est monté sur une aile d'un avion, par exemple d'un avion de transport,20 selon l'invention, pour identifier et corriger les erreurs dues à un défaut d'alignement du magnétomètre provoqué par le battement de ladite aile, on introduit de façon avanta- geuse, dans le modèle théorique, le produit des éléments suivants: 0 0 Pz.Hmz [C]: -Cox 0 et [E] = 0 Cox Pz.Hmy dans lesquels: - Cox est le coefficient d'élasticité en battement de ladite aile, - Pz est la portance de ladite aile, qui est définie avanta- geusement par la relation Pz = M.gz, avec M la masse du véhicule et qz le facteur de charge vertical, et - Hmy et Hmz sont les valeurs des composantes latérale et verticale du champ magnétique mesuré Hm. Outre la prise en compte de perturbations mécaniques, telles que celles présentées ci-dessus, la présente invention permet également de remédier aux conséquences dommageables d'autres types de perturbations, en particulier des pertur- bations d'origine magnétique comme précisé ci-après. Ainsi, par exemple, dans le cas o ledit véhicule est un hélicoptère dont le rotor présente un résidu de magnétisa- tion, selon un mode de réalisation particulièrement avanta- geux de l'invention, pour identifier et corriger les erreurs dues à des variations des perturbations magnétiques provo- quées par ledit rotor aimanté soumis à des contraintes mécaniques élevées, on introduit de façon avantageuse dans le modèle théorique le produit des éléments suivants: Hcx [C] = Hcy et [E] = Cr, Hcz Cr étant le couple rotor et Hcx, Hcy et Hcz des coefficients de proportionnalité à déterminer. Par ailleurs, avantageusement, pour identifier et corriger les erreurs dues à des perturbations magnétiques engendrées par la présence d'au moins un équipement consommateur de courant continu à proximité dudit magnétomètre, on introduit5 dans le modèle théorique le produit des éléments suivants: Hix [C] = Hiy et [E] = I, Hiz I étant le courant continu consommé par ledit équipement et Hix, Hiy et Hiz des coefficients de proportionnalité à déterminer. En outre, un mode de réalisation avantageux de l'invention permet d'identifier et de corriger un défaut de linéarité des mesures du magnétomètre selon au moins l'un des trois axes de mesure possibles dudit magnétomètre. A cet effet, selon l'invention, on détermine une approxima- tion polynomiale de la fonction non linéaire reproduisant ce défaut de linéarité et, pour identifier le terme correcteur d'une harmonique n de ladite approximation polynomiale, on introduit dans le modèle théorique le produit des éléments suivants: axn 0 0 Hmxn [C] = 0 ayn 0 et [E] = Hmyn 0 0 azn Hmzn avec: - axn, ayn et azn, des coefficients à déterminer représen- tant les amplitudes respectives, selon les trois axes de mesure du magnétomètre, des corrections de non-linéarité d'ordre n; et - Hmxn, Hmyn et Hmzn, les puissances nièmes des composantes respectives du champ mesuré Hm selon les trois axes de mesure du magnétomètre. Lors de la mise en oeuvre du procédé conforme à l'invention, on a supposé que le champ magnétique est constant pour toutes les mesures de champ magnétique effectuées. Ceci n'est toutefois pas toujours vérifié en pratique. Aussi, un mode de réalisation avantageux de l'invention permet de déterminer les variations locales des composantes du champ magnétique terrestre, ce qui permet d'améliorer la précision du modèle correcteur obtenu par la présente invention. A cet effet, selon l'invention, on introduit dans le modèle théorique le produit d'une matrice [C], dont les coeffi- cients représentent les différents gradients desdites composantes du champ magnétique terrestre, et d'une matrice [E], dont les coefficients sont calculés à partir des coefficients de la matrice de passage [M] et à partir de la position dans un repère déterminé dudit véhicule lors des mesures du magnétomètre. Finalement, un dernier mode de réalisation de l'invention permet d'identifier et de corriger les biais de gyromètres, dont les mesures sont utilisées notamment pour déterminer les coefficients de la matrice [M], ce qui met bien en évidence la diversité de mises en oeuvre possibles de l'invention qui n'est pas limitée à la seule correction d'erreurs dues à des perturbations magnétiques et/ou mécani- ques, telles que celles précisées ci-dessus. Selon ce dernier mode de réalisation, de façon avantageuse, on introduit dans le modèle théorique le produit des élé-30 ments suivants: 0 -Awz Awy Hmx.t [C] = Awz 0 -Awx et [E] = Hmy.t -Awy Awx 0 Hmz.t dans lesquels: Awx, Awy et Awz sont les biais desdits gyromètres à déterminer, - Hmx, Hmy et Hmz sont les composantes respectives du champ mesuré Hm dans le référentiel du véhicule, et - t est l'instant auquel est effectuée la mesure. La présente invention concerne également un dispositif pour la mise en oeuvre du procédé précité. Selon l'invention, ledit dispositif est remarquable en ce qu'il comporte: - un premier module de calcul recevant des données disponi- bles sur ledit véhicule et susceptible de déterminer à partir de ces données les coefficients des matrices [E]i; - un second module de calcul susceptible de calculer les coefficients de la matrice de passage [M], utiles pour la mise en oeuvre de l'invention, à partir des assiettes et si nécessaire du cap du véhicule; et - un module principal de calcul relié auxdits premier et second modules de calcul, ainsi qu'au magnétomètre, et susceptible de déterminer ledit modèle correcteur et ledit champ magnétique effectif. On notera que la présente invention concerne en outre deux systèmes comportant le dispositif conforme à l'invention, à25 savoir, d'une part, un système de référence d'attitude et de cap embarqué à bord d'un aéronef et comportant un magnétomè- tre utilisé pour corriger une éventuelle dérive de la mesure du cap et, d'autre part, une centrale inertielle d'un aéronef, qui est associée à un magnétomètre. Dans le système de référence d'attitude et de cap, le dispositif conforme à l'invention est utilisé pour corriger les erreurs de mesure dudit magnétomètre, tandis que, dans la centrale inertielle, ledit dispositif est utilisé pour 5 déterminer de façon précise, en collaboration avec ledit magnétomètre, le champ magnétique terrestre, ainsi que pour corriger les erreurs de mesure dudit magnétomètre. Les figures du dessin annexé feront bien comprendre comment l'invention peut être réalisée. Sur ces figures, des réfé-10 rences identiques désignent des éléments semblables. La figure 1 montre de façon schématique un véhicule compor- tant un dispositif conforme à l'invention. La figure 2 est le schéma synoptique du dispositif conforme à l'invention. La figure 3 est une vue latérale de la poutre de queue d'un hélicoptère portant un magnétomètre. La figure 4 est une vue transversale de la poutre de queue de la figure 3. La figure 5 est une vue de dessus de la poutre de queue de la figure 3. Le dispositif 1 conforme à l'invention est destiné à corri- ger des erreurs dans les mesures d'un magnétomètre 2 monté à bord d'un véhicule 3 qui est représenté schématiquement sous la forme d'un hélicoptère sur la figure 1. Le magnétomètre 2 est susceptible de mesurer, dans un référentiel R1 qui est associé audit véhicule 3, les compo- santes d'une grandeur physique magnétique orientée et d'amplitude constante localement, représentée par un vecteur, tel le vecteur H orienté dans l'exemple représenté selon le gradient du potentiel magnétique terrestre. Le champ magnétique Hm mesuré par ledit magnétomètre 2 vérifie la relation Hm = [M].H, dans laquelle H représente la valeur5 effective du champ magnétique dans un référentiel de base R2 et [M] représente la matrice de passage dudit référentiel de base R2 audit référentiel du véhicule R1. Ledit référentiel de base R2 peut par exemple être formé par trois axes V, N et E définis comme suit:10 - V est un axe parallèle au gradient du champ de gravité terrestre, orienté vers le centre de la terre, et donnant la verticale locale; - N est un axe horizontal donnant le nord géographique ou directionnel; et - E est un axe donnant l'est, tel que représenté sur la figure 1. De façon connue, le passage de l'un à l'autre de ces réfé- rentiels R1 et R2 nécessite trois rotations successives définies par les angles d'Euler, et À non représentés,20 parmi lesquels: - est l'angle de cap, - G est l'angle d'assiette longitudinale, et - ' est l'angle d'assiette transversale. Ces trois rotations s'écrivent selon la matrice usuelle des angles d'Euler, c'est-à-dire la matrice de passage [M] 1 permettant le passage du référentiel R1 vers le référentiel R2, définie comme suit: -1 [M] cosD.cosG cosp.sin0.sinc - sinp.cosc sinl.sin> + cosP.sinO.cosc sinP.cosG cos.cose + sini.sinO.sinc sinp.sine.cosq - cosi.sinO - sin cose.sine cose.cose Toutefois, un tel magnétomètre 2 est souvent soumis à diverses perturbations, par exemple des perturbations mécaniques ou des perturbations magnétiques, de sorte que les mesures effectuées présentent des erreurs. Le dispositif 1 conforme à l'invention est destiné à identi- fier simultanément toutes les erreurs existantes, quelles que soient les origines de ces erreurs, et à les corriger. On décrit ci-après un premier mode de mise en oeuvre de l'invention permettant de déterminer simultanément: - un modèle correcteur permettant l'identification et la correction de toutes les erreurs apparaissant dans les mesures du magnétomètre 2; et - les trois composantes du champ magnétique H dans le référentiel de base. A cet effet, on définit tout d'abord un modèle théorique correcteur des erreurs de mesure du magnétomètre de la forme: i=n Hc = [A].Hm + Z ([C]i.[E]i) + Hp (1) i=l n étant un nombre entier supérieur ou égal à 1, dans lequel, en plus de l'élément Hm précité: - Hc est le champ corrigé; [A] est une matrice à déterminer, prenant en compte les perturbations dites de "fers doux" et les mésalignements, sans limitation de valeur angulaire; - les éléments [C]i sont des matrices à déterminer de dimensions respectives (3xmi), mi étant un nombre entier supérieur ou égal à 1; - les éléments [E]i sont des matrices de dimensions respec- tives (mixl) dont les coefficients sont des mesures de paramètres perturbateurs ou des combinaisons de telles mesures, comme précisé cidessous; et - Hp est un champ perturbateur à déterminer. Ainsi, le modèle théorique comporte la somme de produits matriciels: i=n E ([C]i.[E]i) = ([C]l.[E]1) + ([C]2.[E]2) + *- + ([C]n.[E]n) i=l avec: cll c12... clmi [C]i = c21 c22... c2mi c31 c32... c3mi El E2 et [E]i = Emi L J i On considère que le champ corrigé Hc est le champ terrestre projeté dans le référentiel du véhicule de sorte que: i=n [A].Hm + E ([C]i.[E]i) + Hp = [M].H (2) i=l Puis on effectue, lors de manoeuvres dudit véhicule 3, des mesures successives du champ magnétique, au moyen du magné- tomètre 2, et des mesures de cap et d'assiettes au moyen d'un système approprié, comme par exemple une centrale de cap et d'attitude en mode gyroscopique libre non représentée et montée sur ledit véhicule 3, ainsi qu'éventuellement, c'est-à-dire si elles sont nécessaires pour la mise en20 oeuvre de l'invention, des mesures permettant de déterminer les coefficients des éléments [E]i, tel que précisé ci-après à partir de modes de réalisation particuliers. Ces manoeuvres peuvent correspondre, par exemple pour un aéronef, à deux virages, à droite puis à gauche, de 360 avec un roulis de l'ordre de 30 . Les mesures d'assiettes et éventuellement de cap sont utilisées pour calculer la matrice de passage [M] telle que définie précédemment. L'expression précédente (2) du modèle théorique peut s'é- crire sous la forme du système suivant: j=p E (aj) = Rll. Hn + R12.He + R13.Hv j=l j=p E (Bj) = R21 Hn + R22. He + R23.Hv j=l1 j=p E (yj) = R31 Hn + R32. He + R33.Hv j=l dans lequel: - les termes Rpq, (p = 1, 2, 3 et q = 1, 2, 3) sont les coefficients de la matrice [M]; - Hn, He et Hv sont les composantes du vecteur H; - les sommes des termes (aj), des termes (Bj) et des termes ( j) sont les expressions de gauche respectivement des trois lignes de ladite expression (2); et i=n - p= E (mi)+4 i=l On remarquera que ce système ne peut être déterminé qu'à une constante multiplicative près. Pour lever l'indétermination, on peut diviser les différentes équations par un coefficient apparaissant dans chacune d'elles, par exemple Hn, He ou Hv. Toutefois, comme les coefficients He et Hv sont susceptibles de présenter des valeurs très faibles dans certaines régions du globe terrestre, ce qui risque de poser des problèmes de5 calcul, on divise de préférence lesdites équations par Hn de manière à obtenir le système suivant: j=p E (aj)/(Hn) = Rll + R12.he + R13.hv ]=1 j=p E (13j)/(Hn) = R21 + R22 he + R23 hv j=1 j=p E (yj)/(Hn) = R31 + R32.he + R33.hv j=33 avec he = He/Hn et hv = Hv/Hn. Puis on considère le vecteur d'erreur E de composantes Ex, Ey et Ez, défini par la relation: i=n E = [M].H - ([A].Hm + E ([C]i.[E]i) + Hp) i=l Toutefois, afin de pouvoir utiliser le système précédent, on détermine un vecteur d'erreur modifié E de composantes ex, cy et cz vérifiant la relation c = E / Hn. A partir desdites composantes, on détermine une erreur globale de la forme E2 = Ex2 + Ey2 + Ez2 ou de la forme C2 = EX2 + Ey2 + CZ2. Puis on définit un système d'équations formé de m équations du type (aEe2)/acl = 0, 1 = 1 à m, ou du type (âEE2)/3cl = 0, (aEc2)/9cl correspondant à la dérivée partielle de la somme des erreurs globales Ec2, pour toutes les mesures, par rapport à un coefficient cl, les différents coefficients cl représentant les coefficients à déterminer dudit modèle théorique de perturbation, à savoir les coefficients de la matrice [A], des n matrices [C]i et du vecteur Hp, ainsi que les composantes du champ magnétique effectif H; et on résout par une méthode numérique ledit système d'équations. Finalement, afin de retrouver des coefficients de dimensions exactes, on divise encore tous les coefficients obtenus précédemment par un diviseur approprié, à savoir par exemple le terme Qll = (All)/Hn. Les coefficients obtenus alors constituent le modèle correc- teur permettant de déterminer les valeurs exactes du champ magnétique à partir de valeurs de champ mesurées par le magnétomètre 2, c'est-à-dire de corriger les erreurs desdi- tes valeurs de champ mesurées. Le mode de réalisation de l'invention décrit ci-dessus permet donc de déterminer simultanément ledit modèle correc- teur et les trois composantes du champ magnétique terrestre. On décrit ci-dessous un mode de réalisation simplifié n'utilisant pas de référence de cap, et ne permettant de déterminer, en plus dudit modèle correcteur, que la compo- sante verticale Hv dudit champ magnétique terrestre. A cet effet, on écrit l'expression (2) du modèle théorique sous la forme: -1 i=n [M].([A].Hm + Z ([C]i.[E]i) + Hp = H i=l c'est-à-dire: ]=p E (craj) j=l -1 i=P Hn [M]. E (13j) He j=1 Hv i=p E (Yij) j=l Dans la présente hypothèse, les composantes Hn et He du champ magnétique effectif H sont inaccessibles du fait de l'absence de référence de cap. Par conséquent, seule la5 troisième équation du système précédent est exploitable à des fins d'identification d'une composante du champ magnéti- que, ladite troisième équation s'écrivant: j=p j=p i=p R31*.( E (aj)) + R32*.( E (Bj)) + R33 *.( E (yj)) = Hv 31 32 331j= j=l j=l j=l R31*, R32* et R33* sont les trois composantes de la troisième ligne de la matrice de passage [M] définie ci-dessus, c'est-à-dire que: R31* = - sin8 R32 * = cosO.sin R33 * = cosO.cosD. Ces coefficients R31*, R32* et R33* sont donc définis uniquement à partir des angles d'assiettes 8 et b, et non à partir de l'angle de cap. Comme dans le mode de réalisation précédent, cette équation admet une infinité de solutions à une constante multiplicative près. Pour lever l'indétermination existante, on divise les deux membres de l'équation par Hv. On notera que, dans ce cas, on détermine uniquement la composante Ez du vecteur d'erreur g. Puis on définit un système d'équations formé d'équations du type (EEz2)/cl = 0, les coefficients cl représentant de nouveau les coeffi- cients à déterminer du modèle théorique (1), et on résout ce système de façon usuelle. Les coefficients ainsi calculés sont divisés par un diviseur approprié pour obtenir les dimensions correctes et le modèle correcteur recherché. La mise en oeuvre est donc identique à celle décrite précé- demment. Pour ces mises en oeuvre du procédé conforme à l'invention, le dispositif 1, qui est relié par une liaison 4 au magnéto- mètre 2, comporte, selon l'invention et tel que représenté sur la figure 2: - un module de calcul 5 relié au magnétomètre 2 par ladite liaison 4 et recevant, par l'intermédiaire d'une liaison multiple 6, des données disponibles sur le véhicule 3 et précisées ci-après pour des modes de réalisation parti- culiers, ledit module de calcul 5 étant susceptible de déterminer les coefficients des matrices [E]i; - un module de calcul 7 recevant d'un dispositif approprié non représenté, d'une part, les valeurs des assiettes du véhicule 3 par l'intermédiaire d'une liaison 8 et, d'autre part, éventuellement, c'est-à-dire si elle est nécessaire pour la mise en oeuvre de l'invention, la valeur du cap du véhicule 3 par l'intermédiaire d'une liaison 9, ledit module de calcul 7 étant susceptible de calculer les coefficients de la matrice de passage [M] utilisés lors de la mise en oeuvre de l'invention; et - un module principal de calcul 10 relié auxdits modules de calcul 5 et 7 et audit magnétomètre 2, par l'intermédiaire respectivement de liaisons 11 et 12 et de la liaison 4, ledit module principal de calcul 10 étant susceptible de déterminer ledit modèle correcteur et ledit champ magnéti- que effectif et pouvant transmettre les résultats obtenus par une liaison 13, par exemple à un dispositif utilisa- teur non représenté. On notera que le module de calcul 5 est relié en particulier à des dispositifs de mesure et/ou à des capteurs qui lui fournissent, par l'intermédiaire de la liaison multiple 6, des données spécifiques précisées ci-dessous et traitées par ledit module de calcul 5. Toutefois, dans certains cas, comme on le verra également ci-après, les données disponi-15 bles correspondent telles quelles aux coefficients recher- chés des matrices [E]i de sorte que le module de calcul 5 ne sert alors qu'à transmettre ces données au module de calcul principal 10. Bien entendu, ledit module principal de calcul 10 et lesdits modules de calcul 5 et 7 peuvent être implantés soit dans un calculateur unique, soit dans des calculateurs différents. Avant de préciser différents modes de réalisation parti- culièrement avantageux de l'invention, il convient encore de noter que le dispositif 1 conforme à l'invention peut être utilisé avantageusement par différents types de systèmes non représentés et embarqués à bord dudit véhicule 3. Par exemple, un système de référence d'attitude et de cap d'un aéronef, comportant entre autres un magnétomètre pour corriger une éventuelle dérive du cap, peut utiliser le dispositif 1 pour corriger les erreurs dues à des perturba- tions multiples et variées dans les mesures dudit magnéto- mètre. Par ailleurs, une centrale inertielle d'un aéronef, associée à un magnétomètre, peut utiliser le dispositif 1 conforme à l'invention pour déterminer de façon précise, en collabora- tion avec ledit magnétomètre, la déclinaison et l'inclinai- 5 son du champ magnétique terrestre, ainsi que pour corriger les erreurs de mesure dudit magnétomètre. On présente à présent différents modes de réalisation avantageux permettant de bien mettre en évidence la multi- plicité d'utilisations possibles de l'invention. Bien que10 ces différents modes de réalisation soient précisés les uns après les autres ci-dessous, ils peuvent bien entendu être mis en oeuvre simultanément. Ainsi, la présente invention est particulièrement bien adaptée à la correction d'erreurs de mesure du magnétomètre,15 dues à des déformations de l'élément porteur sur lequel est monté ledit magnétomètre. De telles déformations sont susceptibles d'entraîner de légers déplacements du magnéto- mètre, provoquant en particulier un défaut d'alignement du trièdre de référence de ce dernier par rapport au trièdre de20 référence du véhicule, ce qui est source d'erreurs. On sait que, dans le cas d'un hélicoptère 3, le magnétomètre 2 est généralement monté dans la poutre de queue 14, tel que représenté sur les figures 3 à 5, de manière à être le plus éloigné possible des principaux équipements électriques et électroniques et des organes massifs de l'hélicoptère 3, susceptibles de perturber ledit magnétomètre 2. Toutefois, comme la poutre de queue 14 n'est pas parfaitement rigide, elle se déforme sous l'effet des contraintes mécaniques auxquelles elle est soumise, ce qui peut entraîner un défaut d'alignement dudit magnétomètre 2. On sait en outre qu'une erreur d'alignement en assiette se répercute au double sur l'erreur de cap, dans le cas d'une inclinaison de 60 du champ terrestre, telle qu'existant par exemple en France. Comme les déformations observées sont de l'ordre de 0,2 , l'erreur résultante en cap est de l'ordre de 0,4 , cette valeur correspondant environ à la performance5 recherchée sur les centrales de référence d'attitude et de cap actuelles. Une telle erreur n'est donc pas négligeable. Généralement, la poutre de queue 14 est soumise à un effort latéral Fv. Le couple rotor principal de l'hélicoptère 3 provoque alors, sur ladite poutre de queue 14:10 - d'une part, une flexion en lacet A& représentée sur la figure 5 et vérifiant la relation A; = Cez.Fv, Cez étant le coefficient d'élasticité de flexion en lacet de la poutre de queue 14; et -d'autre part, une torsion en roulis A' représentée sur la figure 4 et vérifiant la relation AP = Cex.Fv, Cex étant le coefficient d'élasticité de torsion de la poutre de queue 14. De plus, dans le cas d'un hélicoptère dont l'empennage 15 présente une symétrie latérale, c'est-à-dire avec des20 envergures identiques à gauche et à droite, de longueur D, comme représenté sur la figure 5, l'effort vertical (ou la portance) Fh de cet empennage provoque une flexion en tangage AO représentée sur la figure 3 et vérifiant la relation AL = Cey.Fh, Cey étant le coefficient d'élasticité25 de flexion en tangage de la poutre de queue 14. Par conséquent, les valeurs A$, AL, AL représentent les variations angulaires de la position du magnétomètre, respectivement en lacet, tangage et roulis, dans le cas o la poutre de queue 14 est soumise aux contraintes précitées, par rapport à la position existant lorsque la poutre de queue 14 est hors contrainte. Selon l'invention, pour identifier et corriger les erreurs dues à de telles déformations dans les mesures du magnétomè- tre, on introduit dans le modèle théorique (1) le produit des matrices [C] de dimension (3,5) et [E] de dimension5 (5,1) telles que: 0 -Cez 0 0 Cey [C] = Cez 0 -Cex 0 0 I 0 Cex 0 -Cey 0 et Fv. Hmx Fv.Hmy [E] = Fv.Hmz Fh.Hmx Fh.Hmz avec Hmx, Hmy et Hmz, les composantes respectives du champ mesuré Hm dans le référentiel du véhicule R1. Par conséquent, les coefficients de la matrice [C] sont déterminés par le procédé conforme à l'invention et ceux de la matrice [E] par le module de calcul 5 à partir de mesures effectuées de manière connue. A cet effet, ledit module de15 calcul 5 est relié par l'intermédiaire d'une liaison non représentée au magnétomètre 2 pour recevoir lesdites compo- santes du champ mesuré Hm. Le produit [C].[E] constitue donc un vecteur de correction des erreurs de mesure du champ magnétique, provoquées par20 les contraintes citées précédemment. Ce vecteur de correction peut être défini, pour bien expli- citer la précision de la correction obtenue, à partir d'une matrice [Msd] permettant de passer de la position du magné- tomètre 2, lorsque la poutre de queue 14 est hors contrainte, à la position dudit magnétomètre 2, lorsque ladite poutre de queue 14 est sous contrainte. Cette matrice vaut: 1 - AO [Msd] A 1 -a -A83 As 1 de sorte que le champ corrigé Hc s'écrit alors, selon l'expression (2) précisée ci- dessus du modèle théorique, Hc = [Msd].[M].H En négligeant les infiniment petits du second ordre, il subsiste un terme: O -A0 AO A[ 0 -A.H -As AO 0 qui correspond à: 0 -Cez.Fv Cey.Fh Hmx Cez.Fv 0 -Cex.Fv Hmy = [C].[E] l-Cey.Fh Cex.Fv 0 Hmz Ces explications permettent de préciser ce que représente le vecteur de correction obtenu par le procédé conforme à l'invention et susceptible d'être appliqué essentiellement dans le cas d'un hélicoptère dont la poutre de queue porte le magnétomètre. En revanche, dans le cas d'un avion non représenté, par exemple un avion de transport, le magnétomètre est générale- ment installé à l'extrémité d'une aile également pour être éloigné le plus possible d'éléments perturbateurs. Or, en raison de l'élasticité de la structure de ladite aile, celle-ci est susceptible d'être déformée, ce qui peut provoquer également un défaut d'alignement du magnétomètre. La déformation principale existante et prise en compte dans la présente invention est la déformation en battement de. ladite aile At Cette déformation en battement A, qui est due à la por- tance de l'aile Pz, s'écrit Ab = Cox.Pz, Cox étant le coefficient d'élasticité en battement de l'aile. Selon l'invention, pour identifier et corriger les erreurs dues à un défaut d'alignement du magnétomètre provoqué par le battement de l'aile portant ledit magnétomètre, on introduit dans le modèle théorique (1) le produit des éléments suivants: 0 0 Pz.Hmz [C] = -Cox 0 et [E] = 0 Cox Pz.Hmy La portance Pz vérifie la relation Pz = M.gz dans laquelle: - M représente la masse de l'avion, pouvant être déterminée * à partir de la masse initiale dudit avion, de la consomma- tion de carburant et, le cas échéant, de largages d'em- ports; et - gz est le facteur de charge vertical susceptible d'être mesuré par un accéléromètre du système de référence d'attitude et de cap ou de la centrale inertielle de l'avion. Ces données sont bien entendu transmises par la liaison 6 au module de calcul 5 qui est susceptible de calculer les coefficients de ladite matrice [E]. On notera que le produit [C].[E] vérifie la relation: 0 0 0 Hmx [C].[E] = 0 0 -Cox.M.gz Hmy 0 Cox.M.gz 0 Hmz ce qui permet de faire le rapprochement avec la solution précisée précédemment dans le cas d'un hélicoptère. On vient de décrire ci-dessus deux modes de réalisation avantageux permettant de corriger des erreurs de mesure du magnétomètre dues à des perturbations mécaniques, à savoir des déformations de l'élément porteur. On décrit à présent ci-dessous plusieurs modes de réalisation permettant de10 corriger des erreurs dues à différentes perturbations d'origine magnétique. Un mode de réalisation, particulièrement avantageux et s'appliquant à un hélicoptère, permet de remédier aux instabilités magnétiques provoquées par le rotor de l'héli-15 coptère qui, d'une part, présente un résidu de magnétisation et, d'autre part, est soumis à des contraintes mécaniques élevées. Comme le rotor d'un hélicoptère subit généralement en phase de fabrication des contrôles non destructifs de type magné-20 toscopie, il subsiste en effet souvent sur celui-ci un résidu de magnétisation. Or, de façon connue, un matériau ferromagnétique aimanté soumis à de fortes contraintes voit son état magnétique varier de façon réversible et proportionnelle aux contrain-25 tes appliquées. Ce phénomène transposé au cas d'un rotor d'hélicoptère peut provoquer des perturbations magnétiques susceptibles d'engendrer des erreurs de mesure non négligeables lorsqu'il s'agit d'un petit hélicoptère sur lequel le magnétomètre ne peut être éloigné suffisamment du mât rotor. Selon l'invention, pour identifier et corriger de telles erreurs, on introduit dans le modèle théorique (1) du procédé conforme à l'invention le produit des éléments suivants: Hcx [C] = Hcy et [E] = Cr, Hcz Cr étant le couple rotor et Hcx, Hcy et Hcz des coefficients de proportionnalité à déterminer par ledit procédé, de manière à obtenir un vecteur de correction défini par ledit produit [C].[E]. Un autre mode de réalisation de l'invention permet d'identi- fier et de corriger des erreurs dues à des perturbations15 magnétiques engendrées par la présence d'équipements consom- mateurs de courant continu à proximité dudit magnétomètre. On sait que, généralement, les deux câbles d'une alimenta- tion électrique cheminent au plus proche l'un de l'autre ou sont torsadés de manière que la surface du circuit ainsi20 constitué soit la plus faible possible et que le passage de courant n'engendre pas de champ magnétique. Toutefois, de telles précautions de montage ne sont pas toujours possi- bles, de sorte qu'il apparaît alors un champ magnétique perturbateur. Selon l'invention, pour corriger un tel champ perturbateur, on détermine un terme correcteur [C].[E] avec: Hix [C] Hiy et [E] - I, Hiz I étant le courant continu consommé par ledit équipement non représenté et Hix, Hiy et Hiz étant des coefficients de proportionnalité déterminés par l'invention. Bien entendu, dans ce cas, un dispositif de mesure du courant I, relié au module de calcul 5, est nécessaire pour la mise en oeuvre de l'invention. Les erreurs de mesure peuvent, par ailleurs, provenir d'un défaut de linéarité des voies de mesure et d'acquisition du10 champ magnétique du magnétomètre. On sait que le fonctionne- ment d'un magnétomètre est généralement basé sur la compen- sation du champ à mesurer par un courant continu de sorte que le matériau constituant le noyau des capteurs du magné- tomètre fonctionne toujours par rapport au même état de15 magnétisation, ce qui garantit en principe une assez bonne linéarité. Or, il arrive en pratique, pour certaines condi- tions de fonctionnement particulières, que cette linéarité ne soit pas vérifiée. Aussi, pour corriger les erreurs de mesure susceptibles d'être provoquées par un tel défaut de linéarité, on déter- mine, selon l'invention, une approximation polynomiale de la fonction non linéaire reproduisant ce défaut de linéarité et, pour obtenir le terme correcteur d'une harmonique n de ladite approximation polynomiale, on introduit dans le25 modèle théorique (1) du procédé conforme à l'invention le produit des éléments suivants: axn 0 0 Hmxn [C] =1 0 ayn 0 et [E] = Hmyn 0 0 azn Hmzn avec: - axn, ayn et azn des coefficients à déterminer, représen- tant les amplitudes respectives, selon les trois axes de mesure possibles du magnétomètre 2, des corrections de non-linéarité d'ordre n; et nn n ièmes - Hmxn, Hmy et Hmz les puissances n des composantes respectives du champ mesuré Hm selon les trois axes de mesure du magnétomètre. Bien entendu, dans le procédé conforme à l'invention et précisé cidessus, on a supposé que le champ magnétique mesuré est constant tout au long des différentes manoeuvres du véhicule effectuées pour la mise en oeuvre dudit procédé. Ceci est généralement vérifié en pratique, puisque:15 - les manoeuvres effectuées sont limitées, par exemple à deux virages à 30 de roulis, à gauche et à droite, de sorte que le territoire parcouru par le véhicule n'est pas trop important; et - la variation du champ terrestre est, dans la plupart des régions du globe terrestre, négligeable pour les varia- tions de position obtenues par de telles manoeuvres. Toutefois, il peut arriver que, dans certaines régions du globe terrestre, les variations du champ magnétique ne soient pas négligeables, lorsque le territoire parcouru au25 cours des mesures est relativement important, ce qui peut entraîner des erreurs dans le modèle correcteur obtenu. Dans la présente invention, on tient compte uniquement des variations linéaires dudit champ magnétique. Ainsi, selon l'invention, pour déterminer les variations locales linéaires des composantes du champ magnétique terrestre, permettant d'améliorer la précision du modèle correcteur, on introduit dans le modèle théorique le produit d'une matrice [C], dont les coefficients représentent les différents gradients desdites composantes du champ magnéti- que, et d'une matrice [E], dont les coefficients sont définis à partir des coefficients de la matrice de passage [M] et à partir de la position dans un repère déterminé du véhicule 3 lors des mesures du magnétomètre 2. En effet, le champ magnétique terrestre H, dont les va- riations sont localement linéaires, peut s'écrire sous la forme: gnn gne gnv pn HnO H = gen gee gev pe + HeO gvn gve gvv pv HvO dans laquelle: (pn, pe, pv) représente le vecteur position du véhicule 3 dans un repère déterminé défini par rapport à une position initiale (O, O, O), par exemple la position au début des mesures; - les vecteurs lignes de la matrice (3x3) sont les gradients respectivement des trois composantes du champ magnétique terrestre; et - (HnO, HeO, HvO) est le champ terrestre à ladite position initiale (O, O, O), de sorte que: gnn gne gnv pn HnO [M].H = [M]. ( gen gee gev pe + HeO) gvn gve gvv pv HvO gnn gne gnv pn HnO = [M]. ( gen gee gev pe) +[M]. HeO gvn gve gvv pv HvO De cette expression, on déduit le terme correctif: Rll R12 R13 gnn gne gnv pn R21 R22 R23 gen gee gev. pe R31 R32 R33 gvn gve gvv pv qui permet de retrouver par simple transcription matricielle le produit [C].[E] défini précédemment: Rll R12 R13 gnn gne gnv pn [C].[E] = R21 R22 R23 gen gee gev. pe R31 32 R33 gvn gve gvv pv La matrice [C] apparait sous la forme donnée ci-après, ayant 3 lignes et 27 colonnes. AA6 eA6 UA5 A96 aa65 Ua AU6 aU5 UU5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 AA5 aA5 UA5 A aa6 UaU AUo au oUU 0 0 0 0 0 0 0 0 0O 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 AA5 aA6 UA6 Aa5 aaB Ua AU6 aU6 UU5 Quant à elle, la matrice [E] présente 1 colonne et 27 lignes, telle que ci-dessous: R11 pn R11 pe R pv R11 P R12 pn R12 pe R12 Pv 13 pn R13 pe R13 pv 21 pn R21 pe R21 Pv R22 pn [E] =R22 pe R22 Pv R23 pn R23 pe R23 Pv R31 pn R31 pe R31 P R32 pn R32 pe R32 P R33 pn R33 pe R33 P Par ailleurs, la présente invention permet d'identifier et de corriger les biais de gyromètres utilisés en particulier pour déterminer les coefficients de la matrice de passage [M] du référentiel de base R2 au référentiel du véhicule R1. Le mode de réalisation concerné s'applique plus précisément aux biais gyrométriques qui sont fixes dans le référentiel10 du véhicule R1, comme par exemple les biais gyrométriques susceptibles d'exister dans les mesures effectuées par une centrale inertielle à composants liés. A cet effet, selon l'invention, pour identifier et corriger de tels biais gyrométriques, on introduit dans le modèle théorique (1) le produit des éléments suivants: 0 -Awz Awy Hmx.t [C] - Awz 0 -Awx et [E] = Hmy.t -Awy Awx 0 Hmz.t dans lesquels: - Awx, Awy, Awz sont les biais desdits gyromètres-à détermi- ner; et - t est l'instant auquel est effectuée la mesure. Les corrections de biais gyrométriques effectuées ainsi o10 permettent en particulier d'améliorer le modèle correcteur déterminé par le procédé conforme à l'invention. Elles permettent également d'utiliser sur le véhicule 3 des gyromètres de performance moyenne et de coût réduit, pour lesquels la précision des mesures est limitée en raison de15 l'existence de biais, puisque grâce à l'invention on est en mesure de corriger de tels biais gyrométriques. Les différents modes de réalisation avantageux décrits précédemment permettent de bien mettre en évidence la multiplicité de corrections possibles et les nombreux20 avantages de la présente invention. Ces différents modes de réalisation peuvent bien entendu être réalisés simultané- ment, à partir d'une seule mise en oeuvre du procédé confor- me à l'invention, dans le cas o le magnétomètre est soumis à toutes les perturbations citées précédemment. On notera en outre que l'invention n'est évidemment pas limitée à ces modes de réalisation particuliers, mais qu'elle englobe toutes les identifications et corrections d'erreurs dans les mesures d'un magnétomètre, susceptibles d'être réalisées par le procédé conforme à l'invention. REVENDICATIONS 1. Procédé pour identifier et corriger des erreurs dues à des perturbations déterminées dans les mesures d'un magnéto- mètre (2) monté à bord d'un véhicule (3), caractérisé en ce que: - on définit un modèle théorique correcteur des erreurs de mesure du magnétomètre de la forme: i=n Hc = [A].Hm + E ([C]i.[E]i) + Hp, i=l n étant un nombre entier supérieur ou égal à 1, dans lequel: 10. Hc est le champ corrigé, [A] est une matrice à déterminer, Hm est la valeur mesurée par le magnétomètre (2) du champ magnétique dans un référentiel (Ri) lié audit véhicule (3), 15. les éléments [C]i sont des matrices à déterminer de dimensions respectives (3xmi), mi étant un nombre entier supérieur ou égal à 1, les éléments [E]i sont des matrices de dimensions respectives (mixl) dont les coefficients sont des mesures de paramètres perturbateurs ou des combinaisons de telles mesures, et Hp est un champ perturbateur à déterminer; - on considère que le champ corrigé Hc est le champ terres- tre défini dans le référentiel du véhicule de sorte que: i=n [A].Hm + E ([C]i.[E]i) + Hp = [M].H i=l avec: H la valeur effective du champ magnétique dans un référentiel de base (R2), et [M] une matrice de passage dudit référentiel de base (R2) audit référentiel du véhicule (R1); - on effectue des mesures Hm successives du champ magnétique au moyen du magnétomètre lors d'évolutions du véhicule; - on détermine les coefficients des différentes matrices [E]i, à partir des mesures des paramètres perturbateurs, effectuées simultanément aux mesures Hm; - on détermine les coefficients de la matrice de passage [M]; - on définit un vecteur d'erreur E de composantes Ex, Ey et Ez, défini par la relation: i=n E = [M].H - ([A].Hm + E ([C]i.[E]i) + Hp) i=1 - on définit une erreur globale E2 = Ex2 + Ey2 + Ez2; - on définit un système d'équations formé de m équations du type (EE2)/âcl = 0, 1 = 1 à m, (aEE2)/acl correspondant à la dérivée partielle de la somme des erreurs globales EE2, pour toutes les mesures, par rapport à un coefficient cl, les différents coefficients cl représentant les coeffi- cients à déterminer dudit modèle théorique, à savoir les coefficients de la matrice [A], des matrices [C]i et du vecteur Hp, ainsi qu'au moins une composante du champ magnétique effectif H; - on résout ledit système d'équations de manière à obtenir lesdits coefficients cl, dont ladite composante du champ magnétique effectif H; et - on utilise le modèle correcteur obtenu à partir desdits coefficients pour corriger les erreurs dues à des pertur- bations déterminées dans les mesures dudit magnétomètre (2). 2. Procédé selon la revendication 1, pour déterminer unique- ment la composante verticale du champ magnétique, caractérisé en ce qu'on utilise comme erreur globale E2 = Ez2, de sorte que seuls les coefficients de la troisième ligne de la matrice [M] 1 sont utilisés, ladite matrice [M]- 1 étant la matrice de passage dudit référentiel du véhicule (R1) audit référentiel de base (R2), lesdits coefficients étant déterminés à partir des assiettes dudit véhicule (3). * 3. Procédé selon la revendication 1, pour déterminer les trois composantes du champ magnétique, caractérisé en ce que les coefficients de la matrice [M] sont déterminés à partir des assiettes, ainsi qu'à partir du cap dudit véhicule (3). 4. Procédé selon la revendication 3, caractérisé en ce que ledit cap est un cap géographique. 5. Procédé selon la revendication 3, caractérisé en ce que ledit cap est un cap délivré par une centrale d'attitude et de cap en mode gyroscopique libre. 6. Procédé selon l'une quelconque des revendications 1 à 5, le magnétomètre (2) étant monté sur un élément porteur (14) susceptible d'être déformé lors de la variation de contrain- tes mécaniques agissant sur ledit élément porteur, caractérisé en ce que, pour identifier et corriger les erreurs dues à de telles déformations dans les mesures du magnétomètre (2), on introduit dans le modèle théorique le produit d'une matrice [C] de dimension (3,5) et d'une matrice [E] de dimension (5,1), telles que: 0 -Cez 0 0 Cey [C]x Cez 0 -Cex 0 0 0 Cex 0 -Cey 0 et Fv.Hmx Fv.Hmy [E] = Fv.Hmz Fh.Hmx Fh.Hmz avec: - Fv: l'effort latéral engendré sur l'élément porteur (14), - Fh: l'effort vertical engendré sur l'élément porteur (14), - Cex: le coefficient d'élasticité de torsion de l'élément porteur (14), - Cey: le coefficient d'élasticité de flexion en tangage de l'élément porteur (14), - Cez: le coefficient d'élasticité de flexion en lacet de l'élément porteur (14), et - Hmx, Hmy et Hmz: les composantes respectives du champ mesuré Hm dans le référentiel du véhicule (R1). -7. Procédé selon la revendication 6, ledit véhicule étant un hélicoptère (3), caractérisé en ce que l'élément porteur du magnétomètre (2) est la poutre de queue (14) dudit hélicoptère (3). 8. Procédé selon l'une quelconque des revendications 1 à 5, ledit véhicule étant un avion et le magnétomètre étant monté sur une aile dudit avion, caractérisé en ce que, pour identifier et corriger les5 erreurs dues à un défaut d'alignement du magnétomètre provoqué par le battement de ladite aile portant ledit magnétomètre, on introduit dans le modèle théorique le produit des éléments suivants: 0 0 Pz.Hmz [C] -Cox 0 et [E] = 0 Cox Pz.Hmy dans lesquels: - Cox est le coefficient d'élasticité en battement de ladite aile, - Pz est la portance de ladite aile, et - Hmy et Hmz sont les valeurs des composantes latérale et verticale du champ magnétique mesuré Hm. 9. Procédé selon la revendication 8, caractérisé en ce que l'on détermine la portance Pz à partir de la relation Pz = M.gz, avec M la masse du véhicule et qz le facteur de charge vertical. *10. Procédé selon l'une des revendications 1 à 7, ledit véhicule étant un hélicoptère dont le rotor présente un résidu de magnétisation, caractérisé en ce que, pour identifier et corriger les erreurs dues à des variations des perturbations magnétiques25 provoquées par ledit rotor aimanté soumis à des contraintes mécaniques élevées, on introduit dans le modèle théorique le produit des éléments suivants: [ Hcx l [C] = Hcy et [E] = Cr, Hcz Cr étant le couple rotor et Hcx, Hcy et Hcz des coefficients de proportionnalité à déterminer. e 11. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que, pour identifier et corriger les erreurs dues à des perturbations magnétiques engendrées par la présence d'au moins un équipement consommateur de courant continu à proximité dudit magnétomètre (2), on introduit dans le modèle théorique le produit des éléments suivants: Hix [C] = Hiy et [E] = I, Hiz I étant le courant continu consommé par ledit équipement et Hix, Hiy et Hiz des coefficients de proportionnalité à déterminer. ' 12. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que, pour identifier et corriger un défaut de linéarité des mesures du magnétomètre (2) selon au moins l'un des trois axes de mesure possibles dudit magnétomètre20 (2), on détermine une approximation polynomiale de la fonction non linéaire reproduisant ce défaut de linéarité, et en ce que, pour identifier le terme correcteur d'une harmonique n de ladite approximation polynomiale, on intro- duit dans le modèle théorique le produit des éléments25 suivants: a i ni axn 0 0 Hmx [C] = 0 ayn 0 et [E] = Hmy 0 0 azn Hmzn avec: - axn, ayn et azn, des coefficients à déterminer représen- tant les amplitudes respectives, selon les trois axes de mesure du magnétomètre (2), des corrections de non- linéarité d'ordre n; et - Hmxn, Hmyn et Hmzn, les puissances nimes des composantes respectives du champ mesuré Hm selon les trois axes de mesure du magnétomètre (2). 10. 13. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que, pour déterminer les variations locales des composantes du champ magnétique terrestre, permettant d'améliorer la précision du modèle correcteur, on15 introduit dans le modèle théorique le produit d'une matrice [C], dont les coefficients représentent les différents gradients desdites composantes du champ magnétique terres- tre, et d'une matrice [E], dont les coefficients sont calculés à partir des coefficients de la matrice de passage [M] et à partir de la position dans un repère déterminé dudit véhicule (3) lors des mesures du magnétomètre (2). 14. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que, pour identifier et corriger les biais25 de gyromètres, dont les mesures sont utilisées notamment pour déterminer les coefficients de la matrice [M], on introduit dans le modèle théorique le produit des éléments suivants: 0 -Awz &wy F Hmx.t [C] = Awz 0 -Awx et [E] = Hmy.t -Awy Awx 0 Hmz.t dans lesquels: - Awx, Awy et &wz sont les biais desdits gyromètres à déterminer, - Hmx, Hmy et Hmz sont les composantes respectives du champ mesuré Hm dans le référentiel du véhicule (R1), et - t est l'instant auquel est effectuée la mesure. 15. Dispositif pour la mise en oeuvre du procédé spécifié sous l'une quelconque des revendications 1 à 14, caractérisé en ce qu'il comporte: - un premier module de calcul (5) recevant des données disponibles sur ledit véhicule (3) et susceptible de déterminer à partir de ces données les coefficients des matrices [E]i; un second module de calcul (7) susceptible de calculer les coefficients de la matrice de passage [M], utiles pour la mise en oeuvre de l'invention, à partir des assiettes et si nécessaire du cap du véhicule (3); et - un module principal de calcul (10) relié auxdits premier et second modules de calcul (5, 7), ainsi qu'au magnétomè- tre (2), et susceptible de déterminer ledit modèle correc- teur et au moins une composante du champ magnétique effectif. t 16. Dispositif selon la revendication 15, caractérisé en ce que ledit module principal de calcul (10) et lesdits premier et second modules de calcul (5, 7) sont implantés dans des calculateurs différents. 17. Système de référence d'attitude et de cap embarqué à bord d'un aéronef et comportant un magnétomètre pour corri- ger une éventuelle dérive de la mesure du cap, caractérisé en ce qu'il comporte le dispositif (1) spécifié5 sous l'une des revendications 15 ou 16, ledit dispositif (1) étant utilisé pour corriger les erreurs de mesure dudit magnétomètre. 18. Centrale inertielle d'un aéronef, associée à un magnéto- mètre, o10 caractérisée en ce qu'elle comporte le dispositif (1) spécifié sous l'une des revendications 15 ou 16, ledit dispositif (1) étant utilisé pour déterminer de façon précise, en collaboration avec ledit magnétomètre, le champ magnétique terrestre, ainsi que pour corriger les erreurs de15 mesure dudit magnétomètre.