Let f be a function defined by a power series at the point z=∞:
f(z)=∞∑k=0ckzk.
Given the first N coefficients c0,…,cN−1 of the series in (1), we can construct from them a few constructive approximations to the original function f (see [10]; here we mean ‘constructive’ in the sense specified in [2], § 2). It is natural to wonder if such approximations are optimal. In what follows we assume that N=2n+1=3m+1=4ℓ+1, где n,m,ℓ∈N. This is not a substantative condition and has a purely technical nature (see (2)–(4)).
The best known, popular, and well investigated constructive approximations to a power series are the Padé approximations [n/n]f=Pn/Qn, where the polynomials Pn,Qn≢0, degPn,degQn⩽, are specified by the relation
The convergence of Padé approximations in the class of multivalued analytic functions with a finite number of singular points follows from Stahl’s theory of 1985–1986 [7]. On the other hand there also exist some other methods, which are not as well known, for approximating a power series (1) given the set of its first N coefficients. Here we discuss rational approximations constructed on the basis of Hermite–Padé polynomials of the second type for the systems f, f^2 and f, f^2, f^3. It turns out that, from the standpoint of the optimal use of the N coefficients of the series (1), Hermite–Padé rational approximations have certain advantages over Padé approximations (see (10); also see [4] and [3]). Their convergence has so far been considered only in special cases (see [1], [6], [5], and [4]). The paper [3] presents some numerical results related to van der Pol’s equation. They show that the use of Hermite–Padé polynomials for the systems f, f^2 and f, f^2, f^3 extends significantly one’s abilities to analyse the properties of the function f numerically on the basis of the coefficients of the series (1). In particular, one can distinguish in this way quadratic branch points of a function given by a series (see [9]). The aim of this note is to present theoretical results showing that Hermite–Pade rational approximations are more optimal than Padé approximations.
For the system f, f^2 we define the Hermite–Padé polynomials P^{(2)}_{2m,0}\not\equiv 0, P^{(2)}_{2m,1}, and P^{(2)}_{2m,2}, \deg{P^{(2)}_{2m,j}}\leqslant 2m, by the relations
In a similar way, for the system f, f^2, f^3 we define the Hermite–Padé polynomials P^{(3)}_{3\ell,0}\not\equiv0, P^{(3)}_{3\ell,1}, P^{(3)}_{3\ell,2}, and P^{(3)}_{3\ell,3}, \deg{P^{(3)}_{3\ell,j}}\leqslant 3\ell, by the relations
Let f\in\mathcal{H}(\infty) be explicitly defined by f(z):=[(A-1/{\varphi(z)})\times(B-1/{\varphi(z)})]^{-1/2}, z\notin E:=[-1,1], where 1<A<B, \varphi(z)=z+(z^2-1)^{1/2} and the branch of (\,\cdot\,)^{1/2} is selected so that \varphi(z)\sim 2z as z\to\infty. We denote the class of these functions by \mathcal Z(E). The function f is algebraic of order four, with four branch points \pm1, a, and b, where a=(A+1/A)/2 and b=(B+1/B)/2, 1<a<b. All branch points of f are of the second order. It was shown in [8] that f, f^2, f^3 is a Nikishin system:
where \operatorname{supp}\sigma=\operatorname{supp}{s}_1= \operatorname{supp}{s}_2=E, \operatorname{supp}{\sigma}_2=F:=[a,b], s_1:=\langle\sigma,\sigma_2\rangle, and s_2:=\langle\sigma,\sigma_2,\sigma\rangle.
Let M_1(F) denote the class of unit (Borel) measures with support on F. Let g_E(t,z) be the Green’s function for the domain D:=\widehat{\mathbb{C}}\setminus{E} with singularity at t=z, and let V^\mu(z) be the logarithmic potential and G^\mu_E(z) be the Green’s potential of \mu\in M_1(F):
For each \theta\in(0,\infty) there exists (see [6]) a unique measure \lambda=\lambda(\theta)\in M_1(F) such that the equilibrium condition \theta V^\lambda(x)+G^\lambda_E(x)+ \theta g_E(x,\infty)\equiv\operatorname{const}, x\in F, is satisfied. We have the following result.
Theorem 1. Let f\in\mathcal Z(E). Then as N\to\infty, locally uniformly in D,
We can show that G^{\lambda(1)}_E(z)>2G^{\lambda(3)}_E(z)/3, z\in D. From (9), (7), and (8) we obtain
\begin{equation}
\delta_3(z)<\delta_2(z)<\delta_1(z)<1,\qquad z\in D.
\end{equation}
\tag{10}
Bibliography
1.
A. I. Aptekarev, Dokl. Math., 78:2 (2008), 717–719
2.
P. Henrici, SIAM J. Numer. Anal., 3:1 (1966), 67–78
3.
E. O. Dobrolyubov, N. R. Ikonomov, L. A. Knizhnerman, and S. P. Suetin, Rational Hermite–Padé approximants vs Padé approximants. Numerical results, 2023, 53 pp., arXiv: 2306.07063
4.
A. V. Komlov, Sb. Math., 212:12 (2021), 1694–1729
5.
E. A. Rakhmanov, Russian Math. Surveys, 73:3 (2018), 457–518
6.
E. A. Rakhmanov and S. P. Suetin, Sb. Math., 204:9 (2013), 1347–1390
7.
H. Stahl, J. Approx. Theory, 91:2 (1997), 139–204
8.
S. P. Suetin, Math. Notes, 104:6 (2018), 905–914
9.
S. P. Suetin, Russian Math. Surveys, 77:6 (2022), 1149–1151
10.
E. J. Weniger, Comput. Phys. Rep., 10:5-6 (1989), 189–371
Citation:
S. P. Suetin, “Convergence of Hermite–Padé rational approximations”, Russian Math. Surveys, 78:5 (2023), 967–969
This publication is cited in the following 4 articles:
S. P. Suetin, “O skalyarnykh podkhodakh k izucheniyu predelnogo raspredeleniya nulei mnogochlenov Ermita–Pade dlya sistemy Nikishina”, UMN, 80:1(481) (2025), 85–152
S. P. Suetin, “Maximum principle and asymptotic properties of Hermite–Padé polynomials”, Russian Math. Surveys, 79:3 (2024), 547–549
N. R. Ikonomov, S. P. Suetin, “On some potential-theoretic problems related to the asymptotics of Hermite–Padé polynomials”, Sb. Math., 215:8 (2024), 1053–1064
V. G. Lysov, “Distribution of zeros of polynomials of multiple discrete orthogonality in the Angelesco case”, Russian Math. Surveys, 79:6 (2024), 1101–1103