Skip to content

DSM for solving ill-conditioned linear algebraic systems

A. G. Ramm

math.NAarXiv:math/0601299

Abstract

A standard way to solve linear algebraic systems Au=f,\,\,(*) with ill-conditioned matrices A is to use variational regularization. This leads to solving the equation (A*A+aI)u=A*f, where a is a regularization parameter, and f are noisy data, ||f-f||≤ . Numerically it requires to calculate products of matrices A*A and inversion of the matrix A*A+aI which is also ill-conditioned if a>0 is small. We propose a new method for solving (*) stably, given noisy data f. This method, the DSM (Dynamical Systems Method) is developed in this paper for selfadjoint A. It consists in solving a Cauchy problem for systems of ordinary differential equations.

Create a lesson