Numerical solution of first-kind Volterra equations by sequential Tikhonov regularization

Patricia K. Lamm
Department of Mathematics
Michigan State University
E. Lansing, MI 48824-1027

Lars Eldén
Department of Mathematics
Linkoping University
S-581 83 Linkoping, Sweden

SIAM J. Numerical Analysis 34 (1997) pp 1432-1450.


Abstract:

We consider the problem of finding regularized solutions to ill-posed Volterra integral equations. The method we consider is a sequential form of Tikhonov regularization that is particularly suited to problems of Volterra type. We prove that when this sequential regularization method is coupled with several standard discretizations of the integral equation (collocation, rectangular and midpoint quadrature), one obtains convergence of the method at an optimal rate with respect to noise in the data. In addition we describe a fast algorithm for the implementation of sequential Tikhonov regularization, and show that for small values of the regularization parameter, the method is only slightly more expensive computationally than the numerical solution of the original unregularized integral equation. Finally, numerical results are presented to show that the performance of sequential Tikhonov regularization nearly matches that of standard Tikhonov regularization in practice, but at considerable savings in cost.


Text of paper:
Contact: lamm@math.msu.edu