
- •Введение
- •1. Метод прогонки
- •1.1. Метод прогонки для трехдиагональных матриц
- •1.2. Метод прогонки для пятидиагональных матриц
- •2. Сплайн – интерполяция
- •2.1. Интерполяционный линейный сплайн
- •Пример программы, которая по интерполяционной таблице строит интерполяционный линейный сплайн, приводится в примере выполнения задания 1 расчетно-графического задания 1.
- •Нахождение коэффициентов интерполяционного линейного сплайна
- •Нахождение значений интерполяционного линейного сплайна
- •2.2. Интерполяционный параболический сплайн
- •Нахождение коэффициентов интерполяционного параболического сплайна
- •Нахождение значений интерполяционного параболического сплайна и его производной
- •Результат работы программы
- •2.3. Интерполяционный кубический сплайн
- •Нахождение коэффициентов естественного интерполяционного кубического сплайна
- •Нахождение значений естественного интерполяционного кубического сплайна и его производных
- •2.5. Построение интерполяционных сплайновых кривых при помощи сплайн - функций
- •2.6. Примеры решения задач
- •3. Сглаживание кубическими сплайнами
- •3.1. Постановка задачи сглаживания
- •Определение естественного сглаживающего кубического сплайна
- •3.2. Построение естественного сглаживающего кубического сплайна
- •Нахождение значений естественного сглаживающего кубического сплайна и его производных
- •4. Аппроксимация
Нахождение значений интерполяционного линейного сплайна
Если нам уже известны
все коэффициенты
и
,
интерполяционного линейного сплайна
,
то для нахождения его значения в любой
точке
из отрезка
нужно знать еще номер
отрезка
,
такого что
.
Воспользуемся простой
формулой для нахождения этого номера:
.
Следовательно,
.
Недостатком этой формулы является то,
что при
мы получаем
,
а индекс
меняется в пределах от
до
.
В этом случае можно поступить следующим
образом. Введем дополнительные
коэффициенты
и
.
Определим их следующим образом
,
.
Сложность
вычислительного алгоритма построения
интерполяционного линейного сплайна.
Число арифметических действий, необходимых
для построения интерполяционного
линейного сплайна, пропорционально
числу отрезков (),
объем памяти также пропорционален числу
отрезков (
).
2.2. Интерполяционный параболический сплайн
Из
названия сплайна понятно, что параболический
сплайн – это функция, состоящая из
«кусочков» парабол. Эти «кусочки»
состыкованы таким образом, чтобы
параболический сплайн являлся непрерывно
дифференцируемой функцией. Обозначим
параболический сплайн. Для каждого
из отрезка
,
где
,
и
- числовые коэффициенты. Таким образом,
на каждом отрезке
,
три неизвестных коэффициента:
,
и
,
а на всем отрезке
число неизвестных коэффициентов равно
.
Для того чтобы однозначно определить
интерполяционный параболический сплайн
на отрезке
,
нам требуется
уравнений относительно
неизвестных
,
и
.
Так
как
является интерполяционной функцией,
то первые
уравнения получаем из условий
,
.
Во всех внутренних узлах сетки функция
должна быть непрерывной и непрерывно
дифференцируемой функцией, следовательно,
получаем еще
уравнения. В сумме получаем
линейных уравнений. Нам необходимо еще
одно условие. Естественно ввести краевое
условие, то есть условие либо в точке
,
либо в точке
.
Так как условия на функцию
в этих точках уже заданы, введем условие
на
.
Для того чтобы однозначно
определить интерполяционный параболический
сплайн, требуется одно краевое (граничное)
условие. Либо задается условие в точке
,
либо в точке
:
.
В дальнейшем мы будем предполагать, что
задано левое краевое условие
.
Определение
интерполяционного параболического
сплайна. Пусть
на отрезке
задана сетка
,
в узлах которой заданы значения
,
.
Интерполяционным параболическим
сплайном называется функция
,
удовлетворяющая следующим условиям:
-
функция
– непрерывно дифференцируемая функция на
;
-
на каждом из отрезков
функция
является полиномом второй степени вида
,
;
-
функция
– интерполяционная функция, то есть:
,
;
-
краевому условию
.
Отметим, что, отбросив условие три, мы получаем определение параболического сплайна.
Пример. Функция
является
параболическим сплайном, определенным
на отрезке
и удовлетворяющим краевому условию
.
Эта же функция является интерполяционным
параболическим сплайном, удовлетворяющим
следующей интерполяционной таблице:
|
0 |
1 |
2 |
|
4 |
2 |
1 |
Теорема
существования и единственности. Пусть
на отрезке
заданы интерполяционная таблица
,
,
причем все узлы сетки различны (
при
)
и краевое условие
.
Тогда существует единственный
параболический сплайн
такой, что
и краевое условие
.
Другими словами, если
задана интерполяционная таблица
,
в которой все узлы сетки различны, и
краевое условие, то существует единственный
интерполяционный параболический сплайн,
удовлетворяющий этой таблице и краевому
условию.