Пусть функция f задана степенным рядом в бесконечно удаленной точке z=∞:
f(z)=∞∑k=0ckzk.
Если известны N первых коэффициентов c0,…,cN−1 ряда (1), то на их основе можно определить несколько конструктивных аппроксимаций исходной функции f (см. [10]; “конструктивность” понимается здесь в том смысле, как это сформулировано в [2; § 2]). Естественно поставить вопрос об оптимальности таких аппроксимаций. В дальнейшем предполагается, что N=2n+1=3m+1=4ℓ+1, где n,m,ℓ∈N. Это условие не связано с существом дела, а носит чисто технический характер (см. (2)–(4)).
Наиболее известные, востребованные и исследованные конструктивные аппроксимации степенного ряда – это аппроксимации Паде (АП) [n/n]f=Pn/Qn, где полиномы Pn,Qn≢0, degPn,degQn⩽, определяются из соотношения
Сходимость АП в классе многозначных аналитических функций с конечным числом особых точек вытекает из теории Шталя 1985–1986 гг. [7]. Вместе с тем существуют и другие способы конструктивной аппроксимации степенного ряда (1) по заданному набору из N его первых коэффициентов, которые не так известны, как АП. Здесь мы обсудим рациональные аппроксимации, построенные по полиномам Эрмита–Паде 2-го типа для систем f, f^2 и f, f^2, f^3. Оказывается, с точки зрения оптимального использования N коэффициентов ряда (1) рациональные аппроксимации Эрмита–Паде имеют определенные преимущества перед АП (см. (10), а также [4], [3]). Сходимость таких аппроксимаций исследована пока только для частных случаев (см. [1], [6], [5], [4]). В работе [3] представлены некоторые численные результаты, связанные с уравнением Ван дер Поля. Эти результаты показывают, что использование полиномов Эрмита–Паде для систем f, f^2 и f, f^2, f^3 существенно расширяет возможности численного анализа свойств функции f на основе коэффициентов ряда (1). В частности, таким образом можно распознать квадратичные точки ветвления функции, заданной рядом (см. [9]). Цель настоящей работы – представить теоретические результаты, показывающие, что рациональные аппроксимации Эрмита–Паде оказываются более оптимальными, чем АП.
Для системы f, f^2 определим полиномы Эрмита–Паде P^{(2)}_{2m,0}\not\equiv 0, P^{(2)}_{2m,1}, P^{(2)}_{2m,2}, \deg{P^{(2)}_{2m,j}}\leqslant 2m, из соотношений
Пусть функция f\in\mathcal{H}(\infty) задана явным представлением: f(z):=[(A-1/{\varphi(z)})\times(B-1/{\varphi(z)})]^{-1/2}, z\notin E:=[-1,1], где 1<A<B, \varphi(z)=z+(z^2-1)^{1/2} и выбрана такая ветвь функции (\,\cdot\,)^{1/2}, что \varphi(z)\sim 2z при z\to\infty. Класс таких функций обозначим через \mathcal Z(E). Функция f – алгебраическая функция 4-го порядка с четырьмя точками ветвления \pm1, a, b, где a=(A+1/A)/2, b=(B+1/B)/2, 1<a<b. Все точки ветвления f – 2-го порядка. В [8] показано, что f, f^2, f^3 – система Никишина:
где \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, s_2:=\langle\sigma,\sigma_2,\sigma\rangle.
Обозначим через M_1(F) класс всех единичных (борелевских) мер с носителем на F. Пусть g_E(t,z) – функция Грина для области D:=\widehat{\mathbb{C}}\setminus{E} с особенностью в точке t=z, и пусть V^\mu(z) – логарифмический, а G^\mu_E(z) – гринов потенциалы меры \mu\in M_1(F):
Для любого \theta\in(0,\infty) существует [6] единственная мера \lambda=\lambda(\theta)\in M_1(F) такая, что выполняется условие равновесия \theta V^\lambda(x)+G^\lambda_E(x)+ \theta g_E(x,\infty)\equiv\operatorname{const}, x\in F. Справедливо следующее утверждение.
Теорема 1. Пусть f\in\mathcal Z(E). Тогда при N\to\infty равномерно внутри D
Можно показать, что G^{\lambda(1)}_E(z)>2G^{\lambda(3)}_E(z)/3, z\in D. Из (9), (7) и (8) получаем, что
\begin{equation}
\delta_3(z)<\delta_2(z)<\delta_1(z)<1,\qquad z\in D.
\end{equation}
\tag{10}
Список литературы
1.
А. И. Аптекарев, Докл. РАН, 422:4 (2008), 443–445
2.
P. Henrici, SIAM J. Numer. Anal., 3:1 (1966), 67–78
3.
N. R. Ikonomov, L. A. Knizhnerman, S. P. Suetin, Rational Hermite–Padé approximants vs Padé approximants. I. Numerical results, 2023, 36 pp., arXiv: 2306.07063
4.
А. В. Комлов, Матем. сб., 212:12 (2021), 40–76
5.
Е. А. Рахманов, УМН, 73:3(441) (2018), 89–156
6.
Е. А. Рахманов, С. П. Суетин, Матем. сб., 204:9 (2013), 115–160
7.
H. Stahl, J. Approx. Theory, 91:2 (1997), 139–204
8.
С. П. Суетин, Матем. заметки, 104:6 (2018), 918–929
9.
С. П. Суетин, УМН, 77:6(468) (2022), 203–204
10.
E. J. Weniger, Comput. Phys. Rep., 10:5-6 (1989), 189–371
Образец цитирования:
С. П. Суетин, “О сходимости рациональных аппроксимаций Эрмита–Паде”, УМН, 78:5(473) (2023), 185–186; Russian Math. Surveys, 78:5 (2023), 967–969
С. П. Суетин, “О скалярных подходах к изучению предельного распределения нулей многочленов Эрмита–Паде для системы Никишина”, УМН, 80:1(481) (2025), 85–152
С. П. Суетин, “Принцип максимума и асимптотические свойства многочленов Эрмита–Паде”, УМН, 79:3(477) (2024), 181–182; S. P. Suetin, “Maximum principle and asymptotic properties of Hermite–Padé polynomials”, Russian Math. Surveys, 79:3 (2024), 547–549
Н. Р. Икономов, С. П. Суетин, “О теоретико-потенциальных задачах, связанных с асимптотикой многочленов Эрмита–Паде”, Матем. сб., 215:8 (2024), 52–65; 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
В. Г. Лысов, “Распределение нулей многочленов совместной дискретной ортогональности в случае Анжелеско”, УМН, 79:6(480) (2024), 165–166; V. G. Lysov, “Distribution of zeros of polynomials of multiple discrete orthogonality in the Angelesco case”, Russian Math. Surveys, 79:6 (2024), 1101–1103