![](/user_photo/66067_jeDxE.png)
- •Тема.3. Интерполяция функций
- •3.1. Постановка задачи
- •3.2. Интерполяционная формула Лагранжа
- •3.3. Интерполяционные формулы Ньютона
- •3.3.1. Конечные разности
- •3.3.2. Первая интерполяционная формула Ньютона
- •3.3.3. Вторая интерполяционная формула Ньютона
- •1.1.13.4. Сплайн – интерполяция
- •3.5. Сравнение интерполяционных многочленов по применению
- •3.6. Технология интерполяции функций в среде математических пакетов
- •3.6.1. Технология интерполяции функций в среде системы MathCad
Тема.3. Интерполяция функций
3.1. Постановка задачи
Вычисление значений функции y = f(x) – одна из тех задач, с которой приходится постоянно сталкиваться в инженерной практике. Однако сделать это не всегда возможно. Примером тому следующие типичные ситуации:
функция задана таблицей значений (нет аналитического выражения)
, (i = 0, 1, 2,…, n), необходимо вычислить значения функции в точках, не совпадающих с табличными;
аналитическое выражение f(x) есть, но получение ее значений затруднено громоздкими и сложными вычислениями;
значения функции в требуемых точках могут быть получены только экспериментально.
В этих и ряде других случаев возникает необходимость приближенного вычисления функции y = f(x).
Задача аппроксимации состоит в следующем. Функцию f(x), заданную таблично, требуется приближенно заменить (аппроксимировать) некоторой функцией j(х) так, чтобы отклонение j(х) от f(x) в некоторой области удовлетворяло заданному условию. Функция j(х) называется аппроксимирующей функцией.
В качестве аппроксимирующей функции часто используют алгебраический многочлен вида:
jm(x) = a 0 + a 1 x + a 2 x2 + … + a m xm . (3.1-1)
В этом случае говорят о параболической аппроксимации.
Частным случаем
задачи аппроксимации таблично заданной
функции является интерполирование.
Интерполирование состоит в следующем.
Для функции y
= f(x),
заданной в (n
+ 1) точке
,
найти функцию j(х),
принимающую в этих точках заданные
значения, то есть
,
i
= 0, 1, 2, … n.
(3.1-2)
Будем называть (3.1-2) условием интерполяции, точки – узлами интерполяции, а функцию j(х) – интерполирующей функцией.
При интерполяции критерием приближения аппроксимирующей функции к заданной является совпадение их значений в узлах интерполяции.
Геометрической
интерпретацией задачи интерполяции
является нахождение функции, график
которой проходит через заданную систему
точек
,
i
= 0, 1, …, n
(рис. 3.1-1). Если в качестве интерполирующей
функции используется алгебраический
многочлен
(3.1-1) степени не выше n,
то задача имеет единственное решение.
|
___ интерполируемая функция
----- интерполирующая функция
|
Рис.3.1-1
Применяя интерполирующую функцию (3.1-1), запишем условие (3.1-2) для каждого из (n + 1) узлов. В результате получим следующую систему (n + 1) линейных уравнений:
Эта система однозначно разрешима, так как ее определитель (определитель Вандермонда) отличен от нуля, если узлы интерполяции различны. Решение полученной системы n+1 линейных уравнений относительно неизвестных а0, а1, …, аn позволяет найти коэффициенты интерполирующего многочлена (3.1-1).
Пример 3.1-1. Пусть функция y = f(x) задана таблично:
xi |
1 |
1.2 |
1.4 |
1.6 |
1.8 |
y i |
0 |
-0.16 |
-0.24 |
-0.24 |
-0.16 |
Требуется построить интерполяционный многочлен, позволяющий вычислить значение f(x) в точке x = 1.43.
Полагая x0 = 1.2 , x1 = 1.4 , x2 = 1.6,
y0 =-0.16, y1 = -0.24, y2 = -0.24, получим систему уравнений
Решая систему уравнений, получим следующие значения а0 = 2, а1 = -3, а2 = 1. Тогда интерполяционный многочлен имеет следующий вид: P2(x) = 2 – 3x + x2, а значение многочлена в точке 1.43 равно P2(1.43) = - 0.2451.