CC BY
16
УДК 517.5, 519.6 Вестник СПбГУ. Математика. Механика. Астрономия. Т. 4 (62). 2017. Вып. 1
МЯО 68'№10
НЕОБХОДИМЫЕ И ДОСТАТОЧНЫЕ УСЛОВИЯ НЕОТРИЦАТЕЛЬНОСТИ КООРДИНАТНЫХ ТРИГОНОМЕТРИЧЕСКИХ СПЛАЙНОВ ВТОРОГО ПОРЯДКА*
Ю. К. Демьянович, А. А. Макаров
Санкт-Петербургский государственный университет,
Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9
В работе получены необходимые и достаточные условия неотрицательности непрерывно дифференцируемых координатных тригонометрических сплайнов второго порядка и установлены промежутки выпуклости и вогнутости этих сплайнов. Метод исследования состоит в определении вогнутости на промежутках, прилегающих к концам носителя рассматриваемого координатного сплайна, и в применении рассуждений, связанных с числом нулей решения соответствующей краевой задачи для дифференциального уравнения второго порядка. Библиогр. 16 назв. Ил. 2.
Ключевые слова: неотрицательность, координатные сплайны, тригонометрические сплайны.
1. Введение. В работах [1-6] рассматривались разные варианты построения аппроксимаций. Среди них фигурировали полиномиальные сплайны (кусочно-полиномиальные функции) различной гладкости, а также неполиномиальные, а именно: тригонометрические и экспоненциальные сплайны; сплайны, определяемые тем или иным оператором (так называемые Ь-сплайны); сплайны, определяемые семейством линейных функционалов (С-сплайны); сплайны, получаемые минимизацией квадратичного функционала и т. п. (см., например, библиографию в монографиях [23]). Сплайны, которые получаются из аппроксимационных соотношений с использованием полной цепочки векторов и генерирующей вектор-функции, рассматривались в ряде работ (см., например, [6-8]); подобные сплайны с минимальным носителем были названы минимальными. Определенный способ выбора упомянутой цепочки векторов позволил рассмотреть минимальные сплайны максимальной гладкости и установить единственность пространства таких сплайнов среди множества пространств минимальных сплайнов (определяемых произвольным выбором полной цепочки векторов при фиксированной сетке и фиксированной генерирующей вектор-функции). Оказалось, что такие пространства образуют цепочку вложенных пространств на последовательности вложенных сеток, а это, в свою очередь, позволило построить вэйвлетное разложение таких цепочек (см., например, [9]).
С другой стороны, минимальные сплайны предоставляют прекрасный аппарат аппроксимации (ибо они получаются из аппроксимационных соотношений). Аппроксимация сплайнами является линейной комбинацией координатных сплайнов, так что наблюдение за качеством аппроксимации и разработка путей ее улучшения значительно упрощаются в случае, когда координатные сплайны имеют определенный знак.
Координатный сплайн будем называть положительным, если он положителен во всех внутренних точках своего носителя.
* Работа частично выполнена при финансовой поддержке РФФИ (гранты №15-01-08847, №1631-60060 мол_а_дк).
(¡5 Санкт-Петербургский государственный университет, 2017
Положительность координатных сплайнов важна для практического применения в задачах геометрического моделирования (см., например, [10-14]), и в задачах изо-геометрического анализа (подробности см. в монографии [15]). Знакоопределенность координатных сплайнов важна в использующих эти сплайны методах для минимизации ошибок округления и для поддержания устойчивости вычислительного процесса.
Давно известно (см., например, [3]), что координатные сплайны, генерируемые вектор-функцией y(t) =f (1,t,t2), положительны: это непрерывно дифференцируемые квадратичные сплайны; в этом случае никаких ограничений на сетку не налагается.
Тригонометрические сплайны, получаемые из аппроксимационных соотношений, обладают асимптотически оптимальной аппроксимацией (по N-поперечнику стандартных компактов). Вложенность пространств таких сплайнов, построенных на вложенных сетках, позволяет получать сплайн-всплесковое (вэйвлетное) разложение [8]. Вопросы, связанные с поведением тригонометрических сплайнов и с положительностью координатных тригонометрических сплайнов, рассматривались в [7], где анонсировалась положительность упомянутых сплайнов для сетки с шагом меньше п, и в работе [16], где их положительность доказана для сетки шага меньше п/2.
Цель данной работы — получить необходимые и достаточные условия положительности непрерывно дифференцируемых координатных тригонометрических сплайнов второго порядка. Метод исследования состоит в определении вогнутости на промежутках, прилегающих к концам носителя рассматриваемого координатного сплайна, а также в применении рассуждений, связанных с числом нулей решения соответствующей краевой задачи для дифференциального уравнения второго порядка.
2. Общая схема построений. Пусть z — множество всех целых чисел, а r3 — множество трехмерных векторов с компонентами из множества r1 всех вещественных чисел. На интервале (a,ß) С r1 вещественной оси r1 введем сетку {xj}jez,
. . . < Х-2 < X_1 < Xo < X1 < X2 < ...,
lim = a, lim = ß.
j^—w j^+w
Введем множества Md==f Ujez (xj, Xj+1) и Sj =f [xj_ 1, Xj+2]. Пусть {aj}jez — множество трехмерных вектор-столбцов из r3; предположим, что квадратные матрицы третьего порядка (ak—1, ak, ak+1), составленные из вектор-столбцов ak—1, ak, ak+1, неособенные:
det(afc_1, afc, a^) =0, Vk e z. (2.1)
Рассмотрим трехкомпонентную вектор-функцию (столбец) y(t) с компонентами из пространства C2(a,ß) и с отличным от нуля вронскианом
de%,^ ',<£ '')(t)=0, Vt e (a,ß).
Линейное пространство вещественнозначных функций, заданных на множестве M, обозначим через X.
Пусть функции Wj e X, j e z, удовлетворяют тождествам
ak_1^k—1 (t) + afcwk(t) + afc+1Wfc+1 (t) = y>(t), Vt e (Xfc, Xk+1), Vk e z, (2.2) W(t) = 0, Vt e m\sj, Vj e z. (2.3)
При каждом фиксированном £ € , ), Ук € Ж, соотношения (2.2) можно рассмотреть как систему линейных алгебраических уравнений относительно неизвестных ^ (£). Ввиду предположения (2.1) система (2.2) однозначно разрешима. Из (2.2) и (2.3) видно, что справедливо эирр ^ С ^ и
/(í) =
/(í) =
det a j _ 2, aj _ i, a j det(aj_i,y(t), aj+ij
/(t) =
det aj_i, aj, aj+i det( y(t), aj+i, aj+2
det aj, aj+b aj+2
при t G (xj_i, xj)
при t G (xj, Xj+i),
при t G (xj+i,Xj+2).
(2.4)
(2.5)
(2.6)
Введем обозначения у == у (ж),у^ == у'(ж).
Как известно (см. [10]), если а^- = е^-Ъ, У] € Ж, с произвольными ненулевыми константами е^- и векторами Ъ, определяемыми по формулам
bj tfdet(^j, ^j+i, ^ j+i^ j - det(^ j, ^j+i, ^ j+i,
(2.7)
то w G Ci(a,e), Vj G z.
3. Гладкие тригонометрические сплайны второго порядка. Пусть имеем y(t) = (1, sin t, cos t)T. В этом случае справедливо равенство det(^, <р', <р '')(t) = -1. Для отыскания векторов aj найдем векторы bj по формуле (2.7). Имеем
1
1 0
xj+i cos xj
^cos Xj cos Xj+i -sin Xj+i,
det(^j, fj+i, tp j+i) = det I sin Xj sin Xj+i cos Xj+i I = — 2sin2 Жд+1——.
а кроме того,
0 1 0
det(<^> j, y>j+i, ^ j+i) = det I cos Xj sin Xj+i cos Xj+i I = sin (Xj+i — Xj).
4—sin Xj cos Xj+i —sin Xj+i,
Итак, получаем
b, = —2sin2 Xj+1 Xj j-
0
cos Xj sin Xj
— sin (Xj+i — Xj)
cos Xj
Дальнейшие преобразования дают
bj = —2sin Xj+1 Xj
2
cos sin
cos
xj+i xó 2,
Xj + l^Xj 2
2
Предположим, что выполнено условие (А):
— X' < п, У] € ж.
1
sin X
j
Полагая а,- = — ^бш
ры а, в виде
, после подстановки а, = е, Ъ, определяем векто-
а, = I соэ —-эш —-соэ —--
(3.1)
Подставляя (3.1) в равенства (2.4)—(2.6) находим тригонометрический В^-сплайн шТ второго порядка с носителем эирр шТ = [х,—1 , х,+2]:
ш Т (£) = вт2
£ , — 1 \ , _1 { х, ~~+1 х, — 1 \ . _1 { х, х, — 1
вт
вт
при £ € [х,—1, х,), (3.2)
4 И - 2
• I х,+2 \ х, - х,—1
Э1П ^———-——--^сов——^--ь
| ( ^ + 1 I +
х вт
+ эт-
2
вт
2
х вт
при £ € [х,, х,+1),
(3.3)
+ 2 ].
(3.4)
4. Свойства тригонометрических сплайнов второго порядка. Из формул (3.2) и (3.4) видно, что при условии (А) на промежутках (х,—1 ,х,) и (х,+1, х,+2) функция шТ(£) положительна.
Найдем первую и вторую производные функции шТ(£) на промежутке [х,, х,+1 ]:
У® =
х,+1 I х,+2 л х, - х,—1 / х,—1 | хл х,+2 - х,+1
—СОЭ ( — ---—--¿]С08 —-^--Ьсов и------ ) СОЭ■
2 ) 2 V 2 ) 2
X 81П
2
Бт
2
81П
, (4.1)
х,+1 I х,+2 Л х, - х,—1 . / ——-——--Ь СОв —---81Пи-
, , — 1 » ( , — 1 I , \ ^х, —+ 2 ^х, —+1
X 81П
^лг,;^-! [х3+2 [х3+1 | (4 2)
2
2
Отыскивая первую и вторую производные функции шТ(£) на промежутках
(х, —1, х,] и [х, + 1,х, + 2 ), получим
а
1
^ V = 28Ч* " < > (^Т^К1 > О,
(4.3)
2
2
2
X
X
X
—8111
2
при £ € , х-] и
- ^ - (З^)™-1 С+2 ~ Хз+1) < 0, (4.5)
!2 т, , 1 ТЮ !:
при £ € [ж,-+1, х^+2).
Рассмотрим функцию и(Ь) = на промежутке [ж^-, . Используя (4.3),
(4.5) и учитывая непрывность производной функции щТ(£), придем к неравенствам
и(ж^) > 0, и(ж^+1) < 0. (4.7)
Из соотношений (4.7) следует, что (ввиду аналитичности) функция и(£) может иметь лишь нечетное количество нулей на промежутке [ж^-, ж^+1 ]. Легко видеть, что и(£) удовлетворяет уравнению и'' + и = 0, так что и(£) имеет представление
и(£) = К бш^ + в), (4.8)
где К и в определяются из граничных условий
и{хо) = =
с использованием формул (4.3) и (4.5) соответственно. Представление (4.8) показывает, что, для того чтобы на промежутке [с, !] функция и(£) имела менее трех нулей, необходимо и достаточно, чтобы выполнялось неравенство ! — с < 2п. Применяя это соображение к функции =£ на рассматриваемом промежутке [ж^,
видим, что справедливо следующее утверждение.
Лемма 1. При условии (А) функция (см. формулу (4-2)) на промежутке
Х, х^+1 ] имеет один нуль.
Следствие 1. При условии (А) функция на промежутке (ж^ж^+х)
неположительна.
Теорема 1. 1) Условие (А) необходимо и достаточно для того, чтобы все функции щТ(£) были положительны У] € Ж.
2) При условии (А) на промежутке (ж^-) функция щТ(£) положительна и выпукла. Если ж^- — < п/2, то на промежутке (ж^-_1, ж^) функция щТ(£) положительна и вогнута, а если п/2 < ж^- — < п, то эта функция (оставаясь положительной) вогнута на промежутке (ж^-_1 + п/2) и выпукла на промежутке (ж^-_1 + п/2, ж^-). Аналогичное утверждение относится к промежутку (ж^-+1, ж^+2): если ж^+2 — ж^+1 < п/2, то рассматриваемая функция положительна и вогнута на всем упомянутом промежутке, если же п/2 < ж^+2 — < п, то на промежутке (ж^+1, ж^+2 — п/2) эта функция (оставаясь положительной) выпукла, а на промежутке (ж^+2 — п/2, ж^+2) она вогнута.
Доказательство. Для доказательства необходимости, указанной в пункте 1), покажем, что в классе сеток, нарушающих условие (А), существует такая сетка, на которой хотя бы один координатный сплайн меняет знак.
Рассмотрим сетку, для всех сеточных интервалов которой выполняется условие (А), кроме одного, где это условие нарушено: сумма длин двух соседних интервалов больше 2п. Тогда найдется незнакоопределенная функция шТ(£). Действительно, пусть Хк+2 — Хк > 2п, а все промежутки [х^, х^+1), ] = к, удовлетворяют условию (А). Тогда функция (¿) меняет знак: на промежутке (хк+1 ,хк+2) она отрицательна, а на промежутке (хк_1 ,хк) положительна.
Достаточность, утверждаемая в пункте 1 теоремы, следует из формул (3.2)-(3.4), (4.1)-(4.8) и леммы 1; из этих же формул вытекают утверждения пункта 2. ■
На рис. 1, 2 приведены графики двух вариантов рассматриваемого координатного тригонометрического сплайна на равномерных сетках вида х^- = к > 0, для значений к = 1 и к = 5 соответственно. Эти графики иллюстрируют доказанный выше результат: при к < п упомянутый сплайн положителен, а при к > п он не положителен.
0,8 0,6 0,4
0,2
.........о] |
-2 -1 0
Рис. 1. при h = 1.
Замечание. В работе [9] допущена погрешность, замеченная А. А. Макаровым: в качестве основного условия фигурирует условие (A), но приведенные там доказательства справедливы лишь при более жестком предположении:
Xj+i — Xj < п/2, Vj g z.
Результаты данной работы устраняют эту погрешность и позволяют установить положительность тригонометрических сплайнов при выполнении условия (A).
Литература
1. Bezier P. E. Example of an Existing System in the Motor Industry: The Unisurf System // Proc. Roy. Soc. London. 1971. Vol.A321. P. 207-218.
2. Варга Р. Функциональный анализ и теория аппроксимации в численном анализе. М.: Мир,
1974.
3. .Завьялов Ю. С., Квасов Б. И., Мирошниченко В. Л. Методы сплайн-функций. М., 1980. 352 с.
4. Piegl L., Tiller W. Curve and Surface Constructions Using Rational B-splines // Сотр. Aid. Des. 1987. Vol. 19. P. 485-498.
5. Rogers D. F., Fog N. G. Constrained B-spline Curve and Surface Fitting // CADJ. 1989. Vol. 21. P.641-648.
6. Демьянович Ю. К. Локальная аппроксимация на многообразии и минимальные сплайны. СПб.: Изд-во С.-Петерб. ун-та, 1994. 356 с.
1 2 3 4 5
Тригонометрический сплайн
Рис. 2. Тригонометрический сплайн при h = 5.
7. Бурова И. Г., Демьянович Ю. К. О гладкости сплайнов // Математическое моделирование. 2004. Т. 16, №12. С. 40-43.
8. Демьянович Ю. К. Вложенные пространства тригонометрических сплайнов и их всплесковое разложение // Математические заметки. 2005. Т. 78, вып. 5. С. 658-675.
9. Демьянович Ю. К. Теория сплайн-всплесков. СПб.: Изд-во С.-Петерб. ун-та, 2013. 526 с.
10. Pena J. M. Shape preserving representations for trigonometric polynomial curves // Computer Aided Geometric Design. 1997. Vol.14, N1. P. 5-11.
11. Роджерс Д., Адамс Дж. Математические основы машинной графики. М.: Мир, 2001. 604 с.
12. Buchwald B., Muhlbach G. Construction of B-splines for generalized spline spaces from local ECT-systems // Journal of Computational and Applied Mathematics. 2003. Vol. 159. P. 249-267.
13. Квасов Б. И. Методы изогеометрической аппроксимации сплайнами. М.; Ижевск: НИЦ «Регулярная и хаотическая динамика», Институт компьютерных исследований, 2006. 416 с.
14. Xu G., Wang G.-Z. AHT Bezier Curves and NUAHT B-Spline Curves // Journal of Computer Science and Technology. 2007. Vol. 22, N 4. P. 597-607.
15. Cottrell J. A., Hughes Th. J. R., Bazilevs Yu. Isogeometric Analysis: Toward Integration of CAD and FEA. John Wiley & Sons. 2009. 360 p.
16. Демьянович Ю. К., Лебединский Д. М., Лебединская Н. А. Двусторонние оценки некоторых координатных сплайнов // Записки научн. семинаров ПОМИ. 2015. Т. 439. C. 74-92.
Статья поступила в редакцию 21 апреля 2016 г.; рекомендована в печать 6 октября 2016 г.
Сведения об авторах
Демьянович Юрий Казимирович — доктор физико-математических наук, профессор; [email protected]
Макаров Антон Александрович — доктор физико-математических наук, доцент; [email protected]
NECESSARY AND SUFFICIENT CONDITIONS OF NONNEGATIVITY FOR COORDINATE TRIGONOMETRICAL SPLINES OF THE SECOND ORDER
Yuri K. Dem'yanovich, Anton A. Makarov
St. Petersburg State University, Universitetskaya nab., 7—9, St. Petersburg, 199034, Russian Federation; [email protected], [email protected]
There are many different ways to define coordinate splines. The splines, which are defined with approximation relations, have the best approximation properties as to order of N-width for the standard compact sets. The suitable choice of the approximation relations gives the maximal smoothness under the condition of minimal support for coordinate splines. On the other hand, the choice of different generating function gives opportunity to receive splines of different types (polynomial, exponential, trigonometrical ones and so on).
The paper is discussed the positivity of coordinate trigonometrical splines of the second order with maximal smoothness.
Let us consider the real grid {xj}jgz,
... < X—2 < x-i < xo < xi < X2 < ...
By definition we put M d=2f Ujgz (xj, xj+i), Sj d=2f [xj—i,xj+2], and consider three-component vectors: generating vector function p(i) and family of vectors aj, defined with formulas
a, = ucu^ j V-j + 1, ^ j+i)^ j — jflt^ j + i j >
, def / \ / def < / \
where pj = p(xj),pj = p' (xj).
The coordinate trigonometric splines Uj € C1 are defined by approximation relations
afc—iwfc—i(t)+ afcwfc(i) + afc+iwfc+i(i) = p(t) Vi € (xfc,xfc+i), Vfc € Z, Uj(i) = 0 Vi € M, Vj € Z. Theorem. The inequality xj+i — xj < n Vj € Z is necessary and sufficient condition for nonnegativity of the functions Uj (i) Vj € Z.
Furthermore, the sufficient conditions of convexity on the intervals (xj—i,xj), (xj+i,xj+2) and concavity on the interval (xj, xj+i) have been obtained in the paper. Refs 16. Figs 2. Keywords: nonnegativity, coordinate splines, trigonometrical splines.
<^(i) =f (1, sin t, cos i)T, aj d= det(^j, <^j+i, ^ j+iV j - det(^ j, ^j+i, ^ j+i
References
1. Bezier P. E., "Example of an Existing System in the Motor Industry: The Unisurf System", Proc. Roy. Soc. A321, 207-218 (London, 1971).
2. Varga R. S., Functional Analysis and Approximations in Numerical Analysis (Kent State University, Kent, Ohio, 1971).
3. Zavjalov Yu.S., Kvasov B.I., Miroshnichenko V. L., Methods of Spline Functions (Moscow, 1980) [in Russian].
4. Piegl L., Tiller W., "Curve and Surface Constructions Using Rational B-splines", Соmp. Aid. Des. 19, 485-498 (1987).
5. Rogers D.F., Fog N.G., "Constrained B-spline Curve and Surface Fitting", CADJ 21, 641-648 (1989).
6. Dem'yanovich Yu. K., Local Approximation on Manifold and Minimal Splines. (St. Petersburg University Press, St. Petersburg, 1994) [in Russian].
7. Burova I. G., Dem'yanovich Yu. K., "On smoothness of Splines", J. Mathematical Modeling 16(12), 40-43 (2004) [in Russian].
8. Dem'yanovich Yu. K., "Embedded Spaces of Trigonometrical Splines and Their Wavelet Expansion", Mathematical Notes 78, issue 5, 615-630 (2005).
9. Dem'yanovich Yu.K., Theory of Spline-Wavelets (St. Petersburg University Press, St. Petersburg, 2013, 526 p.) [in Russian].
10. Pena J.M., "Shape preserving representations for trigonometric polynomial curves", Computer Aided Geometric Design 14(1), 5-11 (1997).
11. Rogers D.F., Adams J. A., Mathematical Elements for Computer Graphics (2nd ed., McGraw-Hill, Inc., New York, USA, 1990).
12. Buchwald B., Muhlbach G., "Construction of B-splines for generalized spline spaces from local ECT-systems", Journal of Computational and Applied Mathematics 159, 249-267 (2003).
13. Kvasov B. I., Methods of Iso-geometrical Spline Approximation (Regular and Chaotic Dynamics, Izhevsk, 2006) [in Russian].
14. Xu G., Wang G.-Z., "AHT Bezier Curves and NUAHT B-Spline Curves", Journal of Computer Science and Technology 22(4), 597-607 (2007).
15. Cottrell J. A., Hughes Th. J. R., Bazilevs Yu., Isogeometric Analysis: Toward Integration of CAD and FEA (John Wiley & Sons, 2009, 360 p.).
16. Dem'yanovich Yu.K., Lebedinskii D.M., Lebedinskaja N.A., "Two-Sides Estimates of Some Coordinate Splines", J. of Math. Sci. 216, issue 6, 770-782 (2016).
Для цитирования: Демьянович Ю. К., Макаров А. А. Необходимые и достаточные условия неотрицательности координатных тригонометрических сплайнов второго порядка // Вестник СПбГУ. Математика. Механика. Астрономия. 2017. Т.4 (62). Вып. 1. С. 9-16. DOI: 10.21638/11701/spbu01.2017.102
For citation: Dem'yanovich Yu. K., Makarov A. A. Necessary and sufficient conditions of nonnegativity for coordinate trigonometrical splines of the second order. Vestnik SPbSU. Mathematics. Mechanics. Astronomy, 2017, vol. 4(62), issue 1, pp. 9-16. DOI: 10.21638/11701/spbu01.2017.102
