Numerical Solution of a Coefficient Inverse Problem with Multi-Frequency Experimental Raw Data by a Globally Convergent Algorithm

Dinh-Liem Nguyen,Michael V. Klibanov,Loc H. Nguyen,Aleksandr E. Kolesov,Michael A. Fiddy,Hui Liu
DOI: https://doi.org/10.1016/j.jcp.2017.05.015
2016-12-04
Abstract:We analyze in this paper the performance of a newly developed globally convergent numerical method for a coefficient inverse problem for the case of multi-frequency experimental backscatter data associated to a single incident wave. These data were collected using a microwave scattering facility at the University of North Carolina at Charlotte. The challenges for the inverse problem under the consideration are not only from its high nonlinearity and severe ill-posedness but also from the facts that the amount of the measured data is minimal and that these raw data are contaminated by a significant amount of noise, due to a non-ideal experimental setup. This setup is motivated by our target application in detecting and identifying explosives. We show in this paper how the raw data can be preprocessed and successfully inverted using our inversion method. More precisely, we are able to reconstruct the dielectric constants and the locations of the scattering objects with a good accuracy, without using any advanced \emph{a priori} knowledge of their physical and geometrical properties.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: **How to reconstruct the dielectric constant and position of scattering objects from multi - frequency experimental raw data without using detailed prior knowledge of these objects**. Specifically, the paper focuses on a Coefficient Inverse Problem (CIP), that is, reconstructing the physical and geometric properties of three - dimensional objects given multi - frequency experimental backscattering data, without relying on any detailed prior information about these objects. ### Problem Background 1. **High Non - linearity and Severe Ill - posedness**: The coefficient inverse problem is highly non - linear and severely ill - posed, which means that the solving process is very complex and easily affected by noise. 2. **Limited Data Volume and High Noise**: The experimental data is collected in a non - ideal experimental environment, with a limited data volume and a large amount of noise. 3. **Application Scenarios**: The main application motivation of this research is to detect and identify explosives, such as anti - personnel mines and Improvised Explosive Devices (IEDs). In addition, this method can also be applied in the fields of non - destructive testing and material characterization. ### Core Contributions of the Paper 1. **Global Convergence Algorithm**: The paper proposes a new global convergence numerical method (GCM), which can provide a good approximate solution to the coefficient inverse problem without using any advanced prior knowledge. 2. **Multi - frequency Experimental Data Processing**: For the first time, the performance of the new GCM on multi - frequency experimental raw data is studied, and a new data pre - processing procedure is developed to make the raw data closer to the computational simulation data. 3. **No Special Experimental Conditions Required**: The experimental data is collected in an ordinary room, containing many interference signals (such as office furniture, computers, Wi - Fi signals, etc.), which increases the difficulty of the problem. ### Mathematical Model The mathematical model in the paper is based on the Helmholtz equation: \[ \Delta u + k^2\varepsilon_r(x)u = 0, \quad x\in\mathbb{R}^3 \] where \(u(x,k)\) is the total field, \(k=\frac{\omega}{c}\) is the wave number, and \(\varepsilon_r(x)\) is the relative dielectric constant. ### Data Pre - processing To deal with the noise problem, the paper develops a new data pre - processing method, aiming to extract the reflection signals of the target objects from the raw data, making it closer to the computational simulation data. This process includes: - Calibrating the internal reflection signals - Removing the interference of irrelevant signals ### Results Through this method, the paper shows that the dielectric constant and position of scattering objects can be accurately reconstructed from multi - frequency experimental raw data without using any detailed prior knowledge. This provides important reference information for future classification algorithms, especially in the radar community, where combining the dielectric constant, shape and other parameters can better distinguish explosives from clutter. In conclusion, this paper solves the challenge of reconstructing the characteristics of scattering objects from multi - frequency experimental raw data, especially for the application scenarios of detecting and identifying explosives.