Введение
Интерполяция - это определение промежуточных значений функции по известному дискретному набору значений функции.
Этот способ применяется во многих сферах жизни человека, в частности, в научных и инженерных расчетах, где часто приходится оперировать наборами значений с функциями заданными таблично, полученных опытным путём (при условии, что мы не учитываем различные посторонние шумы при снятии показаний приборов) или методом случайной выборки.
На основании этих наборов требуется построить функцию, на которую могли бы с высокой точностью попадать другие получаемые значения. Такая задача называется аппроксимацией. В итоге можно сказать, что интерполяцией называют такую разновидность аппроксимации, при которой кривая построенной функции проходит точно через имеющиеся точки данных.
Интерполяционные методы применяются во многих сферах жизни человека, в частности, в научных и инженерных расчетах, где часто приходится оперировать наборами значений, полученных опытным путем, не учитывая различны шумы, или методом случайной выборки.
Задача интерполяции
Пусть функция 
задана таблицей своих значений на интервале [ 
:
, 
n - количество узлов
(1.1)
Задача интерполяции - найти функцию 
, принимающую в точках 
те же значения 
.
При этом предполагается, что среди значений 
нет одинаковых. Точки 
называют узлами интерполяции. Узлы интерполяции не обязательно должны располагаться равномерно на отрезке [ .
Функция ƒ(x) называется интерполянтом функции g(x)
Если 
ищется на интервале [ , то эту задачу принято называть задачей интерполяции, а если за пределами этого интервала, то это задачей экстраполяции.
Задача имеет много решений, т.к. через заданные точки, i=0, 1,..., n, можно провести бесконечно много кривых, каждая из которых будет графиком функции, для которой выполнены все условия (1.2).
В зависимости от цели приближения используют либо интерполяцию (точечную аппроксимацию), либо аппроксимацию. Аппроксимация – это замена таблично заданной функции g(x) функцией f(x), которая на рассматриваемом отрезке имеет ограниченное отклонение от функции g(x).
Условие интерполяции:
Где а – вектор неизвестных коэффициентов.
Решить задачу интерполяции значит найти при заданных и.
В общем виде система представляет систему нелинейных уравнений и при больших n часто не имеет решений.
Первым методом решений задачи интерполяции является метод Лагранжа.
Интерполирование состоит в приближённой замене функции g(x), заданной таблично, функцией f(x), которая принимает те же значения, что и функцияg(x).
Простейшим и наиболее часто применяемым функцией f(x) является полином:
где а0, а1, а2, …., аm – коэффициенты полинома,
m – степень аппроксимирующего многочлена.
Все методы интерполяции можно разделить на локальные и глобальные. В случае глобальной интерполяции отыскивается единый полином на всем интервале [ . Методы глобальной интерполяции обычно применяют для функций, заданных небольшим количеством точек, т. к. при увеличении количества точек увеличивается порядок интерполирующего многочлена, что отрицательно сказывается на гладкости получаемой функции. Многочленная аппроксимация, использующая сразу все узлы таблицы (глобальная интерполяция) имеет существенный недостаток – возможность появления больших экстремумов в промежутках между узлами сетки. Т.е. интерполяционный полином может иметь колебания, не свойственные исходным данным. Кроме того, с ростом степени полинома происходит быстрое накопление ошибок округления. Чтобы избежать этих нежелательных эффектов, на практике применяют локальную интерполяцию.. В случае локальной интерполяции на каждом интервале [ 
] строится отдельный полином. Для локальной интерполяции количество узлов большого значения не имеет.
Рассмотрим некоторые виды локальной и глобальной интерполяции.
Локальная интерполяция:
1. Кусочно–линейная интерполяция
2. Интерполяция сплайнами
Глобальная интерполяция:
1. Полином Лагранжа
2. Многочлен Ньютона
Глобальная интерполяция
1. Интерполяция полиномом Лагранжа
При глобальной интерполяции на всем интервале [a,b] строится единый многочлен. Одной из форм записи интерполяционного многочлена для глобальной интерполяции является многочлен Лагранжа:
Интерполяционный полином Лагранжа n-ой степени есть линейная комбинация базисных полиномов Лагранжа:
(1.7)
где 
– базисные многочлены степени n:
(1.8)
То есть многочлен Лагранжа:
(1.9)
Многочлен 
удовлетворяет условию.
Это условие означает, что многочлен равен нулю при каждом кроме, то есть – корни этого многочлена. Таким образом, степень многочлена 
равна n и при в сумме обращаются в нуль все слагаемые, кроме слагаемого с номером i=j, равного.
принимает значение 1 в точке 
и 0 в остальных узлах интерполяции. Следовательно в точке исходный полином принимает значение 
Выражение () применимо как для равноотстоящих, так и для не равноотстоящих узлов. Погрешность интерполяции методом Лагранжа зависит от свойств функции g (x), от расположения узлов интерполяции и точки x. Полином Лагранжа имеет малую погрешность при небольших значениях n (n <20). При больших n погрешность начинает расти, что свидетельствует о том, что метод Лагранжа не сходится (т.е. его погрешность не убывает с ростом n).
Многочлен Лагранжа в явном виде содержит значения функций в узлах интерполяции, поэтому он удобен, когда значения функций меняются, а узлы интерполяции неизменны. Число арифметических операции, необходимых для построения многочлена Лагранжа, пропорционально и является наименьшим для всех форм записи. К недостаткам этой формы записи можно отнести то, что с изменением числа узлов приходится все вычисление проводить заново.
2. Многочлен Ньютона
Пусть функция g(x) задана с произвольным шагом и точки таблицы значений занумерованы в произвольном порядке.
Многочлен Ньютона во многом опирается на понятие разделенных разностей.
Разделенные разности нулевого порядка совпадают со значениями функции в узлах. Разделенные разности первого порядка определяются через разделенные разности нулевого порядка:
(1.3)
Разделенные разности второго порядка определяются через разделенные разности первого порядка:
(1.4.)
Разделенные разности k -го порядка определяются через разделенную разность порядка 
:
(1.5)
Используя понятие разделенной разности интерполяционный многочлен Ньютона можно записать в следующем виде:
(1.6)
Для повышения точности интерполяции в сумму могут быть добавлены новые члены, что требует подключения дополнительных узлов. При этом для формулы Ньютона безразлично, в каком порядке подключаются новые узлы, в то время как для многочлена Лагранжа при добавлении новых узлов все расчеты надо производить заново.
Предположим, что необходимо увеличить степень многочлена на единицу, добавив в таблицу еще один узел 
. Для вычисления достаточно добавить к лишь одно слагаемое.
Локальная интерполяция.
1. Кусочно-линейная интерполяция.
Одним из самых используемых и простейших видов локальной интерполяции, является кусочно-линейная интерполяция, при которой каждые две точки 
и 
табличной функции 
соединяются отрезками прямой(т.е. проводится полином первой степени)
. (1.10)
А сама функция приближается к ломаной с вершинами в данных точках(рис.1).
Рис.1
Коэффициенты уравнений прямых 
и 
разные на каждом интервале 
легко находятся из условий интерполяции на концах отрезка:
(1.11)
Кусочно-линейная интерполяция является самой простой, и поэтому довольно часто применяется для расчета значений между узлами интерполяции. Для построения интерполирующей зависимости, используемой в дальнейших научных и инженерных расчетах, обычно используются более сложные методы интерполяции.
2. Интерполяция сплайнами
Иногда требуется обеспечить непрерывность не только интерполирующей функции, но и нужного количества её производных для этого прибегают к интерполяции сплайнами.
Сплайн – функция, область определения которой разбита на конечное число отрезков, на каждом из которых сплайн совпадает с алгебраическим многочленом. Максимальная степень из использованных полиномов называется степенью сплайна.
Преимущества интерполяции сплайнами по сравнению с обычными методами интерполяции – в сходимости и устойчивости вычислительного процесса. На практике наиболее часто используются кубические сплайны - сплайны третьей степени с непрерывной, по крайней мере, первой производной. При этом величина, называется наклоном сплайна в точке (узле).
Разобьём отрезок [a,b] на N равных отрезков [, ], где, i=0,1,…,N-1,
Если в узлах заданы значения, которые принимает кубический сплайн, то на частичном отрезке [ ] он принимает вид:
. (1)
В самом деле, это легко проверить, рассчитав и в точках
Можно доказать, что если многочлен третьей степени, принимает в точках значения и имеет в этих точках производные, соответственно,, то он совпадает с многочленом (1).
Таким образом, для того, чтобы задать кубический сплайн на отрезке, необходимо задать значения i=0,1…,N в N+1 в узле.
Ошибка интерполяции и многочлены Чебышева.
При интерполяции функции всегда получают ошибку состоящую из погрешности самого метода и ошибок округления.
Ошибка приближения функции интерполяционным полиномом n -й степени в точке x определяется разностью
. (1.12)
Можно показать, что погрешность определяется следующим выражением
. (1.13)
Здесь – производная (n +1) порядка функции в некоторой точке, а функция определена как
(1.14)
Если максимальное значение производной равно
,
то для погрешности интерполяции следует оценка
. (1.15)
Конкретная величина погрешности в точке x зависит, очевидно, от значения функции в этой точке. Качественный характер зависимости показан на рис..
Из рисунка видно, что погрешность интерполяции тем выше, чем ближе точка x лежит к концам отрезка. За пределами отрезка интерполяции (т.е. при экстраполяции) быстро растет, поэтому погрешность возрастает существенно.
Рис. 2
Вследствие описанного поведения погрешности, глобальная интерполяция в некоторых случаях может давать совершенно неудовлетворительный результат.