Subspace linear inverse method

Publisher Website
Google Scholar
Abstract
Abslract. This paper presents a robust, flexible and efficient algorithm to solve large scale linear inverse problems. The method is iterative and at each iteration a perturbation in a q-dimensional subspace of an M-dimensional model space is sought. The basis vectors for the subspace are primarily steepest descent vectors obtained from segmenting the data misfit and model objective functions. The efficiency of the algorithm is realized because only a y x y matrix needs to be inverted at each iteration instead of a matrix of order M. As M becomes large the number of computations per iteration is of order yNM where N is the number of data. An important feature of our algorithm is that positivity can easily be incorporated into the solution. We do this by introducing n two-segment mapping which transforms positive parameters to parameters defined on the rea1 line. The nonlinear mapping requires that a line semh involving forward modclling is implemented so that at each iteration we obtain a model which misfits the data to a predetermined level. This obviates the need to carry out additional inversions with trial and error selection of a Lagrange multiplier. In this paper we present the details of the subspace algorithm and explore the effect on convergence of using different suategies for selecting basis vectors and altering adjustable panmeters which control the rate of decrease in the misfit and rate of increase in the model norm as a function of iteration number.

This publication has 4 references indexed in Scilit: