Пришвидшення ряду

У математиці пришвидшення ряду є одним із Шаблон:Нп для покращення швидкості збіжності ряду. Методи пришвидшення рядів часто застосовуються в чисельному аналізі, де вони використовуються для підвищення швидкості чисельного інтегрування. Методи пришвидшення ряду також можуть бути використані, наприклад, для отримання різноманітних тотожностей для спеціальних функцій. Наприклад, перетворення Ейлера, застосоване до гіпергеометричного ряду, дає деякі з класичних, добре відомих тотожностей для гіпергеометричних рядів.

Означення

Нехай задано послідовність

$S={\left\{s_n\right\}}_{n\in \mathbb{N}},$

що має границю

$\underset{n\to \mathrm{\infty }}{\lim }{s}_n=\ell .$

Пришвидшенним рядом є друга послідовність

$S^{\prime }={\left\{s_n^{\prime }\right\}}_{n\in \mathbb{N}},$

яка збігається швидше до $\ell$, ніж початкова послідовність, в сенсі, що

$\underset{n\to \mathrm{\infty }}{\lim }\frac{s_n^{\prime }-\ell }{s_n-\ell }=0.$

Якщо початкова послідовність є розбіжною, то Шаблон:Нп діє як метод екстраполяції на антиграницю $\ell$.

Відображення від початкового до перетвореного ряду можуть бути лінійними (як визначено в статті про Шаблон:Нп) або нелінійними. У загальному випадку, нелінійні перетворення послідовності мають тенденцію бути більш потужними.

Огляд

Двома класичними методами пришвидшення ряду є перетворення рядів Ейлера і Шаблон:Нп. У XX столітті було розроблено низку більш швидших та спеціальних методів, включаючи Шаблон:Нп, яку запропонував Льюїс Фрай Річардсон на початку XX століття, але який був відомий і використовувався Шаблон:Нп в 1722 році; Шаблон:Нп, запропонований Шаблон:Нп у 1926 році, який також був відомий і використовувався Такакадзу Секі у XVIII столітті; епсилон метод, запропонований Шаблон:Нп у 1956 році; $u$-перетворення Левіна; метод Вілфа—Цейлбергера—Ехада; Шаблон:Нп.

Для знакопереміжних рядів існує декілька потужних методів, описаних Коеном зі співавторами , які забезпечують швидкість збіжності від $5,828^{-n}$ до $17,93^{-n}$ для підсумовування $n$ членів.

Перетворення Ейлера

Основним прикладом Шаблон:Нп, що забезпечує покращену збіжність, є перетворення Ейлера. Воно застосовується до знакопереміжних рядів і має вигляд

$\sum _{n=0}^{\mathrm{\infty }}(-1{)}^n{a}_n=\sum _{n=0}^{\mathrm{\infty }}(-1{)}^n\frac{({\mathrm{\Delta }}^{n}a{)}_0}{2^{n+1}},$

де $\mathrm{\Delta }$ — це оператор правих різниць:

$({\mathrm{\Delta }}^{n}a{)}_0=\sum _{k=0}^n(-1{)}^k(\genfrac{}{}{0ex}{}{n}{k})a_{n-k}.$

Якщо початковий ряд у лівій частині лише повільно збігається, то праві різниці будуть мати тенденцію до швидкого зменшення, а додаткова степінь двійки іще більше збільшує швидкість збіжності правої частини.

Особливо ефективною чисельною реалізацією перетворення Ейлера є перетворення Шаблон:Нп.

Конформні відображення

Ряд

$S=\sum _{n=0}^{\mathrm{\infty }}{a}_n$

можна записати як $f(1)$, де функція $f$ визначається як

$f(z)=\sum _{n=0}^{\mathrm{\infty }}{a}_n{z}^n.$

Функція $f(z)$ може мати особливості в комплексній площині точки розгалуження, полюси чи суттєві особливості), які обмежують радіус збіжності ряду. Якщо точка $z=1$ знаходиться поблизу або на межі круга збіжності, то ряд для $S$ буде збігатися дуже повільно. Можна покращити збіжність ряду за допомогою конформного відображення, яке переміщує сингулярності таким чином, що точка, яка відображається в $z=1$, опиняється глибше в новому крузі збіжності.

Конформне перетворення $z=\mathrm{\Phi }(w)$ потрібно вибрати таким, що $\mathrm{\Phi }(0)=0$, і зазвичай вибирають функцію, яка має скінченну похідну при $w=0$. Без втрати загальності можна вважати, що $\mathrm{\Phi }(1)=1$, оскільки завжди можна перемасштабувати $w$ і перевизначити $\mathrm{\Phi }$. Далі розглянемо функцію

$g(w)=f(\mathrm{\Phi }(w)).$

Оскільки $\mathrm{\Phi }(1)=1$, то $f(1)=g(1)$. Можемо отримати розклад в ряд функції $g(w)$, поклавши $z=\mathrm{\Phi }(w)$ у розкладі в ряд функції $f(z)$, оскільки $\mathrm{\Phi }(0)=0$. Перші $n$ членів розкладу ряду для функції $f(z)$ дадуть перші $n$ членів розкладу ряду для функції $g(w)$, якщо ${\mathrm{\Phi }}^{\prime }(0)\ne 0$. Таким чином, поклавши $w=1$ у цьому розкладі в ряд, отримаємо ряд, який у випадку його збіжності, збігається до того ж значення, що й початковий ряд.

Нелінійні перетворення послідовності

Прикладами таких нелінійних перетворень послідовностей є апроксимації Паде, Шаблон:Нп та перетворення послідовностей типу Левіна.

Особливо нелінійні перетворення послідовностей часто забезпечують потужні числові методи для підсумовування розбіжних рядів або асимптотичних рядів, які виникають, наприклад, у теорії збурень, і можуть бути використані як високоефективні методи екстраполяції.

Метод Ейткена

Основна стаття: Шаблон:Нп

Просте нелінійне перетворення послідовності (екстраполяція Ейткена або дельта-квадратичний метод):

$𝔸:S\to S^{\prime }=𝔸(S)={(s_n^{\prime })}_{n\in \mathbb{N}},$

яке визначається наступним чином:

$s_n^{\prime }=s_{n+2}-\frac{(s_{n+2}-s_{n+1}{)}^2}{s_{n+2}-2s_{n+1}+s_n}.$

Це перетворення дуже часто використовується для покращення швидкості збіжності послідовності, що повільно збігається; евристично, таке перетворення усуває більшу частину абсолютної похибки.

Див. також

• Шаблон:Нп
• Шаблон:Нп
• Шаблон:Нп

Література

• C. Brezinski and M. Redivo Zaglia, Extrapolation Methods. Theory and Practice, North-Holland, 1991.
• G. A. Baker Jr. and P. Graves-Morris, Padé Approximants, Cambridge U.P., 1996.
• Шаблон:Mathworld
• Herbert H. H. Homeier: Scalar Levin-Type Sequence Transformations, Journal of Computational and Applied Mathematics, vol. 122, no. 1–2, p 81 (2000). Шаблон:Cite journal, Шаблон:Arxiv.
• Brezinski Claude and Redivo-Zaglia Michela : "The genesis and early developments of Aitken's process, Shanks transformation, the $ϵ$-algorithm, and related fixed point methods", Numerical Algorithms, Vol.80, No.1, (2019), pp.11-133.
• Delahaye J. P. : "Sequence Transformations", Springer-Verlag, Berlin, ISBN 978-3540152835 (1988).
• Sidi Avram : "Vector Extrapolation Methods with Applications", SIAM, ISBN 978-1-61197-495-9 (2017).
• Brezinski Claude, Redivo-Zaglia Michela and Saad Yousef : "Shanks Sequence Transformations and Anderson Acceleration", SIAM Review, Vol.60, No.3 (2018), pp.646–669. doi:10.1137/17M1120725 .
• Brezinski Claude : "Reminiscences of Шаблон:Iw", Numerical Algorithms, Vol.80(2019), pp.5-10.
• Brezinski Claude and Redivo-Zaglia Michela : "Extrapolation and Rational Approximation", Springer, ISBN 978-3-030-58417-7 (2020).

Зовнішні лінки

• Convergence acceleration of series
• GNU Scientific Library, Series Acceleration
• Digital Library of Mathematical Functions