An inverse random source problem for the Helmholtz equation

This paper is concerned with an inverse random source problem for the one-dimensional stochastic Helmholtz equation, which is to reconstruct the statistical properties of the random source function from boundary measurements of the radiating random electric field. Although the emphasis of the paper is on the inverse problem, we adapt a computationally more efficient approach to study the solution of the direct problem in the context of the scattering model. Specifically, the direct model problem is equivalently formulated into a two-point spatially stochastic boundary value problem, for which the existence and uniqueness of the pathwise solution is proved. In particular, an explicit formula is deduced for the solution from an integral representation by solving the two-point boundary value problem. Based on this formula, a novel and efficient strategy, which is entirely done by using the fast Fourier transform, is proposed to reconstruct the mean and the variance of the random source function from measurements at one boundary point, where the measurements are assumed to be available for many realizations of the source term. Numerical examples are presented to demonstrate the validity and effectiveness of the proposed method.

[1]  A. Devaney,et al.  Nonuniqueness in inverse source and scattering problems , 1982 .

[2]  Nicolas Valdivia,et al.  Acoustic source identification using multiple frequency information , 2009 .

[3]  Daniel Ocone,et al.  Linear stochastic differential equations with boundary conditions , 1989 .

[4]  Vladimir Rokhlin,et al.  On the inverse scattering problem for the Helmholtz equation in one dimension , 1992 .

[5]  T. Faniran.  Numerical Solution of Stochastic Differential Equations , 2015 .

[6]  Thomas William MacDowell.  Boundary value problems for stochastic differential equations , 1968 .

[7]  George Papanicolaou,et al.  Wave Propagation in a One-Dimensional Random Medium , 1971 .

[8]  Akira Ishimaru,et al.  Wave propagation and scattering in random media , 1997 .

[9]  Rafael Yuste,et al.  Fluorescence microscopy today , 2005, Nature Methods.

[10]  M. P. Levin,et al.  Numerical Recipes In Fortran 90: The Art Of Parallel Scientific Computing , 1998, IEEE Concurrency.

[11]  Junshan Lin,et al.  A multi-frequency inverse source problem , 2010 .

[12]  George Papanicolaou,et al.  Fine scale uncertainty in parameter estimation for elliptic equations , 2009 .

[13]  H. Keller.  Numerical Solution of Two Point Boundary Value Problems , 1976 .

[14]  Guillaume Bal,et al.  Central Limits and Homogenization in Random Media , 2007, Multiscale Model. Simul..

[15]  R. Kress,et al.  Inverse Acoustic and Electromagnetic Scattering Theory , 1992 .

[16]  A. Devaney,et al.  The inverse source problem of electromagnetics: linear inversion formulation and minimum energy solution , 1999 .

[17]  Faming Liang,et al.  Statistical and Computational Inverse Problems , 2006, Technometrics.

[18]  Guillaume Bal,et al.  Physics-based models for measurement correlations: application to an inverse Sturm–Liouville problem , 2009 .

[19]  Werner Römisch,et al.  Numerical Solution of Stochastic Differential Equations (Peter E. Kloeden and Eckhard Platen) , 1995, SIAM Rev..

[20]  Anthony J. Devaney,et al.  Inverse problem for random sources (A) , 1978 .

[21]  William H. Press,et al.  Numerical recipes in Fortran 90: the art of parallel scientific computing, 2nd Edition , 1996, Fortran numerical recipes.

[22]  GoldfarbDonald,et al.  Bregman Iterative Algorithms for $\ell_1$-Minimization with Applications to Compressed Sensing , 2008 .

[23]  Shui-Nee Chow,et al.  Numerical solution of an inverse medium scattering problem with a stochastic source , 2010 .

[24]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

[25]  Cao,et al.  FINITE ELEMENT AND DISCONTINUOUS GALERKIN METHOD FOR STOCHASTIC HELMHOLTZ EQUATION IN TWO- AND THREE-DIMENSIONS , 2008 .

[26]  Anthony J. Devaney,et al.  Inverse Source Problem in Nonhomogeneous Background Media , 2007, SIAM J. Appl. Math..

[27]  Josselin Garnier,et al.  Mathematical and Statistical Methods for Imaging , 2011 .

[28]  Desmond J. Higham,et al.  An Algorithmic Introduction to Numerical Simulation of Stochastic Differential Equations , 2001, SIAM Rev..

[29]  Ali Adibi,et al.  Model for efficient simulation of spatially incoherent light using the Wiener chaos expansion method. , 2007, Optics letters.

[30]  G. M.,et al.  Partial Differential Equations I , 2023, Applied Mathematical Sciences.

[31]  G. Choe.  Numerical Solution of Stochastic Differential Equations , 2016 .

[32]  I. Gel'fand,et al.  On the determination of a differential equation from its spectral function , 1955 .

[33]  Peijun Li,et al.  An inverse random source scattering problem in inhomogeneous media , 2011 .

[34]  Steven A. Orszag,et al.  CBMS-NSF REGIONAL CONFERENCE SERIES IN APPLIED MATHEMATICS , 1978 .

[35]  S. Shreve,et al.  Stochastic differential equations , 1955, Mathematical Proceedings of the Cambridge Philosophical Society.

[36]  J. Craggs.  Applied Mathematical Sciences , 1973 .

[37]  V. Naroditsky.  Wave Propagation in Random Media , 1989 .

[38]  R. Dobarzić,et al.  [Fluorescence microscopy]. , 1975, Plucne bolesti i tuberkuloza.

[39]  Wotao Yin,et al.  Bregman Iterative Algorithms for (cid:2) 1 -Minimization with Applications to Compressed Sensing ∗ , 2008 .

[40]  K. Abdella.  Numerical Solution of two-point boundary value problems using Sinc interpolation , 2012 .

[41]  Junshan Lin,et al.  An inverse source problem with multiple frequency data , 2011 .