Когда кубическая сплайн-интерполяция лучше, чем интерполяционный полином?

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

Однако недавно я нашел много примеров сигналов с ограниченной полосой пропускания, в которых интерполяционный полином высокого порядка дает меньшую погрешность аппроксимации, чем интерполяция кубического сплайна. Обычно интерполяционный полином более точен во всем интервале интерполяции, когда частота дискретизации достаточно высока. Кажется, это имеет место, когда выборки равномерно распределены с частотой выборки, по крайней мере, в 3 раза большей, чем частота Найквиста сигнала. Кроме того, преимущество по сравнению с кубической сплайн-интерполяцией улучшается с увеличением (частоты дискретизации) / (частоты Найквиста).

В качестве примера я сравниваю интерполяцию кубического сплайна с интерполяционным полиномом для синусоидальной волны с частотой Найквиста 2 Гц и частотой дискретизации 6,5 Гц. Между точками выборки интерполирующий полином выглядит точно так же, как и фактический сигнал.

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

Обсуждаемый феномен - это феномен Рунге .

Максимальное абсолютное значение й производной от равно . Для функции Рунге максимальное абсолютное значениегрешить ( ω т ) ω $n$$\sin(\omega t)$$\omega^n$ $\frac{1}{25t^2+1}$$n$ й (четной) производной равно гдеобозначает факториал. Это намного быстрее роста. Только если производные растут слишком быстро при увеличении , возможно, что ошибка интерполяции расходится при увеличении порядка интерполяции. Экспонента по еще не слишком быстра. Взгляните на: Джеймс Ф. Эпперсон, На примере Рунге , Американский математический месяц , том. 94, 1987, с. 329-341.$5^nn!,$$n!$$n$$n$

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

Если оба метода сходятся, то (не кусочная) полиномиальная интерполяция имеет преимущество более высокой степени полинома, если используется много выборок, и может обеспечить лучшее приближение, как вы видели в своем примере с синусоидой. Вас также может заинтересовать Л. Н. Трефетен, Два результата о полиномиальной интерполяции в равномерно распределенных точках , Журнал теории приближений, том 65, выпуск 3, июнь 1991 г., стр. 247-260. Quote:

[...] в ограниченной по полосе интерполяции сложных показательных функций ошибка уменьшается до как если и только если достаточно мала, чтобы обеспечить как минимум шесть точек на длину волны.$e^{i\alpha x}\, (\alpha \in \mathbb{R}),$$0$$n \to \infty$$\alpha$

У вас есть 6,5 образцов на длину волны.