Maximum likelihood estimation with side information of a 1-D discrete layered medium from its noisy impulse reflection response

We consider the problem of computing the maximum likelihood estimates of th e reflection coefficients of a discrete 1-D layered medium from noisy obser vations of its impulse reflection response. We have side information in tha t a known subset of the reflection coefficients are known to be zero; this knowledge could come from either a priori knowledge of a homogeneous subreg ion inside the scattering medium or from a thresholding operation in which noisy reconstructed reflection coefficients with absolute values below a th reshold are known to be zero. Our procedure converges in one or two iterati ons, each of which requires only setting up and solving a small system of l inear equations and running the Levinson algorithm, Numerical examples are provided that demonstrate not only the operation of the algorithm but also that the side information improves the reconstruction of unconstrained refl ection coefficients as well as constrained ones due to the nonlinearity of the inverse scattering problem.