Abstract:We propose the generalized quadrature methods for numerical solution of singular integral equation of Abel type. We overcome the singularity using the analytical calculation of the singular integral expression. The problem of solution of singular integral equation is reduced to nonsingular system of linear algebraic equations without shift meshes techniques employment. We also propose generalized quadrature method for solution of Abel equation using the singular integral. Relaxed errors bounds are derived. In order to improve the accuracy we use Tikhonov regularization method. We demonstrate the efficiency of proposed techniques on infrared tomography problem. Numerical experiments show that it make sense to apply regularization in case of highly noisy sources only. That is due to the fact that singular integral equations enjoy selfregularization property.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to propose a new numerical method for solving Abel - type singular integral equations (SIEs) and apply it to infrared tomography problems. Specifically, the author aims to overcome the singularity in the singular integral equation through the generalized quadrature formula and transform the problem of solving the singular integral equation into a non - singular linear algebraic equation system. In addition, the author also explores the effectiveness of using the Tikhonov regularization method to improve the solution accuracy in the case of highly noisy data.
### Background and Problem Description of the Paper
Abel - type singular integral equations are widely used in practical models, including plasma diagnosis, thermal tomography, X - ray CT, spectroscopy, galaxy cluster astrophysics, and other fields. Such equations usually have the following form:
\[
k(r)=-\frac{1}{\pi}\int_{r}^{R}\frac{q^{\prime}(x)}{\sqrt{x^{2}-r^{2}}}\,dx,\quad0\leq r\leq R
\]
Here, \(k(r)\) is the function to be determined, and \(q(x)\) is the source function. Since experimental data usually contains noise, directly solving the above equation will encounter difficulties, especially when it comes to singularity. Therefore, new numerical methods need to be developed to effectively handle these problems.
### Main Contributions
1. **Generalized Quadrature Formula**: The author proposes a numerical method based on the generalized left - rectangle formula for solving Abel - type singular integral equations. This method avoids the problem of introducing grid offset in traditional methods by analytically calculating the singular integral.
2. **Error Estimation**: The author derives the estimation formula for the quadrature error and takes into account the measurement error of the right - hand side term \(q(x)\). This helps to evaluate the accuracy and stability of the numerical solution.
3. **Tikhonov Regularization**: To improve the stability of the solution, especially in the case of highly noisy data, the author applies the Tikhonov regularization method. Research shows that in the case of a relatively high noise level (about 10%), the regularization method can significantly improve the solution results.
4. **Infrared Tomography Application**: The author applies the proposed method to infrared tomography problems, demonstrating its effectiveness in practical applications. In particular, by using spline smoothing techniques to reduce the grid step size and moderately smooth the fluctuations in the measured data, the solution accuracy is further improved.
### Conclusion
This paper proposes a new generalized quadrature method that can effectively solve Abel - type singular integral equations and shows good performance in practical applications such as infrared tomography. The research results show that for cases with a low noise level and a small grid step size, the singular integral equation has self - regularization characteristics and does not require additional regularization processing; while in the case of high noise and a coarse grid, the Tikhonov regularization method can significantly improve the stability of the solution.
### References
[1] E.O. Åkesson, K.J. Daun, Parameter selection methods for axisymmetric flame tomography through Tikhonov regularization, Appl. Opt., 47 (3) (2008) 407–416.
[2] R.S. Anderssen, F.R. de Hoog, Abel integral equations, in Numerical Solution of Integral Equations, edt. M.A.Goldberg, Plenum Press, 1990, 373–410.
[3] S.M. Belotserkovskii, I.K. Lifanov, Method of Discrete Vortices, CRC Press, 1993.
[4] I.V. Boikov, N.Yu. Kudryashova, Approximate methods for singular integral