In this paper, we study the inverse scattering problem for a class of signals that have a compactly supported reflection coefficient. The problem boils down to the solution of the Gelfand–Levitan–Marchenko (GLM) integral equations with a kernel that is bandlimited. By adopting a sampling theory approach to the associated Hankel operators in the Bernstein spaces, a constructive proof of existence of a solution of the GLM equations is obtained under various restrictions on the nonlinear impulse response (NIR). The formalism developed in this article also lends itself well to numerical computations yielding algorithms that are shown to have algebraic rates of convergence. In particular, the use Whittaker–Kotelnikov–Shannon sampling series yields an algorithm that converges as O(N − 1 / 2 ) whereas the use of Helms and Thomas (HT) version of the sampling expansion yields an algorithm that converges as O(N − m − 1 / 2 ) for any m > 0 provided the regularity conditions are fulfilled. The complexity of the algorithms depend on the linear solver used. The use of conjugate-gradient (CG) method yields an algorithm of complexity O(N iter. N 2 ) per sample of the signal where N is the number of sampling basis functions used and N iter. is the number of CG iterations involved. The HT version of the sampling expansions facilitates the development of algorithms of complexity O(N iter. N log N) (per sample of the signal) by exploiting the special structure as well as the (approximate) sparsity of the matrices involved. The algorithms are numerically validated using Schwartz class functions as NIRs that are either bandlimited or effectively bandlimited. The results suggest that the HT variant of our algorithm is spectrally convergent for an input of the aforementioned class.
|Number of pages||33|
|Journal||Journal of Physics A: Mathematical and Theoretical|
|Publication status||Published - 2019|
- Bandlimited signals
- Gelfand–Levitan–Marchenko integral equations
- Inverse scattering