Linear and Nonlinear Inverse Scattering

In this paper we discuss one-dimensional scattering and inverse scattering for the Helmholtz equation on the half-line from the point of view of the layer stripping. By full or nonlinear scattering, we mean the transformation between the index of refraction (actually half of its logarithmic derivative) and the reflection coefficient. We refer to this mapping as nonlinear scattering, because the mapping itself is nonlinear. Another appropriate name is multiple scattering, as this model includes the effects of multiple reflections. By linear scattering we mean the Born, or single scattering, approximation. This is the Frechet derivative of the full scattering transform at the constant index of refraction, which can be calculated to be exactly the Fourier transform. In [J. Sylvester, D. P. Winebrenner, and F. Gylys-Colwell, SIAM J. Appl. Math., 56 (1996), pp. 736-754], we introduced a variant of layer stripping based on causality and the Riesz transform, rather than on trace formulas-see [A. Brickstein and T. Kailath, SIAM Rev., 29 (1987), pp. 359-389], [Y. Chen and V. Rokhlin, Inverse Problems, 8 (1992), pp. 365-390], [J. P. Corones, R. J. Krueger, and M. E. Davison, J. Acoust. Soc. Amer., 74 (1983), pp. 1535-1541], or [W. W. Symes, J. Math. Anal. Appl., 94 (1983), pp. 435-453] for other approaches to layer stripping. A by-product of our layer stripping formalism was the discovery of a nonlinear Plancherel equality, which plays a role in the analysis of the inverse scattering problem, analogous to that played by the linear Plancherel equality in developing the theory of the Fourier transform. That linear-nonlinear analogy sets the theme for this work. In the next section, we review the main results of our article cited above, including a brief derivation of the nonlinear Plancherel equality. We show how this equality suggests a natural metric for measuring the distance between reflection coefficients and show that the scattering and inverse scattering maps become homeomorphisms when we use this metric. In section 2, we exhibit a nonlinear Riesz transform, which plays for the nonlinear inverse scattering problem the same role in enforcing causality as the linear Riesz transform plays in signal processing. (The Riesz transform is perhaps less well known to the signal processing community than the Hilbert transform, which is the principal value integral Hf(ω ) = 1/π∫ f(η)/η -ωdη . The Riesz transform is a projection built from the Hilbert transform via the form