Оглавление
Введение 5
1 О минимальных тригонометрических сплайнах 10
1.1 Построение непрерывных тригонометрическихтплайнов . . 10
1.2 Построение непрерывно дифференцируемых тригонометрических сплайнов первого порядка............................. 12
1.3 Построение гладких тригонометрических сплайнов второго
порядка................................................. 14
1.4 Построение гладких тригонометрических сплайнов третьего
порядка................................................. 19
1.5 О постановках интерполяционных задач.................... 27
1.5.1 Интерполяционные задачи для непрерывно дифференцируемых тригонометрических сплайнов первого порядка .................................................. 27
1.5.2 Интерполяционные задачи для гладких тригонометрических сплайнов второго порядка...................... 27
1.5.3 Интерполяционные задачи для гладких тригонометрических сплайнов третьего порядка..................... 28
1.6 О свойствах точности аппроксимации...................... 30
1.7 Оценки погрешностей..................................... 40
2 О В-сплайнах второй, четвертой и шестой степеней 43
2.1 Аппроксимационные соотношения и формулы базисных
сплайнов................................................ 43
2.1.1 В-сплайны второй степени.......................... 43
2
2.1.2 В-сплайны четвертой степени........................ 44
2.1.3 В-сплайны шестой степени........................... 46
2.2 О постановках интерполяционных задач..................... 49
2.2.1 Интерполяционные задачи для В-сплайнов второй степени ..................................................... 49
2.2.2 Интерполяционные задачи для В-сплайнов четвертой
степени............................................ 50
2.2.3 Интерполяционные задачи для В-сплайнов шестой
степени............................................ 51
2.3 Оценки погрешностей...................................... 52
2.3.1 Для решения интерполяционных задач с помощью В-
сплайнов второй степени............................ 52
2.3.2 Для решения интерполяционных задач с помощью В-
сплайнов четвертой степени......................... 54
2.3.3 Для решения интерполяционных задач с помощью В-
сплайнов шестой степени............................ 72
3 О свойствах гладких тригонометрических сплайнов 75
3.1 Построение гладких приближений, не использующих значения производных............................................... 75
3.2 Построение квадратурных формул, согласованных с гладкими тригонометрическими и полиномиальными сплайнами 98
3.2.1 Тригонометрический случай.......................... 99
3.2.2 Погрешность. Тригонометрический случай.............105
3.2.3 Полиномиальный случай..............................106
3.2.4 Погрешность. Полиномиальный случай.................110
3
3.3 Об устойчивости приближений Приложение Список литературы
Введение
Сплайны в математике имеют достаточно давнюю историю. Так, ломаную Эйлера можно считать простейшей сплайновой аппроксимацией. Дженкинс (W. A. Jenkins) фактически рассматривал сплайны, когда исследовал оскуляторную интерполяцию в 192G г.; изучал ее также и Гревилль (Th.N. E.Grevill) в 1944 г. Слово “сплайн” английского происхождения: так называлась гибкая линейка, применявшаяся английскими инженерами в конце XIX века для проектирования закруглений железных дорог. Математика получила этот термин благодаря работам Шонберга в 1946 г. (см. [54]), который назвал так рассмотренные им функции с “кусочными” свойствами. К настоящему времени имеется большая серия статей и ряд монографий, освещающих многие стороны теоретических исследований и практического применения сплайнов (см. [17, 19, 32, 40, 42, 43, 48] и библиографию в них). Отметим в связи с этим глубокие связи сплайнов с конечно-элементной аппроксимацией [1, 42, 47, 49, 50, 56] и их фундаментальную роль в бурном развитии теории вейвлетов [об].
Аппроксимация и интерполяция широко используются при построении конечно-элементного и сеточного методов приближенного решения задач математической физики, при сжатии и последующем восстановлении с заданным порядком потоков числовой информации, при их фильтрации и статистической обработке. Значительное место среди средств приближения и интерполяции занимают пространства сплайнов. Стремление к разработке экономичных методов приводит к использованию функций с малым носителем (локальных функций — см. [21]) с теми или иными минимальными свойствами. Наиболее часто применяются пространства с
5
базисом из функций, носители которых имеют минимальную кратность перекрытия (обычно именуемую кратностью накрытия). К ним относятся В-сплайны [33, 35, 48], сравнительно недавно разработанные интерполяционные минимальные сплайны [28, 42, 46], а также ряд конечно-элементных аппроксимаций (см., например, [50]). Подобные аппроксимации называются минимальными [22, 26]. Они позволяют значительно снизить трудоемкость вычислений [23].
В монографиях [19, 33, 41] отмечено, что в ряде случаев предпочтение отдается сплайнам с локальным интерполяционным базисом. Отличительная особенность таких сплайнов заключается в том, что решение интерполяционной задачи в точке не зависит от информации о поведении функции в достаточно удаленных узлах сетки и, таким образом, интерполирующая функция строится достаточно просто.
Свойство минимальности рассматриваемых аппроксимаций позволяет экономить ресурсы вычислительного устройства. Например, при решении краевых задач упомянутое свойство приводит к минимизации ширины ленты у матрицы соответствующей конечномерной задачи, а в итоге уменьшается число арифметических операций при отыскании решения.
Свойство минимальности особенно важно при аппроксимации функций с быстро растущими производными, где приходится использовать сильно неравномерную сетку. Характер неравномерности сетки можно определить в соответствии с классом приближаемых функций [25], при этом сохраняются свойства оптимальности приближения (по порядку) и оценка числа арифметических действий.. К рассматриваемому семейству сеток относятся конечные сетки, а также определенные бесконечные сетки с конечными или бесконечными точками сгущения.
6
Сплайн В. С. Рябенького [46] был, по-видимому, первым интерполяционным минимальным сплайном. С помощью сплайнов с локальным интерполяционным базисом, как было отмечено выше, удобно решать краевые задачи вариационно-разностным методом, так как в этом случае в результате решения системы линейных алгебраических уравнений с матрицей ленточного типа мы сразу получим приближенное решение исходной задачи. Преимущества применения базисных фунций с конечным малым носителем при решении задач вариационно-разностным и проекционным методом отмечались в монографиях [41, 42].
Полиномиальные В-сплайны с момента возникновения до настоящего времени привлекают пристальное внимание специалистов, как в плане применения так и с целью развития теории дальнейшего исследования, в частности, получения констант в оценках погрешностей.
В работе [34] построено несложное явное выражение для В-сплайнов произвольной степени на неравномерной сетке. Найдены скалярные произведения этих сплайнов вплоть до 7-ой степени в случае равномерной сетки. Разработан алгоритм нахождения наилучших среднеквадратических аппроксимаций функции такими сплайнами. Предложен тест, пригодный для проверки и сравнения разных классов аппроксимаций. Расчеты на этом тесте показали хорошую обусловленность алгоритма и возможность достижения высокой точности аппроксимации.
В работе [36] предложены простые универсальные алгоритмы нахождения естественных и периодических интерполяционных сплайнов произвольной степени. Расчеты показали высокую устойчивость алгоритма. Тригонометрические сплайны максимальной гладкости первого порядка предложены в работе И.Г.Буровой [7].
7
В данной работе проведено исследование свойств тригонометрических сплайнов первого порядка и предложены тригонометрические локальные сплайны максимальной гладкости второго и третьего порядков, с носителем, занимающем пять и семь сеточных интервалов. Предложены новые интерполяционные задачи использующие гладкие тригонометрические базисные сплайны, для решения которых не требуется решение систем линейных алгебраических уравнений. Даны оценки погрешности приближений.
Диссертация содержит — 3 главы (13 параграфов) и Приложение. Глава 1 посвящена построению тригонометрических сплайнов максимальной гладкости и рассмотрению связанных с ними интерполяционных задач. В первом параграфе приведены основные моменты построения непрерывных тригонометрических базисных сплайнов. Во втором параграфе — основные результаты построения непрерывно дифференцируемых тригонометрические сплайны первого порядка согласно работе [7]. В третьем и четвертом параграфах строятся гладкие тригонометрические сплайны второго и третьего порядков соответственно. В пятом параграфе для построенных сплайнов предложены интерполяционные задачи, аналогичные сформулированным в работе [29]. В шестом параграфе доказывается, что предлагаемые аппроксимации обладают точностью на соответствующих тригонометрических полиномах. В седьмом параграфе получены оценки погрешностей аппроксимации тригонометрическими сплайнами первого и второго порядков.
В главе 2 рассматриваются полиномиальные В-сплайны второй, четвертой и шестой степеней, а также некоторые новые интерполяционные задачи. В первом параграфе приведены формулы указанных сплайнов и ап-проксимационные соотношения, из которых они могут быть получены. Во
8
втором параграфе для полиномиальных В-сплайнов поставлены интерполяционные задачи, аналогичные тригонометрическому случаю. В третьем параграфе получены оценки погрешностей для решений интерполяционных задач, поставленных в параграфе 2.
Глава 3 посвящена рассмотрению прикладных моментов. В первом параграфе рассматриваются аппроксимации функций, аналогичные приведенным в [7], сохраняющие гладкость приближения, но не использующие значения производных функции в узлах сетки. Получены оценки аппроксимации, использующей величины функции в узлах, и построены новые гладкие базисные функции. Во втором параграфе строятся квадратурные формулы, согласованные с гладкими тригонометрическими сплайнами первого порядка и полиномиальными В-сплайнам второй степени. В третьем параграфе рассмотрены вопросы устойчивости приближений непрерывно дифференцируемыми полиномиальными и тригонометрическими сплайнами. Для сравнения приведены аналогичные результаты для В-сплайнов второй степени [33]. Результаты численных экспериментов, подтверждающие тероретические оценки погрешностей, приведены в Приложении.
9
1 О минимальных тригонометрических сплайнах
В этом параграфе напомним основные моменты построения непрерывных тригонометрических базисных сплайнов.
Пусть функция / € С(М.]) задана в узлах равномерной сетки {#/}, х= = 3 — О» • • •) ^ ^ 0, • • • ^ ^'-1 ^ Ху Ху+\ «С ... .
Решение интерполяционной задачи Лагранжа с помощью непрерыв-
ных минимальных интерполяционных тригонометрических сплайнов на [xj1 #;+1) (см., например, [8]) имеет вид
]+М-1
/(х)= 5Г (1.1л)
В случае М2 > 1, М\ > О, М\ 4- М2 = 2т, тбМ, базисные сплайны (1.1.2) являются непрерывными функциями. Нетрудно видеть, что иу{хк) = $зк, где бук — символ Кронекера. Данный базисный сплайн при М\ = М2 = 1 и М\ = М2 = 2, 7 = 2, Л = 1, иллюстрируют рисунки 1 и 2 соответственно.
1.1 Построение непрерывных тригонометрических сплайнов
к = у- М2,... ,3 +
О, х [xj—м21Xj+м^+l]•
о
5 хРис. 1.
10
Рис. 2.
Рассматриваемые базисные сплайны на промежутке [х;*, 2^+1) могут быть получены как решение линейной алгебраической системы уравнений
>+л/2
Wy(x) cos(k,Xj>) = cos(&2;), к = 0,1,... , m,
j+Mj
^ а;;/(т)sin(fcx<;/) = sin(A;x), A; = 1,2,... ,m.
Определитель этой системы
(1.1.3)
VPcs {Xj—Mi i • • • 1 ) —
1 1
sin Xj-Mi sin Xj+\fy
COS Xj—Mj .. COS Xj+A/3
sin(mxj-M,) ■■ sin(m2:;4A/2)
cos {mxj-Mi) cos(mxJ+A/2)
есть обобщение определителя Вандермонда. Известна формула для его вычисления [20]:
Wa(xj.Ml,... , я;+д,2) = 22ш’ П sin Xq 2 Жр~
>-A/i <p<q<j+Mj
Очевидно, что при хр — хдф 2пк, к £ Z, матрица системы невырождена и поэтому решение системы уравнений единственно.
Аппроксимация
j+Мг
f(x) = /(**)"*(*)» xe[Xj,Xj+1)
k=j-Mx
11
с приведенными выше тригонометрическими базисными сплайнами обладает свойством f(x) — f(x) = 0, если функция f(x) равна 1, sinx, cosx,... , sin rax, cos mx.
В дальнейшем, для краткости, при выполнении данного свойства будем говорить, что аппроксимация точна на тригонометрических полиномах не выше ш-го порядка.
1.2 Построение непрерывно дифференцируемых
тригонометрических сплайнов первого порядка
Непрерывно дифференцируемые тригонометрические сплайны первого порядка были построены в работе [7]. Для полноты изложения приведем основные результаты.
Пусть т=1. Система (1.1.3) определяет непрерывные тригонометрические базисные сплайны. Рассмотрим теперь систему уравнений более сложного вида относительно а;*(х) с параметрами соь С1(Ь ^02> С20> считая,
что X е [Xj, Xj+i), suppwj(ar) = [xj_b xj+2]:
'
Wj-i(z) + Uj(x) + Vj+l(x) = 1,
< sin(x;_i)u/;_i(x) + sin(x^)cj;(x) + sin(xj+i)u;;+i(x) = сю sin(x) + c0i cos(x), COs(Xj-\)<Jj-i(x) + C0s(xj)cjj(x) + COS(Xj+i)a^'+i(x) = C02 sin(x) 4* C20 cos(x).
(1.2.1)
На промежутке [xj-i, xj) для определения базисного сплайна имеем сис-
12
тему
4-^_і(х) + и;>(х) = 1,
8Іп(ху_2)и;,_2(х) +5Іп(аг;_і)а7у-і(х) 4-8Іп(х^)и;,(х) = сю 8Іп(х) 4- с0і соб(х), С05(х^_2)и,';-2(т) +соз(а:;_1)а;;_1(х) + сов^иДд;) = с0 2 віп(а:) 4-с20со8(х),
(1.2.2)
а на промежутке [хя;+і) — и>,(х) 4-с^+і (х) +^у+2 (аг) = 1,
8Іп(х;)а;;(х) + 8Іп(ху+і)^-+і(х) 4- 8Іп(ху+2)о^+2(х) = сю 8Іп(х) + с0і сов(х), СОЬ{хі)иі(х) 4- С08(х,+і)а;;-+і(х) 4- С08(х^+2)а^+2(х) = Со28Іп(х) 4- с20 С08(х).
(1.2.3)
Значения параметров сої, сю, со2, с2о находим из условия о>;(х) Є С^К1). При со2 = —сої = сов(/г/2) зіп(Л/2), сю = с2о = со82(Л/2), на промежутке [х^-,х;+і) из (1.2.1) получаем
сов(Л) — сов(х — ЗЬ, — Л./2) соз(/і/2)
а/,-(х) =
соб(Н) — 1
, ч сов(х — — Л) — 1 . . сов(х — У/г) — 1
= 2(соз(А)-1) » ^1(1) = ВДЛ) -Т
Нетрудно видеть, что формула базисного сплайна имеет вид
' С08(х — рь 4- Ь) — 1
(1.2.4)
щ(х) = <
2(сов(Л)-1) ■’ *єІ*і-ь*Д
соз(Л) — со8(х — у/г — /і/2) соз(/і/2)
С08(/і) — 1
соя(х - ІН — 2/г) — 1 г .
х Є [ху+і,хі+2].
2(соз(/і) - 1)
График этой функции при У = 1, Л = 1 приведен на рис. 3.
13
- Київ+380960830922