Кубический интерполяционный сплайн

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

Итак, на каждом -м отрезке , , решение будем искать в виде полинома третьей степени:

Неизвестные коэффициенты , , находим из:

· условий интерполяции: , , ;

· непрерывности функции , ;

· непрерывности первой и второй производной:

, , .

Учитывая, что ,

для определения неизвестных получаем систему уравнений:

, ;

, ;

, ;

, ,

где _. Недостающие два уравнения выводятся из дополнительных условий: . Можно показать, что при этом . Из системы можно исключить неизвестные , , получив систему линейных уравнений (СЛАУ) для определения коэффициентов :

, ,

, (3.1)

Матрица СЛАУ (3.1) является ленточной, поэтому система решается с помощью прямого экономичного метода прогонки. После этого вычисляются коэффициенты , :

, , . (3.2)

В случае постоянной сетки , и тогда система уравнений упрощается:

Для вычисления значения в произвольной точке отрезка необходимо решить систему уравнений на коэффициенты , , затем найти все коэффициенты , . Далее, необходимо определить, на какой интервал попадает эта точка, и, зная номер , вычислить значение сплайна и его производных в точке

,

,

.

ПРИМЕР 3.2. Задана таблица значений функции

		1/4	1/2	3/4
f

Требуется вычислить значения функции в точках 0, 2 и 0, 8, используя сплайн-интерполяцию.

В нашем случае: , .

Выпишем систему уравнений для определения :

.

Решая эту систему линейных уравнений, получим:

.

, ,

, .

,

,

,

,

.

Рассмотрим точку 0, 2, которая принадлежит первому отрезку, т.е. . Следовательно, получим,

Итак, значение функции в точке 0, 2 равняется 1, 944.

Рассмотрим точку 0.8, которая принадлежит четвертому отрезку, т.е. .

Следовательно,