Свойство 1. Пусть 5-сплайн построен на неубывающей по
следовательности т л-2 узлов th tM, ..., ti+m+i. Тогда 5-сплайн N/"(f) отличен от нуля только в области ?, < / < ti+m+i и Nf"(f) = 0 при t < t,- и при t > ti+m+l. Действительно, как только t < между точками /, и /,+га+1 усеченная функция примет вид ст,„(г) = (z ~ ?)"' и в силу того, что (z - t)m является полиномом степени, меньшей т + 1, любая разделенная разность (т + 1)-го порядка, построенная по этому полиному, равна нулю. Как только t > ti+m+u разделенная разность равна нулю в силу того, что между узлами t, и /,+т+] ст„,(г) = 0. Только при t, < t < Г,+т+| разделенная разность ст,„[/,, tM,..., /,+т+| | отлична от нуля. Таким образом, 5-сплайны являются локальными функциями, принимающими ненулевые значения на отрезке, содержащем узлы построения 5-сплайна.
Свойство 2. N,m(t) > 0, т.е. 5-сплайны при любых значениях параметра принимают неотрицательные значения. Это следует из свойства 1и рекуррентного соотношения (1.48) при неубывающем расположении узлов.
Свойство 3. Пусть задана бесконечная (или достаточно длинная) неубывающая последовательность узлов г, и на каждых т + 2 подряд расположенных узлах построен 5-сплайн т - го порядка. Покажем, что для любого /, < 1 < tM отличны от нуля только т + 1 5-сплайнов, а именно, Njmm(t), Nn\ m(0, ■■■, и их сумма равна единице. Действительно, все 5-сплайны Nf(t), где j < i-m, равны нулю, так как построены на узлах, где am(z) = 0; все 5-сплайны Nj”(t), где j > i + 1, равны нулю, так как построены на узлах, где am(z) = (z - t)m является полиномом степени меньшей т + 1. Таким образом, сумма всех 5-сплайнов при t,< t < tM зависит от конечного числа 5-сплайнов, а именно:
XA74/)= х N7«)- (1-51)
j J-i-m
Если воспользоваться определением разделенной разности (1.37), то 5-сплайн (1.44) можно записать в виде
Njm(t) = <3m[tj+\, tj+2, -., tJ+m+1] - ojtj, tJ+b ..., tJ+m\. (1.52)
Используя равенства (1.51) и (1.52), найдем сумму всех 5-сплайнов при t, M
ГЕОМЕТРИЧЕСКОЕ МОДЕЛИРОВАНИЕ 2
fHib 13
«)-£■ ^ р„ 15
... -h 18
+ q,+i(-^2 + w3)(thi - h) = a0(^)p, + «i(^)p,+i + Po(^)(^+i - Oq, + 18
г, = 27
<{ 29
В (1.53) было использовано, что om\tM, ti+2, ..., ti+m+,] = 1 как коэффициент при наивысшей степени аргумента полинома (z - t)m, которым описывается усеченная степенная функция на открытом интервале tM < i+m+1. При этом на интервале /,•_ т усеченная степенная функция равна нулю, следовательно, ее разделенная разность о,,,!/, т, tM. ..., ?,] = 0. Так как участок t,< tM был выбран произвольно, аналогичное равенство выполняется для любого другого участка. Таким образом, при любом t для суммы всех 5-сплайнов справедливо равенство
2>"(0 = 1. (1-54)
j
Соотношение (1.54) 5-сплайнов аналогично соотношению (1.18) для базиса Бернштейна.
Таким образом, при любом фиксированном параметре t совокупность всех ненулевых функций 7V/"(0 представляет собой разложение единицы.
Свойство 4. Для неубывающей последовательности узлов th tM, ..., ti+m+1 5-сплайн вычисляется с помощью рекуррентного соотношения
Nm (?) = ДГ тЧ (7) + ДГШ-1 (/)> (j 55)
^i+m+\ ~Ч+\ h+m~*i
начиная с ^-сплайнов нулевого порядка
(156)
[О, в противном случае.
Соотношение (1.55) следует из рекуррентной формулы Кокса—Де Бура (1.48) при подстановке в нее соотношений между 5-сплайнами и соответствующими ненормированными 5-сплайнами
Л/"(/) = (ti+m„ - /,) (1.57)
Свойство 5. Производная 5-сплайна т-го порядка выражается через 5-сплайны (т - I )-го порядка.
Для доказательства этого утверждения продифференцируем 5-сплайн, построенный на неубывающей последовательности узлов th tM, ..., /I+m+| и представленный в виде (1.52):
dN?(t) = da m\ti+uti+2, -,/<ЧвЫ] dajt
iу ti+1> '''~ 0• n: 1 dt dt dt
Подставив в последнее выражение производные по t усеченной степенной функции
dt
получим
dN™(t)
—±— = [tM,ti+2,..., ti+m^\ + rnam_\\tn tM,..., ti+m\ =
dt
= (t) + тМ™~х (t) =
= iV”r'(0 + ———Ar,m-|(0- (1-58)
/ _ / / _ t
li+m+1 /+1 */+w */
Аналогично получим, что производная ненормированного 5-сплайна /и-го порядка выражается через ненормированные 5-сплайны (т -
dM™(t) _ МГ'(*)~
tj+rя-t-l _ ti
Do'stlaringiz bilan baham: |