CC BY
16
ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
2019 Управление, вычислительная техника и информатика № 47
УДК 519.6
DOI: 10.17223/19988605/47/10
Б.М. Шумилов
МУЛЬТИВЕЙВЛЕТЫ НЕЧЕТНОЙ СТЕПЕНИ, ОРТОГОНАЛЬНЫЕ МНОГОЧЛЕНАМ
Работа выполнена при поддержке фонда РФФИ и Администрации Томской области, грант № 16-41-700400 ра.
Представлена конструкция мультивейвлетов на основе эрмитовых сплайнов произвольной нечетной степени на конечном отрезке. Мультивейвлеты ортогональны многочленам той же степени и имеют наименьший носитель. Предложен алгоритм мультивейвлет преобразования с применением метода блочной матричной прогонки. Исследована устойчивость на больших сетках. Дан численный пример аппроксимации и сжатия данных для случая мультивейвлетов седьмой степени.
Ключевае слова: эрмитовы сплайны; мультивейвлеты; ортогональность многочленам.
Вейвлетом называется короткая или быстро затухающая волновая функция (всплеск), множество сжатий и смещений которой порождает пространство измеримых функций на всей числовой оси [1-3]. Основой для построения вейвлетов является наличие набора аппроксимирующих пространств ... Уь-\ с Уь с Уь+\ ..., таких, что каждая базисная функция в У- может быть выражена в виде линейной комбинации базисных функций в Уь. В частности, таким свойством обладают сплайны - гладкие функции, склеенные из кусков многочленов степени т, на вложенной последовательности сеток. Если порядок склейки равен т - 1, то классические полуортогональные вейвлеты (элементы пространства Уь, ортогональные пространству Уь-0 имеют довольно большой носитель из 2т + 1 шагов сетки. В случае аппроксимации на конечном интервале ситуация усугубляется, так как приходится учитывать краевые эффекты. Это препятствует их широкому использованию для решения задач математического моделирования. В отличие от этого эрмитовы сплайны нечетной степени т = 2г + 1 (соответствующие склейке порядка г) приводят к вейвлетам с носителем из трех шагов сетки, что, несомненно, предпочтительнее. Поскольку в базисе таких функций несколько, они называются мультивейвлетами [4-8]. В [9, 10] были предложены методы построения вейвлетов, ортогональных многочленам, с уменьшенными носителями. Идея уменьшения носителя вейвлета за счет замены свойства ортогональности пространству сплайнов на прореженной сетке ортогональностью многочленам представляется привлекательной. Действительно, с точки зрения скорости приближения гладких функций [3] данные типы вейвлетов эквивалентны, а ортогональность многочленам обеспечивает локально максимальную «похожесть» на наилучшее среднеквадратическое приближение.
В данной работе мы излагаем результаты исследования мультивейвлетов произвольной нечетной степени с носителем [0, 2], ортогональных многочленам той же степени, на конечном отрезке.
1. Построение системы базисных мультивейвлетов произвольной нечетной степени на конечном отрезке
Пусть на отрезке [а, Ь] задана вложенная последовательность равномерных сеток Дь: х, = а + ¡•И, / = 0, 1, ..., 2ь, И = (Ь - а)/2ь, Ь > 0. Мультивейвлеты возникают тогда, когда с каждым узлом сетки связано более одной базисной функции [4]. В частности, если базисные функции ~№Ьг,к(х) = фк(у - ¡), к = 0, 1, ..., г V/, где V = (х - а)/И, с центрами в целых числах, порождены сжатиями и сдвигами г + 1 функций вида [11. С. 82]:
Ф* (г) =
(-1)* ю* (-г) при -1 < г < 0, ю* (г) при 0 < г < 1, 0, г г [-1,1],
г+у г-к (г + В)'
где со^ (V) = (1 — £ )г+ X -— , к = 0,1.....г, то при условии отсечения выступающих за концы от-
р=о £!|3!г!
резка половинок функций ф0(0, ..., фг(0 полученное пространство Уь является пространством эрмитовых сплайнов степени 2г + 1 гладкости Сг.
На любой сетке Дь, Ь > 0, интерполяционный эрмитов сплайн (2г + 1)-й степени может быть представлен как
г 21
Ъ1 (X) = £ £ с^ к;?,* (х), а < х < Ь, (1)
*=01 =0
где коэффициенты С/ь,к, к = 0, ..., г, являются значениями и соответствующими производными аппроксимируемой функции/(х), помноженными на к в соответствующей степени: {/к)(/к)кк, к = 0, ..., г, / = 0, ..., 2Ь}, в узлах сетки.
Поскольку сетка Дь-1, Ь > 1, получена из Дь посредством удаления каждого второго узла, то соответствующее пространство УЬ—1 с базисными функциями ^-1г,к, с носителями, в два раза большими по ширине, и центрами в четных целых числах, вложено в Уь. Остаток Жь-1 от разности пространств Уь и Уь-1, размерности (г + 1)(2ь + 1 - 2ь-1 - 1) = (г + 1)2ь-1, называется пространством вейвлетов.
В классической теории вейвлетов базисные функции ортогональны всем базисным сплайнам на прореженной сетке Уь-1 по отношению к некоторому скалярному произведению. Это означает, что пространство Уь представляет прямую сумму Уь-1 и Жь-1: Уь = Уь-1 © Wь-l. В отличие от этого, будем искать вейвлеты Мьг,к(х), к = 0, 1, ..., г V/, как линейные комбинации базисных эрмитовых сплайнов на сетке Дь+1, удовлетворяющие условиям ортогональности многочленам (2г + 1)-го порядка [12]
ь
|М^(х)хтЛ = 0, Аг = 0,1,...,/■ У/(т = 0,1,...,2г+1) (2)
а
и имеющие минимальные из возможных носителей.
Теорема. Система функций Мьг,к(х), к = 0, 1, ..., г, / = 1, 2, ..., 2ь, удовлетворяющих условиям (2) с носителями не более чем из двух шагов сетки Дь, существует и образует базис в пространстве Wь. Доказательство.
а) Пусть ь > 1, и носители вейвлетов [х2/-1, Х2г+р+1], р > 1, полностью располагаются внутри отрезка аппроксимации [а, Ъ]. Прир > 1 для всех / = 1, 2, ..., 2ь - 1 - [р/2] имеем разложение
М^ (х) = £ £ а.фг (г - Д г = (х - х2;)/И, -1 < г < р +1. (3)
I=0 .=0
Разложение (3) подставим в (2) и вычислим все необходимые интегралы, учитывая при этом, что подынтегральные выражения обращаются в нуль вне носителей вейвлетов. Условия ортогональности (2) определяют однородную систему 2г + 2 уравнений относительно коэффициентов вейвлета
М^*(х). В силу линейной независимости на интервалах [х2/-1, Х2г+р+1] базисных функций N2+/*(х),
I = 0, ..., г,у = 0, 1, ..., р, и поверочных функций хт, т = 0, 1, ..., 2г + 1 [1], ранг полученной системы равняется шт[2г + 2, (г + 1)(р + 1)]. Если предположить, что носитель вейвлета равен трем шагам сетки Дь+1, т.е. р = 1, то однородная система становится квадратной и, будучи невырожденной, может иметь только тривиальное решение. Поэтому будем далее считать, что носитель вейвлета равен четырем шагам сетки Дь+1, те вейвлет построен из 3г + 3 базисных сплайнов. В этом случае ранг матрицы равен 2г + 2. Поэтому однородная система 2г + 2 уравнений с 3г + 3 неизвестными имеет г + 1 линейно независимых частных решений. Длина носителей, полученных для каждого номера / вейвлетов, равна длине носителей базисных сплайнов на сетке Ль. Их линейная независимость вытекает
из того факта, что они представляют собой множество сдвигов ненулевых функций с компактными носителями. При этом количество вейвлетов, полностью лежащих внутри отрезка аппроксимации, равно (г + 1)(2Ь - 2) = (г + 1)2Ь - 2(г + 1), что на 2(г + 1) меньше разности между размерностями пространств У+1 и У: (г + 1)(2Ь+1 + 1) - (г + 1)(2Ь + 1) = (г + 1)2Ь.
Вблизи концов отрезка интервалы интегрирования не выходят за пределы отрезка. Поэтому с учетом того, что по краям отрезка [а, Ь] от выступающих за концы отрезка функций фо(0, • ••, фг(0 остается по половинке, граничные вейвлеты отличаются от вейвлетов, предложенных выше, при условии ортогональности многочленам (2г + 1)-й степени. Аналогичные рассуждения приводят к двум системам линейно независимых функций ^о(0, • • •, ^>(0 отдельно на левом и на правом концах, дополняющих до базиса построенную выше вейвлет-систему.
б) Пусть Ь = 0. Вычисления снова дают г + 1 линейно независимых частных решений, которые дают базис в пространстве У - У о с размерностью (г + 1)3 - (г + 1)2 = (г + 1). Теорема доказана.
В соответствии с правилами линейной алгебры сформулируем следующий алгоритм построения эрмитовых мультивейвлетов нечетной степени.
Следствие. Введем в рассмотрение три матрицы размерности (2г + 2) х (г + 1) с элементами У+1
| 7=0,1,2,/ = 0,1,...,Г, щ = 0,1,...,2г + 1
У-1
и составим из них блочную квадратную матрицу (2г + 2)-го порядка вида [Л0 | Л2]. Тогда при условии, что коэффициенты разложения (3) а! = {1,1 = к; 0,1 Ф к], каждый из г + 1 столбцов матрицы
дшпег у^дшпег = _[/<>0 |] 1/<>' дает значения коэффициентов а'-. / = 0.2. / = 0.1... ..г. соответствующего к-го базисного вейвлета, полностью лежащего внутри отрезка аппроксимации. На левом конце отрезка элементы матрицы Л0 вычисляются по укороченному интервалу:
1
(20^, I = 0,1, • • •, г, т = 0,1,..., 2г +1,
о
а коэффициенты разложения (3) а', у = 1,2, / = 0,1,...,г, соответствующего к-го базисного вейвлета даются значениями столбцов матрицы Га1^/а26Й ] = -[Я11 Я2 ]-1 я0 при условии, что коэффициенты
а.
= {1,1 = к; 0,1 Фк]. Аналогично на правом конце отрезка аппроксимации элементы матрицы Я2 вы-
числяются по укороченному интервалу:
2
= IФ/ (2^ - 2утЛ, I = 0,1,..., г, т = 0,1,..., 2г +1, 1
а коэффициенты разложения (3) а', у = 0,1, / = 0,1,.. .,г, соответствующего к-го базисного вейвлета даются значениями столбцов матрицы Ад8*11/А"®111 = -[Я0 | Я1]-1 Я2 при условии, что коэффициенты а2 = {1,1 = к; 0,1 Ф к] . При Ь = 0 все интегралы вычисляются по интервалу [0, 2], а граничных вейвлетов не возникает: [ А^ /А^ ] = -[Я0 | Я2 ]-1 Я1 при условии, что а1 = {1,1 = к; 0,1 Ф к}.
2. Построение и обращение блока фильтров
Если упорядочить базисные сплайн-функции в виде единой матрицы-строки,
Ф- =
0, N1,1,.", г г то можно записать функции ф в виде линейных ком-
бинаций функций фь: фЬ-1 = фЬРЬ. Здесь блоки матрицы Рь составлены из коэффициентов масштабных
(калибровочных) соотношений [13]:
Ф0 ) Ф0 — -к)'
Фх(0 = ЪНк -к)
к=0
_ф г(0. _фД2г- -к)_
где H2 = Ц~1ЛЦ, Л = (2
—г—1 <-\—2г-
1матрица С/ размерности (г +1) х (г +1) задана элементами
Ц
к' ]
= = ^аЦ^-1,...^), 5=^(1,-1,...,(-1Г), я(
0 =
Аналогично запишем базисные вейвлет-функции (2г + 1)-й степени на уровне масштабирования Ь
в виде матрицы-строки: у =
,1 '1,0'-
1,1
Г1
1 г' •••' -
£
. Тогда для уровня Ь - 1 можно выразить
функции ^ 1 в виде линейных комбинаций функций фЬ, т.е. 1 = ф1^, где блоки матрицы составлены из коэффициентов разложения (3). Соответствующие коэффициенты сплайна будем соби-
т
рать в вектор С =
Гь,0 Гь,1 С0 ' С0 '
• • •> Ч) > м
,Ь,0
С
а соответствующие вейвлет-коэффициенты -
в вектор О =
О 'О1 '
.А
.О
(здесь символ Т обозначает операцию транспонирова-
ния). С использованием обозначений для блочных матриц процесс получения СЬ из СЬ 1 и О 1 может быть записан как [14]
С1 = [ Рь
Я1 ]
С
■Ь—1
Ниже представлены примеры матриц [Р 12
[ Р1121] =
нх
н Т
о
О А ее^ег А0
Н о I
Нх Ас®п1ег
,[ Р2
122] =
1
соответствующих Ь = 1, 2, 3:
н
Н Т2 о
о о
(4)
[ Р3!б3] =
Н\ О О О О I
н Т Н Т О О О А|ей
о нх О О О АlГft
О н Т н Т О О О
О О Нх О О О
О О н Т н Т О О
О О О нх О О
о О О н Т н Т О
О О О О нх О
О О I О
н Т О аГ О
нх О А*,6* А^1
нТ нТ О
О нх О I
О О О
О О О
л тпег А0 О О
I О О
дтпег л тпег А0 О
О I О
О дтпег Ад8Ь1
О О
О О I
Здесь О обозначает матрицу (г + 1)-го порядка с нулевыми коэффициентами, тогда как I - единичная матрица (г + 1)-го порядка. Обратный процесс разбиения коэффициентов СЬ на более грубую версию и уточняющие коэффициенты О 1 состоит в решении системы линейных уравнений (4). Разрешимость данной системы гарантирована линейной независимостью базисных функций.
Процедуру разбиения СЬ на часть СЬ-1, соответствующую низшему масштабу, и уточняющие коэффициенты О1-1 можно применить рекурсивно и к самой этой части СЬ-1. Следовательно, исход-
т
ные значения можно представить в виде иерархии грубых версий с масштабами С0, С1, •.., СЬ-1 и уточняющих деталей D0, D1, • .., DL-1. Результирующее вейвлет-разложение эрмитового сплайна (2г + 1)-й степени SL(x) может быть записано в виде:
Ь-1 2]
(х) = £ £С0,кК0,к(х)+ЕХД],км]к(х) к=01 1=0 ]=0 1=0
, а < х < Ь.
При этом по величине вейвлет-коэффициентов D, ] = 0, 1, — , Ь - 1, можно судить о значимости соответствующих уточняющих деталей. Незначимые убираются с целью сжатия информации.
3. Проблема устойчивости мультивейвлет-преобразования
При больших Ь матрицу [Рь | следует сделать блочной трехдиагональной, изменив порядок неизвестных так, чтобы блоки матрици (У перемежались (см. [15]):
Нх I 0 ••• I
н1 41ей н1 о 4ей ■•. ■•. (
[р1\д1 ] =
А^ Я0Г
0
0
I
н,
£ > 1.
1
Тогда можно применить для решения системы (4) алгоритм блочной матричной прогонки [16]. Введем для Ь > 1 обозначения Ь = Ь1оскп^{И2Т, 0, 0; Л21ей, 0, 0; И2Т, 0, 0; Л21Ш1ег, 0, 0; •.. ; I, 0, 0}; и = Ь1оскиМав{0, 0, I; 0, 0, Нт; 0, 0,Л01тег; • .. ; 0, 0,Л0ПёМ; 0, 0, Нт}; Т = Ь1оска1а8{Тг, I = 0, 1,— , 2Ь}. В результате получим
[Рь | ^] = (ь + Т)Т-1(и + Т), (5)
где для блоков Т справедливы выражения [17]
•/, =Я,:/2 =4ей-Нт2Т~1-,Т3 =Я,-А^Т-'Н1,; Тг =1-Нт2Т~\АГес-,Т1+1 =Н,(г =4,6,...^ -2);
т21 т21+1 =нх-Т-Ы1.
Аналогично для Ь = 1 имеют место Т = Н; Т2 = I - НТ-1-Т = Н - А^^Т-нГ.
Если ввести обозначения Т^СХЛ-Т)^1, ..С^] Т = г, С? = С^1,..., С, ])[' = /)/"',../.)/"''то процесс решения системы (4) разбивается на два этапа:
(6)
•Ь,г
(Ь + Т) 2 = &; (I + Т-1и)
П-1¥
гЬ-1 т^Ь-1
п£-1 ГЬ-1
= 2.
Покомпонентно каждый из них можно записать в следующем виде: для решения (ь + Т) 2 = С1 последовательно выполняются действия
2 = Т
1 = Г1-1 С ) ;22 = Г2-1(С ) - НТ21); 23 = Тз-1((С2^ ) - А^Ч);
= ТТ\(с1$ -Я2Т2г_1);2г+1 )Т -4ППеЧ) (/ = 4,6,...,2* - 2);
Т2^ ((С21 -1 ) - Н2 22х-1); 22ь +1 = Т2^1+1 ((С2Х ) - 22х ),
(7)
22ь =
а для вычисления (I + Т"1и)[с^"1, ..1)'^], формулами:
= 2 можно воспользоваться рекуррентными
Г
Г
(С2Х —1 ) = V +1'(°21 -1 ) = 221 Т2 'н0Т (С2Х —1 ) ' (с^"1 )Т = - т;:114ппег )Т (г = 2^ -1,2Ь~1 - 2,...,2);
(А-1)1 =*2 -г2-чт = ~Г1-1 (А"1)1 •
Здесь Zi суть подвекторы порядка г + 1.
Нетрудно убедиться, что при Ь > 1, кроме первой и последней строк, условия даже слабого диагонального преобладания [16] не выполняются. Объективное представление об устойчивости или неустойчивости алгоритма вычисления мультивейвлет-преобразования дает наблюдение за поведением чисел обусловленности в эвклидовой норме прогоночных матриц Т, (таблица).
Значения стандартной функции MathCad'а 00^2(7;), ; = 1, 2, ..., 33, г = 1, ..., 5
1 2 3 4 5 6
1 \ еопё2(Г,) г = 1 г = 2 г = 3 г = 4 г = 5
1 2 4 8 16 32
2 64,759 2,695е+004 2,799е+007 2,142е+011 8,333е+014
3 2,513 157,2 5,832е+004 8,799е+007 3,717е+010
4 2,953 60,32 7 805 2,468е+006 8,122е+008
5 3,493 478,7 1,817е+005 7,345е+007 7,162е+010
6 2,835 52,05 7 124 1,999е+006 7,353е+008
7 3,56 492,9 1,864е+005 7,264е+007 8,126е+010
8 2,831 51,79 7 102 1,991е+006 7,328е+008
9 3,563 493,4 1,865е+005 7,263е+007 8,158е+010
10 2,831 51,79 7 102 1,991е+006 7,327е+008
11 3,563 493,4 1,865е+005 7,263е+007 8,159е+010
32 533,592 3,61е+005 1,379е+009 5,382е+011 3,713е+017
33 84,003 2,933е+004 1,965е+007 7,101е+009 1,308е+015
Таким образом, можно полагать, что для любого уровня разрешения Ь они являются невырожденными, что на практике означает корректность представленного выше монотонного варианта алгоритма матричной прогонки (6)-(8). Для г = 1 (случай приближения кубическими мультивейвлетами) это согласуется с тем фактом, что на бесконечной числовой оси Ь2-устойчивость базиса гарантирована существованием оценок констант Ритца 0,414076 < А < В < 0,414265 [10], таких что
2
АЕ ЕЕ
к=0 / I
О'
<
Е ЕЕ О/,к М/'к (х)
к=0 / I
О'
<вЕ ЕЕ
¿2(Я) к=0 / 1
В случае приближения на конечном интервале обусловленность мультивейвлет-базиса может быть сделана весьма близкой к обусловленности базиса на бесконечной оси наложением дополнительно однородных краевых условий Дирихле [18, 19]. Однако, как следует из таблицы, с повышением степени мультивейвлета численная устойчивость будет неизбежно ухудшаться.
г
4. Пример для г = 3 (случай мультивейвлетов седьмой степени)
Мультивейвлеты седьмой степени можно применять в тех случаях, когда требуется аппроксимировать функцию или ее производные до пятого порядка с более высокой точностью по сравнению с вейвлетами степеней 3, 5. Кроме того, с помощью вейвлетов седьмой степени можно аппроксимировать на промежутках между узлами х, е ДЬ производные шестого и седьмого порядков - для вейвлетов степеней 3, 5 они равны нулю. Это может быть использовано при решении дифференциальных уравнений высокого порядка, например методом коллокации.
В частности, при Ь > 1 внутри отрезка аппроксимации получаем (все вычисления проводились в системе MathCad, дробные числа там, где они затрудняют восприятие результата, приводятся округленно в виде конечной десятичной дроби)
1,1168 0,082 0,0159 0,0004" —7,894 —0,2968 — 0,0915 — 0,001 —67,9076 —6,9996 —1,1658 — 0,039 ' 1792,4185 129,5384 26,3043 0,6684
На левом конце отрезка аппроксимации имеем
15,9089 1,2448 0,0667 0,0019 —12,3588 —0,7126 — 0,0257 — 0,0004 —1984,0775 —158,6934 —8,679 — 0,2531 ' 30335,1818 2401,4609 130,1543 3,7654
На правом конце отрезка аппроксимации значения совпадают с точностью до знака и перестановки.
При этом н1 = ^ (8,4,2,1)/8,
f 384 840 0 -5040> f384 - 840 0 5040^
1 -132 -228 360 2520 1 132 -228 -360 2520
> H 2 =-
768 18 24 -84 -360 2 768 18 -24 -84 360
v -1 -1 6 18 V v 1 -1 -6 18 V
Для xe[0, 1], полагая на верхнем уровне разрешения L = 5, находим длину шага сетки h = 2-5=1/32. В соотношениях уточнения (4) при выполнении вейвлет-преобразования в качестве исходных используются значения функции и трех производных, всего 132 числа. Рассматривая в качестве тестовой функции многочлен седьмой степени fx) = (2 - x)3(x2 - 1)2, находим на последнем этапе рекуррентного алгоритма вейвлет-разложения восемь значений сплайна <S°(x) и трех его производных в концах отрезка C0 = [8, -12, -20, 138, 0, 0, 8, -48]T. При этом все вейвлет-коэффициенты пренебрежимо малы и могут быть обнулены, обеспечивая в данном случае коэффициент сжатия K = 132/8 = 16,5.
Заключение
В работе представлена общая схема построения и вычисления мультивейвлетов на основе эрмитовых сплайнов и свойства ортогональности многочленам. Исследована устойчивость и даны результаты численного эксперимента по восстановлению функции, заданной многочленом. Дополнительные возможности для оптимизации методов обработки численной информации состоят в построении сглаживающей сплайн-аппроксимации, альтернативной методу наименьших квадратов, и в сжатии дискретных данных на участках гладкости функций.
ЛИТЕРАТУРА
1. Добеши И. Десять лекций по вейвлетам : пер. с англ. Ижевск : Регулярная и хаотическая динамика, 2001. 332 с.
2. Чуи Ч. Введение в вейвлеты : пер. с англ. М. : Мир, 2001. 412 с.
3. Новиков И.Я., Протасов В.Ю., Скопина М.А. Теория всплесков. М. : Физматлит, 2006. 616 с.
4. Strela V. Multiwavelets: regularity, orthogonality and symmetry via two-scale similarity transform // Stud. Appl. Math. 1997.
V. 98, No. 4. P. 335-354.
5. Strela V., Heller P.N., Strang G., Topivala P., Heil C. The application of multiwavelet filterbanks to image processing // IEEE
Trans. Signal Processing. 1999. V. 8, No. 4. P. 548-563.
6. Warming R., Beam R. Discrete multiresolution analysis using Hermite interpolation: Biorthogonal multiwavelets // SIAM J. Sci.
Comp. 2000. V. 22, No. 1. P. 269-317.
7. Dahmen W., Han B., Jia R.-Q., Kunoth A. Biorthogonal multiwavelets on the interval: cubic Hermite splines // Constr. Approx.
2000. V. 16. P. 221-259.
A i
1,1168 -0,082 0,0159 - 0,0004 7,894 -0,2968 0,0915 - 0,001 -67,9076 6,9996 -1,1658 0,039 -1792,4185 129,5384 -26,3043 0,6684
Af =
-27,3194 -2,1593 - 0,115 - 0,0032 -202,3878 -18,5595 -1,1415 -0,0368 2442,3182 192,6016 10,2045 0,2832 61774,985 5405,3691 322,4923 10,2122
Af =
8. Han B. Approximation properties and construction of Hermite interpolants and biorthogonal multiwavelets // J. Approxim. Theory.
2001. V. 110. P. 18-53.
9. Koro K., Abe K. Non-orthogonal spline wavelets for boundary element analysis // Engineering Analysis with Boundary Elements.
2001. V. 25. P. 149-164.
10. Han B., Kwon S.-G., Park S.S. Riesz multiwavelet bases // Appl. Comput. Harmon. Anal. 2006. V. 20. P. 161-183.
11. Завьялов Ю.С., Квасов Б.И., Мирошниченко В. Л. Методы сплайн-функций. М. : Наука, 1980. 352 с.
12. Шумилов Б.М. Мультивейвлеты эрмитовых сплайнов третьей степени, ортогональные кубическим многочленам // Математическое моделирование. 2013. № 4. С. 17-28.
13. Strang G., Strela V. Short wavelets and matrix dilation equations // IEEE Trans. Signal Processing. 1995. V. 43, No. 1. P. 108-115.
14. Столниц Э., ДеРоуз Т., Салезин Д. Вейвлеты в компьютерной графике : пер. с англ. Ижевск : Регулярная и хаотическая динамика, 2002. 272 с.
15. Шумилов Б.М. Алгоритм матричной прогонки вычисления мультивейвлетов нечетной степени, ортогональных многочленам // Автометрия. 2015. Т. 51, № 2. С. 83-92.
16. Самарский А.А., Николаев Е.С. Методы решения сеточных уравнений. М. : Наука, 1978. 591 с.
17. Васильева Е.А. Достаточные условия корректности метода матричной прогонки // Известия КГТУ. 2011. № 20. С. 103-108.
18. Schneider A. Biorthogonal cubic Hermite spline multiwavelets on the interval with complementary boundary conditions // Results in Mathematics. 2009. V. 53, is. 3. P. 407-416.
19. Cerna D., Finek V. Construction of optimally conditioned cubic spline wavelets on the interval // Adv. Comput. Math. 2011. V. 34. P. 519-552.
Поступила в редакцию 1 ноября 2018 г.
Shumilov B.M. (2019) MULTIWAVELETS OF ODD DEGREE, ORTHOGONAL TO POLYNOMIALS. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie vychislitelnaja tehnika i informatika [Tomsk State University Journal of Control and Computer Science]. 47. pp. 84-92
DOI: 10.17223/19988605/47/10
r 2L
For the space of Hermitian splines of 2r + 1-st degree of a kind SL (x) = £ hk £ CL'k NLk (x), a < x < b, where the coefficients
k=0 i=0 '
C,L,k, k = 0, ..., r, are values and corresponding derivatives of the approximated function in the knots of uniform net AL: ui = a + (b - a) i / 2L, i = 0, 1,., 2L, L > 0, and basic functions are NLk (u}-) = 5/5k,l = 0,1,...,r , with the centers in integers, it is offered to use as wavelets the functions MLi,k (x), k = 0, 1, r Vi, with the centers in odd integers, the linear combinations of basic Hermitian splines with the grid AL+1, that are orthogonal to polynomials of2r+ 1-st order \baM^k(x)xmck = 0, £ = 0,l,...,rVz (от = 0,1,. ..,2r + 1) .
If the corresponding spline-coefficients are collected in the vector, CL = [Cq'0, Cq д,..., Cq'r, C{"°,..., C^f] , and the corresponding wavelet-coefficients - in the vector, DL ^D^'0, D^'1,..., ,..., ' ^^ use °f designations for block matrices the
formulas for evaluation of spline-coefficients CL-1 on the thinned grid AL-1 and wavelet-coefficients D-1 in the form of the solution of sparse system of linear algebraic equations are proved:
[PL | Q ]
CL-1
= CL.
Here blocks of the matrix PL are composed from coefficients of the scale relations for basic splines and blocks of the matrix Q are composed from coefficients of the decomposition for basic wavelets MLi,k (x). For the purpose of using of the rarefied structure of a matrix [PL | QL] there is offered to make it block tri-diagonal, having changed an order of unknowns so that blocks of matrixes PL and QL been alternated, to be able to apply an algorithm of a block matrix sweeping to the solution of the received system. The problem of stability of algorithm of multiwavelets-transformation on big grids by means of observation of behavior of condition numbers in Euclidean norm of sweeping matrixes is investigated. The numerical example of approximation and compression of data for a case of Hermitian splines of the 7th degree is given.
Keywords: Hermitian splines; multiwavelets; orthogonal to polynomials.
SHUMILOV Boris Mikhailovich (Doctor of Physical and Mathematical Sciences, Professor, Tomsk State University of Architecture and Building, Tomsk, Russian Federation). E-mail: sbm05@yandex.ru
REFERENCES
1. Daubechies, I. (2001 ) Desyat' lektsiy po veyvletam [Ten lectures on wavelets]. Translated from English. Izhevsk: RKhD.
2. Chui, C.K. (2001) Vvedenie v veyvlety [An introduction to wavelets]. Translated from English. Moscow: Mir.
3. Novikov, I.Ya., Protasov, V.Yu. & Skopina, M.A. (2006) Teoriya vspleskov [Wavelet theory]. Moscow: Fizmat.
4. Strela, V. (1997) Multiwavelets: regularity, orthogonality and symmetry via two-scale similarity transform. Studies of Applied
Mathematics. 98(4). pp. 335-354.
5. Strela, V., Heller, P.N., Strang, G., Topivala, P. & Heil, C. (1999) The application of multiwavelet filterbanks to image processing.
IEEE Trans. Signal Processing. 8(4). pp. 548-563. DOI: 10.1109/83.753742
6. Warming, R. & Beam, R. (2000) Discrete multiresolution analysis using Hermite interpolation: Biorthogonal multiwavelets. SIAM
J. Sci. Comp. 22(1). pp. 269-317. DOI: 10.1137/S1064827597315236
7. Dahmen, W., Han, B., Jia, R.-Q. & Kunoth, A. (2000) Biorthogonal multiwavelets on the interval: cubic Hermite splines.
Constructive Approximation. 16. pp. 221-259. DOI: 10.1007/s003659910010
8. Han, B. (2001) Approximation properties and construction of Hermite interpolants and biorthogonal multiwavelets. Journal
of Approximation Theory. 110. pp. 18-53. DOI: 10.1006/jath.2000.3545
9. Koro, K. & Abe, K. (2001) Non-orthogonal spline wavelets for boundary element analysis. Engineering Analysis with Boundary
Elements. 25. pp. 149-164. DOI: 10.1016/S0955-7997(01)00036-4
10. Han, B., Kwon, S.-G. & Park, S.S. (2006) Riesz multiwavelet bases. Appl. Comput. Harmon. Anal. 20. pp. 161-183. DOI: 10.1016/j.acha.2005.10.002
11. Zavyalov, Yu.S., Kvasov, B.I. & Miroshnichenko, V.L. (1980) Metody splayn-funktsiy [Spline function methods]. Moscow: Nauka.
12. Shumilov, B.M. (2013) Multiwavelets of the third-degree hermitian splines orthogonal to cubic polynomials. Mathematical Models and Computer Simulations. 5(6). pp. 17-28. DOI: 10.1134/S2070048213060100
13. Strang, G. & Strela, V. (1995) Short wavelets and matrix dilation equations. IEEE Trans. Signal Processing. 43(1). pp. 108115. DOI: 10.1109/78.365291.
14. Stollnitz, E.J., DeRose, T.D. & Salesin, D.H. (2002) Veyvlety v komp'yuternoy grafike [Wavelets for Computer Graphics: Theory and Applications]. Translated from English. Izhevsk : Regulyarnaya i khaoticheskaya dinamika.
15. Shumilov, B.M. (2015) Matrix sweep algorithm for computing multiwavelets of odd order orthogonal to polynomials. Optoelectronics, Instrumentation and Data Processing. 51(2). pp. 174-183. DOI: 10.3103/S8756699015020119
16. Samarsky, A.A. & Nikolaev, E.S. (1978) Metody resheniya setochnykh uravneniy [Methods of solving grid equations]. Moscow: Nauka.
17. Vasilieva, E.A. (2011) Dostatochnye usloviya korrektnosti metoda matrichnoy progonki [Ample conditions of existence of the complete block decomposition]. IzvestiyaKGTU. 20. pp. 103-108.
18. Schneider, A. (2006) Biorthogonal cubic Hermite spline multiwavelets on the interval with complementary boundary conditions. Results in Mathematics. 53(3). pp. 407-416. DOI: 10.1007/s00025-008-0352-y
19. Cerna, D. & Finék, V. (2011) Construction of optimally conditioned cubic spline wavelets on the interval. Advances in Computational Mathematics. 34. pp. 519-552. DOI: 10.1007/s10444-010-9152-5
