logo
Шпоры ALL

22.2 Кубические сплайны дефекта 1

Это сплайны, которые произошли от гибкой линейки чертёжника при этом в узлах непрерывна должна быть 1-ая и 2-ая производная.

уравнение изогнутой оси балки, когда заданы изгибные моменты в каждой точке. По аналогии 2-ые производные будем называть моментами.

Пусть в n попарно различных точках отрезка [a,b] a<x1<x2<x3….xn=b

y1=f(x1)…. Yn=f(xn)

непрерывен на [a,b] вместе с 1-ой и 2-ой производной, совпадает с кубическим полиномом на каждом из отрезков и удовлетворяет условиям i=1..n

Поскольку 2ая производная по условию в каждом узле непрерывна, то введя обозначения i=1,2…n и предположив, что меняется линейно на интервале можно записать (1)

Проинтегрируем дважды (1). Интегрирование будем выполнять со скобкой. Тогда после 1-го интегрирования получим (2)

А после второго интегрирования получим

(3)

Это выражение мы получили для отрезка , аналогичное выражение мы можем получить для предыдущего отрезка

(4)

Используя условие непрерывности 1-ой производной для узла xi

Для этого приравняем (2) при x=xi к выражению (4) при x=xi . В результате получим

(5)

(6)

(7) i=2,3…n-1

Т.е. мы получили n-2 линейных алгебраических уравнения относительно Mi

Для замыкания этой системы нужно задать 2 дополнительных уравнения. Эти уравнения можно получить используя условие на концах интервала [a,b] (краевые условия)

1.Периодический сплайн T=b-a.

При этом предполагается, что интерполируемая функция периодическая с T=b-a и поэтому сам сплайн тоже должен быть периодическим с таким же периодом. Запишем условия периодичности:

Тогда можно потребовать, чтобы уравнение (7) было справедливо и для i=n, а условие периодичности примет вид

Тогда уравнение (7) примет вид

(8)

2.Непериодический сплайн с заданными наклонами в т a и b

заданное число

заданное число

Запишем уравнение (2) при i=1

Если записать уравнение (3) при i=n

Тогда замыкающая система примет вид

(9)

Решив системы 8 и 9 можно найти неизвестные M1…..Mn и можем построить сплайн в соответствии с формулой (3)

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

, где

i=1,2..n-1