Методические указания: Интерполяция
7 ИНТЕРПОЛЯЦИЯ, ЭКСТРАПОЛЯЦИЯ
7.1 Постановка задачи
Предположим, что задано различных точек плоскости:
(7.1)
Требуется найти функцию
, значения которой при данных значениях абсциссы
в точности равны соответствующим ординатам заданных точек:
Т.е. нужно найти линию, описываемую уравнением
, проходящую через
данную точку (рис.7.1).
Рис.7.1
Заметим, что здесь приходится различать два случая:
1) интерполяцию (от лат. interpolar Ч подновлять) Ч
восстановление промежуточных значений функции внутри интервала
по ряду известных ее значений;
2) экстраполяцию (лат. приставка extra означает лвне) Ч когда
не вошедшее в исследование значение
лежит вне интервала
.
Очевидно, интерполяция более надежна, чем экстраполяция.
Вообще говоря, существует бесконечное число линий, проходящих через
заданную точку. Потребуем, чтобы искомая линия была простейшей, т.е. значения
функции, задающие эту линию, должны находиться при помощи простейших операций
(сложения, умножения). Этому требованию отвечают многочлены (полиномы), т.е.
выражения вида:
(7.2)
Зная численные значения коэффициентов
многочлена, мы можем найти его ординату при любом значении переменной
. Наконец, из двух многочленов условимся считать простейшим тот, степень
которого ниже.
Итак, приходим к задаче о полиномиальной интерполяции: пусть даны
различных чисел и
соответствующих им чисел
, требуется найти многочлен
наименьшей возможной степени, удовлетворяющий
условиям:
7.2 Интерполяционный многочлен Лагранжа для произвольных узлов
Для решения предложенной задачи зафиксируем одну ординату
, а остальные будем считать равными нулю (рис.7.2), т.е. заданным значениям
абсцисс ставятся
в соответствие значения ординат
Из свойств многочленов следует, что многочлен, обращающийся в нуль в
разных точках, т.е. имеющий
различных корней, должен
делиться на каждую из разностей:
Рис.7.2
а следовательно, и на произведение этих разностей, т.е. его степень не может
быть ниже . В таком
случае многочлен должен иметь вид
(7.3)
Из условия определим значение const
,
таким образом находим
(7.4)
В полученном выражении никакого особого преимущества
не имеет, мы можем приписать эту особую роль любому
, т.е. если абсциссам
поставить в соответствие значения
, указанные в любой из следующих строк:
то выражение для многочлена, принимающего при соответствующих значениях
абсцисс численные значения, выписанные в одной из строк, будет аналогично
рассмотренному, т.е.
(7.5)
Общее решение является суперпозицией (суммой) частных решений (7.5)
(7.6)
Это и есть интерполяционный многочлен Лагранжа. По наборам исходных пар (7.1)
формула (7.6) позволяет достаточно просто составить лвнешний вид многочлена.
Используя обозначение
, (7.7)
формуле Лагранжа можно придать более сжатый вид. Продифференцируем по
;
при имеем:
. (7.8)
Формула Лагранжа с учетом (7.7) и (7.8) примет вид:
или (7.9)
В рассмотренном случае предполагалось, что точки
расположены на отрезке
произвольно. Рассмотрим формулу Лагранжа, для равноотстоящих значений абсцисс.
7.3 Интерполяционный многочлен Лагранжа для равностоящих узлов
Пусть на отрезке
задана система равноотстоящих узлов
которыми отрезок делится на
равных частей
где
В этом случае интерполяционный многочлен Лагранжа строится на равноотстоящих
узлах и имеет более удобный вид.
Обозначим , где . Отсюда:
..................................................
Т.е. в общем случае:
(7.10)
Используя (7.10) и принятое обозначение получим:
(7.11)
Учитывая, что найдем:
(7.12)
Заметим, что в (7.12) ровно
строк (-я строка
отсутствует); причем численные значения
первых строк положительны, а остальные Ч отрицательны. Используя (7.12),
получим:
т.е.
(7.13)
С учетом (7.11) и (7.13) формула Лагранжа для равноотстоящих
узлов примет вид:
(7.14)
7.4 Интерполяционный многочлен Ньютона для равноотстоящих узлов
На практике часто встречается случай, когда интерполяционная функция подбирается
для таблиц с равноотстоящими значениями аргумента
Рассмотрим метод построения интерполирующей функции, основанный на вычислении
конечных разностей.
7.4.1 Конечные разности
Назовем конечными разностями разности между значениями функции в
соседних узлах интерполяции:
где Полученные
конечные разности будем называть разностями первого порядка. Из разностей
первого порядка получим разности второго порядка:
где
Повторяя процедуру, получим конечные разности третьего порядка:
Для конечных разностей -го порядка:
В результате получим таблицу конечных разностей:
| 0,41 | 2,63 |
| 1,55 | 3,75 |
| 2,67 | 4,87 |
| 3,84 | 5,03 |
| Вариант | |||
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 | -1 0 1 2 2 3 4 5 0 2 4 6 7 9 11 13 -3 -1 1 3 1 2 3 4 -2 -1 0 1 2 4 6 8 -4 -2 0 2 -1 1 3 5 2 4 6 8 -9 -7 -5 -3 0 1 2 3 -8 -5 -2 1 -7 -5 -3 -1 1 4 7 10 7 8 9 10 -4 0 4 8 -3 -1 1 3 0 3 6 9 0,7 0,8 0,9 1,0 2,7 2,75 2,8 2,85 3 3,5 4 4,5 10 14 18 22 2 4 6 8 | -3 5 2 3 4 1 7 6 -1 -4 2 -2 2 -2 3 2 7 -1 4 3 -3 -7 2 3 4 9 1 1 9 -3 6 4 2 8 5 3 4 -7 1 -2 -1 -6 3 4 3 -3 4 -4 7 -1 8 0 9 -2 4 2 4 -4 5 2 -2 9 3 0,5 6 -2 7 0 4 8 -2 2 11 -1 6 4 1 5 -4 -2 0,7 1 6 11 0,7 0,8 0,95 1,2 1,7 1,8 1,6 1,4 9,8 9,7 9,4 8,5 2,2 4,2 5,1 1,9 | 2 3,5 1 8 0 0,5 0,1 3 -1 2 5 -3 1,5 -3,2 -6,3 3,5 9,8 -2 0,4 5 0,83 2,72 3,9 23 8 |