Abstract:
We study approximations on the interval $[-1,1]$ of singular integrals of the form
$$
\widehat{f}(x)=\int_{-1}^{1}\frac{f(t)}{t-x}\sqrt{1-t^2}\,dt,
\qquad x \in [-1,1],
$$
by two rational integral operators related to each other in a certain sense. The first is the Fourier–Chebyshev integral operator associated with the Chebyshev–Markov system of rational functions. The second operator is its image under the transformation by the singular integral under consideration.
Approximative properties of the corresponding polynomial analogues of both operators are studied in the case where the density of the singular integral satisfies a Hölder condition of exponent $\alpha \in (0,1]$ on $[-1,1]$.
Rational approximations on $[-1,1]$ of the singular integral with power-law singular density are investigated. In the two cases under consideration the approximating rational functions have arbitrary many fixed geometrically different poles or the
parameters of the approximating rational functions are modifications of the ‘Newman’ parameters.
Bibliography: 34 titles.
Keywords:singular integral on an interval, Fourier–Chebyshev rational integral operators, uniform estimate, Laplace method, strong asymptotics.
Various problems in mathematics and its applications involve singular integrals with Cauchy-type kernel of the form
$$
\begin{equation}
\widehat{f}(x)=\int_{-1}^{1}\frac{f(t)}{t-x}\sqrt{1-t^2}\,dt, \qquad x \in [-1,1],
\end{equation}
\tag{1.1}
$$
which are understood in the sense of the Cauchy principal value. A sufficient condition for such an integral to exist is that the density $f(t)$ obeys a Hölder condition with some exponent (see [1] and [2]). The above integrals can be evaluated precisely (in closed form) only in quite rare particular cases. This calls for the development of approximate methods.
At present, the technique of numerical integration of singular integrals of the form (1.1) is well developed (see, for example, [3]–[8]). For a sufficiently complete account of the results in this direction, see Gabdulkhaev [9].
In 1993 Rusak [10] proposed a method of rational approximation of singular integrals of the form (1.1) with density $f(t)$ in certain function classes on an interval. His studies were continued by his student Boksha [11] and also in the paper [12]. Pointwise approximations by algebraic polynomials of function classes defined by singular integrals of the form (1.1) were studied by Motornyi [13], who obtained asymptotically sharp estimates for approximations. Note that the above studies did not use the classical methods based on Fourier series.
In 1925–1926 Takenaka [14] and Malmquist [15] introduced an orthogonal system of rational functions on the unit circle which generalizes the trigonometric system. Dzhrbashyan [16] constructed rational Fourier series in this system, evaluated the Dirichlet integral, and established analogues of the Jordan–Dirichlet and Dini–Lipschitz tests under the assumption that the poles of the rational functions have no accumulation points on the unit circle. Dzhrbashyan and Kitbalyan [17] constructed systems of rational functions which are orthogonal on $[-1,1]$ and generalize the classical Chebyshev polynomials of the first and second kind.
In 1979 Rovba [18] introduced an integral operator associated with the system of Chebyshev–Markov rational functions and generalizing partial sums of the polynomial Fourier–Chebyshev series. Let $\{a_k\}_{k=1}^n$ be an arbitrary set of numbers, where either the $a_k$, $|a_k|<1$, are real, or they are complex conjugates of one another. On the set of functions $f(x)$ integrable with weight $(1-x^2)^{-1/2}$ on $[-1,1]$ consider the rational integral operator (see [18])
and $A_n$ is the set of parameters $(a_1,\dots,a_{n})$. In particular, if $a_k=0$, $k=1,2,\dots,n$, then $s_n(f,x)$ is a partial sum of the polynomial Fourier–Chebyshev series.
Along with the integral representation (1.2), upper estimates for approximations by the operator $s_n(\,\cdot\,{,}\,\cdot\,)$ of some function classes on an interval were also obtained in [18]. Rational integral operators (1.2) are widely useful in rational approximation [19]–[21].
On the basis of Fourier–Chebyshev rational integral operator (1.2), we introduce the operator
$$
\begin{equation}
\widehat{s}_{n+1}(f,x)=\int_{-1}^{1}\frac{s_{n}(f,t)}{t-x}\sqrt{1-t^2}\,dt, \qquad x \in (-1,1).
\end{equation}
\tag{1.4}
$$
It is known (see, for example, [10]) that $s_{n+1}(f,x)$ is also a rational function of order $\leqslant n+1$ and has the same poles as ${s}_{n}(f,x)$.
It would be interesting to study approximations of singular integrals (1.1) on $[-1,1]$ by the Fourier–Chebyshev rational integral operator (1.2) and the above operator (1.4). In our paper we consider both approaches. As a separate problem, we study approximations of individual singular integrals (1.1) in which the density has a power-law singularity. We consider the cases where the approximating rational function has constraints on the number of geometrically different parameters and where the parameters of the approximating function are certain modifications of the ‘Newman’ parameters.
$$
\begin{equation*}
I_n(x,t)=\int_{-\pi}^{\pi}\frac{\cos u - \cos v}{\cos \tau - \cos v}D_n(v,u)\,dv, \qquad x=\cos u, \quad t=\cos \tau.
\end{equation*}
\notag
$$
Note that the integral $I_n(x,t)$ is singular and is understood in the sense of the Cauchy principal value. Let us evaluate this integral. To do this we change the integration variable to $\zeta=\mathrm{e}^{ i v}$ and set $\xi=\mathrm{e}^{ i u}$ and $z=\mathrm{e}^{ i \tau}$. Then we have
$\Gamma=\{\zeta\colon \zeta=\mathrm{e}^{i v}, -\pi \leqslant v \leqslant \pi \}$. The integral $I_n(x,t)$ contains the points of singularity $\zeta=z$ and $\zeta=1/z$. To evaluate this integral we split it into two. We have
Here we choose the branch of the root function so that $\sqrt{z} \big|_{z=1}=1$ and $\sqrt{\xi} \big|_{\xi=1}=1$. Substituting $I_n(x,t)$ into the integral representation (2.2) we obtain
Note that integral $I_{n+1}(x,v)$ is singular and is understood in the sense of the Cauchy principal value. Putting $t=\cos \tau$ and $x=\cos u$ we can write this integral as
where $\Gamma=\{z\colon z=\mathrm{e}^{i \tau}, -\pi \leqslant v \leqslant \pi \}$.
The integral $I_{n+1}(x,v)$ contains the points of singularity $z=\xi$ and $z=1/\xi$. Our further analysis depends on the method proposed in [22], p. 115. From (1.2) it follows that $I_{n+1}(x,v)$ is a rational function of $\xi=\mathrm{e}^{i u}$ and $\zeta=\mathrm{e}^{i v}$ with first-order poles at the points $z_k$ and $\overline{z}_k$, $k=1,2,\dots,n-1$ (see (1.3)). So it suffices to evaluate the integral $I_{n+1}(x,v)$ for $\zeta=\delta \mathrm{e}^{i v}$, $\delta \in (0,1)$, $v \in (0,\pi)$, and then let $\delta \to 1$. In view of the above we write the integral $I_{n+1}(x,v)$ as
The integrals $J^{(1)}(x,v)$ and $J^{(2)}(x,v)$ contain the points of singularity $z=\xi$ and $z=1/\xi$ on the boundary of the unit disc. To evaluate these integrals we use Sokhotskii’s formulae. Proceeding as with the expression (2.4), we find that
where $\lambda_{n}(v,u)$ was defined in (1.3) and $z_n=0$.
Inserting the last equality into (2.9), splitting the resulting expression into three integrals and changing the variable $v$ to $ -v$ in the integral corresponding to the third term, we obtain
are the approximations of the singular integral (1.1) by the polynomial integral operator that is the image under the transformation (1.1) of partial sums of the Fourier–Chebyshev series. As a result, we have the following statement.
Corollary 2. The following integral representation holds:
§ 3. Approximations of a singular integral with density satisfying a Hölder condition (the polynomial case)
Consider the classes $H^{(\alpha)}[-1,1]$, $\alpha \in (0,1]$, of functions $f(x)$ satisfying the Hölder condition of exponent $\alpha$ with constant 1, that is,
Let us study the approximations (2.7) and (2.14) in the case where $f \in H^{(\alpha)}[-1,1]$, $\alpha \in (0,1]$.
Theorem 3. The approximations on $[-1,1]$ of the singular integral (1.1) with density $f \in H^{(\alpha)} [-1,1]$, $\alpha \in (0,1]$, by partial sums of the Fourier–Chebyshev series satisfy the upper estimate
$$
\begin{equation}
\varepsilon_n^{(0)}(\widehat{f},x) =-\frac{1}{4}\bigl(\sin u I_1 +2 \cos u I_2\bigr), \qquad x=\cos u, \quad x \in [-1,1],
\end{equation}
\tag{3.2}
$$
Consider the integral $I_{11}$ (the arguments for the integral $I_{12}$ are the same). We use the method proposed by Besov in [23] and [24]. The integrand is $2\pi$-periodic, and we can write the integral $I_{11}$ as
Consider the integral $I_{112}$. As before, the integrand is $2\pi$-periodic, and so it suffices to estimate the absolute value of the integral over $[\pi/\lambda_1,\pi]$ (the estimate for the rest of the integral is the same). We have
where $n_0(\alpha)$ is a natural number depending only on $\alpha$.
The proof of the above lemma is similar to the proof of estimate (3.14), and we omit it.
Now we return to the proof of Theorem 3. To derive estimates (3.1) it suffices to substitute (3.3) and (3.15) into (3.2), collect the resulting estimates, and perform some algebraic transformations.
The following theorem gives an estimate for the approximations (2.14) under the condition that $f \in H^{(\alpha)}[-1,1]$, $\alpha \in (0,1]$. We formulate this result without proof because it is similar to the proof of Theorem 3.
Theorem 4. The approximations (2.14) on $[-1,1]$ of the singular integral (1.1) with density $f \in H^{(\alpha)} [-1,1]$, $\alpha \in (0,1]$, satisfy the upper estimate
where $n_0(\alpha)$ is a natural number depending only on $\alpha$.
Remark 1. Estimates (3.1) and (3.16) depend substantially on the position of the point $x$ on $[-1,1]$ (the approximations at the endpoints of the interval have a greater rate of decay than the ones on the whole interval).
§ 4. Approximation of singular integrals with density $|x|^s$ by the Fourier–Chebyshev rational integral operator
4.1. An integral representation and estimates of approximations
Setting $f_s(t)=|t|^s$, $s \in (0,+\infty)$, in representation (1.1) we obtain
We look at the properties of the approximations (2.1) in this case. Since $f_s(t)=|t|^s$, we have
$$
\begin{equation*}
\varepsilon_n(\widehat{f}_{s},x,A_n) =-\frac{1}{2}\int_{-\pi}^{\pi}|{\cos \tau}|^s \sin \tau \frac{\cos \lambda_n(\tau,u)}{\sin \frac{\tau-u}{2}}\,d\tau, \qquad x=\cos u, \quad x \in [-1,1].
\end{equation*}
\notag
$$
Splitting the integral into the three integrals over the intervals $[-\pi,-\pi/2]$, $[-\pi/2,\pi/2]$ and $[\pi/2,\pi]$ and making the changes $\tau+\pi \mapsto \tau$ in the first integral and ${\tau-\pi \mapsto \tau}$ in the third we have
Note that the integral on the right is singular and is understood in the sense of the Cauchy principal value. For further purposes, we need to choose the parameters of the approximating rational function in a special way. We set $n \mapsto 2n-1$. Assume that the $2n-1$ parameters $\{z_k\}_{k=1}^{2n-1}$ have the form
Theorem 5. Under conditions (4.2), for the approximations of the function $\widehat{f}_s(x)$ (see (4.1)) on the interval $[-1,1]$ by the Fourier–Chebyshev rational integral operator (1.2) we have:
where $C=\{z\colon z=\mathrm{e}^{i \tau},\,-{\pi}/{2} \leqslant \tau \leqslant {\pi}/{2}\}$ is the right semicircle of the unit circle which is traversed anticlockwise. Note that in the last integral the integrand has branch points at $z=0$, $z=\infty$ and $z=\pm i$. It it also worth noting that $z=\pm \xi$ are simple poles of the integrand, but only one of these poles lies on the contour of integration. We can assume without loss of generality that this is the pole $z=\xi$. We split the integral (4.9) into two so that
Let us study each integral separately. Consider the domain described by the contour $\Gamma=C\cup C_{\delta_1}^{-}\cup l_1^{-}\cup C_\delta^{-}\cup l_2^{-}\cup C_{\delta_2}^{-}$ (see Figure 1), where
In the above domain $g(z,s)$ splits into the regular branches specified by the conditions $g(1,s)=\mathrm{e}^{2\pi k s i}$, $k\in \mathbb{Z}$. Let $g_0(z,s)$ be the branch satisfying $g_0(1,s)=1$. Then the function
is regular inside this domain and has a singular point on the boundary, which is a simple pole at $z=\xi$. Let us employ Sokhotskii’s formulae. First assume that $\xi \in D$. By Cauchy’s residue theorem we have
Since $2+2p-s>0$ (see (4.2)), the integral over the semicircle $C_\delta$ tends to zero as $C_\delta$ shrinks to a point. Proceeding similarly with the integrals over the arcs $C_{\delta_1}$ and $C_{\delta_2}$ we obtain the asymptotic equalities
Hence the corresponding integrals tend to zero as the arcs $C_{\delta_1}$ and $C_{\delta_2}$ shrink to points. Thus, as $\delta, \delta_1, \delta_2 \to 0$, (4.11) assumes the form
Now we deal with the integral $J_2$. Consider the domain bounded by the contour $\Gamma= l_3^{-}\cup C_{\delta_3}^{-}\cup C^{-}\cup C_{\delta_4}^{-}\cup l_4^{-}\cup C_R$ (see Figure 2), where
In the above domain the function $g(z,s)$ splits into the regular branches specified by the conditions $g(1,s)=\mathrm{e}^{2\pi k s i}$, $k\in \mathbb{Z}$. Let $g_0(z,s)$ be the branch for which $g_0(1,s)=1$. Inside the domain the function
is regular, and on its boundary it has a singular point, which is a simple pole at ${z=\xi}$. We use Sokhotskii’s formulae again. We first set $\xi \in D$, where $D=\{z\colon |z|<1, \mathrm{Re}\,z >0\}$. By Cauchy’s integral formula we have
Proceeding as with the integral $J_1$, we conclude that the integrals over the arcs $C_{\delta_3}$ and $C_{\delta_4}$ tend to zero as $\delta_k \to 0$, $k=3,4$. Now consider the integral over $C_R$. Setting $z=R \mathrm{e}^{i \tau}$ yields
$$
\begin{equation*}
\begin{aligned} \, &\int_{C_R}\varphi_2(z,\xi)\,dz \\ &\quad=\int_{-{\pi}/{2}}^{{\pi}/{2}} \frac{((R \mathrm{e}^{ i \tau})^2+1)^s((R \mathrm{e}^{ i \tau})^2-1)} {(R \mathrm{e}^{ i \tau})^{s+3}((R \mathrm{e}^{ i \tau})^2 - \xi^2)} (R \mathrm{e}^{ i \tau})^{-2p}\prod_{k=p+1}^{n-1} \frac{1 + \alpha_k^2 (R \mathrm{e}^{ i \tau})^2}{(R \mathrm{e}^{ i \tau})^2 + \alpha_k^2} Ri \mathrm{e}^{ i \tau}\,d\tau. \end{aligned}
\end{equation*}
\notag
$$
We have $2-s+2p>0$ (see (4.2)), and so the integral over the semicircle $C_R$ tends to zero as $R \to \infty$. From representation (4.13), as $\delta_k \to 0$, $k=3,4$, and $R \to \infty$, we obtain
where the integrals on the right are taken along the corresponding rays on the imaginary axis. The further change of the variable $z \mapsto {z^{-1}}$ gives
Now, to deduce the integral representation (4.5) it suffices to note that the terms in square brackets in the integrand are complex conjugates of each other and then make the corresponding transformations.
From the integral representation (4.5) we obtain estimates (4.6) and (4.7).
are pointwise and uniform approximations, respectively, of the function $\widehat{f}_s(x)$, ${s \in (0,+\infty)}$, by the polynomial Fourier–Chebyshev sums.
Corollary 3. The approximations of the function $\widehat{f}_s(x)$ on $[-1,1]$ by the Fourier–Chebyshev sums of order $n$, $n>s/2$, satisfy:
We consider the asymptotic behaviour of (4.8) as $n \to \infty$. To do this we change the variable by the formula $t=\sqrt{(1-y)/(1+y)}$, $dt=-dy/((1+y)^{3/2} (1-y)^{1/2})$, on the right-hand side of the integral representation (4.8). As a result,
4.2. An asymptotic expression for the majorant when the approximating function has a fixed number of poles
Let $n>p$, $p=[s/2]$, $n_1=n-p-1$, let $q$ be an arbitrary integer such that $0 < q < n_1$, and let $A_q$ be a set of parameters $(\beta_1,\beta_2,\dots,\beta_{n_1})$ precisely $q$ among which are different, with the multiplicity of each parameter equal to $m$: $n_1=m q$. So we are concerned with approximation by rational functions with a pole of order $2p+2$ at infinity and $2q$ geometrically different poles of multiplicity $m$ each in the extended complex plane.
Note that Lungu [25], [26] was the first to consider the problem of the approximation of continuous functions with typical singularities by rational functions with a fixed number of geometrically different poles.
Under the above assumptions the integral representation (4.15) can be written as
Note that in this case, for each $n \in \mathbb{N}$ the corresponding set of parameters $(\beta_1,\beta_2,\dots,\beta_q)$ can depend on $n$, that is, $\beta_k=\beta_k(n)$, $k=1,2,\dots,q$. We also assume that the following conditions are satisfied:
and $b_j$ is the unique root of the equation $\sum_{k=1}^{q}\beta_{k}/(u^2 - \beta_{k}^2)=0 $ in the interval $(\beta_{j+1},\beta_j)$, $j=1,2,\dots,q-1$.
Proof. We write the integral representation of the approximations (4.17) in the form
Let us study the asymptotic behaviour of each of the above three expressions as $m\to \infty $. To do this we employ the Laplace method [27], [28]. We require the following three lemmas (for a proof, see [19]).
is an asymptotic estimate for the majorant of uniform approximations of the function $\widehat{f}_s(x)$ by partial sums of the polynomial Fourier–Chebyshev series.
Corollary 4. The uniform approximations of the function $\widehat{f}_s(x)$ on the interval $[-1,1]$ by the Fourier–Chebyshev sums of order $n$, $ n>s/2$, satisfy the upper estimate
$$
\begin{equation*}
\varepsilon_{2n-1}^{(0)}(\widehat{f}_{s}) \leqslant2^{1-s}\biggl|\sin\frac{\pi s}{2}\biggr|\frac{\Gamma(s)}{n^s}(1+o(1)), \qquad s \in (0,+\infty), \quad n \to \infty.
\end{equation*}
\notag
$$
4.3. The best majorant of uniform approximations in the case of a fixed number of poles in the approximating function
It would be interesting to minimize the right-hand side of the asymptotic equality (4.19) by choosing a set of parameters $(\beta_1^{*},\beta_2^{*},\dots,\beta_{q}^{*})$ optimal for this problem; that is, it would be interesting to find the best majorant of uniform approximations of the function $\widehat{f}_s(x)$ by Fourier–Chebyshev rational integral operators (1.2) in the case of an arbitrary fixed number of geometrically different poles.
Theorem 7. The uniform approximations of the function $\widehat{f}_s(x)$ by the rational integral operator (1.2) with $2q$ geometrically different poles on $[-1,1]$ satisfy the upper estimate
where the coefficients $c_k$, $k=1,2,\dots,q$, are independent of $m$ (the precise values of the $c_k$ are specified below). Let us find the asymptotic behaviour of the terms $S_n^{(1)}(A_q)$, $S_n^{(2)}(A_q)$ and $S_n^{(3)}(A_q)$ from (4.19) in this case. It is easily seen that $S_n^{(1)}(A_q)$ satisfies the asymptotic equality
(see [19], Lemma 6), which holds for the above parameters $\beta_k$, $k=1,2,\dots,q$. From the last asymptotic equality and Theorem 6 in [19] we find that
Let us choose $c_k$, $k=1,2,\dots,q$, so as to equate the exponents of $m$ in each of the asymptotic equalities for $S_n^{(1)}(A_q)$, $S_n^{(2)}(A_q)$ and $S_n^{(3)}(A_q)$. That is, the $c_k$, $k=1,2,\dots,q$, must satisfy
To show that for the $c_k$ as above, $k=1,2,\dots,q$, the family of parameters $\beta_k$, $k=1,2,\dots,q$, is optimal in the sense that $\varepsilon_{2n-1,2q}^{*}(\widehat{f}_{s},A_q)$ attains its asymptotically minimum value for these $\beta_k$, it suffices to follow the arguments of [29]. Using now the formula $m=(n-1-p)/q$ to get back to the parameter $n$, we have
$$
\begin{equation*}
\varepsilon_{2n-1,2q}^{*}(\widehat{f}_{s}) \sim c(q,s)\biggl(\frac{\log^{2q-1}n}{n^{2q}}\biggr)^{s}, \qquad s \in (0,+\infty), \quad n \to \infty,
\end{equation*}
\notag
$$
where $c(q,s)$ was defined in Theorem 7. Now estimate (4.21) is secured by (4.7). This proves Theorem 7.
It would be interesting to compare the uniform estimate from Theorem 7 with the one for approximations of the function $\widehat{f}_s(x)$ by partial sums of the Fourier–Chebyshev series from Corollary 4. The decay rate of uniform approximations in the polynomial case is $O(1/n^{s})$, whereas the corresponding rational approximations with $2q$ geometrically different poles in the approximating function in the open complex plane decay as $O(({\log^{2q-1}n}/n^{2q})^s)$.
4.4. Asymptotic expression for the majorant in the general case
Let us study the asymptotic behaviour of (4.8) as $n \to \infty$ without any constraints on the number of poles. To do this we use representation (4.15) again. Let the parameters $\beta_k$, $k=1,2,\dots,n_1$, $n_1=n - 1 - p$, be of the form $\beta_{k}=\xi^k$, $\xi \in (0,1)$, $k=1,2,\dots,n_1$. In addition, for each $n$, we can choose its own family of parameters $\beta_k$, that is, $\beta_k=\beta_k(n)$. So we assume that
Conditions (4.22)–(4.24) are consistent. Below we consider a sequence of parameters $\{\beta_{k}\}_{k=1}^{n_1}$ of the approximating function that obey these conditions.
Theorem 8. Under conditions (4.22)–(4.24) the uniform approximations of the function $\widehat{f}_s(x)$ on $[-1,1]$ by the Fourier–Chebyshev rational integral operator satisfy the upper estimate
Let us use the Laplace method [27], [28]. The function $S(y)$ decreases on the interval $[0, \beta_{n_1}]$ since ${S}'(y)<0$, and therefore it attains its maximum value at $y=0$. Using the expansion
where $c_1(s)$ was defined in Theorem 8. It is known (see, for example, [30]) that if $\beta_{k}=\xi^k$, $\xi \in (0,1)$, $k=1,2,\dots,n_1$, then the products on the right-hand side of (4.28) satisfy the estimate
Again, we use the Laplace method [27], [28]. The function $S(y)$ is increasing for $y \in (\beta_{1},1)$ since ${S}'(y)=2\sum_{k=1}^{n_1}\beta_k/(y^2-\beta_{k}^{2})>0$. Hence it assumes its maximum value at $y=1$. From the asymptotic equalities
where $\varphi(n,\varepsilon)=\varepsilon \sqrt{\sum_{k=1}^{n_1}\beta_k/(1-\beta_k^2)} \to \infty $ as $n \to \infty $ in view of condition (4.24). Next,
Substituting the asymptotic estimates (4.27), (4.29) and (4.30) for the integrals $I_n^{(4)}$, $I_n^{(5)}$ and $I_n^{(6)}$ into representation (4.26) we arrive at (4.25).
Let us consider the right-hand side of estimate (4.25) in the case where the parameters $\beta_{k}$, $k=1,2,\dots,n_1$, are modifications of those introduced by Newman in [31] in his famous paper on rational approximation of the function $|x|$. This problem is investigated below.
Let $A_N$ be a family of parameters $\beta_k$, $k=1,2,\dots,n_1$, such that
where $n_1\in \mathbb{N}$ is fixed and $c$ is some positive constant independent of $n$. Note that this family of parameters satisfies conditions (4.22).
Theorem 9. The uniform approximations of the function $\widehat{f}_s(x)$ on $[-1,1]$ by the Fourier–Chebyshev rational integral operator (1.2) with parameters (4.31) satisfy the upper estimate
Proof. We examine the asymptotic behaviour of the right-hand side of (4.25) when the parameters $\beta_{k}$, $k=1,2,\dots,n_1$, are given by (4.31). The proof of the following two lemmas can be found in [32].
Lemma 6. Under conditions (4.31) the following asymptotic equality holds:
Now estimate (4.32) is a consequence of the last relation and (4.7). This proves Theorem 9.
It would be interesting to minimize the majorant in the upper estimate (4.32) by choosing an optimal parameter $c$ for this problem; in other words we are interested in a sharp upper estimate for the uniform approximations with parameters (4.31).
Theorem 10. For the uniform approximations of the function $\widehat{f}_s(x)$ on $[-1,1]$ by the Fourier–Chebyshev rational integral operator (1.2) there exists a family of parameters $A_N^{*}$ of the form (4.31) such that
Proof. The constant $c$ on the right-hand side (4.32) is chosen from the condition of the equality of the arguments of the exponential in each term. Therefore,
To verify that for this $c^{*}$ the parameters of the approximating rational function are indeed optimal, namely, they deliver the best uniform approximations by the rational integral operator (1.2) with parameters satisfying condition (4.31), it suffices to employ the method proposed in [29]. Since $n_1=n - 1 - p$ and $n - 1 - p \sim n$ as ${n \to \infty}$, we have (4.36). This completes the proof of Theorem 10.
Corollary 5. For the uniform approximations of functions (1.1) with density $|x|^s$, $s \in (0,2)$, on $[-1,1]$ by the Fourier–Chebyshev rational integral operator (1.2) there exists a family of parameters of the approximating function such that
(see Stahl [33]), where the $R_{n,n}(|x|^s,[-1,1])$ are best uniform approximations of the function $|x|^s$ on $[-1,1]$ by rational functions of degree $\leqslant n$.
§ 5. Approximation of singular integrals with density $|x|^s$ by the rational integral operator $\widehat{s}_{n+1}$
5.1. Integral representations and estimates for approximations
We consider approximations (2.8) of the function $\widehat{f}_s(x)$ (see (4.1)) on $[-1,1]$ by the rational integral operator (1.4).
For further purposes, we need to choose in a special way the parameters of the approximating rational function. Let the family of parameters $\{ z_k\}_{k=1}^{2n}$ (see (1.3)) be defined by
where the ${\varepsilon}_{2n}(f_s,t,A_n)$ are approximations of the function $|x|^s$, $s \in (0,+\infty)$, on the interval $[-1,1]$ by the rational integral operator (1.2).
Note that the integrals $J^{(1)}(y,x)$ and $J^{(2)}(y,x)$, which contain points of singularity $z=\xi$ and $z=1/\xi$ on the boundary of the unit disc, are understood in the sense of the Cauchy principal value. To evaluate these integrals, we use the Sokhotskii’s formulae. After some calculations we obtain
are, respectively, the pointwise and uniform approximations of the function $\widehat{f}_s(x)$ on the interval $[-1,1]$ by the polynomial analogue of operator (1.4).
Corollary 6. For the approximations of the function $\widehat{f}_s(x)$ on $[-1,1]$ by the operator which is the image under the transformation (1.1) of partial sums of the polynomial Fourier–Chebyshev series of order $n$, $n>s/2$, we have:
Let us study the asymptotic behaviour of (5.6) as $n \to \infty$. To do this we change the variable to $t=\sqrt{(1-y)/(1+y)}$, $dt=-dy/((1+y)^{3/2} (1-y)^{1/2})$, on the right-hand side of (5.6). Then we have
5.2. An asymptotic expression for the majorant in the case of an approximating function with a fixed number of poles
Let $n>p$, $p=[s/2]$, $n_1=n-p$, let $q$ be an arbitrary fixed natural number, and let $A_q$ be a set of parameters $(\beta_{1},\beta_2, \dots, \beta_{n_1})$ that contains precisely $q$ different parameters and the multiplicity of each parameter is $m$, $n_1=m q$. So we are concentred here with approximations by rational functions with a pole of order $2p$ at infinity and with $2q$ geometrically different poles (each of multiplicity $m$) in the extended complex plane.
In view of the above assumptions, (5.10) can be written as
In the next theorem we find an asymptotic expression for the majorant of uniform approximations of the function $\widehat{f}(x)$ on the interval $[-1,1]$ by the rational integral operator (1.4) with an arbitrary fixed number of geometrically different poles. The theorem is stated without a proof, because it is based on the method described in detail in [19] and used already for the verification of Theorem 6.
Theorem 12. For each natural number $q$, $0 < q<n_1$, $n_1=n-p$, $n_1=mq$, and any $s \in (0,+\infty)$, the following asymptotic equality holds:
and $b_j$ is the unique root of the equation $\sum_{k=1}^{q}\beta_{k}/(u^2 - \beta_{k}^2)=0$ on the interval $(\beta_{j+1},\beta_j)$, $j=1,2,\dots,q-1$.
Putting $\beta_j=1$ in Theorem 12, $j=1,2,\dots,q$, we find that the quantity
is an upper bound for uniform approximations of the function $\widehat{f}_s(x)$ on the interval $[-1,1]$ by the polynomial operator equal to the image of the partial sums of the Fourier–Chebyshev series under the transformation (1.1).
Corollary 7. The uniform approximations of the function $\widehat{f}_s(x)$ on $[-1,1]$ by the polynomial operator that is the image of the partial sums of the Fourier–Chebyshev series of order $n$, $n>s/2$, under the transformation (1.1), satisfy the estimate
$$
\begin{equation*}
\widehat{\varepsilon}_{2n+1}^{\,(0)}(f_s) \leqslant 2^{1-s}\biggl|\sin\frac{\pi s}{2}\biggr| \frac{\Gamma(s)}{n^{s}}(1+o(1)), \qquad s \in (0,+\infty), \quad n \to \infty.
\end{equation*}
\notag
$$
It would be interesting to compare the estimate in Corollary 7 with the one for uniform polynomial approximations from Corollary 4.
5.3. The best majorant of uniform approximations for an approximating function with a fixed number of poles
Our aim here is to minimize the right-hand side of the asymptotic equality (5.11) by choosing a family of parameters $(\beta_1^{*},\beta_2^{*},\dots,\beta_q^{*})$ optimal for this problem, that is, we find a best estimate for approximations of the singular integrals under consideration by the rational integral operator (1.4) with a fixed number of poles.
Theorem 13. The uniform approximations of the function $\widehat{f}_s(x)$ on $[-1,1]$ by the rational integral operator (1.4) with $q$ geometrically different poles satisfy the upper estimate
Next we choose numbers $c_k$, $k=1,2,\dots,q$, to equate the exponents of $m$ in the expressions for $S_n^{(1)}(A_{q})$, $S_n^{(2)}(A_{q})$ and $S_n^{(3)}(A_{q})$. It is clear that the $c_k$, $k=1,2,\dots,q$, must satisfy
Now, to prove that for the above values of the $c_k$, $k=1,2,\dots,q$, the family of parameters $\beta_k$, $k=1,\dots,q$, is optimal, that is, the quantity $\varepsilon_{n}^{*}(|\cdot|^s,A_{q})$ assumes its asymptotically minimum values for these $c_k$, it suffices to follow the arguments in [29]. Consequently, we have
Let us compare the estimate from Theorem 13 with the one from Corollary 7 for uniform polynomial approximations of the function $\widehat{f}_s(x)$. The rate of decay of polynomial analogues of the above method of rational approximation is $O(1/n^{s})$ uniformly with respect to $x \in [-1,1]$, whereas in the rational case, when the approximating function has $2q$ geometrically different poles in the open complex plane the rate of decay can be increased to $O(({\log^{2q-1}n}/{n^{2q}})^{s})$.
Remark 2. From Theorems 7 and 13 it follows that, for the asymptotic expressions for best majorants of uniform approximations of the function $\widehat{f}_s(x)$ by the above rational integral operators, in the case of an arbitrary fixed number of geometrically different poles we have the equality
$$
\begin{equation*}
\widehat{\varepsilon}_{2n+1,2q}^{\,*}(f_s)=\varepsilon_{2n-1,2q}^{*}(\widehat{f}_s), \qquad n \to \infty.
\end{equation*}
\notag
$$
5.4. The majorant of uniform approximations in the general case
Let us investigate (5.10) as $n \to \infty$ without any constraints on the number of parameters of the approximating function. As above, we assume that the parameters $\beta_k$, $k=1,2,\dots,n_1$, $n_1=n-p$, are ordered as follows:
and, moreover, they have the form $\beta_{k}=\xi^k$, $\xi \in (0,1)$, $k=1,2,\dots,n_1$.
Theorem 14. The uniform rational approximations of the function $\widehat{f}_s(x)$ on $[-1,1]$ by the integral operator (1.4) with parameters (4.22)–(4.24) satisfy the upper estimate
For a proof it suffices to find the asymptotic behaviour of the integral on the right-hand side of (5.10) in the case where the parameters of the approximating function satisfy (4.22)–(4.24). The proof repeats verbatim the argument in Theorem 8 and is omitted.
5.5. The case of ‘Newman’ parameters
Let us investigate the right-hand side of estimate (5.13) in the case where the parameters $\beta_{k}$, $k=1,2,\dots,n_1$, $n_1=n-p$, are modifications of those introduced by Newman in [31]. For each fixed $n\in \mathbb{N}$ let $A_N$ be a family of parameters $\beta_k$ of the form
where $c$ is some positive quantity independent of $n$. Note that some modifications of the parameters $\beta_k$, $k=1,2,\dots,n_1$, were used by Bulanov (see, for example, [34]) for the solution of some approximation problems.
The theorems establishing the corresponding estimates for uniform approximations with parameters (5.14) of the approximating functions are presented given without proofs, because the corresponding arguments are quite similar to the ones used for Theorems 9 and 10.
Theorem 15. The uniform approximations of the function $\widehat{f}_s(x)$ on $[-1,1]$ by the rational integral operator (1.4) with parameters (5.14) satisfy the upper estimate
where $c_1(s)$ is defined in Theorem 11 and $\Gamma (\,\cdot\,)$ is the Euler gamma function.
Next we minimize (5.15) by choosing an optimal parameter $c$ for this problem, that is, we find a best estimate for uniform approximations with parameters (5.14). We set
Theorem 16. For the uniform approximations of the function $\widehat{f}_s(x)$ on $[-1,1]$ by the rational integral operator (1.4) there exists a family of parameters $A_N^{*}$ of the form (5.14) such that the following upper estimate holds:
Corollary 8. For the uniform approximations of functions (1.1) with density $|x|^s$, $s \in (0,2)$, on $[-1,1]$ by the rational integral operator (1.4) there exists a family of parameters of the approximating function such that the upper estimate
Remark 3. It follows from Theorems 10 and 16 that for sharp estimates of uniform approximations of the function $\widehat{f}_s(x)$ by the above rational integral operators in the case of ‘Newman’ poles we have
$$
\begin{equation*}
\widehat{\varepsilon}_{2n+1}^{\,*}(f_s)=\varepsilon_{2n-1}^{*}(\widehat{f}_s), \qquad n \to \infty.
\end{equation*}
\notag
$$
§ 6. Conclusions
In the present paper we have studied approximations on the interval $[-1,1]$ of singular integrals of the form (1.1) by two rational integral operators. The first operator here is the Fourier–Chebyshev rational integral operator [18] associated with the system of Chebyshev–Markov rational functions and which is a natural generalization of partial sums of the polynomial Fourier–Chebyshev series. The second operator is the image of the above integral operator under the transformation (1.1). For each of these operators an integral representation of approximations has been presented.
In the polynomial case, for both operators we have established estimates for approximations on $[-1,1]$ of singular integrals of the form (1.1) with density satisfying a Hölder condition of exponent $\alpha \in (0,1]$.
As a separate problem for both methods we have studied rational approximations in the case where the density of the singular integrals (1.1) has a power-law singularity on the interval $[-1,1]$. Two radically different approaches have been considered. The first approach involves constraints on the number of geometrically different poles in the approximating functions. In the second approach the parameters of the approximating functions are modifications of the ‘Newman’ parameters. In both cases we have established estimates for pointwise approximations and for uniform approximations with a certain majorant, found an asymptotic expression for this estimate and provided the optimal values of the parameters under which the majorant has the greatest decay rate.
As a corollary, we have obtained estimates for approximations by polynomial analogues of both integral operators of singular integrals of the form (1.1) with density having a power-law singularity on $[-1,1]$.
The above findings show that the classes of singular integrals of the form (1.1) with density having a power-law singularity on the interval $[-1,1]$ reflect the peculiarities of the rational approximation by the integral operators under consideration in the following sense: for certain values of the parameters of the approximating functions these function provide approximations of higher order in comparison to their polynomial analogues.
It is worth pointing out that the estimates of best uniform approximations of the function $\widehat{f}_s(x)$ on $[-1,1]$ by the above methods are asymptotically equal.
Bibliography
1.
F. D. Gakhov, Boundary value problems, Pergamon Press, Oxford; Addison-Wesley Publishing Co., Inc., Reading, MA–London, 1966, xix+561 pp.
2.
N. I. Muskhelishvili, Singular integral equations, 3d revised and augmented ed., Nauka, Moscow, 1968, 511 pp. ; English transl. of 2nd ed., P. Noordhoff N. V., Groningen, 1953, vi+447 pp.
3.
F. Erdogan and G. D. Gupta, “On the numerical solution of singular integral equations”, Quart. Appl. Math., 29 (1972), 525–534
4.
D. Elliott and D. F. Paget, “On the convergence of a quadrature rule for evaluating certain Cauchy principal value integrals”, Numer. Math., 23 (1975), 311–319
5.
M. A. Sheshko, “Convergence of quadrature processes for a singular integral”, Soviet Math. (Iz. VUZ), 20:12 (1976), 86–94
6.
A. V. Saakyan, “Quadrature formulae of Gauss's type for singular integrals”, Problems in mechanics of thin deformed bodies, GITUTIUN of the National Academy of Sciences of Armenia, Erevan, 2002, 259–265 (Russian)
7.
Sh. S. Khubezhty, “Quadrature formulas for singular integrals with a Cauchy kernel”, Vladikavkaz. Mat. Zh., 10:4 (2008), 61–75 (Russian)
8.
Sh. S. Khubezhty and A. O. Tsutsaev, “Quadrature formulae for singular integrals which have almost Gaussian accuracy order”, Izv. Vyssh. Uchebn. Zaved. Sev.-Kavkaz. Region Estestv. Nauki, 2015, no. 2, 53–57 (Russian)
9.
B. G. Gabdulkhaev, “Finite-dimensional approximations of singular integrals and direct methods of solution of singular integral and integrodifferential equations”, J. Soviet Math., 18:4 (1982), 593–627
10.
V. N. Rusak, “Uniform rational approximation of singular integrals”, Izv. Akad. Nauk Belarusi Ser. Fiz.-Mat. Nauk, 1993, no. 2, 22–26 (Russian)
11.
A. N. Boksha, “Approximation of singular integrals by rational functions in uniform metric”, Vestn. Beloruss. Gos. Univ. Ser. 1, Fiz. Mat. Inform., 1997, no. 3, 68–71 (Russian)
12.
V. N. Rusak and A. Kh. Uasis, “Rational approximation of singular integral with differentiable density”, Izv. Belarus. Gos. Ped. Univ. Ser. 1 Fiz. Mat. Inform. Biol. Geograf., 59:1 (2009), 8–11 (Russian)
13.
V. P. Motornyi, “Approximation of certain classes of singular integrals by algebraic polynomials”, Ukrainian Math. J., 53:3 (2001), 377–394
14.
S. Takenaka, “On the orthogonal functions and a new formula of interpolation”, Japan. J. Math., 2 (1925), 129–145
15.
F. Malmquist, “Sur la détermination d'une classe de fonctions analytiques par leurs dans un ensemble donné de points”, Comptes rendus du 6ème congrès des mathématiciens scandinaves (Kopenhagen 1925), Det Hoffenbergske Etablissement, Kopenhagen, 1926, 253–259
16.
M. M. Dzhrbashyan, “On the theory of Fourier series in rational functions”, Izv. Akad. Nauk Arm. SSR Ser. Fiz.-Mat. Estestv. Tekhn. Nauki, 9:7 (1956), 3–28 (Russian)
17.
M. M. Dzhrbashyan and A. A. Kitbalyan, “One generalization of Chebyshev polynomials”, Dokl. Akad. Nauk Arm. SSR, 38:5 (1964), 263–270 (Russian)
18.
E. A. Rovba, “One direct method in rational approximation”, Dokl. Akad. Nauk BSSR, 23:11 (1979), 968–971 (Russian)
19.
P. G. Patseika (Potseiko), Y. A. Rouba (E. A. Rovba) and K. A. Smatrytski, “On one rational integral operator of Fourier–Chebyshev type and approximation of Markov functions”, Zh. Beloruss. Gos. Univ. Mat. Inform., 2 (2020), 6–27
20.
P. G. Potseiko and E. A. Rovba, “Approximations on classes of Poisson integrals by Fourier–Chebyshev rational integral operators”, Siberian Math. J., 62:2 (2021), 292–312
21.
P. G. Potseiko and Ye. A. Rovba, “Conjugate rational Fourier–Chebyshev operator and its approximation properties”, Russian Math. (Iz. VUZ), 66:3 (2022), 35–49
22.
V. N. Rusak, Rational functions as approximation apparatus, Belarusian State University, Minsk, 1979, 174 pp. (Russian)
23.
O. V. Besov, “Estimate of the approximation of periodic functions by Fourier series”, Math. Notes, 79:5 (2006), 726–728
24.
O. V. Besov, Lectures on mathematical analysis, 4th revised and augmented ed., Fizmatlit, Moscow, 2020, 476 pp. (Russian)
25.
K. N. Lungu, “On best approximations by rational functions with a fixed number of poles”, Math. USSR-Sb., 15:2 (1971), 313–324
26.
K. N. Lungu, “Best approximations by rational functions with a fixed number of poles”, Siberian Math. J., 25:2 (1984), 289–296
27.
M. A. Evgrafov, Asymptotic estimates and entire functions, 3d revised and augmented ed., Nauka, Moscow, 1979, 320 pp. ; English transl. of 1st ed., Gordon and Breach, Inc., New York, 1961, x+181 pp.
28.
M. V. Fedoryuk, Asymptotics: integrals and series, Nauka, Moscow, 1987, 544 pp. (Russian)
29.
E. A. Rovba and E. G. Mikulich, “Constants in the approximation of $|x|$
using rational interpolation processes”, Dokl. Nats. Akad. Nauk Belarusi, 53:6 (2009), 11–15 (Russian)
30.
A. A. Gončar (Gonchar), “On the rapidity of rational approximation of continuous functions with characteristic singularities”, Math. USSR-Sb., 2:4 (1967), 561–568
31.
D. J. Newman, “Rational approximation to $|x|$”, Michigan Math. J., 11:1 (1964), 11–14
32.
P. G. Potseiko and Y. A. Rovba, “On estimates of uniform approximations by rational Fourier–Chebyshev integral operators for a certain choice of poles”, Math. Notes, 113:6 (2023), 815–830
33.
H. R. Stahl, “Best uniform rational approximation of $x^\alpha$ on $[0,1]$”, Acta Math., 190:2 (2003), 241–306
34.
A. P. Bulanov, “Asymptotics for least deviation of $|x|$ from rational functions”, Math. USSR-Sb., 5:2 (1968), 275–290
Citation:
P. G. Potseiko, E. A. Rovba, “Approximations of one singular integral on an interval by Fourier–Chebyshev rational integral operators”, Sb. Math., 215:7 (2024), 953–992
\Bibitem{PotRov24}
\by P.~G.~Potseiko, E.~A.~Rovba
\paper Approximations of one singular integral on an interval by Fourier--Chebyshev rational integral operators
\jour Sb. Math.
\yr 2024
\vol 215
\issue 7
\pages 953--992
\mathnet{http://mi.mathnet.ru//eng/sm10030}
\crossref{https://doi.org/10.4213/sm10030e}
\mathscinet{http://mathscinet.ams.org/mathscinet-getitem?mr=4813936}
\zmath{https://zbmath.org/?q=an:07945704}
\adsnasa{https://adsabs.harvard.edu/cgi-bin/bib_query?2024SbMat.215..953P}
\isi{https://gateway.webofknowledge.com/gateway/Gateway.cgi?GWVersion=2&SrcApp=Publons&SrcAuth=Publons_CEL&DestLinkType=FullRecord&DestApp=WOS_CPL&KeyUT=001346292600001}
\scopus{https://www.scopus.com/record/display.url?origin=inward&eid=2-s2.0-85208431361}