Sampling and reconstruction of solutions to the Helmholtz equation

We consider the inverse problem of reconstructing general solutions to the Helmholtz equation on some domain $\Omega$ from their values at scattered points $x_1,\dots,x_n\subset \Omega$. This problem typically arises when sampling acoustic fields with $n$ microphones for the purpose of reconstructing this field over a region of interest $\Omega$ contained in a larger domain $D$ in which the acoustic field propagates. In many applied settings, the shape of $D$ and the boundary conditions on its border are unknown. Our reconstruction method is based on the approximation of a general solution $u$ by linear combinations of Fourier-Bessel functions or plane waves. We analyze the convergence of the least-squares estimates to $u$ using these families of functions based on the samples $(u(x_i))_{i=1,\dots,n}$. Our analysis describes the amount of regularization needed to guarantee the convergence of the least squares estimate towards $u$, in terms of a condition that depends on the dimension of the approximation subspace, the sample size $n$ and the distribution of the samples. It reveals the advantage of using non-uniform distributions that have more points on the boundary of $\Omega$. Numerical illustrations show that our approach compares favorably with reconstruction methods using other basis functions, and other types of regularization.

[1]  Ralf Hiptmair,et al.  Error analysis of Trefftz-discontinuous Galerkin methods for the time-harmonic Maxwell equations , 2011, Math. Comput..

[2]  Albert Cohen,et al.  On the Stability and Accuracy of Least Squares Approximations , 2011, Foundations of Computational Mathematics.

[3]  Michael P. Friedlander,et al.  Probing the Pareto Frontier for Basis Pursuit Solutions , 2008, SIAM J. Sci. Comput..

[4]  L. Daudet,et al.  Low-complexity computation of plate eigenmodes with Vekua approximations and the method of particular solutions , 2013, 1301.0908.

[5]  Ralf Hiptmair,et al.  Vekua theory for the Helmholtz operator , 2011 .

[6]  M. Abramowitz,et al.  Handbook of Mathematical Functions With Formulas, Graphs and Mathematical Tables (National Bureau of Standards Applied Mathematics Series No. 55) , 1965 .

[7]  Laurent Daudet,et al.  Plate impulse response spatial interpolation with sub-Nyquist sampling , 2011 .

[8]  Alan B. Tayler,et al.  New methods for solving elliptic equations , 1969 .

[9]  A. Moiola Plane wave approximation in linear elasticity , 2013 .

[10]  Ralf Hiptmair,et al.  Plane wave approximation of homogeneous Helmholtz solutions , 2011 .

[11]  Y. C. Pati,et al.  Orthogonal matching pursuit: recursive function approximation with applications to wavelet decomposition , 1993, Proceedings of 27th Asilomar Conference on Signals, Systems and Computers.

[12]  Irene A. Stegun,et al.  Handbook of Mathematical Functions. , 1966 .

[13]  H. Rauhut,et al.  Interpolation via weighted $l_1$ minimization , 2013, 1308.0759.

[14]  Emmanuel Perrey-Debain,et al.  Plane wave decomposition in the unit disc: convergence estimates and computational aspects , 2006 .