Brevet FR2735224A1 — publié le 19961213

La présente invention concerne un capteur électromagnétique et un procédé de détermination simultanée de la perméabilité magnétique et la conductivité d'une cible métallique ainsi que la distance de cette cible par rapport à un capteur électromagnétique. On sait qu'un capteur électromagnétique comporte une bobine associée à un noyau à forte perméabilité magnétique. On sait également que si on dispose le capteur électromagnétique en face d'une cible, c'est-à-dire un objet métallique, et on excite la bobine par un courant alternatif, l'impédance du capteur électromagnétique varie en fonction de paramètres qui sont la fréquence du courant d'excitation, la perméabilité magnétique de la cible, la conductivité électrique de la cible et la distance du capteur par rapport à la cible. Lorsque la perméabilité magnétique et la conductivité électrique d'une cible sont connues il est donc possible d'utiliser un capteur électromagnétique pour déterminer la distance du capteur électromagnétique par rapport à la cible, en particulier pour effectuer des mesures dimensionnelles de pièces métalliques.Toutefois, il s'avère que lors de la fabrication d'une pièce métallique la perméabilité magnétique et la conductivité électrique dépendent non seulement du métal utilisé, mais également des différents traitements auquel le métal est soumis lors de la fabrication de sorte qu'il n'est pas possible de considérer que différentes pièces fabriquées successivement à partir d'un même métal auront la même perméabilité magnétique et la même conductivité électrique. En outre la perméabilité magnétique varie avec la température. Lorsque le capteur électromagnétique est utilisé pour effectuer des mesures dimensionnelles, il n'est donc pas possible d'étalonner le capteur préalablement à la mesure d'une série de pièces. Par ailleurs, les capteurs électromagnétiques connus actuellement ont généralement un circuit magnétique en forme de U, la bobine de mesure étant alors disposée sur une branche du U, ou en forme de pot comprenant un plot central sur lequel est disposé la bobine. Dans un cas comme dans l'autre le nombre important d'interfaces entre les différents milieux qui composent le capteur électromagnétique rend inexploitable en pratique la réalisation d'un modèle permettant de déterminer simultanément les différents paramètres qui influent sur l'impédance du capteur électromagnétique. Le but de l'invention est de proposer un capteur électromagnétique pouvant être modélisé, et un procédé de détermination simultanée de la perméabilité magnétique et la conductivité d'une cible métallique ainsi que la distance de cette cible par rapport à un capteur électromagnétique. Selon l'invention on propose un capteur électromagnétique comportant au moins une bobine associée à un noyau à forte perméabilité magnétique, dans lequel le noyau est une plaque plane, la bobine est fixée sur une face du noyau et comporte un axe de symétrie s'étendant perpendiculairement à la face du noyau sur laquelle elle est fixée, et le noyau a une épaisseur supérieure à la dimension de la bobine prise selon son axe de symétrie de façon que les lignes de champs magnétique résultant de l'excitation de la bobine soient sensiblement dans leur totalité incluses dans le noyau. On obtient ainsi un capteur électromagnétique qui tout en ayant une sensibilité suffisante pour effectuer des mesures présente un nombre d'interfaces réduit et des fuites magnétiques suffisamment faibles pour permettre une modélisation. Selon un aspect avantageux de l'invention, le noyau et la bobine sont des cylindres circulaires coaxiaux, le noyau ayant un diamètre supérieur à un diamètre externe de la bobine de façon que les lignes de champ magnétique se referment sensiblement en totalité dans l'espace délimité par le noyau et la cible en regard. On obtient ainsi un capteur ayant des fuites magnétiques très faibles même lorsque celui-ci est disposé en face d'une cible très conductrice et ayant une perméabilité magnétique faible telle que du laiton. Selon l'invention on prévoit également un procédé de détermination simultanée de la perméabilité magnétique et la conductivité d'une cible métallique ainsi que la distance de cette cible par rapport à un capteur électromagnétique, ce procédé comportant les étapes de a) déterminer l'impédance normalisée complexe, ZN = RN + j X N du capteur électromagnétique à différentes fréquences d'excitation (f1, f2, f,... fQ), et pour différents exemples d'étalonnage chacun défini par une perméabilité magnétique et une conductibilité de la cible ainsi que par une distance du capteur électromagnétique par rapport à la cible, avec RN = R-Ro et XN = X X0 X0 où R et X sont respectivement la composante résistive et la composante réactive du capteur électromagnétique en présence de la cible, et R0 et X0 sont respectivement la composante résistive et la composante réactive du capteur électromagnétique en l'absence de cible. b) établir la matrice d'un plan d'expérience 1 RN1(f1) XN1(f1) RN1(f2) XN1(f2) ... RN1(fk) XN1(fk) ... RN1(fM) XN1(fM) 1 RN2(f,) XN2(f,) RN2(f2) XN2(2) ... RN2(fk) XN2(fk) ... RN2(fM) XN2(fM) # = 1 ...... ....... ....... ....... ....... ....... ....... ....... #1 RNi(f1) XNi(f1) RNi(f2) XNi(f2) ... RNi(fk) XNi(fk) ... RNi(fM) XNi(fM) 1 ...... ....... ....... ....... ....... ....... ....... ....... 1 RNQ(f1) XNQ(f1) RNQ(f2) XNQ(f2) ... RNQ(fk) XNQ(fk) ... RNQ(FM) XNQ(fM) dans laquelle RNj(fk) est la composante résistive normalisée de l'impédance pour l'exemple d'étalonnage d'ordre i à la fréquence d'excitation d'ordre k, et XNi(fk) est la composante réactive normalisée de l'impédance pour l'exemple d'étalonnage d'ordre i à la fréquence d'excitation d'ordre k. c) établir la matrice de coefficients linéaires modélisés par la formule C = (#t x #)-1 x #t x Pe où P,est la matrice des valeurs d'étalonnage r1 #1 d1 r2 #2 d2 Pe -- ri #i di ... ... ... rQ #Q dQ où ri, #i et di sont respectivement la perméabilité magnétique relative, la conductivité de la cible et la distance de la cible par rapport au capteur dans l'exemple d'étalonnage d'ordre i. d) effectuer aux fréquences f1, f2, ... fk,... f une mesure d'impédance normalisée réelle du capteur électromagnétique en regard de la cible dont les paramètres sont à déterminer et établir le vecteur de mesure &gamma; = [1 RN(f1) XN(f1).... RN(fk) XN(fk)... RN(fM) XN(fM)]t e) calculer le vecteur des paramètres inconnus par la formule e = x C avec # = [ r, # , d]t où Zrv a et d sont les paramètres recherchés. D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de la description qui suit d'un exemple de réalisation particulier non limitatif de l'invention, en relation avec les figures ci-jointes parmi lesquelles: - la figure 1 est une vue en coupe schématique d'un capteur électromagnétique selon l'invention en regard d'une cible - la figure 2 est une vue en perspective illustrant la réalisation de la bobine dans le mode de réalisation particulier de l'invention - la figure 3 est une vue partielle agrandie en coupe selon la ligne III-III de la figure 2 d'un capteur électromagnétique selon le mode de réalisation de l'inven tion - la figure 4 donne les formules de modélisation du capteur selon l'invention. En référence à la figure 1, le capteur électromagnétique selon l'invention comporte un noyau 1 de forte perméabilité magnétique, par exemple un noyau en ferrite, en forme de cylindre plat circulaire dont une face porte une bobine annulaire 2 de forme circulaire fixé coaxialement au noyau 1. L'axe de symétrie de la bobine 2 s'étend perpendiculairement à la face du noyau sur laquelle la bobine est fixée. Sur la figure 1, le capteur électromagnétique ainsi formé est disposé en regard d'une cible métallique 3. Le noyau 1 a une perméabilité magnétique absolue yaet une conductivité q. L'air environnant a une perméabilité magnétique absolue 'La et une conductivité a et la cible a une perméabilité absolue 'Lc et une conductivité ae. On a noté par D le diamètre du noyau 1, L l'épaisseur du noyau 1, r1 le rayon interne de la bobine 2, r2 le rayon externe du bobinage 2, c la distance de la face du capteur électromagnétique tournée vers la cible par rapport à la cible, l la dimension du bobinage 2 prise selon son axe de symétrie, et d la distance de la cible par rapport à la face en regard de la bobine 2. On a également figuré en trait mixte 4 la représentation d'une ligne du champ magnétique créé par la bobine 2 lorsque celle-ci est soumise à un courant d'excitation alternatif. Selon l'invention l'épaisseur L du noyau est très supérieure à la dimension l de la bobine 2 prise selon son axe longitudinal de façon que les lignes de champs magnétiques 4 résultant de l'excitation de la bobine soient sensiblement dans leur totalité incluses dans le noyau, c'est-à-dire qu'elles se referment à l'intérieur du noyau sans sortir de la face du noyau opposée à la face portant la bobine 2. On notera à ce propos que les lignes de champs magnétique ont tendance à pénétrer d'autant plus profondément dans le noyau magnétique que la cible en regard a une perméabilité magnétique plus élevée. On pourra donc adapter l'épaisseur du noyau à la gamme de cible à laquelle le capteur électromagnétique est plus particulièrement destiné.En pratique, une épaisseur du noyau dix fois supérieure à la dimension l de la bobine s'est avérée satisfaisante quel que soit le type de cible. Dans la suite de l'exposé de l'invention, on supposera également que la cible a elle-même une épaisseur suffisante pour que les lignes de champ magnétique 4 se referment à l'intérieur de la cible. Par ailleurs, le diamètre du noyau est de préférence supérieur au diamètre externe de la bobine de façon que les lignes de champ magnétique se referment dans l'espace délimité par le noyau et la cible en regard, c'est-à-dire ne s'échappent pas par la paroi latérale du noyau 1. On rappellera à ce propos que les lignes de champs magnétique 4 se referment à une distance d'autant plus grande de l'axe de la bobine que la cible a une conducti vité plus élevée. En pratique, un diamètre du noyau égal au double du diamètre externe de la bobine s'est avéré suffisant même pour des cibles très conductrices telles que du laiton. Les figures 2 et 3 illustrent un mode de réalisation particulier de l'invention dans lequel la bobine 2 est constituée d'un enroulement d'excitation 5 et d'un enroulement de mesure 6 imbriqués l'un dans l'autre et chacun formé par une nappe de fils 7 disposés parallèlement les uns aux autres et réunis entre eux, par exemple par une résine. Après réalisation des nappes de fils, celles-ci sont disposées sur la tranche et enroulées en spirale comme illustré par la figure 2 où les nappes ont été représentées écartées les unes des autres pour une meilleure compréhension de la structure de la bobine 2. La figure 3 illustre de façon très agrandie la structure finale de la bobine 2 dans laquelle les différents tours de la bobine d'excitation 5 sont intercalés avec les tours de la bobine de mesure 6.A chacune des extrémités des enroulements, les fils 7 sont reliés en parallèle par une barrette 8 convenablement reliée à un fil de raccordement 9. Lorsque l'enroulement d'excitation 5 est soumis à un courant alternatif I, il apparaît aux bornes de l'enroulement de mesure une tension V qui est caractérisée par l'impédance complexe du capteur électromagnétique Z = V/I. Dans la suite de l'exposé, il sera fait référence à l'impédance normalisée ZN = RN + j X N avec RN = R-Ro et XN = X Xo où R et X sont respectivement la composante résistive et la composante réactive du capteur électromagnétique en présence de la cible, et Ro et X0 sont respectivement la composante résistive et la composante réactive du capteur électromagnétique en l'absence de cible. Selon l'invention on établit une matrice de plan d'expérience en déterminant l'impédance ZN à différentes fréquences d'excitation f1 f2.. fk...fM pour différents exemples d'étalonnage E1, E2, . E, EQ, chacun défini par une perméabilité magnétique relative Hi et une conductibilité ai de la cible, ainsi que par une distance di du capteur électromagnétique par rapport à la cible. On obtient ainsi la matrice de plan d'expérience 1 RN1(f1) XN1(f1) RN1(f2) XN1(f2) ... RN1(fk) XN1(fk) ... RN1(fM) XN1(fM) 1 RN2(f1) XN2(f1) RN2(f2) XN2(f2) ... RN2(fk) XN2(fk) ... RN2(fM) XN2(fM) # = #1 ...... ....... ....... ....... ....... ....... ....... ....... 1 RNi(f1) XNi(f1) XNi(f2) XNi(f2) ... RNi(fk) XNi(fk) ... RNi(fM) XNi(fM) 1 ...... ....... ....... ....... ....... ....... ....... ....... 1 RNQ(f1) XNQ(f1) RNQ(f2) XNQ(f2) ... RNQ(fk) XNQ(fk) ... RNQ(fM) XNQ(fM) Cette matrice de plan d'expérience est associée à une matrice des valeurs d'étalonnage Pe avec r1 #1 d1 r2 #2 d2 Pe = # ... ... ... ri #i di ... ... ... rQ #Q dQ où ri, #i et di sont respectivement la perméabilité magnétique relative, la conductivité de la cible et la distance de la cible par rapport au capteur dans l'exemple d'étalonnage d'ordre i. On notera que les exemples d'étalonnage pourraient être réalisés en utilisant des cibles réelles ayant une perméabilité magnétique relative et une conductivité bien définie, l'impédance normalisée étant alors obtenue par une mesure de la tension obtenue aux bornes de l'enroulement de mesure lors d'une alimentation de l'enroulement d'excitation à chacune des fréquences d'étalonnage. Toutefois, il est extrêmement difficile en pratique de réaliser des cibles ayant des paramètres voisins déterminés de façon précise, et selon un aspect de l'invention, on prévoit de réaliser la matrice de plan d'expérience à partir d'une modélisation de l'impédance du capteur électromagnétique selon l'invention, l'impédance modélisée associée au capteur selon l'invention étant reproduite à la figure 4 où r et a sont des variables d'intégration, N est le nombre de spires de l'enroulement d'excitation, #k = 2 # fk, a = 4 # x 10-7, J1 (&alpha;;r) est la fonction de Bessel de la première espèce d'ordre 2, et r1, r2, et l sont tels que définis précédemment et fL, ai, di, et ci sont pour l'exemple d'étalonnage d'ordre i les valeurs de , a, d et c telles que définies ci-dessus à propos de la description du capteur. A partir de l'impédance Zi(fj) ainsi modélisée on peut extraire la partie résistive Ri(fj) et la partie réactive Xi(fj) et en déduire les valeurs de RNi(fj) et XNi(fj) à insérer dans la matrice #. En pratique, il est souhaitable d'établir une matrice de plan d'expérience en encadrant chacune des valeurs probables des paramètres recherchés par des valeurs voisines donnant pour chaque paramètre une plage de détermination des paramètres. A partir d'une estimation initiale ro, #0, et d0 des paramètres recherchés, on établit la matrice de plan d'expérience en partant d'exemples d'étalonnage définis par différentes combinaisons des estimations initiales ro, #0 et d0, et de valeurs d'encadrement ro+# , ro-# ; #0+##, #0-## ;; d0+#d, d0-#d, où # , ## et Ad sont les écarts définissant la plage d'évolution estimée de chacun des paramètres. Le nombre total des différentes combinaisons des paramètres initiaux et des valeurs d'encadrement est égal à 33 = 27, ce qui conduirait donc à réaliser Q = 27 exemples d'étalonnage pour chacune des fréquences. Sachant que la matrice de plan d'expérience # a un nombre de lignes égal à Q et un nombre de colonnes à 2M + 1, on comprendra que la réalisation d'une matrice de plan d'expérience en utilisant toutes les combinaisons possibles des valeurs mentionnées ci-dessus conduirait à des temps de calcul extrêmement importants incompatibles avec l'utilisation pratique du capteur électromagnétique selon l'invention. Par ailleurs on sait qu'à haute fréquence la réponse du capteur est surtout sensible à la distance d et peu sensible à r et a tandis que pour des fréquences basses la réponse du capteur est sensible à r et a tout en étant faiblement sensible à une variation de la distance.En pratique on limitera donc le nombre d'exemples d'étalonnage en prenant différents exemples de perméabilité et de conductivité seulement en relation avec une fréquence basse du courant d'excitation, par exemple f1 = 1 kHz, et différents exemples de distance seulement en relation avec une fréquence élevée, par exemple f2 = 100 kHz.On obtient alors une matrice simplifiée des valeurs d'étalonnage ne comportant que 12 lignes r0 #0 d0 r0 #0+## d0 r0 #0-## d0 r0+# #0 d0 r0+# #0+## d0 Pe = # r0+# #0-## d0 r0-# #0 d0 r0-# #0+## d0 r0-# #0-## d0 zrO aO do r0 #0 d0+#d r0 #0 d0-#d et une matrice de plan d'expérience correspondante ne comportant que trois colonnes 1 RN1(f1) XN1(f1) 1 RN2(fl) XNt(fi) 1 RN3(fi) XN3(fl) 1 RN4(fl) XN4(fl) 1 RN5(f1) XN5(f1) 1 = RN6(f1) XN6(fl) 1 RN7(fi) XN7(fi) 1 RN8(f1) XN8(f1) 1 RN9(f1) XN9(f1) 1 RN10(f2) XN10(f2) 1 RN11(f2) XN11(f2) 1 RN12(f2) XN12(f2) Quoi qu'il en soit, à partir de la matrice de valeurs d'étalonnage et la matrice de plan d'expérience, on établit une matrice C de coefficients linéaires modélisés: c = (#t x #)-1 x 'Pt x puis on effectue une mesure de l'impédance réalisée réelle du capteur électromagnétique en regard de la cible dont les paramètres sont à déterminer et on établit le vecteur de mesure &gamma; = [1 RN(f1) XN(f1).... RN(fk) XN(fk)... RN(fM) XN(fM)]t RN et XN étant déterminés de façon connue en soi à partir de la mesure de l'impédance. On calcule alors le vecteur des paramètres à déterminer par la formule # = &gamma;t x C avec # = [ r, # , d]t Lorsque l'estimation initiale r0, #0 et do des paramètres à déterminer est très proche de la valeur réelle de ces paramètres, et que la tolérance sur la détermination de ces paramètres est suffisante par rapport aux performances de la modélisation utilisée, on peut considérer que le vecteur e ainsi trouvé est satisfaisant et cesser les calculs.Si l'on souhaite affiner la modélisation et obtenir une valeur des paramètres encore plus proche de la réalité, on considère alors que la valeur obtenue de e est une première estimation réelle, et l'on note cette estimation e (t = a(') d(') | ' puis on établit à nouveau la matrice des valeurs d'étalonnage et la matrice de plan d'expérience en utilisant ces paramètres de façon analogue à celle dont les paramètres initiaux r0, #0 et do avaient été utilisés.Après réalisation d'un nouveau calcul selon l'algorithme décrit ci-dessus, on obtient une nouvelle estimation réelle 0(2) que l'on peut comparer à l'estimation pour vérifier s'il y a convergence et si l'écart par rapport à l'estimation précédente justifie un nouveau calcul ou non. En pratique, il a été constaté que même pour des valeurs des paramètres initiaux très éloignés des valeurs et paramètres réels une convergence était obtenue au bout de trois ou quatre itérations. Bien entendu l'invention n'est pas limitée au mode de réalisation décrit et on peut y apporter des variantes de réalisation sans sortir du cadre de l'invention. En particulier, bien que le capteur selon l'invention ait été illustré avec une bobine 2 comportant deux enroulements imbriqués l'un dans l'autre, on peut prévoir d'utiliser un seul enroulement servant à la fois à l'excitation et à la mesure, l'excitation étant assurée par un générateur de courant, tandis que la mesure est effectuée par une mesure de la tension aux bornes de l'enroulement unique. De même, bien qu'un mode de réalisation ait été décrit qui a pour avantage d'assurer une grande linéarisation en raison du fait que chaque spire est contenue dans un plan, il est possible de réaliser la bobine par un enroulement hélicoïdal conventionnel en s'assurant toutefois que le pas de l'enroulement reste très faible pour éviter des singularités qui risqueraient de perturber la modélisation utilisée. De même, à propos du procédé selon l'invention, on peut utiliser toute modélisation qui sera jugée appropriée en fonction de la structure précise du capteur utilisé et assurer une simplification de la modélisation lorsque certains des paramètres sont constants d'une cible à l'autre. Lorsque RN et XN ne varient plus linéairement en fonction des paramètres influents, on peut également prévoir des changements de variables pour linéariser la relation entrée-sortie du capteur, par exemple en rempla çant RN et XN dans la matrice de plan d'expérience par des valeurs In(RN) et ln Bien qu'un exemple de modélisation simplifié ait été donné avec des valeurs d'encadrement symétriques par rapport à l'estimation initiale, on peut choisir des valeurs d' encadrement différentes. REVENDICATIONS 1. Capteur électromagnétique comportant au moins une bobine (2) associée à un noyau (1) à forte perméabilité magnétique, caractérisé en ce que le noyau (1) est une plaque plane, en ce que la bobine (2) est fixée sur une face du noyau et comporte un axe de symétrie s'étendant perpendiculairement à la face du noyau sur laquelle elle est fixée, et en ce que le noyau (1) a une épaisseur (L) supérieure à la dimension (1) de la bobine prise selon son axe de symétrie de façon que les lignes de champs magnétiques (4) résultant de l'excitation de la bobine soient sensiblement dans leur totalité incluses dans le noyau. 2. Capteur électromagnétique selon la revendication 1, caractérisé en ce que le noyau (1) et la bobine (2) sont des cylindres circulaires coaxiaux, le noyau (1) ayant un diamètre (D) supérieur au diamètre externe de la bobine de façon que les lignes de champs magnétiques (4) se referment sensiblement en totalité dans l'espace délimité par le noyau et une cible (3) en regard. 3. Capteur électromagnétique selon la revendication 1, caractérisé en ce que la bobine (2) est formée par une série de fils (7) reliés en parallèle selon une nappe qui est enroulée en spirale. 4. Capteur électromagnétique selon la revendication 1, caractérisé en ce que la bobine comporte un enroulement d'excitation (5) et un enroulement de mesure (6) séparés. 5. Capteur électromagnétique selon la revendication 4, caractérisé en ce que l'enroulement d'excitation (5) et l'enroulement de mesure (6) sont imbriqués l'un dans l'autre. 6. Procédé de détermination simultanée de la perméabilité magnétique et la conductivité d'une cible métallique ainsi que la distance de cette cible par rapport à un capteur électromagnétique selon la revendication 1, caractérisé en ce qu'il comprend les étapes de a) déterminer l'impédance normalisée complexe, ZN = RN + j X N du capteur électromagnétique à différentes fréquences d'excitation (f1, ... fk, ... fM), et pour différents exemples d'étalonnage (E1 ... Ek, ...EQ) chacun défini par une perméabilité magnétique relative ri et une conductibilité (ai) de la cible ainsi que par une distance (di) du capteur électromagnétique par rapport à la cible, avec RN = R-Ro et XN = X X0 X0 où R et X sont respectivement la composante résistive et la composante réactive du capteur électromagnétique en présence de la cible, et R0 et X0 sont respectivement la composante résistive et la composante réactive du capteur électromagnétique en l'absence de cible. b) établir la matrice d'un plan d'expérience 1 RN1(f1) XN1(f1) RN1(f2) XN1(f2) ... RN1(fk) XN1(fk) ... RN1(fM) XN1(fM) 1 RN2(f1) XN2(f1) RN2(f2) XN2(f2) ... RN2(fk) XN2(fk) ... RN2(fM) XN2(fM) t= 1 1 RNi(f1) XNi(f1) RNi(f2) XNi(f2) ... RNi(fk) XNi(fk) ... RNi(fM) XNi(fM) 1 ...... ....... ....... ....... ....... ....... ....... ....... 1 RNQ(f1) XNQ(f1) RNQ(f2) XNQ(f2) ... RNQ(fk) XNQ(fk) ... RNQ(fM) XNQ(fM) dans laquelle RNi(fk) est la composante résistive normalisée de l'impédance pour l'exemple d'étalonnage d'ordre i a la fréquence d'excitation d'ordre k, et XNj(fk) est la composante réactive normalisée de l'impédance pour l'exemple d'étalonnage d'ordre i à la fréquence d'excitation d'ordre k. c) établir la matrice de coefficient linéaire modélisé par la formule C = (#t x #)-1 x #t x Pe où P. est la matrice des valeurs d'étalonnage r1 #1 d1 zr2 a2 d2 Pe = ri #i di ... ... ... rQ #Q dQ d) effectuer aux différentes fréquences (f1, ... fk, ... fM) une mesure d'impédance normalisée réelle du capteur électromagnétique en regard de la cible dont les paramètres sont à déterminer et établir le vecteur de mesure &gamma; = [1 RN(f1) XN(f1).... RN(fk) XN(fk)... RN(fM) XN(fM)]t e) calculer le vecteur des paramètres inconnus par la formule e = c avec # = [ r, # , d]t où r, a et d sont les paramètres recherchés. 7. Procédé selon la revendication 6, caractérisé en ce que la matrice de plan d'expérience est établie en modélisant chaque exemple d'étalonnage par la formule reproduite à la figure 4. 8. Procédé selon la revendication 7, caractérisé en ce qu'on réitère le calcul des paramètres recherchés en réactualisant la matrice du plan d'expérience et la matrice des valeurs d'étalonnage par une utilisation de la valeur des paramètres trouvés à l'itération précédente. 9. Procédé selon la revendication 6, caractérisé en ce qu'on utilise une matrice simplifiée des valeurs d'étalonnage r0 #0 d0 r0 #0+## d0 r0 #0-## d0 r0+# #0 d0 r0+# #0+## d0 Pe = # r0+# #0-## d0 r0-# #0 d0 r0-# #0+## d0 r0-# #0-## d0 r0 #0 d0 r0 #0 d0+#d r0 #0 d0-#d associé à une matrice de plan d'expérience simplifiée 1 RN1(f1) XN1(f1) 1 RN2(f1) XN2(f1) 1 RN3(f1) XN3(f1) 1 RN4(fi) XN4(fl) 1 RN5(f1) XN5(f1) # = # 1 RN6(f1) XN6(f1) 1 RN7(fi) XN7(fi) 1 RN8(f1) XN8(f1) 1 RN9(f1) XN9(f1) 1 RN10(f2) XN10(f2) 1 RN11(f2) XN11(f2) 1 RN12(f2) XN12(f2) où f1 est une fréquence basse et f2 une fréquence élevée. 10. Procédé selon la revendication 9, caractérisé en ce que f1 est environ égale à 1 kHz et f2 est environ égale à 100 kHz.