A Large-Scale Trust-Region Approach to the Regularization of Discrete Ill-Posed Problems
This work was also published as a Rice University thesis/dissertation: http://hdl.handle.net/1911/19422
We consider the problem of computing the solution of large-scale discrete ill-posed problems when there is noise in the data. These problems arise in important areas such as seismic inversion, medical imaging and signal processing. We pose the problem as a quadratically constrained least squares problem and develop a method for the solution of such problem. Our method does not require factorization of the coefficient matrix, it has very low storage requirements and handles the high degree of singularities arising in discrete ill-posed problems. We present numerical results on test problems and an application of the method to a practical problem with real data.
Citable link to this pagehttp://hdl.handle.net/1911/101904
MetadataShow full item record
- CAAM Technical Reports