Математика / Численный поиск корней

Процесс Эйткена Δ²: ускорение сходящейся последовательности

Процесс Эйткена Δ² строит оценку предела по трём последовательным членам. Он точно устраняет одну геометрическую составляющую ошибки и может ускорить линейную сходимость.

Опубликовано: Обновлено:

Формула

$$\widehat s_n=s_n-\frac{(s_{n+1}-s_n)^2}{s_{n+2}-2s_{n+1}+s_n}$$
Схема Удаление геометрической ошибки
1410s14121110,50123nУскоренная оценка: 10

Синяя линия: sₙ=10+4·2⁻ⁿ. Зелёная линия: оценка 10 после Δ². Значения первых четырёх членов подписаны.

Обозначения

$sₙ,sₙ₊₁,sₙ₊₂$
три соседних члена исходной последовательности, одинаковые единицы
$ŝₙ$
ускоренная оценка предела, единицы последовательности
$Δsₙ$
первая разность sₙ₊₁−sₙ, единицы последовательности
$Δ²sₙ$
вторая разность sₙ₊₂−2sₙ₊₁+sₙ, единицы последовательности

Условия применения

Вторая разность в знаменателе должна отличаться от нуля и быть достаточно большой по сравнению с погрешностью исходных чисел. Для гарантии ускорения используют предпосылки о линейной сходимости к конечному пределу и ненулевом предельном отношении ошибок по модулю меньше единицы.

Ограничения

Преобразование не является универсальным доказательством сходимости. Малый знаменатель усиливает ошибки округления и исходных данных. Несколько сопоставимых колебательных составляющих могут не устраняться одним шагом. При почти постоянных членах лучше остановиться по критерию точности, а не делить две малые разности.

Подробное объяснение

Пусть последовательность имеет вид sₙ=L+Cλⁿ, где L — искомый предел и |λ|<1. Тогда первая разность равна Cλⁿ(λ−1), а вторая — Cλⁿ(λ−1)². Деление квадрата первой разности на вторую возвращает Cλⁿ. После вычитания остаётся ровно L.

Именно эта модель объясняет пользу трёх соседних значений. Метод оценивает ведущую геометрическую ошибку и удаляет её, не зная заранее ни коэффициента C, ни множителя λ. Точность для одной геометрической составляющей не переносится автоматически на произвольные последовательности.

Если исходная ошибка содержит меньшую добавку, преобразование часто устраняет главный член и оставляет эту добавку. Например, сумма двух геометрических ошибок обычно не исчезнет целиком. Проверка нескольких соседних троек показывает, насколько устойчивой стала оценка предела.

Численно опасно вычитать почти равные значения sₙ₊₂ и 2sₙ₊₁−sₙ. Незначительное возмущение входов способно сильно изменить частное. Поэтому знаменатель сравнивают с масштабом данных и их точностью, а результат проверяют по независимому свойству задачи.

Эйткен преобразует уже имеющуюся последовательность. В методе Стеффенсена, напротив, после каждой ускоренной оценки заново выполняют два шага исходной итерации. Эти процедуры используют одну алгебраическую формулу, но создают разные последовательности и требуют разных критериев остановки.

Как пользоваться формулой

  1. Возьмите три соседних приближения одного вычислительного процесса, сохранив достаточно цифр для устойчивого вычисления их разностей.
  2. Вычислите Δsₙ и Δ²sₙ; сравните вторую разность с масштабом входов и погрешностью, чтобы исключить опасное деление.
  3. Вычтите (Δsₙ)²/Δ²sₙ из первого члена тройки, не изменяя порядок индексов и не округляя промежуточное частное.
  4. Повторите преобразование для следующей тройки и проверьте согласие оценок с исходной задачей или известной границей ошибки.

Историческая справка

Имя Александра Эйткена закрепилось за Δ²-преобразованием последовательностей в численном анализе XX века. Современный обзор Acceleration methods for fixed-point iterations относит работу Эйткена к ранним методам экстраполяции. В справочнике NIST DLMF процесс выделен в отдельный подраздел среди методов ускорения сходимости.

Контекст этого направления — получение более точного результата из уже вычисленных приближений. Вместо простого продолжения медленной последовательности исследуют структуру её остатка. Для Δ² ключевой стала геометрическая модель ошибки; последующие преобразования обобщают идею на несколько составляющих. Современная запись через конечные разности делает эту модель явной и одновременно показывает источник численной неустойчивости.

Историческая линия формулы

Название процесса указывает на Эйткена и принятую в численном анализе запись через вторую конечную разность. Приведённый вывод является современной проверкой для геометрической ошибки; он не означает, что ускорение гарантируется для любой последовательности или любых округлённых данных.

Пример

Последовательность приближений задана правилом sₙ=10+4·2⁻ⁿ. Её первые члены: s₀=14, s₁=12, s₂=11; требуется оценить предел только по этой тройке. Все величины безразмерны. Первая разность Δs₀=12−14=−2. Вторая разность Δ²s₀=11−2·12+14=1. Формула даёт ŝ₀=14−(−2)²/1=10. Для следующей тройки 12;11;10,5 получаем Δs₁=−1, Δ²s₁=0,5 и ŝ₁=12−1/0,5=10. Независимая проверка следует из определения: 2⁻ⁿ стремится к нулю, поэтому L=10. Ответ совпал точно, поскольку ошибка является одной геометрической прогрессией. Совпадение двух оценок здесь подтверждает вычисление, но для произвольных данных само по себе ещё не доказывает точность.

Частая ошибка

Не используйте несоседние значения с разными шагами без отдельной модели ошибки. Не теряйте квадрат первой разности: её знак не сохраняется после возведения в квадрат. Не заменяйте вторую разность на разность двух соседних членов. Не считайте нулевой знаменатель доказательством отсутствия предела: постоянная последовательность уже имеет предел.

Практика

Задачи с решением

Знакочередующаяся ошибка

Условие. Даны члены 5;0,5;2,75 последовательности с геометрической ошибкой. Примените Δ².

Решение. Δs=−4,5; Δ²s=2,75−1+5=6,75. Ускоренная оценка 5−20,25/6,75=2.

Ответ. 2.

Вырожденная тройка

Условие. Можно ли применить формулу к значениям 1;2;3? Найдите знаменатель.

Решение. Δ²s=3−2·2+1=0. Формула требует деления на нуль и для этой тройки не определена.

Ответ. Нет, знаменатель равен нулю.

Дополнительные источники

  • https://dlmf.nist.gov/3.9
  • https://www.maths.dur.ac.uk/users/anthony.yeates/NOTES/na17.pdf
  • https://www.cambridge.org/core/services/aop-cambridge-core/content/view/2FABE7BD50568960C8CD86A6F4D77D72/S0962492924000096a.pdf/acceleration-methods-for-fixed-point-iterations.pdf

Связанные формулы

Математика

Метод Стеффенсена: ускорение простой итерации

$y=g(x_n),\quad z=g(y),\quad x_{n+1}=x_n-\frac{(y-x_n)^2}{z-2y+x_n}$

Метод Стеффенсена ускоряет поиск неподвижной точки без вычисления производной. Он выполняет два шага простой итерации и применяет к полученной тройке преобразование Эйткена.

Математика

Метод простой итерации для поиска корня уравнения

$x_{k+1}=g(x_k),\quad x=g(x)$

Метод простой итерации заменяет уравнение f(x)=0 равносильным уравнением x=g(x) и строит последовательность приближений. Условие сжатия отображения даёт проверяемую гарантию сходимости к единственной неподвижной точке.

Математика

Метод секущих: поиск корня без вычисления производной

$x_{k+1}=x_k-f(x_k)\frac{x_k-x_{k-1}}{f(x_k)-f(x_{k-1})}$

Метод секущих строит прямую по двум последним точкам графика и использует её пересечение с осью для нового приближения корня. Производная заменяется разностным наклоном.