An Algorithm for Improving the Condition Number of Matrices and Its Application for Solving the Inverse Problems of Gravimetry and Magnetometry

Alexander Leonov, Dmitry Lukyanenko (), Anatoly Yagola and Yanfei Wang
Additional contact information
Alexander Leonov: Department of High Mathematics, National Research Nuclear University MEPhI, Moscow 115409, Russia
Dmitry Lukyanenko: Department of Mathematics, Faculty of Physics, Lomonosov Moscow State University, Moscow 119991, Russia
Anatoly Yagola: Department of Mathematics, Faculty of Physics, Lomonosov Moscow State University, Moscow 119991, Russia
Yanfei Wang: Innovation Academy for Earth Science, Chinese Academy of Sciences, Beijing 100029, China

Mathematics, 2025, vol. 13, issue 8, 1-15

Abstract: The paper considers one of the possible statements of inverse problems in gravimetric and magnetometric remote sensing, proposes a new approach to its solution and formulates algorithms that implement this approach. The problem under consideration consists of finding hypothetical sources of the corresponding potential fields at a given depth based on these fields measured on the Earth’s surface. The problem is reduced to solving systems of linear algebraic equations (SLAE) with ill-conditioned matrices. The proposed approach to the numerical solution is based on improving the condition number of the SLAE’s matrix. A numerical algorithm implementing the proposed approach that is applicable to the stable solution of degenerate and ill-conditioned SLAEs with an approximately given right-hand side is formulated in general form. The algorithm uses the SVD decomposition of the SLAE’s matrix and constructs a new matrix close to the original one with a better (smaller) condition number. An approximate solution to the original SLAE is calculated using the pseudoinverse of the new matrix. The results of a theoretical study of the algorithm are presented and the main properties of the new matrix are given. In particular, the reduction of its condition number is estimated. Several implementations of this algorithm are considered, in particular, the MPMI method, which is based on the use of so-called minimal pseudoinverse matrices. For the model problem, the advantage of the MPMI method over a number of other common methods is shown. The MPMI method is applied to solve the considered problems of gravity exploration and magnetic exploration both in the separate solution of these inverse problems and in their joint solution when processing geophysical data for the Kathu region, in the Northern Cape area of South Africa.

Keywords: gravity and magnetic exploration; inverse problem; minimal pseudoinverse matrix method; improvement of condition number; gravity–magnetic joint inversion (search for similar items in EconPapers)
JEL-codes: C (search for similar items in EconPapers)
Date: 2025
References: Add references at CitEc
Citations:

Downloads: (external link)
https://www.mdpi.com/2227-7390/13/8/1280/pdf (application/pdf)
https://www.mdpi.com/2227-7390/13/8/1280/ (text/html)

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:gam:jmathe:v:13:y:2025:i:8:p:1280-:d:1633880

Access Statistics for this article

Mathematics is currently edited by Ms. Emma He

More articles in Mathematics from MDPI
Bibliographic data for series maintained by MDPI Indexing Manager ().