+7 (495) 957-77-43
Контакты

T-Comm_Article 8_10_2021

Извините, этот техт доступен только в “Американский Английский”. For the sake of viewer convenience, the content is shown below in the alternative language. You may click the link to switch the active language.

THE CALCULATION OF THE SOLUTION OF MULTIDIMENSIONAL INTEGRAL EQUATIONS WITH METHODS MONTE CARLO AND QUASI-MONTE CARLO

Abas Wisam Mahdi Abas,  South Russian State Polytechnic University
(Novocherkassk Polytechnic Institute), Novocherkassk, Russia, abas.wisam.82@mail.ru

Abstract
The article considers an approach based on the random cubature method for solving both single and multidimensional singular integral equations, Volterra and Fredholm equations of the 1st kind, for ill-posed problems in the theory of integral equations, etc. A variant of the quasi-Monte Carlo method is studied. The integral in an integral equation is approximated using the traditional Monte Carlo method for calculating integrals. Multidimensional interpolation is applied on an arbitrary set of points. Examples of applying the method to a one-dimensional integral equation with a smooth kernel using both random and low-dispersed pseudo-random nodes are considered. A multidimensional linear integral equation with a polynomial kernel and a multidimensional nonlinear problem – the Hammerstein integral equation – are solved using the Newton method. The existence of several solutions is shown. Multidimensional integral equations of the first kind and their solution using regularization are considered. The Monte Carlo and quasi-Monte Carlo methods have not been used to solve such problems in the studied literature. The Lavrentiev regularization method was used, as well as random and pseudo-random nodes obtained using the Halton sequence. The problem of eigenvalues is solved. It is established that one of the best methods considered is the Leverrier-Faddeev method. The results of solving the problem for a different number of quadrature nodes are presented in the table. An approach based on parametric regularization of the core, an interpolation-projection method, and averaged adaptive densities are studied. The considered methods can be successfully applied in solving spatial boundary value problems for areas of complex shape. These approaches allow us to expand the range of problems in the theory of integral equations solved by Monte Carlo and quasi-Monte Carlo methods, since there are no restrictions on the value of the norm of the integral operator. A series of examples demonstrating the effectiveness of the method under study is considered.

Keywords: integral equation, high dimension, Monte Carlo method

References

  1. Ermakov S. M. (2011). The Monte Carlo Method in computational mathematics (introductory course). Saint Petersburg: Binom Publishing House. 192 p.
  2. Ivanov V. M., Kulchitsky O. Yu., Korenevsky M. L. (1998). Combined method for solving integral equations. Differential equations and control processes. No. 1. Pp. 1-40.
  3. Sipin A. S. (2016). Grid-Free Monte Carlo methods for solving boundary value problems for partial differential equations. The author’s abstract Diss. for the degree of doctor of physics.- math. sciences’. Saint Petersburg. 32 p.
  4. Ivanov V. M., Kulchitsky O. Yu. (1990). Method of numerical solution of integral equations on a random grid. Differents. Equations. Vol. 26, No. 2. Pp. 333-341.
  5. Berkovsky N. A. (2006). Modernization of the semi-statistical method for numerical solution of integral equations. The author’s abstract Diss. for the degree of candidate of physical Sciences. — math. sciences’. Saint Petersburg. 15 p.
  6. Korenevsky M. L. (2000). Development of adaptive statistical methods for calculating certain integrals. Diss. for the degree of candidate of physical Sciences. — math. sciences’. Saint Petersburg. 161 p.
  7. Skorikova O. V. (2013). Strong uniform distribution of the system of van der Korput-Hemmersley functions. Proceedings of the Tula state University. Natural science. Issue 3. Pp. 91-102.
  8. Antonov A. A., Ermakov S. M. (2014). Empirical estimation of integration error by the quasi-Monte Carlo method. Vestnik of Saint Petersburg University. Series 1. Mathematics. Mechanics. Astronomy. Issue #1. Vol. 1. Pp. 3-11.
  9. Antonov A. A. (2015). Algorithm for numerical integration by the quasi-Monte Carlo method with a posteriori error estimation. Bulletin of the Saint Petersburg University. Series 1. Mathematics. Mechanics. Astronomy. Issue no. 1. Vol. 2. Pp. 3-11.
  10. Stoyantsev V. T. (1973). Solution of the Cauchy problem for a parabolic equation by the quasi-Monte Carlo method. mod. and math. Phys. Vol. 13. No. 5. Pp. 1153-1160.
  11. Stoyantsev V. T. (1975). Solution of the Dirichlet problem by the quasi-Monte Carlo method. UMN. Vol. 30. no. 1 (181). Pp. 263-264.
  12. Korn G., Korn T. (1973). Handbook of mathematics of engineers. Moscow: Nauka. 832 p.
  13. Nekrasov S. A., Tkachev A. N. (2007). Probability Theory and its applications: Textbook. South-Russian state technical University. — Novocherkassk: YURSTU. 148 p.
  14. Nekrasov S. A. (2008). Methods of accelerated statistical modeling and their application in electrical problems. Izv.vuzov. Electromechanics. No 5. Pp. 13-19.
  15. Nekrasov S. A. (2016). Solution of integral equations by the Monte Carlo method. NovaInfo.Ru (Electronic journal). 2016. No. 53.
  16. Nekrasov S. A. (2016). Solution of the n-dimensional Schrodinger equation by the method of integral equations on a pseudo-random grid. NovaInfo.Ru (Electronic journal). No. 55.

Information about author:

Abas Wisam Mahdi Abas, Post-graduate student of the South Russian State Polytechnic University (Novocherkassk Polytechnic Institute), Department of Applied Mathematics, Novocherkassk, Russia