Ce travail a pour cadre la diffusion des ondes électromagnétiques par une ou deux interfaces rugueuses séparant des milieux homogènes. Il s'intéresse plus particulièrement aux modèles asymptotiques, qui permettent de résoudre le problème posé de manière rapide, au détriment d'un domaine de validité restreint.<br />Pour le cas simple interface, après un panorama des méthodes existantes, une étude approfondie de l'approximation dite de Kirchhoff est menée pour le cas de la diffraction en réflexion et en transmission par une simple interface. Cette méthode est réduite à l'approximation dite de l'optique géométrique, valide pour des interfaces fortement rugueuses comparativement à la longueur d'onde, pour calculer plus simplement et plus rapidement la puissance diffusée. Le phénomène d'ombrage de la surface, bien connu pour le cas de la réflexion, l'est beaucoup moins pour le cas de la transmission ; c'est pourquoi il est étudié en détail dans cette thèse. <br />Pour le cas double interface, une étude bibliographique des méthodes existantes nous permet de constater l'absence de méthode basée sur l'extension de l'approximation de Kirchhoff au cas de deux interfaces fortement rugueuses. Ainsi, la méthode développée dans cette thèse permet de pallier ce manque. Cette méthode est exposée en supposant que les deux surfaces sont décorrélées, afin de pouvoir obtenir une expression de la puissance diffusée simple à mettre en œuvre. Par comparaison avec une méthode numérique de référence, la méthode développée a été validée dans le cas bidimensionnel. Une application à la détection de nappes de pétrole sur la surface de la mer est présentée, et la méthode est étendue au cas tridimensionnel.

Cette thèse porte sur la diffraction d’une onde plane électromagnétique par une surface lisse parfaitement conductrice (PEC). Elle présente l’analyse des propriétés d’une nouvelle formulation des trois principales équations intégrales de frontières de la théorie de la diffraction électromagnétique (EFIE, MFIE et CFIE). L’idée est d’adapter les équations intégrales conventionnelles à la diffraction d’une onde plane en supposant que la fonction de phase de l’onde plane incidente détermine la fonction de phase de la distribution de courant induit sur la surface.L’idée d’utiliser la phase dans la diffraction d’ondes planes a déjà été étudiée pour les hautes fréquences, notamment dans les thèses de Zhou (1995) et Darrigrand (2002) qui adaptèrent les espaces d’approximation des éléments finis. Dans cette thèse, cependant, nous suivons une formulation plus récente, donnée par Herberthson (2008), où la fonction de phase est incorporée dans la distribution du noyau des opérateurs intégraux.En présentant les versions modifiées de l’EFIE et de la MFIE (dénommées HEFIE et HMFIE)dans des espaces fonctionnels appropriés, nous prouvons ici l’existence d’une solution unique à cette formulation spécifique et présentons une mise en oeuvre pratique originale qui tire parti de l’expérience acquise sur l’EFIE/MFIE. Par la suite, nous explorons une propriété importante offerte par ces nouvelles formulations: la possibilité de réduire le nombre de degrés de liberté requis pour obtenir une solution précise du problème. / This thesis is about the scattering of an electromagnetic plane wave incidenton a perfectly conducting smooth surface. It presents the analysis of the properties of a newformulation of the three principal boundary integral equations of electromagnetic scattering theory(EFIE, MFIE and CFIE). The basic idea is to adapt the conventional integral equations toplane-wave scattering by supposing that the phase function of an incident plane wave determinesthe phase function of the induced boundary current distribution.This idea of using the phase in plane wave scattering has previously been studied in highfrequencyscattering, in particular in the theses by Zhou (1995) and Darrigrand (2002) whoadapt the finite element approximation spaces. In this thesis, though, we follow a more recentformulation, given by Herberthson (2008), where the phase function is incorporated in the kerneldistribution of the integral operators.Presenting the modified version of the EFIE and the MFIE (denoted HEFIE and HMFIE) inappropriate function spaces, we prove the existence of a unique solution to this specific formulationand developp an original practical implementation which takes advantage of the gainedexperience on the EFIE/MFIE. Then, we explore another important property provided by thenew formulations: the possibility to reduce the number of degrees of freedom required to get anaccurate solution of the problem.

L'évaluation de la teneur en eau est actuellement un objectif important pour la caractérisation des sols et la compréhension des mécanismes liés au transfert hydrique et à l'absorption racinaire. Il existe un lien entre les variations de la teneur en eau et la permittivité des sols. L'objectif de ces travaux est d'exploiter ce couplage pour définir un protocole d'imagerie active micro-onde d'un monolithe. En effet, lorsqu'une onde électromagnétique rencontre un objet ayant des dimensions caractéristiques de l'ordre de sa longueur d'onde, celle-ci après interactions, va induire un courant qui à son tour générera un champ électromagnétique diffracté. La mesure de ce champ va permettre de revenir aux informations relatives à l'objet, en particulier sa permittivité. Nous nous appuierons ici sur un système de mesure hyperfréquence constitué d'une cavité métallique cylindrique. Il a donc fallu mettre au point des outils numériques de modélisation et d'imagerie de permettant de décrire le phénomène de diffraction dans un tel environnement. L'imagerie micro-onde quantitative est un problème inverse non-linéaire mal posé. De plus, les mesures du champ électromagnétique sont perturbées, ce qui va induire de l'instabilité dans le problème inverse. C'est pourquoi il est nécessaire de mettre en place un algorithme itératif avec ajout d'informations a priori (polynômes de Zernike) pour obtenir une solution stabilisée. Des reconstructions obtenues à partir de champ mesurés sur différents fantômes sont discutées. De plus, les informations fournies par l'étude du comportement de l'opérateur de diffraction renseignent sur la quantité de données accessibles et les informations a priori qui peuvent être légitimes de considérer.

La thèse porte sur la modélisation électromagnétique et la simulation de composites stratifiés plans (laminés), renforcés par des fibres organisées périodiquement. L'objectif est d'acquérir une bonne compréhension du comportement électromagnétique de telles structures, en première et étape de ce que pourrait ultérieurement être la production d’images mettant en évidence la localisation de zones éventuellement endommagées, et fournissant une certaine quantification de celles-ci. La thèse proprement dite se concentre donc sur la construction et l’évaluation de modèles de la diffraction électromagnétique par des composites multicouches tels que chaque couche est renforcée par des fibres disposées périodiquement.Est d’abord investiguée la diffraction par une plaque diélectrique (mono-couche) au sein de laquelle des fibres cylindriques de section circulaire de même rayon sont incorporées périodiquement, ces fibres ayant la même orientation de leurs axes et la même distance de centre à centre. Un cas bidimensionnel impliquant des ondes planes E ou H-polarisées, ainsi que des faisceaux gaussiens, normalement ou obliquement incidents, est d'abord pris en considération afin de mieux comprendre principes et philosophies des méthodes de choix, le couplage de mode et l'expansion multipolaire. Puis le travail est étendu, la diffraction de la plaque sous un éclairement tridimensionnel (conique) étant alors traitée en détail, ce qui montre aussi le potentiel de la méthodologie mise en œuvre si l’on souhaite obtenir la réponse électromagnétique de la structure à une source ponctuelle.Un composite multicouche, plus courant, mais plus complexe, qui est fait d’un empilement de plaques l'une sur l'autre, est alors étudié. Deux différentes espèces de composites sont ici prises en compte. Pour étudier la première, dont les fibres dans les différentes couches possèdent les mêmes orientations, des méthodes à base de matrices dites S ou dites T sont introduites, impliquant entre autre de s’intéresser à une résolution convenable du système linéaire produit selon le couplage de mode à la transition entre deux couches adjacentes. Une investigation de la deuxième espèce de composites suit alors, pour lequel les fibres au sein des différentes couches sont orientées dans des directions différentes quelconques, ce que permet une extension précautionneuse des approches précédentes. Une certaine attention est également portée au problème de l'homogénéisation des composites, de manière à lier les démarches à petite échelle telles que développées dans la thèse à celles à grande échelle souvent les seules prises en compte dans le contrôle non destructif et l’imagerie des composites stratifiés.De nombreux résultats de simulations numériques sont proposés et validés autant que possible par des résultats de référence de la littérature (notamment dans le cas de cristaux photoniques) et l'utilisation de solveurs « brute-force ». L'accent est aussi mis sur des cas particuliers de matériaux composites (ceux à base de fibres de verre et ceux à base de fibres de carbone) qui sont le plus souvent rencontrés dans les applications pratiques, avec des bandes de fréquences appropriées choisies en accord avec le comportement des fibres, principalement diélectrique ou principalement conducteur. / The contribution is about the electromagnetic modeling of fiber-reinforced periodically organized composite laminates. The final goal is to gain a good understanding of their electromagnetic behavior as well as to acquire images that should exhibit the location of possibly damaged zones, and provide some quantification of these zones. The thesis focuses on the scattering of well-organized periodic structures and building up an efficient full-wave computational model for multilayered composites, wherein each layer is reinforced by periodically arranged fibers, which is the first step for further investigation of the disorganized one.The work firstly considered the scattering problem of a slab in which infinite circular fibers, with the same radius, are periodically embedded with the same orientation of their axes and the same center-to-center distance. A 2-dimensional problem with normally and obliquely incident E- and H-polarized plane waves as well as Gaussian beams is firstly considered for understanding the principles and philosophies of the used mode-matching method and multipole expansion. Then the work is extended to the investigation if the scattering of the slab under illumination of a conically incident 3-dimensional electromagnetic wave, which shows the potential of the work for obtaining the response of the structure to a point source.A more practical but complicated multilayered composite, constructed by stacking up the slabs one over the other, is further investigated. Two different composites are taken into account. To study the first composite, with fibers in different layers having the same orientations, T-matrix- and S-matrix-based methods are introduced into the work for solving the linear system produced by mode-matching at the boundaries between two adjacent layers. Then, further investigation of the second kind of composite, wherein the fibers within different layers are orientated into different directions, is carried out by extending the approach properly.Some attention is also given to homogenization issues, so as to link small-scale approaches as developed in the thesis with large-scale ones as often considered in non-destructive testing of composite laminates.Extensive numerical simulations are proposed, validated whenever possible by reference results taken from the literature (notably in the case of photonic crystals) and the use of brute-force solvers. Emphasis is also on special cases of composites (glass-fiber- and graphite-fiber-based ones) as most often faced in practical applications, with appropriate frequency bands in harmony with the dielectric or conductive aspect of the fibers.

Le Contrôle Non Destructif (CND) de matériaux composites multicouches pour des problèmes de qualité, viabilité, sécurité et disponibilité des systèmes qui impliquent des pièces fabriquées dans les industries aéronautiques et de l’automobile est devenu une tâche essentielle aujourd’hui. L'objectif visé par cette thèse est l’imagerie électromagnétique de structures complexes multicouches anisotropes, de plus en plus utilisées dans des applications, et encore source de sérieux défis à l'étape de leur modélisation et encore plus à l'étape souvent en enfance de leur imagerie. En utilisant une vaste gamme de fréquences, qui va des courants de Foucault jusqu’aux micro-ondes, il y a un fort besoin de rendre disponibles des procédures de modélisation et d'imagerie qui sont robustes, rapides, précises et utiles à la décision des utilisateurs finaux sur des défauts potentiels, tant donc en basse fréquence (BF) (matériaux conducteurs, type fibre de carbone) qu’en haute fréquence (HF) (matériaux diélectriques, type fibre de verre). De plus, il est important d'obtenir des résultats en des temps brefs. Cependant, cela nécessite la connaissance d’une réponse précise à des sources externes aux multicouches, en considérant les couches des composites comme non endommagées ou endommagées : on parle donc de solution du problème direct, avec le cas particulier de sources élémentaires conduisant aux dyades de Green (DGF). La modélisation et la simulation numérique du problème direct sont gérés principalement via une solution au premier ordre de la formulation intégrale de contraste de source impliquant le tenseur de dépolarisation des défauts, quand ceux-ci sont assez petits vis-à-vis de l’épaisseur de peau locale (cas BF) ou de la longueur d'onde locale (cas HF). La précision des DGF doit nécessairement être assurée alors, même si les sources se situent loin de l'origine, ce qui donne un spectre de dyades qui oscille très rapidement. La technique d'interpolation-intégration dite de Padua-Domínguez est ainsi introduite dans le but d'évaluer de façon efficace des intégrales fortement oscillantes.Néanmoins, les matériaux composites peuvent souffrir de divers défauts, lors du processus de fabrication ou pendant leurs utilisations. Vides d’air, cavités remplies de liquide, fissures, etc., peuvent affecter le fonctionnement correct des structures composites. Il est donc indispensable de pouvoir détecter la présence des défauts. Ici, l’insistance est sur la méthode bien connue d’imagerie dite MUltiple SIgnal Classification (MUSIC), qui est basée sur la décomposition en valeurs singulières (SVD) des DGF ; celle-ci est développée afin de localiser les positions de multiples petits défauts volumiques en interaction faible enfouis dans des milieux anisotropes uniaxiaux. Le principal inconvénient de la méthode MUSIC est cependant sa sensibilité par rapport au bruit. Par conséquent, des méthodes MUSIC avec une résolution améliorée et la Recursively Applied and Projected (RAP) MUSIC sont introduites afin de surmonter un tel inconvénient de l'algorithme standard et de fournir des résultats de qualité avec une meilleure résolution. De nombreuses simulations numériques illustrent ces investigations. / Non-Destructive Testing/Evaluation (NdT/E) of multi-layered composite materials for problems of quality, viability, safety and availability of systems involving manufactured parts (in aeronautics and in automotive industry, as a good example) has become an interesting and challenging task nowadays. The focus of the PhD thesis is on the electromagnetic (EM) imaging of complex anisotropic multi-slab composite panels as increasingly encountered in applications, yet source of strong challenges at modeling stage and even more at often-in-infancy imaging stage. From eddy-currents to microwaves, there is a strong need to make available modeling and imaging procedures that are robust, fast, accurate and useful to potential end-users’ decision about potential defects both at low-frequency (LF) (conductive materials, carbon-fiber like) and high-frequency (HF) (dielectric materials, glass-fiber like). Moreover, it is important to get the results in close-to-real-time. However, this requires an accurate response to external sources of the multilayers, considering the layers which these composite structures are made of as undamaged or damaged. The modeling at forward stage is managed via a first-order solution involving the dyadic Green’s functions (DGF) of the layers along with the depolarization tensor of the assumed defects when they are small enough vis-à-vis the skin depth (LF case) or the wavelength (HF case). The accuracy of the DGF has to be ensured even if the sources lie far away from the origin, which yields a fast-oscillating spectrum of the dyads. The Padua-Domínguez interpolation-integration technique is introduced herein in order to evaluate in an effective fashion fast-oscillating integrals.Damages or disorders, which these composite structures may suffer from, are of many kinds. One could mention voids, fluid-filled cavities or uniaxial defects with obvious impacts on the electromagnetic and geometric parameters of the multilayers. That is, the task to make available to end-users imaging algorithms tailored to detect the presence of defects. The well-known standard MUltiple SIgnal Classification (MUSIC) algorithm, which is based on the Singular Value Decomposition (SVD) of such DGF, is here applied to localize the positions of small multiple defects with weak interaction embedded in anisotropic uniaxial media. The main drawback of MUSIC is its sensitivity with respect to the noise. Therefore, MUSIC with enhanced resolution and Recursively Applied and Projected (RAP) MUSIC are introduced to overcome such a drawback of the standard algorithm and to provide quality results with better resolution.

Les radars millimétriques en bande W (75-110 GHz) sont en plein essor, grâce notamment aux progrès des circuits intégrés, permettant de réaliser des systèmes compacts à bas coût et haute résolution due à la courte longueur d’onde. Dans un premier temps, ces systèmes ont été utilisés à des fins de détection et de localisation, avec à terme, pour objectif l’identification. Ainsi, des systèmes d’imagerie radar ont été développés, notamment grâce à l’imagerie qualitative, basée par exemple sur l’imagerie radar par synthèse d’ouverture (SAR). Cependant, afin de reconstruire les propriétés électromagnétiques des objets pour une identification complète, il est nécessaire de développer des algorithmes de reconstruction quantitatifs. Le travail présenté dans ce manuscrit est de poser les bases d’un système d’imagerie qualitative et quantitative en gamme millimétrique pour la détection et l’identification des objets sur les pistes d’aéroport par tomographie, tenant compte de la polarisation de l’onde incidente. Au cours de cette thèse, un outil de simulation permettant de la résolution des problèmes direct et inverse, pour les deux types de polarisation à deux dimensions 2D-TE et 2D-TM, basé sur la méthode des moments (MoM) a été développé. La première étape a consisté en la validation du problème direct en effectuant des comparaisons numériques avec des solutions analytiques pour des cibles canoniques. Ensuite, des mesures expérimentales ont été effectuées et comparées aux résultats numériques. Enfin, les résultats des reconstructions obtenus ont permis de valider l’algorithme de reconstruction 2D développé pour l’imagerie quantitative. / Millimeter-wave radar systems in W-band (75-110 GHz) are booming, due to advances in integrated circuits, allowing the fabrication of low-cost and high-resolution compact systems, thanks to the short wavelength. First, these systems were used for detecting and localizing purposes, with the aim of identification. Thus, imaging radar systems have been developed, especially using qualitative imaging, based for example, on Synthetic Aperture Radar (SAR). Nevertheless, in order to reconstruct the electromagnetic properties of objects, for a complete identification, we must develop quantitative reconstruction algorithms. The work presented in this manuscript is to give the basis of a qualitative and quantitative millimeter wave imaging system for detecting and identifying foreign debris on airport runways using tomography, taking into account the polarization of the incident wave. In this thesis, a simulation tool for solving forward and inverse problems, for the two-dimensional polarization cases 2D-TM and 2D-TE, based on the method of moments (MoM) has been developed. The first step was to study the validation of the direct problem by comparing numerical results with analytical solutions for canonical targets. Then, experimental measurements have ben carried out and compared with numerical results. Finally, reconstruction results obtained have validated the reconstruction algorithm developed for quantitative imaging.

La profilométrie optique est une technique de métrologie de surface rapide et non destructive. Dans ce mémoire, nous avons abordé cette problématique par des méthodes inverses de diffraction électromagnétique et dans une configuration de type Microscopie Tomographique Optique par Diffraction (ODTM). La surface est sondée par un éclairement sous plusieurs angles d'incidences ; la mesure en amplitude et en phase du champ lointain diffracté constitue les données du problème. Des profils de surfaces ont été reconstruits en considérant différents modèles de diffraction, parmi lesquelles une méthode approchée fondée sur les approximations de diffusion simple et de paraxialité. La résolution latérale de cette méthode et des techniques classiques de profilométrie est limitée par le critère d'Abbe-Rayleigh, défini sur la base de l'ouverture numérique pour l'éclairement et la détection du champ. Afin de dépasser cette limite de résolution, nous avons développé une méthode itérative de Newton-Kantorovitch régularisée. L'opérateur de diffraction y est rigoureusement modélisé par une méthode des moments, résolution numérique des équations du formalisme intégral de frontière, et l'expression de la dérivée de Fréchet de cet opérateur est obtenue par la méthode des états adjoints, à partir du théorème de réciprocité. Pour les surfaces unidimensionnelles métalliques, notre technique permet d'inverser à partir de données synthétiques des surfaces très rugueuses avec une résolution au delà du critère d'Abbe-Rayleigh. / Optical profilometry is a nondestructive and fast noncontact surface metrology technique. In this thesis, we have tackled this issue with inverse scattering electromagnetic methods and in an Optical Digital Tomographic Microscopy (ODTM) configuration. The surface is probed with illuminations under several incidence angles; the measure of far scattered field amplitude and phase constitutes the problem data. Surface profiles have been reconstructed using different scattering models among which an approximate theory based on single scattering and paraxiality. The lateral resolution of this technique and classical profilometric approaches is limited by the so-called Abbe-Rayleigh's criterion defined out of the numerical aperture for illumination and field detection. In order to overpass this resolution limit, we have developed a regularized iterative Newton-Kantorovitch's method. The scattering operator is rigorously modelized with the method of moments, that is a numerical solution of boundary integral equations, and its Fréchet derivative adjoint states expression is deduced from the reciprocity theorem. For one-dimensional metallic surfaces, our method succeeds in inverting from synthetic data very rough surfaces with the resolutions beyond the Abbe-Rayleigh's criterion. The performance of this technique and inversion conditions clearly differ from one polarization to the other : in the TM case, interactions at longer distance than in the TE case improve yet the resolution. This work includes also an experimental validation of our inverse model on grooves in indium phosphure substrate at 633 nm.

Durant la dernière décade, les radars millimétriques en bande W (75 - 110 GHz) pour les applications civiles que ce soit dans le domaine de l'aide à la conduite ou de la sécurité. La maturité de ces systèmes et les exigences accrues en termes d'application, orientent actuellement les recherches vers l'insertion de fonctions permettant l'identification. Ainsi, des systèmes d'imagerie radar ont été développés, notamment à l'aide d'imagerie qualitative (SAR). Les premiers résultats sont très prometteurs, cependant, afin de reconstruire les propriétés électromagnétiques des objets, il faut travailler de manière quantitative. De nombreux travaux ont déjà été conduits en ondes centimétriques, mais aucun système d'imagerie quantitative n'existe, à notre connaissance, en gamme millimétrique. L'objectif du travail présenté dans ce manuscrit est de poser les bases d'un système d'imagerie quantitative en gamme millimétrique et de le comparer à l'imagerie radar de systèmes développés en collaboration avec l'Université d'Ulm (Allemagne). L'ensemble des résultats obtenus valide le processus développé pour d'imagerie quantitative. Les recherches doivent être poursuivies. D'une part le système de mesure doit évoluer vers un vrai système multi-incidences/multivues. D'autre part, le cas 2D-TE doit être implémenté afin de pouvoir traiter un objet 2D quelconque dans n'importe quelle polarisation. Enfin, les mesures à partir de systèmes radar réels doivent être poursuivies, en particulier pour rendre exploitables les mesures des coefficients de transmission. Ces dernières sont indispensables si l'on veut un jour appliquer les algorithmes d'inversion à des mesures issues de systèmes radar. / In the last decade, microwave and millimeter-wave systems have gained importance in civil and security applications. Due to an increasing maturity and availability of circuits and components, these systems are getting more compact while being less expensive. Furthermore, quantitative imaging has been conducted at lower frequencies using computational intensive inverse problem algorithms. Due to the ill-posed character of the inverse problem, these algorithms are, in general, very sensitive to noise: the key to their successful application to experimental data is the precision of the measurement system. Only a few research teams investigate systems for imaging in the W-band. In this manuscript such a system is presented, designed to provide scattered field data to quantitative reconstruction algorithms. This manuscript is divided into six chapters. Chapter 2 describes the theory to compute numerically the scattered fields of known objects. In Chapter 3, the W-band measurement setup in the anechoic chamber is shown. Preliminary measurement results are analyzed. Relying on the measurement results, the error sources are studied and corrected by post-processing. The final results are used for the qualitative reconstruction of all three targets of interest and to image quantitatively the small cylinder. The reconstructed images are compared in detail in Chapter 4. Close range imaging has been investigated using a vector analyzer and a radar system. This is described in Chapter 5, based on a future application, which is the detection of FOD on airport runways. The conclusion is addressed in Chapter 6 and some future investigations are discussed.

L'objet de la thèse est la détection électromagnétique non-itérative de petits objets enfouis. Le problème direct de diffraction est abordé en utilisant une formule asymptotique rigoureuse du champ diffracté par des inclusions dont la taille caractéristique est petite devant la longueur d'onde de leur illumination dans le milieu d'enfouissement. La prise en compte de la diffraction multiple dans le cas de deux inclusions sphériques est abordée grâce à un tenseur de polarisation spécifique qui est calculé dans un système approprié de coordonnées bisphériques. Le modèle de Foldy-Lax est aussi utilisé afin de prendre en compte le couplage entre plusieurs inclusions. Les simulations numériques montrent que cet effet de couplage ne peut être ressenti qu'en leurs voisinages immédiats. Une configuration d'enfouissement en demi-espace est aussi étudiée en détail. Les dyades de Green alors nécessaires sont calculées de manière exacte par "force brutale" numérique. Puis trois méthodes approchées de calcul des intégrales de Sommerfeld qui sont impliquées sont proposées, les simulations montrant qu'elles font gagner un temps de calcul significatif dans le calcul de ces dyades, tout en étant de précision convenable. La prise en compte du couplage entre une sphère et l'interface est aussi investiguée grâce à un tenseur de polarisation adéquat en coordonnées bisphériques (de facto, une des deux sphères dégénère en cette interface). A chaque fois, les champs diffractés simulés par la méthode asymptotique sont comparés à des champs obtenus par la méthode dite des dipôles couplés (CDM). Les résultats montrent que la méthode asymptotique fournit des valeurs du champ diffracté satisfaisantes tant que les tailles des inclusions restent assez petites devant la longueur d'onde. L'algorithme d'imagerie MUSIC est quant à lui utilisé pour détecter ces inclusions à partir de leur matrice de réponse multistatique (MSR) collectée via un réseau plan d'extension limitée de dipôles émetteurs-récepteurs idéaux. L'analyse des valeurs et des vecteurs singuliers de la matrice MSR montre qu'il existe une différence entre les données calculées par la méthode asymptotique et celles calculées par la méthode CDM. Mais cette différence ne persiste pas si l'on considère des données bruitées, même à relativement faible niveau de bruit. Dans les deux cas, MUSIC permet une estimation fiable de la position des inclusions, la notion de "super-localisation" étant en particulier discutée. Une méthode est par ailleurs proposée afin de détecter l'angle d'inclinaison d'un ellipsoïde incliné enfoui.