CC BY
34
ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
2012 Управление, вычислительная техника и информатика № 4(21)
УДК 519.6
Б.М. Шумилов, А.Ж. Кудуев
НОВЫЙ ТИП МУЛЬТИВЕЙВЛЕТОВ ПЯТОЙ СТЕПЕНИ, ОРТОГОНАЛЬНЫХ МНОГОЧЛЕНАМ ПЯТОЙ СТЕПЕНИ1
Для случая эрмитовых сплайнов пятой степени исследован новый метод построения системы базисных сплайн-мультивейлетов, использующий условия ортогональности многочленам пятой степени. Рассмотрены варианты построения и обращения блока фильтров. Приведены результаты численных экспериментов.
Ключевые слова: эрмитовы сплайны пятой степени, мультивейвлеты, ортогональность многочленам.
Вейвлетами называются короткие или быстро затухающие волновые функции (всплески), такие, что множество их двоичных сжатий и двоичнорациональных сдвигов в совокупности образует базис пространства Ь2(К) [1-3]. В случае, если волновых функций несколько, они образуют ряд мультивейвлетов [4-8]. Одно из основных преимуществ мультивейвлетов над единичными вейвлетами состоит в том, что они ортогональны многочленам одной и той же степени при меньшей длине носителя. Например, мультивейвлеты пятой степени имеют носитель [0, 3], который значительно меньше носителя [0, 11] единичных вейвлетов пятой степени. В работе [9] были построены мультивейвлеты пятой степени, вторые производные которых ортогональны вторым производным базисных эрмитовых сплайнов пятой степени. Данные мультивейвлеты имеют носитель [0, 2]. При численном решении дифференциальных уравнений по способу Галеркина это увеличивает разреженность решаемых систем. Недостаток состоит в том, что они ортогональны многочленам не выше первой степени, что ухудшает сжатие числовых данных по сравнению с обычными вейвлетами и мультивейвлетами. В данной работе, используя подход [10], построим базисные мультивейвлеты пятой степени с носителями [0, 2], которые ортогональны любым многочленам пятой степени.
1. Построение системы базисных сплайн-мультивейлетов на конечном отрезке
Исходным для построения вейвлет-преобразования является наличие набора вложенных пространств ...Уь-1сУьсУь+1.... В данном случае пространство У является пространством сплайнов степени 5 гладкости С2 на отрезке [а,Ь] с равномерной сеткой узлов Дь: Хі = а + (Ь - а) і / 2Ь, і = 0, 1,..., 2Ь, Ь > 0, и базисными функциями Мьк ( х) = фк (у - і), к = 0,1, Ті і, где V = 2Ь(х-а)/(Ь-а) + 1, с центрами в целых
числах, порожденными сжатиями и сдвигами трех масштабирующих функций вида [11] (рис. 1)
1 Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта № 12-0190915 мол_снг_нр
9o(t )" t 6 1 5 + 0) 1 1 2 t 6 1 t + ) 1
9i(t) + £ 1 (N СП СО "Г ,0 < t < 1; 2 t 3 1 t + ) ,1 < t < 2
Ф2О1) _ 1 + <N 1 (N 3 1 1 1 + <N 1 (N СО (2 1
фк(t) = о, k = 0,1,2, t g [0,2].
I I I
Рис. 1. Графики масштабирующих функций ф0(Г), Ф1(Г), ф2(Г)
Поскольку сетка Д]-1, ]>1, получается из сетки Д] посредством удаления каждого второго узла, то базисными функциями для У]-1 будут функции ^-1(х),
с носителями в два раза большими по ширине и центрами в четных целых числах. В соответствии с масштабными соотношениями [12]
~Фо(t)" 2 Ф0(2t - k)
Фі (t) = Е Hk Фі (2t - k)
_Ф2(5) _ k=0 Ф2(25 - k)
" 1 15 0 - " " 1 15 0
2 16 1 0 0 2 16
где H 0 = 5 32 - 3 1 , 2 3 8 , H1 = 0 1 2 0 , H2 = 5 32 - 3 1 , 2 - 3
1 1 1 0 0 1 1 1 1
_ 64 64 16 _ 4 _ _ 64 64 16 _
каждую широкую базисную функцию внутри отрезка аппроксимации можно построить из трех, а по краям отрезка из двух троек узких базисных функций. Обратное, очевидно, возможно только приближенно. Остается некоторое пространство WL-1, содержащее уточняющие подробности, которые и позволяют восстановить У] из У]-1: У], = У]-1+ ^]-1. Следующим этапом является определение вейвлетов как базиса WL-1. В классической теории вейвлетов базисные функции WL-1 ортогональны всем базисным сплайнам на прореженной сетке Д]1 по отношению к некоторому скалярному произведению. Это означает, что пространство У] представляет прямую сумму У]-1 и WL_1: У] = У]_1© WL_1.
В отличие от этого будем искать базисные функции пространства WL как линейные комбинации базисных эрмитовых сплайнов на сетке Дм, удовлетворяющие условиям ортогональности всем многочленам шестого порядка, то есть
[Ь(х)хтй?х = 0, к = 0,1,2 V/ (т = 0,1,...,5). (2)
* а 5
Теорема 1. Пусть
М0,о (х) = Мо (V), МОд (х) = Мі (V), М0,2 (х) = М2 (V),
М1,о (х) = ^о (V), М1,1 (х) = М’і (V), Мі,2 (х) = ^2 (V),
MLL,0(х) = Wо(2L -V), MLLд(х) = -V), М^,2(х) = w2(2L -V), L > 1,
М^;0(х) = v0(v + 3/2-/), М^(х) = + 3/2-/), М^,2(х) = у2(v + 3/2-/),
/ = 2,3,...,2і- -1, L > 2,
где V = 2ь(х - а)/(Ь - а), и
Мо (/) = Фо (2/)- 4фо (2/ +1)- 4фо (2/ -1) +
+ 84(ф1 (2/ +1) - ф1 (2/ -1)) - 828 (ф2 (2/ +1) + ф2 (2/ -1)),
М1 (/) = 35ф1 (2/) - 200(ф1 (2/ +1) + ф (2/ -1)) +
8ф0 (2/ +1) - 8ф0 (2/ -1) + 2280(ф2 (2/ +1) - ф2 (2/ -1)),
М2 (/) = 35ф2 (2/)- 251 (ф2 (2/ +1) + ф2 (2/ -1)) - ф0 (2/ +1) --ф0 (2/-1) + 23( (2/ +1) -ф| (2/-1),
^0 (/) = 1296ф0 (2/ +1) + (7990ф0 (2/) - 538ф0 (2/-1) - 41545ф1 (2/) --41140ф1 (2/ -1) - 924040ф2 (2/) + 438040ф2 (2/ -1)),
^ (/) = 432ф1 (2/ +1) + (283ф0 (2/) - 12ф0 (2/ -1) -1161ф1 (2/) -- 1288ф1 (2/ -1) - 32372ф2 (2/) + 13520ф2 (2/ -1)),
^2 (/) = 2160ф2 (2/ +1) + (77ф0 (2/) - 2ф0 (2/ -1) - 245ф1 (2/) -
- 320ф1 (2/ -1) - 8900ф2 (2/) + 3320ф2 (2/-1)), у0(/) = 164ф0 (2/) + (28ф0 (2/ +1) + 28ф0 (2/ -1) + 115ф1 (2/ +1)- 115ф1 (2/ -1) - 13200ф2 (2/)),
(/) = 80ф] (2/) + (-ф0 (2/ +1) + ф0 (2/-1) + 37ф] (2/ +1) +
+ 37ф1 (2/ -1) + 412ф2 (2/ +1) - 412ф2 (2/ -1)),
у2 (/) = 464ф2 (2/) + (-4ф0 (2/) - 11ф1 (2/ +1) + 11ф1 (2/ -1)-
-112ф2 (2/ +1)- 112ф2 (2/ -1)). (3)
Тогда система функций { М]к (х), / = 1, 2,..., 2]; к = 0,1,2|, удовлетворяет условиям (2) и образует базис в пространстве WL, ]>0.
Доказательство. Проверка условий ортогональности (2) выполняется непосредственным вычислением. Для установления базисности функций М],к (х) в случае конечного отрезка достаточно их линейной независимости, вытекающей из того факта, что они представляют собой множество сдвигов ненулевых функций с компактными носителями. При этом их количество в точности равно разности
между размерностями пространств У]+1 и У] : 3 •(+1 +1)- 3 •( +1) = 3 • 2] . Теорема 1 доказана.
Носители полученных функций равны носителям базисных сплайнов на сетке Д] (рис. 2).
Рис. 2. Вид двух симметричных и одного антисимметричного «материнских» вейвлетов у0(Г), ^(Г), Уг(/)
2. Построение и обращение блока фильтров
На любой сетке Д], ]>0, интерполяционный эрмитов сплайн 5-й степени может быть представлен как
^ (х)_Х £с],к Щк (х), а < х < Ь, (4)
к= 0 I=0
где коэффициенты С],к, к = 0, 1, 2, являются значениями и соответственно первыми и вторыми производными аппроксимируемой функции в узлах сетки.
Если записать базисные сплайн-функции в виде единой матрицы-строки
ф] = [ <0, <д, <2, <0, Л^,..., М] ,2 ]
и упорядочить коэффициенты сплайна в виде вектора
с] _Гс],0 С],1 С],2 с],0 с],1 с],2 ] Т
С _ [С0 , С0 , С0 , С1 , С1 ,..., С2] ] .
то уравнение (4) переписывается как 8](х) = ф](х)С]. Аналогично, запишем базисные вейвлет-функции на уровне разрешения ] в виде матрицы-строки
у] _[М]0,М],,М]2,...,М]],2] .
Соответствующие вейвлет-коэффициенты будем собирать в вектор
Б] _ [Б]’0, Б],1, б],2, ..., ] Т .
Тогда с использованием обозначений для блочных матриц процесс получения С из
С-1
и Б-1 может быть записан как [13]
cl = [pl | gL;
cl-1
dl-i
(5)
Ниже представлены примеры матриц Р, О, соответствующих Ь = 1, 3:
Р3 =
" 1 0 0 0 0 0 0 ■ 0 0 0 0
0 1 2 0 0 0 0 0 ■ 0 0 0 0
0 0 1 4 0 0 0 0 ■ 0 0 0 0
1 2 5 32 1 64 1 2 5 32 1 64 0 ■ 0 0 0 0
15 16 7 32 1 64 15 16 7 32 1 64 0 ■ 0 0 0 0
0 3 8 1 16 0 3 8 1 16 0 ■ 0 0 0 0
0 0 0 1 0 0 0 ■ 0 0 0 0
0 0 0 0 1 2 0 0 ■ 0 0 0 0
0 0 0 0 0 1 4 0 ■ 0 0 0 0
0 0 0 1 2 5 32 1 64 1 2 0 0 0 0
0 0 0 15 16 7 32 1 64 15 16 0 0 0 0
0 0 0 0 _ 3 8 1_ 16 0 ■ 0 0 0 0
0 0 0 0 0 0 1 ■ 0 0 0 0
0 0 0 0 0 0 0 ■ 1 4 0 0 0
0 0 0 0 0 0 0 ■ 1 * 64 1 2 -5 32 1 64
0 0 0 0 0 0 0 ■ , 1 64 15 16 7 32 1 64
0 0 0 0 0 0 0 ■ 1 16 0 3 8 1 16
0 0 0 0 0 0 0 ■ 0 1 0 0
0 0 0 0 0 0 0 ■ 0 0 1 2 0
0 0 0 0 0 0 0 ■ 0 0 0 1 4 _
е3 =
1296 0 0 0 0 0 0 0 0 0 0 0
0 432 0 0 0 0 0 0 0 0 0 0
0 0 2160 0 0 0 0 0 0 0 0 0
7990 283 77 0 0 0 0 0 0 0 0 0
-41545 _1161 -245 0 0 0 0 0 0 0 0 0
-924040 -32372 -8900 0 0 0 0 0 0 0 0 0
_538 _12 -2 28 -1 0 0 0 0 0 0 0
-41140 -1288 -320 115 37 -11 0 0 0 0 0 0
438040 13520 3320 0 412 -112 0 0 0 0 0 0
0 0 0 164 0 -4 0 0 0 0 0 0
0 0 0 0 80 0 0 0 0 0 0 0
0 0 0 -13200 0 464 0 0 0 0 0 0
0 0 0 28 1 0 28 -1 0 0 0 0
0 0 0 -115 37 11 115 37 -11 0 0 0
0 0 0 0 -412 -112 0 412 -112 0 0 0
0 0 0 0 0 0 164 0 -4 0 0 0
0 0 0 0 0 0 0 80 0 0 0 0
0 0 0 0 0 0 -13200 0 464 0 0 0
0 0 0 0 0 0 28 1 0 -538 12 -2
0 0 0 0 0 0 -115 37 11 41140 -1288 320
0 0 0 0 0 0 0 -412 -112 438040 13520 3320
0 0 0 0 0 0 0 0 0 7990 -283 77
0 0 0 0 0 0 0 0 0 41545 -1161 245
0 0 0 0 0 0 0 0 0 -924040 32372 -8900
0 0 0 0 0 0 0 0 0 1296 0 0
0 0 0 0 0 0 0 0 0 0 432 0
0 0 0 0 0 0 0 0 0 0 0 2160
Р1 =
1 0 0 0 0 0 "
0 1 2 0 0 0 0
0 0 1 4 0 0 0
1 5 1 1 -5 1
2 32 64 2 32 64
15 7 1 15 7 1
16 32 64 16 32 64
0 3 1 0 3 1
8 16 8 16
0 0 0 1 0 0
0 0 0 0 1 2 0
0 0 0 0 0 1 4 _
О =
-4 8 -1
84 -200 23
-828 2280 -251
1 0 0
0 35 0
0 0 35
-4 -8 -1
-84 -200 -23
-828 -2280 -251
Здесь многоточия обозначают повторяющиеся блоки матрицы Р3. Блоки матрицы Р составлены из коэффициентов масштабных соотношений для эрмитовых сплайнов 5-й степени (1), тогда как блоки матрицы О - из коэффициентов соотношений (3). Обратный процесс разбиения коэффициентов С на более грубую версию СЬ— и уточняющие коэффициенты Б-1 состоит в решении системы линейных уравнений (5). Разрешимость данной системы гарантирована линейной независимостью базисных функций. При этом для ускорения вычислений матрицу [рь | О] можно сделать ленточной, изменив порядок неизвестных так, чтобы столбцы матриц Р и О перемежались (см. [13]).
Процедуру разбиения С на часть С-1, соответствующую низшему разрешению, и уточняющие коэффициенты Б-1 можно применить рекурсивно и к самой этой части С-1. Следовательно, исходные коэффициенты можно представить в виде иерархии грубых версий с разрешениями С0, С1, ..., С-1 и уточняющих деталей Б0, Б1, ..., Б-1. Результирующее вейвлет-разложение эрмитового сплайна 5-й степени (х) может быть записано в виде
^ (х)= £ С0,к < (х) + £ £ Б/,к М{к
к=0 ^ i=0 ] =0 г=0
При этом по величине вейвлет-коэффициентов Б, } = 0, 1,., Ь-1, можно судить о значимости соответствующих уточняющих деталей. Незначимые убираются с целью сжатия информации.
(х)
а < х < Ь.
Для хе[0,1] положим
и пусть
3. Пример
у0 (х) = V231/272 М°0 (х), у2( х) = 3/8М0,1( х),
^0 (х) = л/ 231/13 М°2 (х), \1/2
уЬ-2 (х) = (21'1695 /30371) /1296 мЬ/30 (х),
уЬ-1(х) = (21'4322/705)1/2 /432МЬ/31 (х), уЬ (х) = (21'505979/208)1/2 /2160МЬ/3 2(х),
г = 3,3• 2Ь, Ь > 1;
уІ_2(х) = (2І8733/1762)1/2 /164МІ/3,0(х), у^(х) = (2І86625/716)1/2 /80МІ/31(х), уІ (х) = (088893/20/ /464М[/3,2(х), і = 6,9,...,3 • 2і _3, і > 2.
Заметим, что у і (х) нормированы так, что ||у і (х)||^ ^01 ^ = 1 для і = 1,2,...,3 • 2і.
Рассмотрим в качестве тестовой функции многочлен пятой степени А(х) = (1/2 - х)5 с числом разбиений п = 25 на интервале 0 < х < 1, длина шага Ах = 1/п. Число используемых значений функции и двух производных на интервале при выполнении вейвлет-преобразования равно 3(п+1) = 99.
Начиная с верхнего уровня разрешения і = 5, находим последовательно при
і = 5:
Б4 = [0.08679; -0,1174; 0,03294; -0,005044; -0,002738; -0,002359; -0,002481; -0,001461; -0,001019; -0,001459; -0,0006021; -0,0006303; -0,0006729; -0,0001976; -0,0002875; -0,0002498; -3,499 10-5; -0,0001086; -5,528 10-5; 2,879 10-6; -2,483 10-5; -2,848 10-6; 1,055-10-6; -1,689 10-6; 2,848-10-6; 1,055-10-6; 1,689-10-6; 5,528 10-5; 2,879 10-6; 2,483 10-5; 0,0002498;
-3,499 10-5; 0,0001086; 0,0006729; -0,0001976; 0,0002875; 0,001459; -0,0006021; 0,0006303; 0,002481; -0,001461; 0,001019; 0,005044;
-0,002738; 0,002359; -0,08679; -0,1174; -0,03294] т;
і = 4:
Б3 = [-0,07943; 0,1773; -0,09874; -0,00107; 0,00016; -0,0005311; -8,338 10-5; 6,975 10-5; -3,565 10-5; -1,852 10-5; 1,13510-5; -1,869 10-5; 1,852 10-5; 1,135-10-5; 1,869-10“5; 8,338 10-5; 6,975 10-5; 3,565 10-5; 0,00107; 0,00016; 0,0005311; 0,07943; 0,1773; 0,09874] т;
і = 3:
Б2 = [-0,07274; 0,1526; -0,08045; -0,0002049; 1,39110-5; -0,0001316;
0,0002049; 1,39110-5; 0,0001316; 0,07274; 0,1526; 0,08045] т;
і = 2:
Б1 = [-0,04817; 0,09931; -0,05134; 0,04817; 0,09931; 0,05134] т;
і = 1: на последнем шаге остается шесть значений многочлена пятой степени и двух его производных в концах отрезка С0 = [0,05701; -0,3206; -0,858; -0,05701; -0,3206; -0,858] т, и в центре три вейвлет-коэффициента разложения функции и производных Б0 = [4,653-1046; -0,001413; -3,571-10-16] т.
После обнуления вейвлет-коэффициентов, по модулю меньших 0,032, остается 30 значений, что дает коэффициент сжатия 99/30 = 3,3.
На рис. 3, 4 представлены результаты реконструкции значений сплайна 56(х) и двух его производных в узлах сетки А5: у, у1і, у2і, і = 0,1,..., 32. Здесь сплошными линиями обозначаются исходная функция и ее производные: Ах), /1(х), /2(х), 0 < х< 1.
Рис. 3. Сравнение результатов вейвлет-реконструкций узловых значений сплайна 5-й степени для случая функции у(х) = (1/2 - х)5
х, г/32 х, г/32
Рис. 4. Сравнение результатов вейвлет-реконструкций узловых значений производных сплайна 5-й степени для случая функции ^х) = (1/2 - х)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. Iss. 4. P. 335-354.
5. Strela V., Heller P.N., Strang G., Topivala P., Heil C. The application of multiwavelet filter-banks 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: Bior-thogonal 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.
8. Han B. Approximation properties and construction of Hermite interpolants and biorthogonal multiwavelets // J. Approxim. Theory. 2001. V. 110. P. 18-53.
9. Турсунов Д.А., Шумилов Б.М., Эшаров Э.А., Турсунов Э.А. Новый тип эрмитовых мультивейвлетов пятой степени // Пятая Сибирская конференция по параллельным и высокопроизводительным вычислениям / под ред. проф. А.В. Старченко. Томск: Изд-во Том. ун-та, 2010. С. 101-105.
10. Koro K., Ade K. Non-orthogonal spline wavelets for boundary element analysis // Engineering Analysis with Boundary Elements. 2001. V. 25. P. 149-164.
11. Завьялов Ю.С., Квасов Б.И., Мирошниченко В.Л. Методы сплайн-функций. М.: Наука, 1980. 352 с.
12. Strang G., Strela V. Short wavelets and matrix dilation equations // IEEE Trans. Signal Processing. 1995. V. 43. No. 1. P. 108-115.
13. Столниц Э., ДеРоуз Т., Салезин Д. Вейвлеты в компьютерной графике: пер. с англ. Ижевск: НИЦ «Регулярная и хаотическая динамика», 2002. 272 с.
Шумилов Борис Михайлович Томский государственный университет,
Кудуев Алтынбек Жалилбекович Ошский государственный университет
E-mail: [email protected], [email protected] Поступила в редакцию 31 августа 2012 г.
Shumilov Boris M., Kuduev Altynbek Z. (Tomsk State University, Osh State University). New type multiwavelets of the fifth degree orthogonal to quintic polynomials.
Keywords: Hermitian splines of the fifth degree, multiwavelets, orthogonality to polynomials.
For the space of Hermitian splines of the fifth degree of a kind
SL W = E SCf,k NLk (x), a < x < b,
k =0 i=0
with a uniform grid of nodes AL : ui = a + (b - a) i / 2L, i = 0, 1,..., 2L, L > 0, and the basic functions NLk(l) (uj ) = Sj • Sk, l = 0,1,2 , with the centers in integers, it is proposed to use as wavelets functions MLi,k(x), satisfying the conditions of orthogonality to all polynomials of sixth order, i.e. I" MLk(x)xmdx = 0, k = 0,1,2 Vi (m = 0,1,.. .,5). For the wavelets centered at even integers,
Ja ’
and the supports equal to the supports of basic splines on the grid AL, received two-scale relations of expansion and formulas for calculating the coefficients in the thinned grid AL_1 from the spline coefficients in a dense grid AL in the form of solution of linear algebraic equations system with band matrix. There are results of numerical experiments presented. Justified the improving of compression of numerical data in comparison with the known quintic wavelets and multiwavelets.
