Составители:
интерполяционные формулы при заданной таблицу значений функции. Это следует из единственности
интерполяционного многочлена заданной степени (при отсутствии совпадающих узлов интерполяции).
3.4 Сплайны
Рисунок 6.3
Сейчас широкое распространение для интерполяции получило использование кубических сплайн-функций —
специальным образом построенных многочленов третьей степени. Они представляют собой некоторую
математическую модель гибкого тонкого стержня из упругого материала. Если закрепить его в двух соседних
узлах интерполяции с заданными углами наклонов A и B (Рисунок 6.3), то между точками закрепления этот
стержень (механический сплайн) примет
некоторую форму, минимизирующую его потенциальную энергию.
Пусть форма этого стержня определяется функцией y=S(x). Из курса сопротивления материалов известно, что
уравнение свободного равновесия имеет вид S
IV
(x)= 0. Отсюда следует, что между каждой парой соседних
узлов интерполяции функция S(x)является многочленом третьей степени. Запишем ее в виде
3
1
2
11
)()()()(
−−−
−+−+−+=
iiiii
xxdxxcxxbaxS
(6.38)
ii
xxx ≤≤
−1
Для определения коэффициентов a
i
, b
i
, c
i
, d
i
на всех n элементарных отрезках необходимо получить 4n
уравнений. Часть из них вытекает из условий прохождения графика функции S(x) через заданные точки, т.е
S(x
i-1
)=y
i-1
. S(x
i
)=y
i
. Эти условия можно записать в виде
11
)(
−−
=
=
iii
yaxS (6.39)
iii
i
iiiii
yhdhchbaxS =+++=
32
)(
(6.40)
1−
−
=
iii
xxh , i=1, 2, , …, n.
Эта система содержит 2n уравнении. Для получения недостающих уравнений зададим условия непрерывности
первых и вторых производных в узлах интерполяции, т. е. условия гладкости кривой во всех точках.
Вычислим производные многочлена (6.38).
2
11
)(3)(2)(
−−
−+−+=
′
iiiii
xxdxxcbxS
,
)(62)(
1−
−
+
=
′
′
iii
xxdcxS .
Приравнивая в каждом внутреннем узле х = x
i
значения этих производных, вычисленные в левом и правом от
узла интервалах, получаем 2n—2 уравнений
Страницы
- « первая
- ‹ предыдущая
- …
- 86
- 87
- 88
- 89
- 90
- …
- следующая ›
- последняя »
