A STOPPING RULE FOR THE CONJUGATE GRADIENT REGULARIZATION METHOD APPLIED TO INVERSE PROBLEMS IN ACOUSTICS

Abstract

We present a novel parameter choice strategy for the conjugate gradient regularization algorithm which does not assume a priori information about the magnitude of the measurement error. Our approach is to regularize within the Krylov subspaces associated with the normal equations. We implement conjugate gradient via the Lanczos bidiagonalization process with reorthogonalization, and then we construct regularized solutions using the SVD of a bidiagonal projection constructed by the Lanczos process. We compare our method with the one proposed by Hanke and Raus and illustrate its performance with numerical experiments, including detection of acoustic boundary vibrations.