The Landweber Iteration for an Inverse Scattering Problem

Author(s):  
Martin Hanke ◽  
Frank Hettlich ◽  
Otmar Scherzer

Abstract A Landweber iteration scheme is presented for the numerical solution of an inverse obstacle problem. The method uses a recently obtained characterization of the Fréchet derivative of the far field operator and its adjoint. The performance of the method is illustrated by some numerical examples. Some theoretical aspects are pointed out to motivate the use of nonlinear Landweber iteration.

2007 ◽  
Author(s):  
Κωνσταντίνος Αναγνωστόπουλος

The scope of this doctoral thesis is, first, to develop an analytical, in principle, method for the solution of the two-dimensional scattering problem of time-harmonic elastic plane waves by a homogeneous orthotropic scatterer, second, to establish the complete theoretical framework, which is necessary for the application of the Linear Sampling Method (LSM) to the problem of reconstructing the support of twodimensional elastic anisotropic inclusions embedded in isotropic media and, third, to derive an extension of the Factorization Method (FM) to the inverse elastic scattering problem by penetrable isotropic bodies for time-harmonic plane wave incidence. Aconcise description of the contents of the thesis is outlined below. Chapter one contains a detailed bibliographical search, which is related to the analytical and numerical methods (with emphasis on the former) usually employed for the solution of the direct scattering problem by anisotropic elastic bodies as well as to those inverse scattering techniques, which are usually referred to as sampling and probe methods and, in particular, the LSM and the FM. Chapter two commences with a brief discussion of some fundamental results from the linearized theory of dynamic elasticity. The problem of a rigorous analysis of the elasticity equation governing the elastic behaviour of an orthotropic material in two dimensions is then addressed. This analysis, which is based on a suitable diagonalization applied to the underlying differential system and a plane wave expansion of the sought field, results in a Fourier series expansion for the displacement field describing the elastic deformations of the orthotropic medium and is complemented by the results of appendix A. A mathematical model for the solution of the associated transmission scattering problem, taking advantage of the aforementioned expansion, is also settled and analyzed. The details of its numerical treatment can be found in appendix B. Finally, numerical results for several inclusion geometries and a system thereof with material properties characterized by the cubic symmetryclass -a special case of the orthotropic class of symmetry- are presented. In chapter three, the LSM is extended to the case of a two-dimensional homogeneous anisotropic inclusion embedded in an isotropic background medium. The concepts of the elastic Herglotz function, the elastic far-field operator and the corresponding far-field equation, on which the formulation of the LSM heavily relies, are first introduced. Then, the proposed inverse scattering scheme is introduced and discussed in detail. By means of an appropriate operator decomposition of the far-field operator,the main theorem of the method, concerning the characterization of the behaviour of an approximate solution to the far-field equation at the boundary of the scatterer, is proved. In the end of the third chapter, the performance of the LSM is examined by applying it to a set of different geometric configurations of the elastic inclusion, filled with a cubic anisotropic material. An investigation of the effect of the various parameters entering the problem, such as the scatterer’s degree of anisotropy, the polarization of the elastic point source located at the sampling point and the noise level in the synthetic far-field data, on the reconstructed geometric profiles’ quality,is carried out. In the fourth chapter, the FM is elaborated for the shape reconstruction of a penetrable isotropic elastic body from the knowledge of the far-field pattern of the scattered fields for plane incident waves. The theoretical analysis is conducted in three dimensions and focuses on deriving a factorization of the far-field operator, which is the cornerstone for the applicability of the particular inversion scheme, and investigating thorougly the properties of the involved operators. This investigation gave birth to a number of interesting by-products and one of them, namely, a regularity estimate for the solution of a particular form of the corresponding interior transmission problem, is the subject matter of appendix C. By means of the proposed factorization, a series of theorems, which finally lead to an explicit characterization of the scattering obstacle, is then proved. In the end of the chapter, the performance of the investigated inverse scattering technique is demonstrated by applying it to specific two-dimensional elastic scatterer reconstruction problems involving different scatterer configurations and various choices for their constitutive parameters. The effect of using different levels of additive random noise in the forward synthetic data and combining results obtained for different polarizations of the elastic point source located at the sampling point, on the quality of the reconstructed profiles, is also examined. Finally, chapter five draws the conclusions that flow from the foregoing chapters and discusses the contribution of this doctoral thesis. A brief discussion about possible future studies is also included.


Author(s):  
Jianli Xiang ◽  
Guozheng Yan

Abstract This paper is concerned with the inverse scattering problem of time-harmonic elastic waves by a mixed-type scatterer, which is given as the union of an impenetrable obstacle and a crack. We develop the modified factorization method to determine the shape of the mixed-type scatterer from the far field data. However, the factorization of the far field operator $F$ is related to the boundary integral matrix operator $A$, which is obtained in the study of the direct scattering problem. So, in the first part, we show the well posedness of the direct scattering problem by the boundary integral equation method. Some numerical examples are presented at the end of the paper to demonstrate the feasibility and effectiveness of the inverse algorithm.


1991 ◽  
Vol 34 (2) ◽  
pp. 295-311 ◽  
Author(s):  
Yongzhi Xu

The inverse scattering problem for acoustic waves in shallow oceans are different from that in the spaces of R2 and R3 in the way that the “propagating” far-field pattern can only carry the information from the N +1 propagating modes. This loss of information leads to the fact that the far-field pattern operator is not injective. In this paper, we will present some properties of the far-field pattern operator and use this information to construct an injective far-field pattern operator in a suitable subspace of L2(∂Ω). Based on this construction an optimal scheme for solving the inverse scattering problem is presented using the minimizing Tikhonov functional.


2015 ◽  
Vol 2015 ◽  
pp. 1-9 ◽  
Author(s):  
M. T. Bevacqua ◽  
L. Di Donato

Total Variation and Compressive Sensing (TV-CS) techniques represent a very attractive approach to inverse scattering problems. In fact, if the unknown is piecewise constant and so has a sparse gradient, TV-CS approaches allow us to achieve optimal reconstructions, reducing considerably the number of measurements and enforcing the sparsity on the gradient of the sought unknowns. In this paper, we introduce two different techniques based on TV-CS that exploit in a different manner the concept of gradient in order to improve the solution of the inverse scattering problems obtained by TV-CS approach. Numerical examples are addressed to show the effectiveness of the method.


Mathematics ◽  
2021 ◽  
Vol 9 (17) ◽  
pp. 2087
Author(s):  
Won-Kwang Park

In this study, we consider a sampling-type algorithm for the fast localization of small electromagnetic inhomogeneities from measured far-field pattern data in the limited-aperture inverse scattering problem. For this purpose, we designed an indicator function based on the structure of left- and right-singular vectors of a multistatic response matrix, the elements of which were measured far-field pattern data. We then rigorously investigated the mathematical structure of the indicator function in terms of purely dielectric permittivity and magnetic permeability contrast cases by establishing a relationship with an infinite series of Bessel functions of an integer order of the first kind and a range of incident and observation directions before exploring various intrinsic properties of the algorithm, including its feasibility and limitations. Simulation results with synthetic data corrupted by random noise are presented to support the theoretical results.


Sign in / Sign up

Export Citation Format

Share Document