Обратная квадратичная интерполяция - Википедия - Inverse quadratic interpolation

В числовой анализ, обратная квадратичная интерполяция это алгоритм поиска корней, что означает, что это алгоритм решения уравнений вида ж(Икс) = 0. Идея состоит в том, чтобы использовать квадратичная интерполяция приблизить обратный из ж. Этот алгоритм редко используется сам по себе, но он важен, потому что он является частью популярных Метод Брента.

Метод

Алгоритм обратной квадратичной интерполяции определяется отношение повторения

{ Displaystyle x_ {n + 1} = { frac {f_ {n-1} f_ {n}} {(f_ {n-2} -f_ {n-1}) (f_ {n-2} -f_ {n})}} x_ {n-2} + { frac {f_ {n-2} f_ {n}} {(f_ {n-1} -f_ {n-2}) (f_ {n-1 } -f_ {n})}} x_ {n-1}}

{ displaystyle {} + { frac {f_ {n-2} f_ {n-1}} {(f_ {n} -f_ {n-2}) (f_ {n} -f_ {n-1}) }} x_ {n},}

где ж_k = ж(Икс_k). Как видно из рекуррентного соотношения, для этого метода требуются три начальных значения: Икс₀, Икс₁ и Икс₂.

Объяснение метода

Мы используем три предыдущих итерации, Икс_п−2, Икс_п−1 и Икс_п, с их значениями функций, ж_п−2, ж_п−1 и ж_п. Применяя Формула интерполяции Лагранжа делать квадратичную интерполяцию на обратном ж дает

{ displaystyle f ^ {- 1} (y) = { frac {(y-f_ {n-1}) (y-f_ {n})} {(f_ {n-2} -f_ {n-1) }) (f_ {n-2} -f_ {n})}} x_ {n-2} + { frac {(y-f_ {n-2}) (y-f_ {n})} {(f_ {n-1} -f_ {n-2}) (f_ {n-1} -f_ {n})}} x_ {n-1}}

{ displaystyle {} + { frac {(y-f_ {n-2}) (y-f_ {n-1})} {(f_ {n} -f_ {n-2}) (f_ {n} -f_ {n-1})}} x_ {n}.}

Ищем корень ж, поэтому подставляем у = ж(Икс) = 0 в приведенном выше уравнении, и это приводит к приведенной выше формуле рекурсии.

Поведение

Асимптотическое поведение очень хорошее: обычно итерация Икс_п быстро сходятся к корню, как только они приблизятся. Однако производительность часто бывает довольно низкой, если начальные значения не близки к фактическому корню. Например, если случайно два из значений функции ж_п−2, ж_п−1 и ж_п совпадают, алгоритм полностью не работает. Таким образом, обратная квадратичная интерполяция редко используется в качестве автономного алгоритма.

Порядок этой сходимости составляет приблизительно 1,84, что доказывает Секущий метод анализ.

Сравнение с другими методами поиска корней

Как отмечалось во введении, обратная квадратичная интерполяция используется в Метод Брента.

Обратная квадратичная интерполяция также тесно связана с некоторыми другими методами поиска корней. линейная интерполяция вместо квадратичной интерполяции дает секущий метод. Интерполяция ж вместо инверсии ж дает Метод Мюллера.

Смотрите также

Последовательная параболическая интерполяция - связанный метод, использующий параболы для поиска экстремумов, а не корней.