Brevet FR2692700A1 — publié le 19931224

"DISPOSITIF DE TRAITEMENT D'UN SIGNAL MESURE CORRESPONDANT A L'INTENSITE DE RAYONS X REFLECHIE PAR UNE STRUCTURE DE COUCHES MULTIPLES SUR UN SUBSTRAT" Description: L'invention concerne un dispositif de traitement d'un signal mesuré correspondant à l'intensité réfléchie, en fonction de l'angle de rasance (E), à une longueur d'onde ( 4 x) de rayons X, par une structure réelle de couches multiples empilées sur un substrat, ces couches présentant à leurs interfaces des différences faibles d'indice de réfraction. Par angle de rasance, il faut entendre l'angle que fait le faisceau incident avec la surface réfléchissante, c'est-à-dire l'angle complémentaire de l'angle d'incidence. L'invention trouve son application dans la réalisation de dispositif de tests de structures multicouches; par exemple, l'invention peut être utilisée pour la détermination automatique de l'épaisseur et de la nature ou composition des couches de miroirs multicouches périodiques, ou bien d'hétérostructures de matériaux semiconducteurs empilés. Il est déjà connu de l'état de la technique, par le brevet EP O 468 582, un appareil de diffractométrie à rayons X comprenant des moyens optiques pour mesurer la réflectivité d'un échantillon solide. Cet appareil trouve son application dans la détermination des paramètres structuraux de miroirs multicouches travaillants dans le domaine des rayons X, formés de l'empilement de bicouches de matériaux de fort et faible indices de réfraction Cet appareil permet d'effectuer une comparaison entre des paramètres calculés à partir de la courbe de réflectivité et des paramètres connus par des Tables de classement. Le document cité divulgue un système de commande pour sélectionner des filtres optiques appropriés parmi des filtres ayant différents coefficients d'absorption, et pour déclencher leur mise en place sur le trajet du faisceau de rayons X réfléchis par une structure à tester Grâce à ce système, la mesure de l'intensité réfléchie peut être effectuée au moyen d'un compteur proportionnel, de manière telle que ce dernier ne sorte jamais de sa région de fonctionnement linéaire quelles que soient les valeurs de l'intensité Par compteur proportionnel, il faut entendre un dispositif comportant des gaz ionisables par le flux de photons à détecter, et fournissant un signal sous la forme d'une tension proportionnelle au nombre de photons reçus. L'appareil connu comprend un système de report du signal de sortie du compteur proportionnel sur un dispositif de traitement de ce signal dans lequel un programme d'ordinateur permet de raccorder entre elles les différentes parties de la courbe de réflectivité obtenues en relation avec chacun des filtres du système commandé On obtient ainsi une courbe de réflectivité lissée. Cette courbe de réflectivité, établie à longueur d'onde x constante d'une source de rayon X, et en fonction des variations de l'angle de rasance O, est fortement décroissante entre la valeur dite critique Oc o la réflectivité est égale à 1, et des valeurs o l'angle de rasance est encore faible, c'est-à-dire inférieur à, ou de l'ordre de 4 (degrés) C'est pourquoi un système commandé de filtres est utilisé afin de pouvoir effectuer toutes les mesures d'intensité sans sortir du domaine de linéarité du compteur proportionnel Cette courbe montre en outre un certain nombre de maxima d'intensité réfléchie assez marqués qui apparaissent pour des angles de rasance spécifiques. Les valeurs des angles de rasance O correspondant aux maxima vérifient la loi de Bragg qui s'énonce: 2 d sin O = k X (avec k = constante = un nombre entier) et permet le calcul de la période D des couches d'un système multicouches périodique. Cependant cette détermination correspondant aux ordres de la réflexion doit aussi prendre en compte la nature des matériaux des couches C'est pourquoi, selon l'enseignement du document cité, le pas D des couches périodiques de miroirs multicouches est déterminé, à partir de la courbe de réflectivité lissée, par comparaison avec des courbes théoriques construites à partir de Tables de classement. Néanmoins, ce dispositif n'est applicable que dans les cas o la courbe de réflectivité montre des maxima lisibles, c'est-à-dire suffisamment marqués pour permettre une mesure précise de leur position en valeur d'angle de rasance Ceci limite cette méthode aux structure multicouches formées d'empilements périodiques de couches dont les compositions des matériaux sont connues et dont les indices de réfraction sont très différents, ainsi qu'aux miroirs multicouches de structure connue n'ayant pas une épaisseur trop importante, et dans lesquels la période n'est pas trop petite. C'est pourquoi un des buts de la présente invention est de fournir un dispositif pour la détermination automatique de l'épaisseur et de la composition des couches de structure multicouches, dans les cas particulièrement difficiles o: les épaisseurs des couches ne sont pas périodiques mais différentes dans toute la structure; les matériaux de la structure ne se retrouvent pas non plus de manière périodique; de plus, les indices de réfraction des matériaux sont très proches d'une couche à l'autre; et en outre certaines couches peuvent être extrêmement fines, par exemple moins d'l nm. Ce but est atteint au moyen du dispositif de traitement défini plus haut comprenant des moyens pour effectuer la transforma- tion de la fonction (f(E)) correspondant au signal d'intensité réflé- chie mesuré en fonction de l'angle de rasance (e) dans un domaine d'angles de rasance limité par des valeurs d'angles minimales et maximales emin max' en une fonction (F(d)) correspondant au signal d'intensité formulé en fonction de la profondeur (d) dans la structure réelle de couches multiples, cette fonction (F(d)) obtenue permettant d'établir une cartographie directe des profondeurs de chaque interfa- ces (di, d 2, dn) de la structure réelle, et des intensités réflé- chies associées. L'invention est décrite ci-après en détail en référence avec les figures annexées dont: la figure la représente une structure de couches réalisée sur un substrat, et testée au moyen du dispositif selon l'invention; la figure lb représente l'intensité réfléchie f(O) en fonction de l'angle de rasance 0, pour la structure de la figure la; la figure lc repré- sente la fonction g(E) obtenue par le filtrage mathématique de la fonction f(e) représentée sur la figure lb; la figure 2 a représente une structure théorique à une couche sur un substrat; la figure 2 b représente l'intensité réfléchie t 1 ( 0) par la structure théorique de la figure 2 a, calculée en fonction de l'angle de rasance e; la figure 2 c représente la fonction h 1 ( O > obtenue après un filtrage de la fonction tl(E) de la figure 2 b au moyen d'un filtre identique à celui qui est utilisé pour obtenir la courbe de la figure lc; la figure 3 a représente une structure théorique à une couche sur un substrat, en outre munie d'une couche supérieure d'oxyde; la figure 3 b représente l'intensité réfléchie t 2 (e) par la structure théorique de la figure 3 a, calculée en fonction de l'angle de rasance E; la figure 3 c représente la fonction h 2 (O) obtenue après un filtrage de la fonction t 2 (e) de la figure 3 b au moyen d'un filtre identique à celui qui est utilisé pour obtenir la courbe de la figure lc: la figure 3 d représente la variation de la longueur d'onde If(O) calculée à partir de la fonction filtrée de la figure 3 c; la figure 3 e représente la fonction CHIRP(E); la figure 4 a représente un algorithme simplifié mis en oeuvre selon l'invention, et la figure 4 b un algorithme pour une mise en oeuvre plus précise; la figure 5 représente schématiquement la fonction du résultat G(d) de l'intégration sans lissage; la figure 6 représente-l'intensité réfléchie F(d) en fonction de la distance d mesurée dans la structure multicouches de la figure la, à partir de sa surface supérieure, cette intensité étant obtenue par le traitement des données d'intensité de la figure lb au moyen du dispo- sitif selon l'invention; la figure 7 représente un diffractomètre à rayons X. En référence avec la figure la et les figures 4, l'inven- tion se propose de fournir un dispositif de traitement d'un signal mesuré correspondant à l'intensité réfléchie I, à une longueur d'onde CRX de rayons X, par une structure réelle de couches multiples empi- lées sur un substrat 100, ces couches présentant à leurs interfaces des différences faibles d'indice de réfraction, ce dispositif de traitement comprenant des moyens pour effectuer la transformation de la fonction f(E) correspondant au signal d'intensité réfléchie mesuré en fonction de l'angle de rasance ( 5 dans un domaine d'angles de rasance limité par des valeurs d'angles minimales et maximales O min max' en une fonction F(d) correspondant au signal d'intensité formulé en fonction de la profondeur d dans la structure réelle de couches multiples, cette fonction F(d) obtenue permettant d'établir une cartographie directe des profondeurs de chaque interfaces dl, d 2, d, de la structure réelle, et des intensités réfléchies associées. Le domaine Emin Omax sera défini ultérieurement. En référence avec la figure 7, l'homme du métier peut utiliser, pour réaliser la mesure de l'intensité réfléchie par un échantillon à tester, en fonction de l'angle de rasance G, un appareil de diffractométrie à rayons X. Les appareils connus, tels que par exemple commercialisés par Philips sous la référence PW 1050, comprennent au moins un système source 1, 2, 3 comportant une source linéaire 1 à rayons X un système collimateur 2 et 3, incluant des fentes de Sollers 2 et une fente de divergence 3; un porte-échantillons goniométrique 9 comprenant des moyens d'orientation (non représentés) pour un échantillon 10: ces moyens comprennent en particulier un système pour sélectionner et enregistrer la valeur de l'angle de rasance E que fait le faisceau provenant du système source 1, 2, 3 vis-à-vis du plan supérieur 4 de l'échantillon 10; et un système capteur 120 du flux réfléchi comportant un système collimateur 5 formé d'une fente dite de réception; un monochromateur 7 dans le domaine des rayons X, par exemple un monochromateur graphite; un détecteur 8 du flux de rayons X, par exemple un compteur proportionnel; éventuellement un système de filtre 62 disposés dans un porte-filtres 6. Tout appareil capable d'effectuer la mesure de l'intensité I = f(G) d'un faisceau rasant de rayons X réfléchi par un échantillon à tester, en fonction de l'angle de rasance O, peut être utilisé au lieu de l'appareil décrit plus haut. Dans le cas o l'appareil de la figure 7 est utilisé, le compteur proportionnel 8 fournit un signal proportionnel au flux des photons reçus et qui est reporté, avec la valeur des angles G correspondant à chaque valeur proportionnelle à l'intensité I, sur un moyen de stockage 130 de ces données Le dispositif connu comprend éventuellement un traceur de courbe 131, qui, à partir des données issues du moyen de stockage 130, est capable de tracer une courbe continue représentant l'intensité f(e) portée en ordonnée en unités arbitraires (UA) comprises entre Q et 1 en fonction de l'angle de rasance e, porté en abscisse en degrés. En référence avec la figure la, on se propose de tester au moyen du dispositif selon l'invention, une structure multicouches 10 a. Cette structure i Qa comprend, à partir de sa surface supérieure 4 une couche d'oxyde 11 b; des couches de matériau semiconducteur d'arséniure de gallium (Ga As) 12, 14, 16, 18; des couches d'un autre matériau semiconducteur d'arséniure de gallium et d'aluminium (Gax All_x As) 13, 15, 17, 19, formant des hétérostructures avec les couches de Ga As alternées; un substrat en arséniure de gallium (Ga As) 100 Dans cet échantillon de test: les distances dl, d 2, d 3,, d,, entre la surface supérieure 4 et chacune des interfaces Ga As/Ga Al As sont inconnues, car les épaisseurs exactes de dépôt par croissance épitaxiale de chacune des couches sont inconnues; la composition exacte des couches de Ga Al As est inconnue, c'est-à-dire que la concentration x est inconnue; l'indice de réfraction du matériau Ga Al As est très proche de celui du Ga As. La figure lb montre la courbe expérimentale de réflectivité de l'échantillon i Qa à tester représentant les valeurs mesurées f(e) de l'intensité réfléchie I des rayons X, en fonction de l'angle de rasance e. On a vu précédemment que dans le diffractomètre de la figure 7, la source de rayons X émet un faisceau de rayons incidents qui font un angle e avec la surface 4 de l'échantillon 10 à tester, et qui sont réfléchis sous le même angle e par les différentes surfaces et interfaces de cet échantillon Chacune des interfaces en commençant poar l'interface air-première couche de matériau, provoque des réflexions Les rayons ainsi réfléchis sont en interférences optiques qui sont constructives pour certaines valeurs de l'angle de rasance E, et qui sont destructives pour d'autre valeurs Il se trouve que dans le cas o la première interface est celle de l'air-semiconducteur, la première réflexion est d'intensité beaucoup plus grande que les autres Il en résulte que cet effet des interférences tantôt constructives tantôt destructives produit une surmodulation qui se superpose à une courbe générale décroissante de l'intensité I en fonction de e, telle que représentée sur la figure lb Si la structure consistait en une seule interface air-matériau, la courbe de réflectivité serait seulement descendante et lisse La surmodulation est un effet caractéristique des structures multicouches. On constate qu'à partir de cette mesure directe de l'intensité réfléchie par un échantillon multicouches, on est absolument incapable de réaliser la mesure qui donne la position dl, d 2, d des interfaces par rapport à la surface supérieure de l'échantillon, ou bien la mesure d'épaisseur d'une couche dans la structure multicouches. -L'invention a pour but de fournir un dispositif pour déterminer directement les positions dl, d 2, d 3, d des interfaces par rapport à la surface supérieure de l'échantillon à partir de l'intensité réfléchie I = f(E) du faisceau de rayons X mesurée en fonction de E. Ce problème est résolu par un dispositif qui inclut des moyens pour réaliser la Transformation de la fonction f(E) fournissant la courbé de la figure lb, en une fonction F(d) de l'intensité en fonction de la profondeur d fournissant la courbe de la figure 6. A partir des seules données connues au départ qui sont la valeur de l'intensité I réfléchie par la structure réelle en fonction de e; la composition approximative des matériaux des couches, et leur nombre; l'indice de réfraction moyen de la structure l'indice de réfraction du substrat; on va réaliser une Transformation de la fonction représentant l'intensité I en fonction de e, en une autre fonction représentant l'intensité I en fonction de la profondeur d à l'intérieur de la structure Ci-après on se referera à cette transformation par le terme "la Transformation" ou le terme "ladite Transformation". La Transformation proposée permet d'arriver à la détermination des profondeurs d'interfaces dl, d 2, dn avec une précision de 0,5 nm même si: les indices de réfraction des diverses couches sont peu différents; les couches ont des épaisseurs très différentes. les couches sont très minces par exemple monoatomique. La courbe d'intensité réfléchie de la figure lb a été réalisée en faisant varier l'angle de rasance e entre O et 40 En effet, pour apprécier des distances d dans le matériau à tester, qui sont de l'ordre de 1 à quelques centaines de nm, il faut que l'angle e soit petit afin que d et e vérifient la loi de Bragg. Par e petit, on entend que l'angle de rasance reste compris entre ec, angle critique proche de O o l'intensité réfléchie I est égale à 1, et t qui est inférieur à, ou de l'ordre de 4 . L'angle E? est celui pour lequel on considère que l'on a recueilli suffisamment d'informations pour identifier toutes les interfaces de la structure. Dans l'exemple décrit ci-après la longueur d'onde de la source de rayons X, 4 x = 4,5 nm et e est mesuré entre O et 2 . Le domaine emin Ema X sur lequel on effectue la Transformation sera défini ultérieurement à l'intérieur de l'intervalle E O OM- La courbe de la figure 6 montre de nombreux pics Chaque pic ou maximum caractérise une interface de la structure de la figure la En coordonnée horizontale est portée la profondeur d dans la structure l Oa et en coordonnée verticale est portée l'intensité I en unités arbitraires, entre O et 1. La position de chaque pic caractérise la profondeur d'une interface dans la structure i Qa par rapport à la surface 4. La hauteur de chaque pic est influencée par de nombreux facteurs En principe, la hauteur des pics dépend essentiellement de la différence des indices de réfraction des matériaux de chaque côtés de l'interface. Mais la hauteur des pics est également influencée par la raideur de l'interface la rugosité de l'interface la résolution en angle du système de mesure l'un des moyens du présent dispositif formé par un système de filtrage. Néanmoins, le dispositif selon l'invention permet d'arriver à une fonction représentée par la courbe de la figure 6, de l'intensité en fonction de d, (F(d)), qui fournit les profondeur d et les intensités associées de manière telle que les performances et précisions citées plus haut sont atteintes. Le dispositif qui permet de réaliser la Transformation du signal représenté sur la figure lb pour arriver à celui qui est représenté sur la figure 6 offre l'avantage de fournir cette image I fonction de d directement en un seul passage. Un autre avantage réside dans le fait que, lors de la détection des interfaces par la position des pics, respectivement en dl, d 2, d 3 etc, sur la courbe de la figure 6, les pics parasites qui apparaissent aux positions 2 d 1, 3 d 1, etc 2 d 2, 3 d 2 etc, 2 d 3, 3 d 3 etc sont de tailles très réduites devant celles des pics principaux qui apparaissent en dl, d 2, d 3 De même que, si l'on détecte deux interfaces à des profondeurs dl et d 2 respectivement, les pics parasites qui apparaissent à dl d 2 et à dl + d 2 sont également petits par rapport aux pics principaux. Donc le dispositif selon l'invention a permis de rendre minimes les contributions des pics correspondant aux multiples, aux sommes et aux différences des profondeurs d, et par conséquent de rendre la courbe de la figure 6 exploitable pour arriver aux buts de l'invention. Il se trouve que, lorsqu'une structure telle que celle de la figure la montre de nombreuses interfaces, la surmodulation qui apparaît sur la courbe de réflectivité de la figure lb est engendrée par la somme de plusieurs fonctions périodiques. On pourrait alors penser à appliquer une simple intégrale de Fourier pour obtenir des valeurs importantes d'une intégrale pour chaque valeur relative à chacun des sinus de la somme, dans le but d'obtenir une courbe telle que celle de la figure 6. Un problème apparaît aussitôt du fait que la longueur d'onde de la surmodulation montrée sur la courbe de la figure lb n'est pas constante Elle varie beaucoup, et notamment dans la région des petites valeurs de l'angle e. Donc, si l'on envisageait d'utiliser une intégrale de Fourier pour obtenir une courbe ayant un pic bien défini, relatif à un sinus donné, et si on ne tenait pas compte de cette variation de longueur d'onde, alors on obtiendrait un pic très large et qui serait inexploitable pour la détermination de la profondeur d'interface. L'expérience a montré que, si une structure à plusieurs couches présentait peu d'interfaces peu profondes, on pourrait appliquer l'intégrale de Fourier, quoique sans beaucoup de précision dans le résultat obtenu. Mais dès que la structure comprend des couches ayant une grande profondeur d'interface, ou bien comprend des interfaces très rapprochées, on n'a aucune chance de détecter ces interfaces en appliquant une intégrale de Fourier simple. Une structure, dont l'analyse en matière de profondeur d'interfaces se révèle difficile, est par exemple celle d'un transistor à haute mobilité électronique (en anglais: HEMT) pseudomorphique qui présente trois interfaces principales, dont la dernière est très éloignée de l'interface air-première couche de matériau semiconducteur. Le but de l'invention est de fournir un dispositif capable de traiter ces cas difficiles Pour tenir compte du fait que la fonction f(G) correspondant au signal d'intensité réfléchie en fonction de l'angle de rasance e présente une surmodulation qui est la somme approximativement de fonctions périodiques à périodes variables de façon monotone en fonction de l'angle de rasance e, selon l'invention, les moyens de Transformation de la fonction mesurée f(e) du signal d'intensité en fonction de l'angle de rasance O, en la fonction F(d) du signal d'intensité en fonction de la profondeur d dans la structure comprennent des moyens d'intégration, sur le domaine il 0 min O émax de la fonction mesurée f(E) multipliée par une fonction Ki(d, O) périodique de la profondeur d et de l'angle O, dont la période varie en fonction de l'angle e avec la même variation monotone que celle des fonctions périodiques qui sont les composantes de la fonction mesurée f(O). En particulier la fonction périodique Ki(d,O) à variation monotone de période peut être une fonction pseudosinusoïdale dont la longueur d'onde X varie de façon monotone en fonction de O, dans un intervalle donné O min o - Le dispositif selon l'invention comprend des moyens pour synthétiser une telle fonction périodique Ki(d,e) de la profondeur d et de l'angle O, dont la période varie de façon monotone en fonction de E. La figure 2 a représente schématiquement une structure l Ob théorique bâtie à partir des données approximatives seules connues relativement à la structure réelle l Oa. La structure théorique de la figure 2 a comporte deux couches: une couche inférieure M 2 d'épaisseur telle que d 2 = pour symboliser l'épaisseur d'un substrat, et d'indice de réfraction N 2 égal à l'indice du substrat 100 de la structure réelle à tester l Oa; une couche supérieure M 1 d'épaisseur choisie arbitrairement telle que dl = 100 à 200 nm, et d'indice de réfraction ni égal à l'indice de réfraction moyen ou approximatif de la structure réelle 10 a. Dans ce modèle théorique, on considère que des rayons incidents 50 sont réfléchis d'une part par la surface supérieure 4 a en 51 et d'autre part par l'interface 4 b des couches M,/M 2 en 52. Dans l'exemple décrit ci-après, on a choisi pour M, di = 100 mn et Partie réelle de ( 1-n>) = 13 ppm pour M 2, d 2 = o et Partie réelle de ( 1-n 2) = 14,5 ppm. Dans le courant de l'expérimentation à l'aide du dispositif selon l'invention, il est apparu que l'utilisation du modèle l Ob théorique de la figure 2 a donnait de bons résultats, mais que ces résultats pouvaient encore être améliorés par l'utilisation d'un autre modèle théorique l Oc représenté sur la figure 3 a. Ce modèle théorique l Oc qui est utilisé de préférence, comporte trois couches et est également bâti à partir des données approximatives relatives à la structure réelle 10 a. Ces couches sont: une couche inférieure M 2 d'épaisseur telle que d 2 = pour symboliser le substrat, et d'indice de réfraction N 2 égal à l'indice du substrat 100 de la structure réelle i Qa; une couche intermédiaire M 1 d'épaisseur choisie arbitrairement telle que dl = 100 à 200 nm, et d'indice de réfraction ni de l'ordre de grandeur de l'indice moyen de la structure 10 a; une couche supérieure Mo d'oxyde, d'épaisseur choisie arbitrairement telle que do 1 nm, et d'indice de réfraction no présentant une petite différence avec celui de la couche M 1 sous- jacente. Dans ce modèle théorique l Oc, on considère que des rayons incident 50 sont réfléchis sous l'angle de rasance E, selon 53 par la surface supérieure 4 c de la structure l Oc, selon 51 par l'interface 4 a des couches M /Ml, et selon 52 par l'interface 4 b des couches M 1/M 2. Dans ce modèle théorique l Oc, on entend maintenant par réflexion principale la réflexion qui est la somme de la réflexion à la surface de la couche d'oxyde M, et de la réflexion à l'interface entre l'oxyde et la première couche de matériau M 1 En effet, ces deux surfaces ont des coefficients de réflexion sensiblement identiques. Dans l'exemple décrit ci-après on a choisi pour Mo, do = 1 in Partie réelle de ( 1-n) = 9 ppm pour M 1, dl = 100 nm et Partie réelle de ( 1-n 1) = 13 ppm pour M 2, d 2 = X et Partie réelle de ( 1-n 2) = 14,5 ppm. A l'aide du dispositif de l'invention, un calcul théorique de la réflectivité en fonction de e est alors effectué, pour l'une des structures l Ob ou l Oc consistant essentiellement en un substrat d'indice donné, surmonté d'une couche d'épaisseur donnée et d'indice de réfraction très proche de celui du substrat. Ce calcul théorique peut être effectué par des moyens simples connues de l'homme du métier, par exemple programmé avec la fonction donnant la réflectivité Puis les résultats sont stockés dans des moyens de stockage 140. Il en résulte respectivement pour les structures l Ob et l Oc, les courbes tj(e) et t 2 ( 0) des figures 2 b et 3 b Ces courbes représentent la réflectivité théorique qui décroît à partir de éc, l'angle critique, lorsque la valeur de e augmente Comme précédemment, en référence avec ce qui a été dit pour la figure lb, les courbes des figures 2 b et 3 b montrent une surmodulation. Dans une variante du dispositif, on peut remplacer le calcul théorique de l'intensité I en fonction de E, par le calcul théorique de Log I en fonction de e. Selon l'invention, pour rendre exploitable les informations relatives à la surmodulation des courbes des figures lb et 2 b ou 3 b, on soumet la fonction f (E) donnant la réflectivité mesurée en fonction de e, ainsi que les fonctions tl(E) et t 2 (E) des réflectivités théoriques à un filtre mathématique F 1 Le filtre F 1 est un filtre passe-bande de type passe-haut (en anglais HIGH-PASS FILTER) qui laisse passer les courtes longueurs d'onde et ne laissent pas passer les grandes longueurs d'onde de la courbe de réflectivité. Le succès des moyens préconisés tient beaucoup aux caractéristiques de filtre F 1 Ce filtre ne doit pas décaler le signal en 0, et la longueur d'onde de coupure doit être soigneusement choisie. Après le filtrage de la fonction f(E) des valeurs mesurées en relation avec la structure de test 10 a, on obtient la fonction g(O) représentée sur la figure lc qui ne montre plus que la surmodulation. Il résulte de ce qu'on a dit précédemment, que ce signal représenté par la fonction g(e) sur la figure lc est un signal composé de différents signaux de base ayant des fréquences (ou longueur d'onde k(E)) différentes s'exprimant sous la forme de somme de sinus, du fait que chaque signal est issu de l'interférence entre le signal correspondant à la réflexionprincipale sur la surface supérieure de la structure et le signal correspondant à la réflexion sur une interface à l'intérieur de la structure. Après le filtrage des fonctions théoriques tl(E") ou t 2 (), on obtient les fonctions théoriques hl(E) ou h 2 (E) respectivement représentées sur les figures 2 c et 3 c qui ne montrent plus que la surmodulation. L'allure de la courbe lc correspondant à la structure réelle est quasipériodique, alors que l'allure des courbes théoriques h 1 (G), h 2 ( 0) de la figure 2 c ou 3 c est maintenant celle d'une fonction quasisinusoïdale dont la longueur d'onde A( 8) varie en fonction de O. Dans tous les cas ces courbes théoriques montrent des passages à zéro. Il est alors possible de définir (ou calculer) cette longueur d'onde A(@) en calculant les distances entre les zéros. On a fait en outre une découverte importante en ce qui concerne cette fonction A(O): elle ne dépend absolument pas de la valeur de la profondeur d Donc elle est identique pour toutes les interfaces, quelle que soit la distance d de l'interface avec la surface supérieure de la structure considérée Tout au moins, la variation de k(O) est identique pour les sytèmes de matériaux semiconducteurs. Donc, il suffit de connaître les paramètres approchés des matériaux de la structure à étudier, par exemple une valeur approchée de son indice de réfraction par rapport à l'air, pour être en mesure de déterminer la variation de la longueur d'onde X(@) de la courbe de réflectivité quasi- sinusoïdale filtrée h,(O) ou h 2 (O) obtenue à partir d'une structure théorique, telle que celle des figures 2 a ou 3 a, construite à l'aide de cette valeur approchée. Or, la Transformation selon l'invention vise à créer un signal ki(d, E) quasi-périodique théorique ayant une période qui a la même dépendance de sa période, vis-à-vis de l'angle O, que le signal mesuré. La Transformation vise ensuite à multiplier le signal expérimental g(e) (ou mesuré) avec le signal théorique quasi- périodique ki(d,E) crée, puis à prendre l'intégrale de ce résultat sur tous les domaines d'angle O. Il se trouve que, à partir de la connaissance de cette variation de longueur d'onde A(o) de la fonction théorique de réflectivité, on peut définir une telle Transformation applicable à la fonction réelle g(O) pour obtenir finalement une courbe donnant l'intensité réfléchie uniquement en fonction de la profondeur d'interface d C'est donc uniquement pour pouvoir obtenir la variation de la longueur d'onde X(E) de la courbe quasi-sinusoïdale que l'on a bâti la structure théorique et réalisé la courbe d'intensité réfléchie théorique t, (G) ou t 2 (E) relative aux figures 2 a, 2 b, 2 c, ou bien 3 a, 3 b, 3 c. Dès qu'on obtient la connaissance de la variation de longueur d'onde k(G) de la courbe quasi-sinusoïdale 2 c ou 3 c, on peut réaliser l'intégration de la courbe mesurée de la figure lc incluant une somme de sinus, en tenant compte du fait que chacun des sinus varie en longueur d'onde avec la même variation de longueur d'onde que la courbe théorique. Donc après le filtrage, dont le résultat h 1 (E) ou h 2 (e) est représenté sur la figure 2 c ou 3 c selon le modèle de structure théorique choisi, on réalise k(O) puis on synthétise à partir de X(G) ledit signal quasi-périodique ki(d,E) ayant une dépendance de sa période en fonction de e identique à celle de chaque sinus dans le signal mesuré. Cette fonction X(E) est croissante à partir de Gr et tend vers un maximum Xmax(e) pour un angle de rasance Emax L'angle O,,. pris en considération sera en principe l'angle maximal utilisé lors d'une mesure de test sur la structure réelle 10 a. On définit ensuite, de préférence, mais non obligatoirement, une longueur d'onde relative Xrel(O) par le rapport k(E) sur k Xax(O) afin de normaliser ladite longueur d'onde X(E). irel (E) =I(Of) / Xmax(O)- La fonction,reî(O) est représentée sur la figure 3 d. Elle a la même allure que la fonction k(O). En fait, on pourrait se contenter d'appliquer les opérations ultérieures sur la fonction X(O) elle-même Cependant au cours des expériences pratiquées pour la mise en oeuvre de l'invention, on a trouvé de meilleurs résultats en appliquant ces opérations sur Xrel(e) ce qui constitue donc une optimisation de la méthode d'utilisation du dispositif selon l'invention. La fonction Xrel(e) ne dépend pas non plus de la valeur de la profondeur d Donc elle est identique pour toutes les interfaces, quelle que soit la distance de l'interface avec la surface supérieure de la structure considérée. Maintenant, à partir de la fonction ?( 1), ou bien Xrei (E)' on définit une nouvelle fonction que l'on a appelé CHIRP (et que l'on écrira CH) Par la fonction CHIRP, on va créer ladite fonction quasi- sinusoïdale à changement monotone de période ki(d,E), c'est-à-dire dont la période varie de façon monotone avec la longueur d'onde k(E) Par changement "monotone", il faut entendre "ayant toujours le même signe". La fonction CHIRP s'écrit selon la formule CH(E 3 k) = CH(Ekl) + 1 Bre I (Ek) et est représentée sur la figure 3 e On voit quelle est croissante. La courbe X( 0) a été déterminée de façon discrète par une série de points reportés sur la figure 3 d pour chaque angle de rasance e 1, e 2, * k-1, k auquel le calcul de l'intensité théorique est fait C'est pourquoi la fonction CH est également calculée de façon discrète La formulation proposée permet de passer de la série de points de la courbe de la figure 3 d, à la série de points de la courbe de la figure 3 e qui représente la fonction CH(@) obtenue en fonction de E. Donc, pour arriver à cette étape de l'utilisation du dispositif selon l'invention, on a jusqu'ici eu besoin seulement de l'indice de réfraction moyen de la structure réelle à tester, pour bâtir une structure théorique l Ob ou bien 1 Oc Et avec cette seule donnée, on a calculé la fonction CH(e). On formule alors la fonction pseudo-sinusoïdale ki(d,O) recherchée au moyen de la relation K (d,) = sin(P d CH(G) + i 2) o 5 est une constante et o i est égal à O ou 1. On a maintenant, par rapport à l'intégrale de Fourier classique, remplacé la fonction sin G par une fonction sin (d,E) dont la longueur d'onde varie comme décrit précédemment On s'est ainsi affranchi des inconvénients que l'on aurait obtenus en utilisant une intégrale de Fourier simple. Pour éviter que le résultat final ne soit trop influencé par un décalage d'angle du signal mesuré, le dispositif selon l'invention prévoit deux opérations dans l'intégration: une première multiplication de g(E) avec le signal théorique quasi-sinusoïdal créé, puis une seconde multiplication avec un signal créé identique, de même quasi-période, mais décalé d'1/4 de période. La transformation comprend donc une première intégrale avec le signal créé non décalé, et une seconde intégrale avec le signal créé décalé, suivi de la sommation de la valeur absolue de ces intégrales. Ce résultat est obtenu en appliquant le signal théorique quasi-sinusoïdale sous la forme d'un sinus dans la première intégrale, puis sous la forme d'un cosinus dans la seconde intégrale, ce qui donne le décalage voulu d'1/4 de période. Si dans la transformation, on se restreignait au calcul d'une seule intégrale, on pourrait se trouver dans des cas o le résultat serait zéro du fait du décalage du signal mesuré. On considère donc le signal g(E) représenté sur la figure lc, obtenu à l'issue du filtrage par le filtre mathématique F 1 de la fonction f(G) représentée sur la figure lb et on lui applique les deux intégrales suivantes e G(d)= f g(>) K (, x) de + f g(E> Kl(e X)dé Une amélioration est encore obtenue en réalisant une nouvelle opération donnant une fonction G'(d) qui est la fonction G(d) divisée par l'intervalle jmin, Emax telle que: G'(d) = G(d)/(emax Omîn)* Les fonctions Ki et K 2 sont données par: Ko = Sin( d CH(O)) Ki = Cos( d CH(e)) dans lesquelles 5 est une constante égale à 4 n/lx et d est la profondeur dans la structure réelle (On rappelle que 4 x est la longueur d'onde de la source de rayons X). On prend donc la valeur absolue de chacune des intégrales et on fait la somme de ces valeurs absolues. Dans ces intégrales, on reconnaît outre la fonction expérimentale filtrée g(E), les deux fonctions mettant en oeuvre l'une le sinus de la fonction pseudo-sinusoïdale synthétisée plus haut, l'autre un cosinus des mêmes paramètres et variables. Pour que l'intégrale préconisée selon l'invention donne de bons résultats, il faut choisir correctement l'intervalle d'intégration: il faut intégrer sur un nombre entier de périodes de la fonction h 1 (E) ou h 2 (E) théorique choisie Il faut en outre intégrer sur un nombre entier de périodes le plus grand possible à l'intérieur du domaine considéré ec % Donc sur la courbe de la figure 2 c, ou celle de la figure 3 c, on choisit d'abord un angle minimum et un angle maximum qui définissent le domaine sur lequel on veut effectuer la Transformation On définit ensuite l'intervalle entier pour l'intégration: Emin Emax L'angle emin est défini comme le premier zéro de la courbe de la figure 2 c ou 3 c, que l'on trouve lorsque e augmente On regarde alors le signe de la dérivée de la fonction h 1 ou h 2 à ce premier zéro, puis on se reporte vers la fin de l'intervalle choisi, c'est-à-dire vers l'angle maximum choisi; et à partir de ce dernier on redescend jusqu'au premier zéro o la fonction h 1 ou h 2 présente une dérivée de même signe On a ainsi défini l'intervalle d'intégration le plus grand possible emin Emax avec un nombre entier de périodes de la fonction h 1 ou h 2 Les angles Emin et Ema,, dépendent évidemment de la variable d et devraient en toute rigueur être écrits Emin(d) et emax(d) On écrira ci-après e i et E par souci de simplicité dans l'écriture. L'angle emin est choisi proche, mais néanmoins pas trop près de l'angle critique E),, afin d'éviter d'éventuelles distorsions. Dans la pratique de l'analyse des structures de matériaux semiconducteurs, on choisit emin de l'ordre de 0,40 L'angle emax quant à lui est choisi aussi éloigné que possible, tout en évitant la région o le bruit devient trop important En effet, si l'on considère la courbe f(e) de la figure lb, on constate que le bruit devient important aux fortes valeurs de e du fait que l'intensité mesurée devient faible Il faut donc parfois couper l'intervalle d'intégration avant l'angle maximum t utilisé pour réaliser la mesure I = f(E). Il faut cependant prendre en compte que plus grand est le domaine de E choisi, c'est-à-dire l'intervalle Emin ma X' plus fins seront les pics de la courbe I = F(d) de la figure 6, et meilleure la définition des profondeurs d'interface d 1, d 2, dn. Le résultat de la Transformation n'est pas encore la fonction G(d) ellemême, telle que donnée plus haut En effet, cette fonction G(d) est très accidentée comme montré sur la figure 5; elle présente de nombreux pics très fins et très rapprochés Si le résultat de la Transformation s'arrêtait à la courbe représentative de la fonction G(d), il serait encombré d'informations gênantes non essentielles. C'est pourquoi la fonction G(d) telle qu'obtenue après l'intégration est passée dans un filtre F 2 qui réalise une diffusion symétrique La diffusion symétrique effectue un lissage des pics fins et multiples sans déplacement de la valeur maximale du pic en e Ce filtre F 2 est un filtre passe-bas. Il en résulte la fonction F(d) recherchée qui est représentée sur la figure 6, sur laquelle la position des pics et leur hauteur relative donnent une cartographie directe des profondeurs dl, d 2, d d'interface de la structure réelle i Qa à tester et les intensités associées D'o l'on peut calculer la concentration x exacte des éléments dans les couches. En conclusion, la Transformation selon l'invention ressemble à une transformée de Fourier simple Ce n'en est cependant pas une. En particulier, les fonctions sin et cos de la transformée de Fourier classique sont considérées entre les limites, ±, alors qu'ici nous avons des fonctions périodiques sur un intervalle limité. D'autre part la longueur d'onde de ces fonctions varie comme la longueur d'onde relative à une structure théorique l Ob ou l Oc, bâtie au moyen d'une couche sur un substrat, couche dont l'indice de réfraction diffère seulement légèrement de celui du substrat C'est-à- dire que l'argument des sinus et cosinus varie. Dans la plupart des applications de l'invention, la fonction CHIRP peut être calculée directement à partir de O, au moyen de la relation suivante: CH(@) = (sin 2 O sin 2 Oc)V/cos E et utilisée sous cette forme dans les intégrales de calcul de la fonction G(d). Cette fonction peut être encore écrite pour les petits angles avec l'approximation suivante CH(@) = (sin 2 G sin 2 O,)l Ce moyen direct d'évaluation de la fonction CH( 8) a été obtenu selon l'invention à partir de la mise en oeuvre plus complexe décrite précédemment, et est donc une approximation de la fonction obtenue point par point à partir de k(G) représentée sur la figure 3 e. Cependant cette fonction calculée est très proche de la précédente La différence est tout à fait minime et négligeable dans la plupart des cas. Dans la description précédente, on a énoncé les moyens pour parvenir aux buts de l'invention par la création d'une fonction quasi-sinusoïdale hi ou h 2, puis par la création de la fonction CH(s), suivies de la synthèse de la fonction quasi-sinusoïdale sin(d,G). Mais les fonctions quasi-sinusoïdales Ki(d,E) de l'intégrale G(d) peuvent très bien être remplacées par des fonctions quasi-périodiques dont la période varie de façon monotone en fonction de O de la même manière Par exemple on peut utiliser des fonctions triangulaires périodiques. On a présenté dans ce qui précède la possibilité de partir d'une structure hypothétique l Oc avec une couche d'oxyde: cette hypothèse est particulièrement valable dans le cas o la structure réelle l Oa est elle-même recouverte d'une couche d'oxyde Ce cas est très fréquent dans la réalité quotidienne du fabricant de dispositifs semiconducteurs. Donc si la structure réelle présente une surface oxydée, et si on utilise pour arriver aux buts de l'invention la filière de calcul passant par la structure théorique sans oxyde (Fig 2 a), ou bien passant par la formulation mathématique directe de CHIRP (e), on arrivera à évaluer les positions des interfaces dl, d 2, dn avec un écart de moins de 1 nm par rapport à leur position exacte Ces positions seront cependant mieux évaluées en passant par la filière de calcul correspondant à la structure avec oxyde (fig 3 a). Néanmoins, tant que l'oxyde de la structure réelle est mince, on a une très bonne approximation en réalisant le calcul direct de CH(@) et on gagne considérablement en temps de calcul et en complexité du dispositif. On peut donc dans la plupart des cas s'en tenir au calcul direct de CH(@), tout en conservant, dans le dispositif, les moyens pour un calcul plus précis éventuellement. L'éventuel décalage entre les résultats obtenus avec et sans tenir compte de la couche d'oxyde vient du fait suivant. Lorsque la surface supérieure d'une structure de couches est couverte d'un oxyde, on trouve deux interfaces proches l'une de l'autre la première air-oxyde, la seconde oxyde-semiconducteur, ces deux interfaces ayant des réflectivités du même ordre de grandeur. On a dit plus haut que l'on allait dans ce cas considérer la couche d'oxyde comme une interface moyenne Mais, le positionnement de cette interface moyenne située entre ladite première et seconde interfaces réelles dépend de l'épaisseur effective de l'oxyde Plus l'oxyde est épais, plus forte est l'erreur induite par l'adoption de l'interface moyenne. Le dispositif selon l'invention s'applique, à la différence des dispositifs qui font appel aux comparaisons avec des Tables de Classement, à des structures dans lesquelles les couches ont de très faibles différences d'indices, hormis la couche d'oxyde supérieure Pour des différences d'indices supérieures à 1,5 ppm, les linéarités sont perdues, et il en résulte que la cartographie directe de d, devient moins précise. L'invention vise une application pratique à des hétérostructures de matériaux semiconducteurs ayant des faibles différences d'indice parmi lesquels par exemple: Ga As, Ga Al As, Ga In As, Al In As, In P, avec de nombreuses variations et combinaisons possibles de composés ternaires et quaternaires. les systèmes Si-Ge, Cd-Hg-Te, Zn-S, Zn-Se, Zn- Cd, etc On peut aussi appliquer cette méthode à des couches de revêtement sur des substrats non semiconducteurs. En référence avec la figure 4 b, le dispositif selon l'invention comporte donc les blocs fonctionnels suivants: un bloc fonctionnel 130 pour le stockage des données mesurées I = f(O) relatives à la structure réelle i Qa à tester; un bloc fonctionnel 131 pour le traçage de la courbe représentative I = f(G) telle que montrée sur la figure lb. Le dispositif comporte d'autre part les blocs fonctionnels suivants: un bloc 140 pour le stockage de données concernant la structure réelle 1 Oa, c'est-à-dire l'indice de réfraction du substrat, l'indice de réfraction moyen de l'ensemble des couches de la structure réelle, la longueur d'onde ?x des rayons X incidents; un bloc 141 pour la définition d'une structure hypothétique l Ob ou l Oc à partir des données du blocs précédent 140 un bloc 142 pour le calcul de la réflectivité I = t(e) ou bien I = t 2 ( 9) relative à l'une ou l'autre des structures hypothétiques choisies l Ob ou 1 Oc. La sortie du bloc 142 est ensuite passée dans le bloc 131 de traçage de courbe représentative pour obtenir les courbes tj(G) et t 2 (O) telles que montrées sur les figures 2 b et 2 c. La sortie du bloc 131 de traçage de courbe est portée dans un bloc fonctionnel 132 réalisant le filtrage par le filtre F 1 de toutes les fonctions à ce stade, c'est-à-dire aussi bien de la fonction mesurée f(O) pour obtenir la fonction g(e) que de la fonction théorique t,(e) pour obtenir la fonction hl(E) * ou de la fonction théorique t 2 (e) pour obtenir la fonction h 2 (O). Les fonctions théoriques h 1 (O) ou h 2 () sont ensuite portées dans le bloc fonctionnel 143 qui effectue le calcul de X(O) et éventuellement de Xrel(e) puis dans le bloc 144 qui effectue le calcul de CH(@) à partir de k(E) ou Xrei(e)* La fonction g(e) issue du bloc de filtrage 132 et la fonction CH(e) issue du bloc 144 sont ensuite introduites dans un bloc de calcul 134 des intégrales pour fournir la fonction G(d). Puis la fonction G(d) est passée dans le filtre F 2 de diffusion matérialisé par le bloc 135 La sortie du bloc 135 fournit alors la fonction F(d) résultat de la Transformation selon l'invention. Enfin la fonction finale F(d) est passée dans le bloc 137 de traçage qui en fournit la courbe représentative de la figure 6. Entre le bloc 135 de filtrage et le bloc 137, un bloc 136 effectue les calculs relatifs à la position exacte des interfaces et à l'interpolation des intensités associées aux pics pour déterminer la composition des couches de part et d'autre de chaque interface. En référence avec la figure 4 a, le dispositif est simplifié Il ne comporte pas de calcul de X(O) mais seulement le calcul direct de CH(E) = (sin 2 e sin 2 EO)I/cos E dans un bloc 144 ' Les autres blocs ont les mêmes fonction que ceux de la figure 4 b. REVENDICATIONS 1 Dispositif de traitement d'un signal mesuré correspondant à l'intensité réfléchie, en fonction de l'angle de rasance ( 9) à une longueur d'onde ( 4 x) de rayons X, par une structure réelle de couches multiples empilées sur un substrat ( 100), ces couches présentant à leurs interfaces des différences faibles d'indice de réfraction, ce dispositif de traitement comprenant des moyens pour effectuer la transformation de la fonction (f(O)) correspondant au signal d'intensité réfléchie mesuré en fonction de l'angle de rasance (O) dans un domaine d'angles de rasance limité par des valeurs d'angles minimales et maximales Emin émax' en une fonction (F(d)) correspondant au signal d'intensité formulé en fonction de la profondeur (d) dans la structure réelle de couches multiples, cette fonction (F(d)) obtenue permettant d'établir une cartographie directe des profondeurs de chaque interfaces (dl, d 2, dj) de la structure réelle, et des intensités réfléchies associées. 2 Dispositif selon la revendication 1, caractérisé en ce que, la fonction (f(O)) correspondant au signal d'intensité réfléchie en fonction de l'angle de rasance (O) présentant une surmodulation qui est la somme approximativement de fonctions périodiques à périodes variables de façon monotone en fonction de l'angle de rasance (O), les moyens de transformation de la fonction mesurée (f(O)) du signal d'intensité en fonction de l'angle de rasance (O), en la fonction (F(d)) du signal d'intensité en fonction de la profondeur (d) dans la structure comprennent des moyens d'intégration, sur le domaine O min - Omax de la fonction mesurée (f(O)) multipliée par une fonction ki(d,O) périodique de la profondeur (d) et de l'angle de rasance (O), dont la période varie en fonction de l'angle de rasance (O) avec la même variation monotone que celle des fonctions périodiques qui sont les composantes de la fonction mesurée (f(e)). 3 Dispositif selon la revendication 2, caractérisé en ce qu'il comprend en outre un premier filtre F, passe-bande de type passe- haut pour couper les grandes longueurs d'onde du signal d'intensité réfléchie mesurée f(O). 4 Dispositif selon la revendication 3, caractérisé en ce qu'il comprend des moyens pour synthétiser ladite fonction périodique ki(d,E) de la profondeur (d) et de l'angle de rasance (O), dont la période varie de façon monotone en fonction de l'angle (O) d'incidence dans le domaine d'angles Emin emax' ces moyens de synthèse comprenant au moins des moyens de calcul appliqués à une structure hypothétique comportant au moins une couche d'indice de réfraction proche de l'incidence de réfraction moyen de la structure à couches multiples réelles, cette couche étant disposée sur un substrat d'indice de réfraction sensiblement égal à celui de la structure réelle et éventuellement une couche supérieure d'oxyde, lesdits moyens de calculs incluant: des premiers moyens de calcul d'une fonction correspondant à un signal d'intensité réfléchie théorique à la longueur d'onde de rayons x, (?x), en fonction de l'angle de rasance (O), dans le domaine d'incidence rasante Gmin emax' des seconds moyens de calcul d'une fonction (?"O)) correspondant à la longueur d'onde calculée pour chacune des valeurs ( 1 D 2, E)k s de l'angle d'incidence O dans le domaine d'angles min - max ' des troisièmes moyens de calcul d'une fonction (CH (G)> présentant la même variation monotone en fonction de l'angle (G) que celle des composantes de la fonction mesurée, ladite fonction CH (O) étant formulée par la relation CH (EK) = CH (EK-1) + 1/X(@), ladite fonction périodique ki(d,E) étant formulée par sin (P d CH O)) + i n/2), avec i = O ou 1, et g = constante. 5 Dispositif selon la revendication 4, caractérisé en ce qu'il comprend aussi un filtre F 1 passe-bande du type passe-haut pour couper les grandes longueurs d'onde du signal d'intensité réfléchie théorique. 6 Dispositif selon la revendication 3, caractérisé en ce qu'il comprend des moyens pour synthétiser ladite fonction périodique de la profondeur (d) et de l'angle ( 0) de rasance, formulée par sin (d, 0) dont la période varie de façon monotone en fonction de l'angle (O) d'incidence dans le domaine Emin Emax, ces moyens de synthèse comprenant au moins des moyens de calcul d'une fonction (CH(E)) présentant la même variation monotone en fonction de l'angle O que celle des composantes de la fonction mesurée, ladite fonction CH étant formulée par la relation: CH( 8) = (sin 28 sin 28 c)1/2/cos E o Ec = la valeur de l'angle critique pour la structure réelle formée de couches multiples, à la longueur d'onde >x de rayons x de travail, ladite fonction périodique étant formulée par Ks(d, 8) = sin(P d CH() + i-) o 5 est une constante égale à 4 n/4 x et o i égal soit O soit 1. 7 Dispositif selon l'une des revendications 5 ou 6, caractérisé en ce que les moyens d'intégration réalisent le calcul des intégrales, et la somme de leur valeur absolue selon la relation: p G(d) = f f( 8) ko(d, 8) d + f f( 8) k 1 (d,8 e) d 8 8 Dispositif selon la revendication 7, caractérisé en ce que les moyens d'intégration réalisent le calcul des intégrales et la somme de leur valeur absolue selon la relation: G'(d) = G(d)/(Emax Emin) 9 Dispositif selon l'une des revendications 7 ou 8, caractérisé en ce qu'il comprend en outre un filtre F 2 de diffusion symétrique pour lisser la fonction G(d) issue de l'intégration, et pour obtenir la fonction de transformation F(d).