Феномен Рунге — проблема, що виникає в обчислювальній математиці при використанні поліноміальної інтерполяції на рівновіддалених вузлах за допомогою поліномів високих порядків (степенів). Була описана Карлом Рунге при вивченні поводження похибок при використанні поліноміальної інтерполяції для апроксимації функцій.
Проблема
Розглянемо функцію:
Спробуємо інтерполювати цю функцію в рівновіддалених точках xi між −1 і 1 таких що:
де
за допомогою полінома порядок якого . Рунге виявив, що отримана інтерполяція коливається і має досить велику похибку поблизу кінців інтервалу, тобто поблизу точок −1 і 1. Можна довести, що верхня межа похибка інтерполяції прямує до нескінченності, коли ступінь полінома зростає:
Цей приклад показує, що поліноміальна інтерполяція високого порядку на рівновіддалених вузлах може бути небезпечною.
Причина
Похибка наближення заданої функції f(x) інтерполяційним поліномом порядку визначається так:
для деякого у (−1, 1). Таким чином,
Позначимо вузлову функцію:
і нехай буде максимумом функції
Тоді, можна довести, що коли використовувати рівновіддалені вузли,[1] тоді:
де це розмір кроку.
Більше того, припустімо, що n-та похідна обмежена, тобто
- .
З цього,
Але для випадку функції Рунге, заданої вище: перші дві похідні
Величина похідних вищого порядку наведеної функції Рунге зростає. Отже, верхня межа похибка (у проміжках між точками інтерполяції) зі збільшенням степеня інтерполюючих поліномів стає ще більшою.
Проте факт збільшення верхньої межі похибки до нескінченності не обов'язково означає, що сама похибка також зростає з n.[джерело?]
Розв'язання проблеми
Феномен Рунге демонструє, що коли обирати рівновіддалені вузли інтерполяції, то високий порядок апроксимуючого полінома не завжди підходить.
Однак, апроксимаційна теорема Вейєрштраса стверджує, що для неперервної на відрізку функції існує послідовність поліноміальних апроксимацій, для яких похибка прямує до нуля. Розбіжність можна мінімізувати, якщо замість рівновіддалених вузлів обирати вузли інтерполяції в коренях поліномів Чебишова першого роду. Такий підхід гарантує, що максимальна похибка зменшуватиметься з підвищенням порядку полінома[2].
Проблему також можна вирішити шляхом застосуванням сплайнів, зокрема, поліноміальних, тобто, замість підвищення степеня поліномів збільшити кількість поліноміальних частин.
Примітки
- ↑ Heath, Michael (2000). Scientific Computing. McGraw-Hill. с. 324. ISBN 0072399104.
- ↑ Lloyd N. Trefethen. Approximation Theory and Approximation Practice : With 163 figures and 210 exercises : [англ.]. — 2013.