The inverse problem about recovering a parameter function by measurements of solutions of the system partial differential equations is considered. A typical formulation of this inverse problem consists minimization of the sum of a data fitting error term and a regularization term, subject to the forward problem being satisfied. The problem is typically ill-posed without regularization and it is ill-conditioned with it, since the regularization term is aimed at removing noise without overshadowing the data.
Let the forward problem be a linear elliptic differential equation where refers a differential operator depending on a parameter vector function , defined on an appropriate domain and equipped with suitable boundary conditions. The discrimination of this problem is studied and for regularization the Tikhonov method with introducing the Lagrangian approach is applied. Finally the problem is numerically solved by the Gauss-Newton method, and a preconditioned conjugate gradient algorithm is applied at each iteration for the resulting reduced Hessian system. Alternatively, a preconditioned Krylov method is applied to arising system.
The considered problem the arises in many applications. The results are illustrated by different computational examples.