Next Article in Journal
Nearly Cosymplectic Manifolds of Constant Type
Next Article in Special Issue
On the Bias in Confirmatory Factor Analysis When Treating Discrete Variables as Ordinal Instead of Continuous
Previous Article in Journal
Towards Predictive Vietnamese Human Resource Migration by Machine Learning: A Case Study in Northeast Asian Countries
Previous Article in Special Issue
Some Notes for Two Generalized Trigonometric Families of Distributions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Approximate Methods for Calculating Singular and Hypersingular Integrals with Rapidly Oscillating Kernels

1
Department of Higher and Applied Mathematics, Penza State University, Penza 440026, Russia
2
Department of Computational Physics, Saint Petersburg State University, 7/9 Universitetskaya Emb., Saint Petersburg 199034, Russia
*
Author to whom correspondence should be addressed.
Axioms 2022, 11(4), 150; https://doi.org/10.3390/axioms11040150
Submission received: 30 January 2022 / Revised: 13 March 2022 / Accepted: 23 March 2022 / Published: 24 March 2022
(This article belongs to the Collection Mathematical Analysis and Applications)

Abstract

:
The article is devoted to the issue of construction of an optimal with respect to order passive algorithms for evaluating Cauchy and Hilbert singular and hypersingular integrals with oscillating kernels. We propose a method for estimating lower bound errors of quadrature formulas for singular and hypersingular integral evaluation. Quadrature formulas were constructed for implementation of the obtained estimates. We constructed quadrature formulas and estimated the errors for hypersingular integrals with oscillating kernels. This method is based on using similar results obtained for singular integrals.

1. Introduction

Recent years have shown the importance of evaluating singular and hypersingular integrals with rapidly oscillating kernels in mathematical modeling of wave processes in many areas of physics and technology: electrodynamics (waveguides, gyrotrons), aerodynamics, geophysics (transformation of gravity and magnetic fields), etc.
Today, there are very few manuscripts devoted to approximate methods for evaluating singular integrals with rapidly oscillating kernels. We are unaware of papers dealing with approximate methods for evaluating hypersingular integrals with rapidly oscillating kernels.
In this paper, we construct an optimal with respect to order quadrature formulas for evaluating singular and hypersingular integrals on Hölder functions and differentiable function classes.
The paper is organized as follows:
Section 1 contains the review of publications evaluating singular integrals with rapidly oscillating functions. In this section, we give the definitions of singular and hypersingular integrals and optimal algorithms for their evaluation.
In Section 2, Levin’s method is extended to singular and hypersingular integrals.
In Section 3, we introduce an optimal with respect to order quadrature formulas for calculating singular integrals with oscillating functions.
In Section 4, we present methods for evaluating the hypersingular integrals with rapidly oscillating functions.
In Section 5, we give the conclusions of our study.

1.1. Literature Review

An extensive literature exists regarding approximate methods for singular and hypersingular integral evaluation. Detailed reviews are given in [1,2,3,4,5,6,7,8]. Very few publications concentrate on approximate methods for evaluating singular integrals with oscillating kernels [9,10,11,12,13,14,15,16,17].
Here, we give a brief overview of the manuscripts considering singular integrals with oscillating kernels.
The chapter “Oscillatory Singular Integrals” in [18] is devoted to the study of oscillatory singular integrals. The authors consider singular integrals of the form
T f ( x ) = p . v . R n e i p ( x , y ) K ( x y ) f ( y ) d y ,
where
(i)
K is a C 1 function away from the origin;
(ii)
K is homogeneous of degree -n;
(iii)
the mean value of K on the unit sphere vanishes;
(iiii)
p ( x , y ) is a real-value polynomial on R n × R n .
The boundedness of the operator T is investigated in a number of function spaces. Integrals of the form (1) are widely used in the Radon transform.
The paper [11] deals with approximate methods for evaluating the integral
I ( f , w ) = 1 1 f ( x ) x cos w x d x ,
where f ( x ) is an analytic function in [ 1 ; 1 ] , w R   \   { 0 } . The integral (2) is converted to the form
I ( f , w ) = 1 1 cos w x 1 x f ( x ) d x + 1 1 f ( x ) x d x = I 0 ( f ) + I ( f ) .
Philon’s method is used for the integral I 0 ( f ) .
The second method proposed in [11] consists in approximating the function f ( x ) in I 0 ( f ) by segment of the Taylor series. The well-known methods are used to approximate the integral I ( f ) [19].
In [20], Levin’s method [21,22] for evaluating integrals with oscillating kernels was extended to weakly singular integrals with logarithmic singularities.
Approximate methods for evaluating singular integrals with oscillating kernels of the form
K ( f , w ) ( t ) = 1 1 e i w x f ( x ) x t d x
have been studied in [9].
In [16], quadrature formulas are constructed for evaluating singular integrals of the form
I ( w , a ) = 1 1 e i w x 1 x 2 f ( x ) x a d x ,
where w is a large positive number, 1 < a < 1 .
To construct a quadrature formula, the function f (x) is approximated by interpolation polynomial P n ( f , x ) on nodes x 1 , x 2 , , x n , a , and a does not match x 1 , x 2 , , x n . As a result, the integral (5) is approximated by the quadrature formula
I ( w , a ) = 1 1 e i w x P n ( f , x ) 1 x 2 ( x a ) d x + R n ( f , a ) .
The estimate for R n ( f , a ) is given in [16].

1.2. Definitions of Singular and Hypersingular Integrals

Recall the definitions of function classes.
Let γ be the unit circle centered at the origin in the plane of the complex variable. Let A = [ 1 , 1 ] or A = γ .
Definition 1.
Hölder function class H α ( M ; A ) ( 0 < α 1 ) consists of functions f ( x ) defined on A satisfying the inequality | f ( x ) f ( x ) | M | x x | α .
Definition 2.
The class W r ( M ; A ) ,   r = 1 , 2 , , consists of functions f C [ a , b ] which have absolutely continuous derivatives of orders j = 0 , 1 , , r 1 and a piecewise continuous derivative f ( r ) satisfying f ( r ) ( x ) M .
Definition 3.
The class W r H α ( M ; A ) consists of functions f ( x ) belonging to the class W r ( M ; A ) and satisfying the additional condition f ( r ) ( x ) H α ( M , A ) .
Consider the integral
a b f ( τ ) τ t d τ , a < t < b .
Definition 4.
The Cauchy principal value of the singular integral (6) is called the limit
lim η 0 a c η f ( τ ) τ c d τ + c + η b f ( τ ) τ c d τ .
Recall the definitions of hypersingular integrals. Hadamard [23] introduced a new type of integral, hypersingular integrals:
Definition 5.
The integral of the type
a b A ( x ) d x ( b x ) p + α
for an integer p and 0 < α < 1 defines a value of the above integral (“finite part”) as the limit of the sum
a x A ( t ) d t ( b t ) p + α + B ( x ) ( b x ) p + α 1 ,
as x b if one assumes that A ( x ) has p derivatives in the neighborhood of point b. Here, B ( x ) is any function that satisfies the following two conditions:
(i)
The above limit exists;
(ii)
B ( x ) has at least p derivatives in the neighborhood of a point x = b .
An arbitrary choice of B ( x ) does not depend on the value of the limit in (i). Condition (ii) defines the values of first ( p 1 ) derivatives of B ( x ) at point b. An arbitrary additional term in the numerator is infinitesimal, of order ( b x ) p .
Notation 1.
Hadamard [24] gave a fascinating report of various aspects of the creative process in solving mathematical problems and, in particular, on his discovery of hypersingular integrals.
Chikin [25] introduced the definition of the Cauchy–Hadamard type integral that generalized a singular integral in the Cauchy principal and Hadamard sense.
Definition 6.
The Cauchy–Hadamard principal sense of the following integral
a b φ ( τ ) d τ ( τ c ) p , a < c < b
is defined as the limit of the expression
a b φ ( τ ) d τ ( τ c ) p = lim v 0 a c v φ ( τ ) d τ ( τ c ) p + c + v b φ ( τ ) d τ ( τ c ) p + ξ ( v ) v p 1 ,
where ξ ( v ) is a function chosen so as to provide the existence of the limit above.
In some cases, it is more convenient to use the following definition of hypersingular integrals, which is equivalent to Definition 6.
Definition 7.
Let f ( t ) W r H α ( M ; [ 1 , 1 ] ) , r = 1 , 2 , , 0 < α 1 . A hypersingular integral with order of p + 1 singularity, p r , is defined by
1 1 f ( τ ) ( τ t ) p + 1 d τ = 1 p ! p t p 1 1 f ( τ ) τ t d τ .

1.3. Optimal Quadrature Formulas for Calculating Singular and Hypersingular Integrals

Formulation of the problem of constructing the best quadrature formula belongs to Kolmogorov. Bakhvalov introduced [26] the concepts of asymptotically optimal and optimal with respect to order passive algorithms for solving problems in numerical analysis. Other approaches to determine optimal passive algorithms are given in [27,28,29].
We give now the definition of optimal quadrature formulas for singular integrals.
Consider the quadrature rule
C φ 1 1 φ ( τ ) τ t d τ = k = 1 N l = 0 ρ p k l ( s ) φ ( l ) ( t k ) + R N ( φ , t , p k l , t k ) .
The error (8) is
R N ( φ , p k l , t k ) = sup 1 < t < 1 | R N ( φ , t , p k l , t k ) | .
The error of (8) on Ψ class is
R N ( Ψ , p k l , t k ) = sup φ Ψ R N ( φ , p k l , t k ) | .
We introduce the functional
ζ N [ Ψ ] = inf p k l , t k R N ( Ψ , p k l , t k ) ,
where the lower bound takes over all the nodes t k ,   1 < t k < 1 , and the coefficients p k l ,   k = 1 , 2 , , N ,   l = 0 , 1 , , ρ .
The quadrature Formula (8) is defined by a set of nodes t k * ,   k = 1 , 2 , , N , and coefficients p k l * ,   k = 1 , 2 , , N ,   l = 0 , 1 , , ρ , called optimal, asymptotically optimal and optimal with respect to order if ζ N [ Ψ ] / R N ( p k l * , t k * , Ψ ) = 1 , 1 , 1 , respectively.
In a similar way, the concept of optimal, asymptotically optimal and optimal with respect to order quadrature formulas for evaluating hypersingular integrals is introduced.
Remark 1.
Let Ψ = W r ( [ 1 , 1 ] , 1 ) . Let the integral C φ be evaluated with quadrature formula
C φ = k = 1 N l = 0 ρ p k l ( t ) φ ( l ) ( t k ) + R N ( φ , t , p k l , t k )
with fixed nodes t k ,   k = 1 , 2 , , N , and fixed coefficients p k l ( t ) ,   k = 1 , 2 , , N ,   l = 0 , 1 , , ρ . In this case, the functional R N ( Ψ , p k l , t k ) is equivalent to the Peano constant. Theory of the Peano constants is a very important part of classical numerical theory (see [30]). Comparing the definitions of the Peano constant and optimal quadrature formulas, one can observe that the Peano constant theory is a special case of optimal algorithms theory.

2. Levin’s Method for Evaluating Singular and Hypersingular Integrals with Rapidly Oscillating Kernels

We present an application of Levin’s method for evaluating hypersingular integrals with rapidly oscillating kernels.
Consider the integral
1 1 f ( τ ) e i ω g ( τ ) ( τ t ) p d τ , 1 < t < 1 , p = 1 , 2 ,
The integral (9) is associated with the differential equation
d d τ [ x ( τ ) e i ω g ( τ ) ] = f ( τ ) e i ω g ( τ ) ( τ t ) p ,
where t is a parameter.
Differentiating the left-hand side, we have
x ( τ ) e i ω g ( τ ) + i ω g ( τ ) x ( τ ) e i ω g ( τ ) = f ( τ ) e i ω g ( τ ) ( τ t ) p ,
moreover, it is enough to consider the equation
x ( τ ) + i ω g ( τ ) x ( τ ) = f ( τ ) ( τ t ) p .
If it is possible to find an analytical solution of Equation (10), then
1 1 f ( τ ) e i ω g ( τ ) ( τ t ) p d τ = 1 1 d [ x ( τ ) e i ω g ( τ ) ] = = x ( 1 ) e i ω g ( 1 ) x ( 1 ) e i ω g ( 1 ) .
Note that when solving the differential Equation (11), the singularity can be avoided for τ = t .
Indeed, by the definition of the hypersingular integral, we have
1 1 f ( τ ) e i ω g ( τ ) d τ ( τ t ) p = lim η 0 1 t η f ( τ ) e i ω g ( τ ) ( τ t ) p d τ + + t + η 1 f ( τ ) e i ω g ( τ ) ( τ t ) p d τ + φ ( τ ) η p 1 .
The function φ ( τ ) has continuous derivatives up to p 1 order in a neighborhood of zero and is chosen such that the limit exists.
Taking the integrals separately on the right-hand side of (12) and applying the formula (11) to each of them, we have
1 1 f ( τ ) e i ω g ( τ ) d τ ( τ t ) p = lim η 0 x ( 1 ) e i ω g ( 1 ) + x ( t η ) e i ω g ( t η ) x ( t + η ) e i ω g ( t + η ) x ( 1 ) e i ω g ( 1 ) + φ ( η ) η p 1 .
The functions x ( t ± η ) e i ω g ( t ± η ) can be represented as a sum:
x ( t ± η ) e i ω g ( t ± η ) = x 1 ( t ± η ) e i ω g ( t ± η ) + x 2 ( t ± η ) e i ω g ( t ± η ) ,
where the first term tends to infinity as η 0 , and the second term tends to the finite limit.
Obviously,
1 1 f ( τ ) e i ω g ( τ ) d τ ( τ t ) p = lim η 0 x ( 1 ) e i ω g ( 1 ) + x 2 ( t η ) e i ω g ( t η ) x 2 ( t + η ) e i ω g ( t + η ) x ( 1 ) e i ω g ( 1 ) .
Assuming that f ( t ) , g ( t ) W p ( [ 1 , 1 ] , M ) , the function ψ ( t ) = 1 1 f ( τ ) e i ω g ( τ ) d τ ( τ t ) p is continuous at t ( 1 , 1 ) [4]. Then, lim η 0 ( x 2 ( t η ) e i ω g ( t η ) x 2 ( t + η ) e i ω g ( t + η ) ) = 0 and from (13) the final formula follows:
1 1 f ( τ ) e i ω g ( τ ) d τ ( τ t ) p = x ( 1 ) e i ω g ( 1 ) x ( 1 ) e i ω g ( 1 ) .
Thus, the analogue of the Newton–Leibniz formula for hypersingular integrals has been obtained. The application of the Newton–Leibniz formula for hypersingular integrals for certain function classes has been shown in [31].
It follows from the above that for evaluating hypersingular integrals with rapidly oscillating kernels, one can use numerical methods for solving ordinary differential equations.

3. Quadrature Formulas for Evaluating Singular Integrals with Rapidly Oscillating Functions

In this section, we study methods for evaluating the following types of singular integrals with rapidly oscillating functions
( C φ ) ( t ) 1 π i γ φ ( τ ) τ m d τ τ t , t γ ,
( H φ ) ( s ) 0 2 π φ ( σ ) { sin m σ cos m σ } c t g σ s 2 d σ , s [ 0 , 2 π ) ,
where m is a natural number.
Note that integral (14) is reduced by the Hilbert transformation to integral (15). Therefore, in this section we can restrict ourselves to considering the integral (15).

3.1. Lower Bound Estimates for Quadrature Formula Errors

First, we find a lower bound estimate for the quadrature formula errors using N values of integrands.
The integral (15) will be evaluated using the quadrature formula
( H φ ) ( s ) = j = 1 N p j ( s ) φ ( w j ) + R N ( s , p j ( s ) , w j , φ ) .
We find a lower bound estimate of the error for (16) provided that φ H 1 ( [ 0 , 2 π ] , 1 ) .
In doing so, we generalize the method for constructing optimal quadrature formulas for evaluating singular integrals proposed in [1,3].
There are two cases: (1) N 2 m ; (2) N 2 m .
First, assume N 2 m .
Introduce the nodes t k = 2 k π / N ,   k = 0 , 1 , , N ,   s l = l π / m ,   l = 0 , 1 , , 2 m . Note t 0 = t N , s 0 = s 2 m .
Let v j form a union of nodes t k ,   k = 0 , 1 , , N ,   w i ,   i = 1 , 2 , , N , and s l ,   l = 0 , 1 , , 2 m . Without loss of generality, assume the sets of nodes t k ,   k = 0 , 1 , , N ,   w i ,   i = 1 , 2 , , N , and s l = 0 , 1 , , 2 m , do not intersect.
It follows from below that when some nodes coincide, the lower bound error of the quadrature formula does not decrease. Thus, we assume that the number of nodes v j , j = 0 , 1 , , n 1 , is equal to n = 2 m + 2 N .
Let Δ k be segments Δ k = [ v k , v k + 1 ] ,   k = 0 , 1 , , n 1 ,   v n = v 0 .
To each node t k ,   k = 0 , 1 , , N 1 , we assign the following function
φ k ( s ) = ( sgn ( sin m s ) ) min ( | s v i | , | v i + 1 s | ) , s [ v i , v i + 1 ] , [ v i , v i + 1 ] [ t j , t j + 1 ] , j = k , k + 1 , , k + [ N / 2 ] 1 ; ( sgn ( sin m s ) ) min ( | s v i | , | v i + 1 s | ) , s [ v i , v i + 1 ] , [ v i , v i + 1 ] [ t j , t j + 1 ] , j = k + [ N / 2 ] + 1 , k + [ N / 2 ] + 2 , , k + N 1 ; 0 , s [ k + t [ N / 2 ] 1 , k + t [ N / 2 ] + 1 ] .
Then,
( H φ k ) ( t k ) = 0 2 π φ k ( σ ) sin m σ c t g σ t k 2 d σ l = 1 [ N / 2 ] 2 t l + k t l + k + 1 φ k ( σ ) sin m σ c t g σ t k 2 d σ + + l = [ N / 2 ] + 1 N 1 t l + k t l + k + 1 φ k ( σ ) sin m σ c t g σ t k 2 l = 1 [ N / 2 ] 2 c t g π ( l + 1 ) N t l + k t l + k + 1 | φ k ( σ ) | | sin m σ | d σ + + l = [ N / 2 ] + 1 N 1 c t g π l N t l + k t l + k + 1 | φ k ( σ ) | | sin m σ | d σ .
Let us average the previous inequality over k ,   k = 0 , 1 , , N 1 . We have
max 0 l N 1 ( H φ l ) ( t l ) 1 N k = 0 N 1 ( H φ k ) ( t k ) 1 N k = 0 N 1 l = 1 [ N / 2 ] 2 c t g π ( l + 1 ) N t l + k t l + k + 1 | φ k ( σ ) | | sin m σ | d σ + + 1 N k = 0 N 1 [ N / 2 ] + 1 N 1 c t g π l N t l + k t l + k + 1 | φ k ( σ ) | | sin m σ | d σ = = 1 N l = 1 [ N / 2 ] 2 c t g π ( l + 1 ) N 0 2 π φ * ( σ ) | sin m σ | d σ + + 1 N l = [ N / 2 ] + 1 N 1 c t g π l N 0 2 π φ * ( σ ) | sin m σ | d σ = = ( 1 + o ( 1 ) ) 2 N 0 2 π φ * ( σ ) | sin m σ | d σ l = 1 [ N / 2 ] 2 c t g π ( l + 1 ) N .
Here, φ * ( σ ) = min ( | σ v i | , | v i + 1 σ | ) ,   σ [ v i , v i + 1 ] ,   i = 0 , 1 , , n .
Estimate from below the integral 0 2 π φ * ( σ ) | sin m σ | d σ under the following conditions: (1) the function φ * ( σ ) H 1 ( [ 0 , 2 π ] , 1 ) ; (2) the function φ * ( σ ) is non-negative; (3) the function φ * ( σ ) vanishes at points v k ,   k = 0 , 1 , , n 1 .
The set of nodes { v k } is a union of three sets { t i } ,   i = 0 , 1 , , N 1 ,   { w i } ,   i = 1 , 2 , , N ,   { s j } ,   j = 0 , 1 , , 2 m 1 .
Let N m / 2 .
Then, there are at least 2 m 2 N segments [ s l , s l + 1 ] , in which there are no nodes from the sets { w j } , j = 1 , 2 , , N , and { t i } ,   i = 0 , 1 , , N 1 . We call such segments marked.
It is easy to see that in the marked segments
s l s l + 1 φ * ( σ ) | sin m σ | d σ = 0 π / 2 m σ sin m σ d σ + π / 2 m π / m ( π m σ ) sin m σ d σ = 2 m 2 .
Thus, for N m / 2 , we have the estimate
0 2 π φ * ( σ ) | sin m σ | d σ ( 2 m 2 N ) C m 2 = C m .
From this estimate and the inequality (18), we have
R N ( H 1 ( [ 1 , 1 ] , 1 ) ) C m N k = 1 N 2 c t g k + 1 N = C m ln N .
Here, and below C, are the constants independent of N and m .
For m 2 N 2 m , we must change the proof. Let { v i } , i = 0 , 1 , , n be a union of node sets { w k } , k = 1 , 2 , , N and { s l } , i = 0 , 1 , , 2 m .
Let
φ * ( s ) = min ( | s v i | , | v i + 1 s | ) , s [ v i , v i + 1 ] , i = 0 , 1 , , n 1 .
Each node s l is associated with the function
φ l ( s ) = ( s g n ( sin m s ) ) min ( | s v i | , | v i + 1 s | ) , s [ v i , v i + 1 ] , [ v i , v i + 1 ] [ s j , s j + 1 ] , j = l , l + 1 , , m 1 ; ( s g n ( sin m s ) ) min ( | s v i | , | v i + 1 s | ) , s [ v i , v i + 1 ] , [ v i , v i + 1 ] [ s j , s j + 1 ] , j = m + 1 , , 2 m 1 .
Then,
( H φ l ) ( s l ) k = 1 m 2 c t g π ( k + 1 ) m s l + k s l + k + 1 | φ l ( σ ) | | sin m σ | d σ + + k = m + 1 2 m 1 c t g k π l m t l + k t l + k + 1 | φ l ( σ ) | | sin m σ | d σ .
Averaging the previous inequality over l , l = 0 , 1 , , 2 m 1 , we have
max 0 l 2 m 1 ( H φ l ) ( s l ) 1 2 m k = 1 m 2 c t g π ( k + 1 ) m 0 2 π φ * ( σ ) | sin m σ | d σ + 1 2 m k = m + 1 2 m 1 c t g k π m 0 2 π φ * ( σ ) | sin m σ | d σ = = ( 1 + o ( 1 ) ) 1 m 0 2 π φ * ( σ ) | sin m σ | d σ k = 1 m 2 c t g π ( k + 1 ) m .
Estimate from below the integral
0 2 π φ * ( σ ) | sin m σ | d σ .
The integral takes the smallest value if in each interval ( s i , s i + 1 ) there is at most one node w j , j = 1 , 2 , , N .
It was shown above (19) that if there are no nodes w j , j = 1 , 2 , , N in ( s i , s i + 1 ) , then
s i s i + 1 φ * ( σ ) | sin m σ | d σ = 2 m 2 .
Next, consider the case when the interval ( s i , s i + 1 ) contains more than one node from { w j } , j = 1 , 2 , , N . Without loss of generality, we assume N = 2 m .
First, let the interval ( s i , s i + 1 ) contain one node from the set { w j } , j = 1 , 2 , , N . Obviously,
s i s i + 1 φ * ( σ ) | sin m σ | d σ 2 s i s i + π 4 m ( σ s i ) | sin m σ | d σ + + 2 s i + π 4 m s i + π 2 m ( s i + π 2 m σ ) | sin m σ | d σ = = 2 m 2 2 sin m 2 s i + π 2 m 2 sin m s i sin m s i + π 2 m = = 2 m 2 2 sin m s i + π 4 m sin m s i sin m s i + π 2 m = = 2 m 2 2 cos i π sin π 4 ( cos i π ) sin π 2 = = 2 m 2 | ( 1 ) i 2 ( 1 ) i | = 2 m 2 | 2 1 | 2 m 2 0.41 .
Assume now that the interval ( s i , s i + 1 ) contains two nodes from { w j } , j = 1 , 2 , , N . This means there is an interval ( s l , s l + 1 ) , which does not contain nodes from { w j } , j = 1 , 2 , , N .
From (22) and (23), we have
s i s i + 1 φ * ( σ ) | sin m σ | d σ + s l s l + 1 φ * ( σ ) | sin m σ | d σ = = 2 m 2 + s l s l + 1 φ * ( σ ) | sin m σ | d σ 2 m 2 0.82 .
Thus, the minimum of the integral 0 2 π φ * ( σ ) | sin m σ | d σ is achieved under the assumption that each interval has at most one node { w j } , j = 1 , 2 , , N .
Other cases for the distribution of nodes w j , j = 1 , 2 , , N over intervals ( s i , s i + 1 ) , i = 0 , 1 , , 2 m 1 are studied similarly.
Thus, from (21)–(23), we have for m 2 N 2 m
R N ( H 1 ( [ 1 , 1 ] , 1 ) ) C m 2 k = 1 m 2 c t g ( k + 1 ) π m = C m ln m = C N ln N .
From this estimate and the inequality (20), we have for N 2 m
R N ( H 1 ( [ 1 , 1 ] , 1 ) ) C m 2 k = 1 m 2 c t g ( k + 1 ) π m = C m ln N .
For the second case, consider N 2 m . It is enough to introduce additional nodes t k = 2 k π / q , k = 0 , 1 , , q , q = 3 N . Let the set v k ,   v k = 0 , 1 , , n ,   n = N + 2 m + q = 4 N + 2 m , be the union of nodes w i , i = 1 , 2 , , N ,   s j , j = 0 , 1 , , 2 m ,   t l , l = 0 , 1 , , q . Each node t k , k = 0 , 1 , , q 1 , is associated with the function ψ k ( s ) , constructed by analogy with the function φ k ( s ) (see (17)). It is easy to see that there are at least N intervals ( t l , t l + 1 ) , in which there are no nodes from the sets s l ,   i = 0 , 1 , , 2 m , and w j , j = 1 , 2 , , N .
Repeating the above arguments yields
R N [ H 1 ( 1 ) ] C N 2 k = 1 q 1 c t g π ( k + 1 ) 2 q = C N ln N .
From the inequalities (24) and (25), the next statement follows.
Theorem 1.
Let φ ( t ) H 1 ( 1 ) . For all possible quadrature formulas of the form (16) using N nodes, the following estimate holds
R N ( H 1 ( [ 0 , 2 π ] , 1 ) ) { C 1 ln N m , N 2 m ; C 2 ln N N , N 2 m ,
where C 1 , C 2 are constants independent of N.
Making the proof more difficult yields the following statement.
Theorem 2.
Let be φ ( t ) H α ( 1 ) , 0 < α 1 . For all possible quadrature formulas of the form (16) using N nodes, the following estimate holds:
R N ( H α ( [ 0 , 2 π ] , 1 ) ) { C 1 ln N m α , N 2 m ; C 2 ln N N α , N 2 m .
Estimate from below the error of quadrature formulas of (16) on the W r ( [ 0 , 2 π ] , 1 ) ,   r = 1 , 2 , class.
In order to simplify the presentation, we give the proof for N / 2 m only.
We introduce the nodes { t k } ,   t k = 2 k π / N , k = 0 , 1 , , N 1 ;   w j , j = 1 , 2 , , N ;   { s k } ,   k = 0 , 1 , , 2 m 1 ;   { v j } ,   j = 1 , 2 , , n ,   n = 2 N + 2 m .
To each node t k ,   k = 0 , 1 , , N 1 we assign the function
φ k ( s ) = A ( s v l ) r ( v l + 1 s ) r h l r sgn sin ms , s [ v l , v l + 1 ] , [ v l , v l + 1 ] [ t j , t j + 1 ] j = k , k + 1 , , k + [ N / 2 ] 1 ; A ( s v l ) r ( v l + 1 s ) r h l r sgn sin ms , s [ v l , v l + 1 ] , [ v l , v l + 1 ] [ t j , t j + 1 ] j = k + [ n / 2 ] + 1 , k + [ N / 2 ] + 2 , , k + N 1 ; 0 , s [ v k + [ N / 2 ] 1 ] , [ v k + [ N / 2 ] + 1 ] ,
h l = | v l + 1 v l | , l = 0 , 1 , , n 1 .
Estimating the integral
0 2 π φ k ( σ ) sin m σ c t g σ t k 2 d σ
and then averaging the result over k ,   k = 0 , 1 , , N 1 , we obtain the required estimate for N / 2 m .
Theorem 3.
Let φ W r ( [ 0 , 2 π ] , 1 ) . For all possible quadrature formulas of the form (16) using N nodes, the following estimate holds:
R N ( W r ( [ 0 , 2 π ] , 1 ) ) { C 1 ln N m r , N 2 m ; C 2 ln N N r , N 2 m .

3.2. Quadrature Formulas

Let us construct quadrature formulas for evaluating integrals of (15).
We start by considering singular integrals with the Hilbert kernel:
1 2 π 0 2 π f ( σ ) cos m σ sin m σ c t g σ s 2 d σ ,
where f W r ( [ 0 , 2 π ] , M ) , M is an integer.
First, we consider the integral
H f = 1 2 π 0 2 π f ( σ ) sin m σ c t g σ s 2 d σ .
The function f ( s ) is approximated by the interpolation polynomial
f n ( s ) = k = 0 2 n f ( s k ) ψ k ( s ) ,
where ψ k ( s ) = 1 2 n + 1 sin 2 n + 1 2 ( s s k ) sin s s k 2 = 1 2 n + 1 ( 1 + 2 cos ( s s k ) + + 2 cos n ( s s k ) ) .
We have for m > n
H f = 1 2 n + 1 k = 0 2 n f ( s k ) 1 2 π 0 2 π sin m σ c t g σ s 2 d σ + + 1 2 n + 1 k = 0 2 n f ( s k ) 1 π 0 2 π l = 1 n cos l ( σ s k ) sin m σ c t g σ s 2 d σ + R n ( f ) = = 1 2 n + 1 k = 0 2 n f ( s k ) cos m s + 1 2 n + 1 k = 0 2 n f ( s k ) 1 2 π 0 2 π l = 1 n [ ( sin ( m + l ) σ ) ( cos l s k ) ( cos ( m + l ) σ ) ( sin l s k ) + + ( sin ( m l ) σ ) ( cos l s k ) ( cos ( m l ) σ ) ( sin l s k ) ] c t g σ s 2 d σ + R n ( f ) = = 1 2 n + 1 k = 0 2 n f ( s k ) cos m s + 1 2 n + 1 k = 0 2 n f ( s k ) l = 1 n [ ( cos ( m + l ) s ) ( cos l s k ) + + ( sin ( m + l ) s ) ( sin l s k ) + ( cos ( m l ) s ) ( cos l s k ) + sin ( m l ) s ) ( sin l s k ) ] + R n ( f ) = = 1 2 n + 1 k = 0 2 n f ( s k ) cos m s + + 2 l = 1 n [ ( cos m s ) ( cos l s ) ( cos l s k ) + ( cos m s ) ( sin l s ) ( sin l s k ) ] = = 1 2 n + 1 cos m s k = 0 2 n f ( s k ) 1 + 2 l = 1 n cos l ( s s k ) = = cos m s k = 0 2 n f ( s k ) ψ k ( s ) = ( cos m s ) f n ( s ) .
Above, we used ([32], p.36)
1 2 π π π e i m σ c t g σ s 2 d σ = 0 , m = 0 i e i m s , m > 0 i e i m s , m < 0
for m integer.
Next, consider the integral
K f = 1 2 π 0 2 π f ( σ ) cos m σ c t g σ s 2 d σ .
As above, the function f ( s ) is approximated by the interpolation polynomial f n ( s ) .
Obviously, for m > n ,
K f = 1 2 n + 1 k = 0 2 n f ( s k ) 1 2 π 0 2 π cos m σ c t g σ s 2 d σ + + 1 2 n + 1 k = 0 2 n f ( s k ) l = 1 n 1 π 0 2 π cos l ( σ s k ) cos m σ c t g σ s 2 d σ + R n ( f ) = = 1 2 n + 1 k = 0 2 n f ( s k ) sin m s cos m s c o s e c s s k 2 cos s s k 4 + cos 2 n + 1 2 ( s s k ) 2 + R n ( f ) .
Now, we study error estimates for constructed quadrature formulas. It is enough to consider the quadrature Formula (29). It is easy to see that the error of (29) is estimated by the inequality
| R n ( f ) | = 1 2 π 0 2 π ψ n ( σ ) sin m σ c t g σ s 2 d σ 1 2 π 0 2 π ( ψ n ( σ ) ψ n ( s ) ) sin m σ c t g σ s 2 d σ + 1 2 π ψ n ( s ) 0 2 π sin m σ c t g σ s 2 d σ = = I 1 ( s ) + I 2 ( s ) .
Here, ψ n ( s ) = f ( s ) f n ( s ) .
Evaluate each term separately,
I 2 ( s ) = | ψ n ( s ) cos m s | C E n ( f ) ln n ,
where E n ( f ) is the best approximation in the uniform metric for the function f by nth-order trigonometric polynomials
I 1 ( s ) 1 2 π 0 2 π | ψ n ( σ ) ψ n ( s ) | 1 β | ψ n ( σ ) ψ n ( s ) | β | sin m σ | | c t g σ s 2 | d σ 1 2 π 0 2 π | ψ n ( σ ) ψ n ( s ) | 1 β | ψ n ( s + θ ( σ s ) ( σ s ) ) | β | c t g σ s 2 | d σ C ( max 0 s 2 π | ψ n ( s ) | ) 1 β ( max 0 s 2 π | ψ n ( s ) | ) β 0 2 π d σ ( σ s ) 1 β d σ C β ( max 0 s 2 π | ψ n ( s ) | ) 1 β ( max 0 s 2 π | ψ n ( s ) | ) β .
The following statement is well known.
Theorem 4
([33]). Let f ( t ) C ˜ [ 0 , 2 π ] , P n ( t ) be an n-order trigonometric polynomial satisfying the inequality | f ( t ) P n ( t ) | η n .
Then, | f ( t ) P n ( t ) | C n η n , where C is a constant independent of n.
Thus
( max 0 s 2 π | ψ n ( s ) | ) β C ( E n ( f ) ln n ) β .
Setting β = 1 ln n , we have
I 2 C E n ( f ) ln 2 n .
From the estimates I 1 and I 2 , we have
| R n ( f ) C E n ( f ) ln 2 n ,
and, therefore, on the function class W r H α ( [ 0 , 2 π ] , M ) .
R n ( W r H α ) C ln 2 n n r + α .
The final estimate is valid for any m 1 .
Now, we consider the following quadrature formula for H f evaluation.
We approximate the function f ( s ) by the polygon f N ( s ) , constructed on the nodes t k = 2 k π / N , k = 0 , 1 , , N .
The integral H f will be evaluated using the quadrature formula
( H f ) ( s ) 0 2 π f N ( σ ) sin m σ c t g σ s 2 d σ + R N ( s , p k ( s ) , t k , f ) .
The error of (32) is estimated by
| R N ( s , s k , f ) | | l = 0 N 1 t l t l + 1 ( f ( σ ) f N ( σ ) ) sin m σ c t g σ s 2 d σ | .
Consider two cases: (1) 2 m N , (2) N < 2 m .
Start with the first one.
Let be s Δ j , Δ k = [ t k , t k + 1 ] , k = 0 , 1 , , N 1 .
Estimate the integral
J 1 = | t l t l + 1 ψ N ( σ ) sin m σ c t g σ s 2 d σ |
for j l 1 , l , l + 1 .
Set ψ N ( s ) = f ( s ) f N ( s ) . Obviously, ψ N ( t l ) = ψ N ( t l + 1 ) = 0 .
Obviously,
| t l t l + 1 ψ N ( σ ) sin m σ c t g σ s 2 d σ | t l t l + 1 | ψ N ( σ ) | | c t g σ s 2 | d σ 2 π N 2 max ( | c t g ( l j 1 ) π N | , | c t g ( j l 1 ) π N | ) .
Let s Δ j , j = l . Estimate the integral
J 2 = t l 1 t l + 2 ψ N ( σ ) ) sin m σ c t g σ s 2 d σ .
Let d = min ( | t l + 2 s | , | s t l 1 | ) . Set d = | s t l 1 | .
Represent the previous integral as
J 2 = t l 1 t l + 2 ψ N ( σ ) ) sin m σ c t g σ s 2 d σ = = t l 1 s + d ψ N ( σ ) ) sin m σ c t g σ s 2 d σ + + s + d t l + 2 ψ N ( σ ) ) sin m σ c t g σ s 2 d σ = J 21 + J 22 .
Estimate each of the integrals J 21 , J 22 separately.
Obviously,
| J 22 | s + d t l + 2 | ψ N ( σ ) | | c t g σ s 2 | d σ 4 π 2 N 2 c t g π N .
Estimate the integral J 21 . We have
| J 21 | = | s l 1 s + d ψ N ( σ ) sin m σ c t g σ s 2 d σ | | t l 1 s + d ( ψ N ( σ ) ψ N ( s ) ) sin m σ c t g σ s 2 d σ | + + | s l 1 s + d ψ N ( s ) ( sin m σ sin m s ) c t g σ s 2 d σ | = J 211 + J 212 .
Estimate J 211 , j 212 .
Obviously,
J 211 t l 1 s + d | ψ N ( σ ) ψ N ( s ) | | c t g σ s 2 | d σ 6 π 2 N ,
J 212 s l 1 s + d | ψ N ( s ) | | sin m σ sin m s | | c t g σ s 2 | d σ 6 π 2 m N 2 .
From inequalities (36)–(40), it follows
| J 2 | C N .
From inequalities (33)–(41), it follows that for 2 m < N the inequality holds
R N ( H 1 ( 1 ) ) C ln N N .
The inequality
R N ( H α ( 1 ) ) C ln N N α
is proved in a similar way.
Consider the second case, N < 2 m . Let s Δ j , Δ k = [ t k , t k + 1 ] , k = 0 , 1 , , N 1 .
Estimating (34), we again consider two cases, j l 1 , l , l + 1 and j = l .
For the first one, after making some calculations, it can be shown that the largest error is yielded for functions of the form ψ N ( s ) = ( min ( s s i , s i + 1 s ) sgn sin ms ,   s [ s i , s i + 1 ] , s i = π i / m , i = 0 , 1 , , 2 m .
Then,
J 1 = | t l t l + 1 ψ N ( σ ) sin m σ c t g σ s 2 d σ | 2 π 2 m 2 max ( | c t g ( l j + 1 ) π N | , | c t g ( j l 1 ) π N | ) .
Now, let j = l , s [ t j , t j + 1 ] and s [ s v , s v + 1 ] .
We represent the integral J 1 as
J 1 | t j 1 s v 1 ψ N ( σ ) sin m σ c t g σ s 2 d σ | + + | s v 1 s v + 2 ψ N ( σ ) sin m σ c t g σ s 2 d σ | + | s v + 2 t j + 1 ψ N ( σ ) sin m σ c t g σ s 2 d σ | = = J 11 + J 12 + J 13 .
Obviously,
J 11 C 1 m 2 l = 1 L 1 c t g l π 2 m C m ln m ; J 13 C 1 m 2 l = 1 L 1 c t g l π 2 m C m ln m ; J 12 = | s v 1 s v + 2 ( ψ N ( σ ) sin m σ ψ N ( s ) sin m s ) c t g σ s 2 d σ | C m ,
where L = 2 m / N .
Making calculations similar to those above, we obtain the estimate
J 1 { C 1 1 N ln N , 2 m N ; C 2 1 m ln m , N 2 m ,
where constants C 1 , C 2 are independent of N.
Theorem 1 and inequalities (33)–(45) yield the following statement.
Theorem 5.
Among all quadrature formulas of the form (16) using N nodes, the optimal with respect to order on the function class H 1 ( 1 ) turns out to be Formula (32). The estimate is valid:
R N ( H 1 ( 1 ) ) { C 1 1 N ln N , 2 m N ; C 2 1 m ln m , N 2 m ,
where constants C 1 , C 2 are independent of N.
Similarly, we can prove the following.
Theorem 6.
Among all quadrature formulas of the form (16) using N nodes, the optimal with respect to order on the function class H α ( 1 ) turns out to be Formula (32). The estimate holds:
R N ( H α ( 1 ) ) { C 1 1 N α ln N , 2 m N ; C 2 1 m α ln m , N 2 m ,
where constants C 1 , C 2 are independent of N.
Using the Hilbert transformation, we obtain
C f = 1 π i γ f ( τ ) τ m τ t d τ = = 1 2 π i 0 2 π f ( e i σ ) cos m σ c t g σ s 2 d σ + 1 2 π 0 2 π f ( e i σ ) sin m σ c t g σ s 2 d σ + + 1 2 π 0 2 π f ( e i σ ) cos m σ d σ + i 2 π 0 2 π f ( e i σ ) sin m σ d σ .
To evaluate the integral C f , we use the quadrature formula
C f = 1 2 π i 0 2 π P n [ f ( e i σ ) ] cos m σ c t g σ s 2 d σ + + 1 2 π 0 2 π P n [ f ( e i σ ) ] sin m σ c t g σ s 2 d σ + + 1 2 π 0 2 π P n [ f ( e i σ ) ] cos m σ d σ + i 2 π 0 2 π P n [ f ( e i σ ) ] sin m σ d σ + R n [ f ] ,
where P n is a projection operator onto a set of interpolating trigonometric polynomials on nodes s k = 2 π k / ( 2 n + 1 ) , k = 0 , 1 , , 2 n .
The error of (46) is estimated by
| R n ( f ) | 1 2 π 0 2 π ψ n ( σ ) cos m σ c t g σ s 2 d σ + + 1 2 π 0 2 π ψ n ( σ ) sin m σ c t g σ s 2 d σ + + 1 2 π 0 2 π ψ n ( σ ) cos m σ d σ + 1 2 π 0 2 π ψ n ( σ ) sin m σ d σ = = I 1 + I 2 + I 3 + I 4 ,
where ψ n ( t ) = f ( t ) P n [ f ] ( t ) .
The estimates I 1 and I 2 have been obtained above (see (31)).
The estimates hold:
I 3 C E n ( f ) ln n ,
I 4 C E n ( f ) ln n .
The inequalities (31), (48) and (49) yield the estimate
| R n ( f ) C E n ( f ) ln 2 n .
Transform the integrals
1 2 π i 0 2 π P n [ f ( e i σ ) ] cos m σ c t g σ s 2 d σ
and
1 2 π i 0 2 π P n [ f ( e i σ ) ] sin m σ c t g σ s 2 d σ .
It was shown above that
1 2 π i 0 2 π P n [ f ( e i σ ) ] sin m σ d σ = P n [ f ( e i s ) ] cos m s = = ( cos m s ) 1 2 n + 1 k = 0 2 n f ( e i s k ) 1 + 2 l = 1 m cos l ( s s k ) = 1 2 e i m s + e i m s 1 2 n + 1 k = 0 2 n f ( e i s k ) 1 + l = 1 m e i l ( s s k ) + e i l ( s s k ) = = 1 2 t m + t m 1 2 n + 1 k = 0 2 n f ( t k ) 1 + l = 1 m t l t k l + t k l t l ;
1 2 π i 0 2 π P n [ f ( e i σ ) ] cos m σ d σ = = 1 2 n + 1 k = 0 2 n f ( e i s k ) sin m s + 2 cos m s l = 1 m sin l ( s s k ) = = 1 2 k = 0 2 n f ( t k ) 1 2 i e i m s + e i m s + e i m s + e i m s × × l = 1 n 1 2 i e i l ( s s k ) + e i l ( s s k ) = = 1 2 n + 1 k = 0 2 n f ( t k ) 1 2 i t m t m + t m t m l = 1 n 1 2 i t l t k l t k l t l .
Thus, the following quadrature formula is valid:
C φ = 1 2 t m t m 1 2 n + 1 k = 0 2 n f ( t k ) 1 + l = 1 m t l t k l t k l t l 1 i 1 2 n + 1 l = 1 2 n f ( t k ) 1 2 i t m t m + t m t m l = 1 n 1 2 i t l t k l t k l t l + + 1 2 π 0 2 π P n f ( e i σ ) cos m σ d σ + 1 2 π 0 2 π P n f ( e i σ ) sin m σ d σ + R n ( f ) .
Estimates (48)–(50) hold:
| R n ( f ) | C E n ( f ) ln 2 n ,
where E n ( f ) is the best uniform approximation of the function f by trigonometric polynomials of order n.
Let us take a look at an illustration of the quadrature formulas we have discussed. Consider an integral
J m ( s ) = 1 2 π π π σ ( σ 2 π ) ( σ π ) 12 sin m σ c t g σ s 2 d σ = cos m s k = 1 sin k s k 3
where m is integer. Let us apply the quadrature Formula (29) for evaluation of the integral:
J m ( s ) = 1 2 n + 1 cos m s k = 0 2 n f ( s k ) 1 + 2 l = 0 n cos l ( s s k ) ,
where f ( s ) = s ( s 2 π ) ( s π ) 12 and s k = 2 k π 2 n + 1 for k = 0 , 1 , , 2 n . We present the results of evaluation of the integral by series summation and by quadrature formula in Figure 1. We observe rapid convergence of the quadrature formula to the exact value of the integral. We also show that the amplitude of the oscillations is determined by the function f ( s ) , as is suggested by Equation (29).

4. Approximate Evaluation of Hypersingular Integrals with Rapidly Oscillating Functions

In this section, we study approximate methods for evaluating hypersingular integrals of the form
( H f ) ( t ) = 1 π i γ f ( τ e i m τ ) d τ ( τ t ) p , t γ , p = 2 , 3 , .
Here, γ = { z : | z | = 1 } is a unit circle centered at the origin in the complex plane, m is a natural number. To obtain a lower bound estimate for the error of the quadrature formula, we use the Hilbert transformation from the integral (56) to the hypersingular integral with a Hilbert kernel. We change the variables in (56): τ = e i σ ,   t = e i s ,   σ [ 0 , 2 π ] .
Now, we have
1 π i γ f ( τ ) e i m τ d τ ( τ t ) p = 1 π 0 2 π f ( e i σ ) ( cos m σ + i sin m σ ) e i σ d σ ( e i σ e i s ) p .
Converting the fraction yields
1 ( e i σ e i s ) p ) = e i σ s 2 p ( e i σ e i s ) p e i σ s 2 p = = e σ s 2 ( e i σ e i s ) e i σ + s 2 p = = cos σ s 2 i sin σ s 2 e i σ s 2 e i σ s 2 p = cos σ s 2 + i sin σ s 2 2 i sin σ s 2 p = = ( 1 ) p 2 p i c t g σ s 2 + 1 p .
Thus, to estimate from below the error for evaluation integrals of the form (56) by quadrature formulas constructed on N nodes, it is enough to study the integrals of the form
1 2 π 0 2 π φ ( σ ) { sin m σ cos m σ } c t g p σ s 2 d σ .
Evaluate the hypersingular integral
F φ = 0 2 π φ ( σ ) sin m σ c t g p σ s 2 d σ
by the quadrature formula
F φ = k = 1 N p k ( s ) φ ( w k ) + R N ( s , p k ( s ) , w k , φ )
on the function class W r ( [ 0 , 2 π ] , 1 ) , r p .
When estimating the error of the quadrature formula from below, two cases should be considered:
(1)
p is an even natural number;
(2)
p is an odd natural number.
Let us first study an integral with a singularity of the even order.
Let t k = 2 k π / N ,   k = 0 , 1 , , N 1 .
The set of nodes { v j } ,   j = 0 , 1 , , n 1 ,   n = 2 N + 2 m is obtained by merging { t k } ,   k = 0 , 1 , , N 1 ,   { w k } ,   k = 1 , , N , and s j = π j / m ,   j = 0 , 1 , , 2 m 1 .
We introduce the function
φ * ( s ) = min 0 i 2 N + 2 m A ( s v i ) r ( v i + 1 s ) r h i r ,
s [ v i , v i + 1 ] , h i = | v i + 1 v i | , i = 0 , 1 , , n 1 .
The constant A is chosen such that φ * W r ( [ 0 , 2 π ] , 1 ) .
To each node t k ,   k = 0 , 1 , , N 1 , we assign the function
φ k ( s ) = φ * ( s ) sgn sin ms .
Then,
( F φ k ) ( s k ) = 0 2 π φ k ( σ ) sin m σ c t g p σ s k 2 d σ = = l = 0 N 1 t l t l + 1 φ k ( σ ) sin m σ c t g p σ s k 2 d σ l = 0 [ N / 2 ] 1 t k + l t k + l + 1 φ k ( σ ) sin m σ c t g p σ s k 2 d σ + + l = [ N / 2 ] + 1 N 2 t k + l t k + l + 1 φ k ( σ ) sin m σ c t g p σ s k 2 d σ l = 0 [ N / 2 ] 1 c t g p ( l + 1 ) π m t k + l t k + l + 1 φ * ( σ ) | sin m σ | d σ + + l = [ N / 2 ] + 1 N 2 c t g p l π m t k + l t k + l + 1 φ * ( σ ) | sin m σ | d σ .
Averaging the previous inequality over k ,   k = 0 , 1 , , N 1 , we have
max 0 j N 1 ( F φ j ) ( s j ) 1 N k = 0 N 1 ( F φ k ) ( s k ) 1 N k = 0 N 1 l = 0 [ N / 2 ] 1 c t g p ( l + 1 ) π N t k + l t k + l + 1 φ * ( σ ) | sin m σ | d σ + + 1 N k = 0 N 1 l = [ N / 2 ] + 1 N 2 c t g p l π N t k + l t k + l + 1 φ * ( σ ) | sin m σ | d σ = = 2 ( 1 + o ( 1 ) ) N l = 0 [ N / 2 ] 1 c t g p ( l + 1 ) π N 0 2 π φ * ( σ ) | sin m σ | d σ .
Now, we obtain the estimate of the integral 0 2 π φ * ( σ ) | sin m σ | d σ under the following assumptions:
(1)
φ * ( σ ) W r ( 1 ) ;
(2)
( φ * ) ( i ) ( v k ) = 0 , i = 0 , 1 , , r 1 , k = 0 , 1 , , n ;
(3)
φ * ( s ) 0 , s [ 0 , 2 π ] .
First, let N 2 m . Repeating the arguments above when studying singular integrals, we have
0 2 π φ * ( σ ) | sin m σ | d σ C m r .
For 2 m < N , we have
0 2 π φ * ( σ ) | sin m σ | d σ C N r .
Thus, for even p, we obtain the estimate
R N ( W r ( 1 ) ) { C 1 1 m r 1 N k = 0 [ N / 2 ] 1 c t g p k π N = 1 m r N 1 p , N 2 m ; C 2 1 N r 1 N k = 0 [ N / 2 ] 1 c t g p k π N = 1 N r + 1 p , 2 m N .
For odd p, the construction is more complicated. Nevertheless, we are sure that the following estimate is valid:
R N ( W r ( 1 ) ) { C 1 1 m r 1 N k = 0 [ N / 2 ] 1 c t g p k π N = 1 m r N 1 p , N 2 m ; C 2 1 N r 1 N k = 0 [ N / 2 ] 1 c t g p k π N = 1 N r + 1 p , 2 m N .
Finally, the theorem is proved.
Theorem 7.
Let φ W r ( 1 ) . Then, for all possible quadrature formulas of the form (57), the estimate
R N ( W r ( 1 ) ) { C 1 1 m r 1 N k = 0 [ N / 2 ] 1 c t g p k π N = 1 m r N 1 p , N 2 m ; C 2 1 N r 1 N k = 0 [ N / 2 ] 1 c t g p k π N = 1 N r + 1 p , 2 m N .
is valid.
When constructing quadrature formulas for hypersingular integral evaluation, we will use Definition 7, which allows us to construct methods for hypersingular integral evaluation based on well-known methods for evaluating singular integrals.
From Definition 7, it follows that
1 π i γ f ( τ ) τ m ( τ t ) p d τ = 1 ( p 1 ) ! d p 1 d t p 1 1 π i γ f ( τ ) τ m τ t d τ = = 1 ( p 1 ) ! d p 1 d t p 1 1 2 t m t m + 1 2 n + 1 k = 0 2 n f ( t k ) 1 + l = 1 m t l t k l t k l t l 1 i 1 2 n + 1 k = 0 2 n f ( t k ) 1 2 i t m t m + t m t m l = 1 n 1 2 i t l t k l t k l t l + R n ( f ) .
To estimate | R n ( f ) | , note that the integrals
1 2 π 0 2 π P n [ f ( e i σ ) ] sin m σ c t g σ s 2 d σ
and
1 2 π 0 2 π P n [ f ( e i σ ) ] cos m σ c t g σ s 2 d σ .
are trigonometric polynomials of the ( n + m ) th order.
As shown above,
1 2 π 0 2 π f ( e m σ ) { cos m σ sin m σ } c t g σ s 2 d σ 1 2 π 0 2 π P n f ( e m σ ) { cos m σ sin m σ } c t g σ s 2 d σ C E n ( f ) ln 2 n .
Using Theorem 5, we obtain the estimate
| R n ( f ) | C ( n + m ) p 1 E n ( f ) ln 2 n .
Setting f W r H α ( γ , 1 ) , we finally have
R n [ W r H α ( [ 0 , 2 π ] ) , 1 ] C ( n + m ) p 1 n r + α ln 2 n .

5. Conclusions

We studied approximate methods for evaluating Cauchy and Hilbert singular and hypersingular integrals with rapidly oscillating kernels. In the case of periodic integrable functions, lower and upper bound quadrature formula estimates have been obtained. Optimals with respect to order quadrature formulas for certain classes of functions have been constructed. We developed a method for constructing and estimating quadrature formulas for hypersingular integrals, based on similar results for singular integrals.
Finally, we point out a few key points of our study presented in this paper:
(1) We introduced a method to estimate below quadrature formulas for evaluating singular and hypersingular integrals with rapidly oscillating kernels (in this paper, a method to obtain lower bound estimates by functional ζ N [ Ψ ] in the class of functions Ψ ). Moreover, these estimates can be obtained from any set of N nodes located in the range of integration and N values of integrand function.
The method can be extended to singular and hypersingular integrals defined on other varieties, to polysingular and polyhypersingular integrals and to many dimensional singular and hypersingular integrals. The existence of lower bound estimates of functional ζ N [ Ψ ] allows us to construct an optimal with respect to order (to accuracy) passive algorithms for evaluating corresponding integrals in the classes of functions Ψ .
(2) We proposed a method to construct quadrature formulas for evaluating hypersingular integrals and their error estimates based on quadrature formulas for evaluating singular integrals.
(3) We proposed optimals with respect to order quadrature formulas, which are apparently the most effective among known methods for evaluating singular and hypersingular integrals with rapidly oscillating kernels. We made a comparison of the efficiency of quadrature Formulas (29), (30) and (58) with well-known rules.
Using Levin’s method, one must analytically solve the equation
x ( τ ) + i w g ( τ ) x ( τ ) = f ( τ ) ( τ t ) p ,
where g and f are known functions.
Applying numerical methods to solve this equation might cause some difficulties due to singularity at point t on the right-hand side of the equation.
Thus, Levin’s method has some application limitations.
In [16], an interpolation quadrature formula was constructed to evaluate integrals of the type
1 1 f ( x ) d x 1 x 2 ( x a )
with fixed singularity a.
When constructing a quadrature formula, function f ( x ) is approximated by an interpolation polynomial with n + 1 nodes. A set of nodes contains a particular point a.
Thus, in [16], it was necessary to construct the quadrature formula for each value of 1 < t < 1 .
Implementation of other methods, constructed in the cited literature here, causes greater smoothness of integrand functions than in our computations.

Author Contributions

I.B. conceived of the presented idea and developed the theory. I.B. and V.R. performed the computations and verified the analytical methods. I.B. wrote the manuscript with support from V.R. and A.B. All authors discussed the results and contributed to the final manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Boykov, I.V. Optimal with Respect to Accuracy Algorithms for Approximate Computation of Singular Integrals; Saratov State University Press: Saratov, Russia, 1983. (In Russian) [Google Scholar]
  2. Boikov, I.V. Numerical methods of computation of singular and hypersingular integrals. Int. J. Math. Math. Sci. 2001, 28, 127–179. [Google Scholar] [CrossRef] [Green Version]
  3. Boykov, I.V. Approximate Methods for Evaluation of Singular and Hypersingular Integrals, Part 1, Singular Integrals; Penza State University Press: Penza, Russia, 2005. (In Russian) [Google Scholar]
  4. Boykov, I.V. Approximate Methods for Evaluation of Singular and Hypersingular Integrals, Part 2, Hypersingular Integrals; Penza State University Press: Penza, Russia, 2009. (In Russian) [Google Scholar]
  5. Boykov, I.V.; Ventsel, E.S.; Boykova, A.I. Accuracy optimal methods for evaluating hypersingular integrals. Appl. Numer. Math. 2009, 59, 1366–1385. [Google Scholar]
  6. Boykov, I.V.; Aikashev, P.V. Approximate methods for calculating hypersingular integrals. Univ. Proc. Volga Reg. Phys. Math. Sci. Math. 2021, 1, 66–84. (In Russian) [Google Scholar] [CrossRef]
  7. Khubezhty, S.S. Quadrature Formulas for Singular Integrals and Some of Their Applications; Southern Mathematical Institute: Vladikavkaz, Russia, 2011; Volume 13. (In Russian) [Google Scholar]
  8. Lifanov, I.K.; Vainikko, G.; Poltavskii, L.N. Integral Equations and Their Applications; CRC Press: Boca Raton, FL, USA, 2003. [Google Scholar]
  9. Capobianco, M.R.; Criscuolo, G. On quadrature for Cauchy principal value integrals of oscillatory functions. J. Comput. Appl. Math. 2003, 156, 471–486. [Google Scholar] [CrossRef] [Green Version]
  10. Chung, K.; Evans, G.A.; Webster, J.R. A method to generate generalized quadrature rules for oscillatory integrals. Appl. Numer. Math. 2000, 34, 85–93. [Google Scholar] [CrossRef]
  11. Hota, M.K.; Soha, A.K.; Ojha, P.; Mohanty, T. On the approximate evolution of oscillatory-singular integrals. Cogent Math. 2017, 4, 1314066. [Google Scholar]
  12. Keller, P. A practical algorithm for computing Cauchy principal value integrals of oscillatory functions. Appl. Math. Comput. 2012, 218, 4988–5001. [Google Scholar]
  13. Milovanovic, G.V. Numerical calculation of integrals involving oscillatory and singular kernels and some applications of quadratures. Comput. Math. Appl. 1998, 36, 19–39. [Google Scholar] [CrossRef] [Green Version]
  14. Okecha, G.E. Quadrature formulae for Cauchy principal value integrals of oscillatory kind. Math. Comp. 1987, 49, 259–268. [Google Scholar]
  15. Okecha, G.E. Hermite Interpolation and a method for evaluating Cauchy principal value integrals of oscillatory kind. Kragujev. J. Math. 2006, 29, 91–98. [Google Scholar]
  16. Ryabov, V.M. Quadrature formulae for evaluating singular integrals of oscillating functions. Comput. Math. Math. Phys. 2006, 36, 1011–1015. [Google Scholar]
  17. Wang, H.; Xiang, S. On the evaluation of Cauchy principal value integrals of oscillatory functions. J. Comput. Appl. Math. 2010, 234, 95–100. [Google Scholar] [CrossRef] [Green Version]
  18. Lu, S.; Ding, Y.; Yan, D. Singular Integrals and Related Topics; World Scientific Publications: Singapore, 2007. [Google Scholar]
  19. Lebedev, V.L.; Baburin, O.V. Calculation of the principal values, weights and nodes of the Gauss quadrature formulae of integrals. USSR Comput. Math. Math. Phys. 1965, 5, 81–92. [Google Scholar] [CrossRef]
  20. Wang, Y.; Xiang, S. A Levin method for logarithmically singular oscillatory integrals. arXiv 2019, arXiv:1901.05192v1. [Google Scholar]
  21. Levin, D. Procedures for computing one- and two-dimensional integrals of functions with rapid irregular oscillations. Math. Comput. 1982, 38, 531–538. [Google Scholar]
  22. Levin, D. Analysis of a collocation method for integrating rapidly oscillatory functions. J. Comput. Appl. Math. 1997, 78, 131–138. [Google Scholar]
  23. Hadamard, J. Lecons sur la Propagation des Ondes et les Equations de l’Hydrodynamique; Reprinted by Chelsea: New York, NY, USA, 1949. [Google Scholar]
  24. Hadamard, J. A Study of the Psychology of the Invention Process in the Field of Mathematics; Trans. from French 1959; Publishing House “Soviet Radio”: Moscow, Russia, 1970. [Google Scholar]
  25. Chikin, L.A. Special cases of the Riemann boundary value problems and singular integral equations. Sci. Notes Kazan State Univ. 1953, 113, 53–105. (In Russian) [Google Scholar]
  26. Bakhvalov, N.S. Properties of optimal methods for the solution of problems of mathematical physics. USSR Comput. Math. Math. Phys. 1970, 10, 555–568. [Google Scholar]
  27. Babenko, K.I. Fundamentals of Numerical Analysis; Research Center “Regular and Chaotic Dynamics”: Moscow-Izhevsk, Russia, 2002. (In Russian) [Google Scholar]
  28. Sukharev, A.G. Minimax Algorithms in Problems of Numerical Analysis; Nauka: Moscow, Russia, 1989. (In Russian) [Google Scholar]
  29. Traub, J.F.; Werschulz, A.G. Complexity and Information; Cambridge University Press: Cambridge, UK, 1998. [Google Scholar]
  30. Brab, H. Bounds for Peano Kernals. In Numerical Integration IV, INSM 112; Brab, H., Hamerlin, G., Eds.; Birkhauser: Basel, Switzerland; Boston, MA, USA; Berlin, Germany, 1993; pp. 39–55. [Google Scholar]
  31. Boikov, I.V.; Syomov, M.A. On a method of calculation of hypersingular integrals. Russ. Math. 2016, 60, 1–13. [Google Scholar]
  32. Mikhlin, S.G. Singular integral equations. Uspekhi Mat. Nauk 1948, 3, 29–112. (In Russian) [Google Scholar]
  33. Gaier, D. Methoden der Konformen Abbildung; Springer: Berlin/Heidelberg, Germany, 1964. [Google Scholar]
Figure 1. Convergence of the quadrature formula (55) to the exact value of the integral (54) for m = 11 and n = 1 , 2 , 8 .
Figure 1. Convergence of the quadrature formula (55) to the exact value of the integral (54) for m = 11 and n = 1 , 2 , 8 .
Axioms 11 00150 g001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Boykov, I.; Roudnev, V.; Boykova, A. Approximate Methods for Calculating Singular and Hypersingular Integrals with Rapidly Oscillating Kernels. Axioms 2022, 11, 150. https://doi.org/10.3390/axioms11040150

AMA Style

Boykov I, Roudnev V, Boykova A. Approximate Methods for Calculating Singular and Hypersingular Integrals with Rapidly Oscillating Kernels. Axioms. 2022; 11(4):150. https://doi.org/10.3390/axioms11040150

Chicago/Turabian Style

Boykov, Ilya, Vladimir Roudnev, and Alla Boykova. 2022. "Approximate Methods for Calculating Singular and Hypersingular Integrals with Rapidly Oscillating Kernels" Axioms 11, no. 4: 150. https://doi.org/10.3390/axioms11040150

APA Style

Boykov, I., Roudnev, V., & Boykova, A. (2022). Approximate Methods for Calculating Singular and Hypersingular Integrals with Rapidly Oscillating Kernels. Axioms, 11(4), 150. https://doi.org/10.3390/axioms11040150

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop