Научная статья на тему 'Мультивейвлет пятой степени'

Мультивейвлет пятой степени Текст научной статьи по специальности «Математика»

CC BY
186
48
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ЭРМИТОВЫ СПЛАЙНЫ / ВЕЙВЛЕТЫ / РАЗЛОЖЕНИЕ / HERMITIAN SPLINES / WAVELETS / DECOMPOSITION

Аннотация научной статьи по математике, автор научной работы — Шумилов Борис Михайлович, Эшаров Элзарбек Асанович, Кудуев Алтынбек Жалилбекович, Ыманов Улукбек Сайдакрамович

Предложены два новых типа мультивейвлетов на основе эрмитовых сплайнов 5-й степени. Получен алгоритм вейвлет-разложения. Представлены результаты численных экспериментов.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

The authors have proposed two new types of multi-wavelets on the basis of quintic Hermitian splines. The algorithm of wavelet decomposition was obtained. The article introduces the results of the numerical experiments.

Текст научной работы на тему «Мультивейвлет пятой степени»

УДК 519.6

МУЛЬТИВЕЙВЛЕТ ПЯТОЙ СТЕПЕНИ

Б.М. Шумилов, Э.А. Эшаров, А.Ж. Кудуев, У.С. Ыманов

Томский государственный архитектурно-строительный университет E-mail: sbm@tsuab.ru; elzare@mail.ru Ошский государственный университет, Кыргызская Республика E-mail: altun_12@rambler.ru; ymanv8106@rambler.ru

Предложены два новых типа мультивейвлетов на основе эрмитовых сплайнов 5-й степени. Получен алгоритм вейвлет-разложения. Представлены результаты численных экспериментов.

Ключевые слова:

Эрмитовы сплайны, вейвлеты, разложение.

Key words:

Hermitian splines, wavelets, decomposition.

Вейвлетом называется маленькая, т. е. короткая или быстро затухающая волна, множество сжатий и смещений которой порождает некоторое пространство ограниченных функций на всей числовой оси [1-4]. Если таких волн несколько, то возникают мультивейвлеты [5, 6].

В данной статье мы построим базисные мультивейвлеты на основе эрмитовых сплайнов пятой степени. При этом, наряду с классическим, рассмотрим неизвестный ранее тип «ленивых» мультивейвлетов и обоснуем новый подход к вычислению вейвлет-преобразования на основе конечных неявных соотношений разложения с расщеплением по четным и нечетным узлам.

Основой для построения вейвлет-преобразования является набор вложенных пространств ...УисУ£сУш.... В данном случае пространство Уь является пространством сплайнов степени 5 гладкости С2 на отрезке [а,Ь] с равномерной сеткой узлов Аь: и=а +(Ь-а) {/2Ь, 1=0,1,...,2Ь, L>0, и базисными функциями Щк^)=щ^Н), &=0,1,2У£, где v=2L(u-a)ДЬ-a)+1, с центрами в целых числах, порожденными сжатиями и сдвигами трех функций вида [7]:

Ро(0 Pi(t) P2(t)

t3(6t2 -ist + іо)

-t 3(3t2 - 7t + 4)

t3

— (t2 - 2t +1)

,о < t < 1;

(2 - t)3(6t2 - 9t + 4) (2 - t)3(3t2 - St + 2)

(2 -1) (t2 - 2t +1)

,1 < t < 2

(рк(г) = 0, к = 0,1,2, г й [0,2].

На любой сетке А\ L>0, интерполяционный эрмитов сплайн 5-й степени может быть представлен как

(u) (u), a <u <b,

(1)

где коэффициенты С^\ к=0,1,2, являются значениями и, соответственно, первыми и вторыми производными аппроксимируемой функции в узлах.

Суть вейвлет-преобразования состоит в том, что оно позволяет иерархически разложить заданную функцию на серию все более грубых приближенных представлений и локальных уточняющих подробностей. При этом «более грубый» уровень представления функции в У^ получается из «более подробного» уровня представления функции в УL посредством прореживания (удаления каждого второго, как правило, отсчета). Здесь необходимо лишь, чтобы каждая базисная функция в У^ могла быть выражена в виде линейной комбинации базисных функций в Уь. В частности, двухмасштабное соотношение для эрмитовых сплайнов 5-й степени можно записать в виде следующей векторной формулы [5]:

(2)

Ро(t) 2 Ро(2 - к)

P(t) - t (N

_Р2(0_ к=о p2(2t - к)

где

Hо и

1 1S о

2 16

S 7 3

32 32 в

1 1 1

64 64 16

- -

1 о о

, H1 и о 1 2 о

о о 1 4 _

H2 и

1 1S о

2 16

S 7 3

32 32 в

1 1 1

64 64 16

В этом случае базисными функциями для У1А будут функции ИУ-1, с носителями в два раза большими по ширине и центрами в четных целых числах. Следующим этапом является определение

її

пространства уточняющих подробностей WL-1. В отличие от классического определения вейвлетов, мы не требуем, чтобы базисные функции из WL-1 были ортогональны базисным функциям в У;-1. Вместо этого просто потребуем, чтобы пространство WL-1 являлось дополнением У;-1 до У;. Следовательно, любая функция в УL может быть записана в виде суммы некоторой функции в У;-1 и некоторой функции в WL-1. Очень простой способ получения базисных функций в WL-1 заключается в использовании функций Щк в УL с центрами в нечетных целых числах [3].

Для облегчения выполнения дальнейших действий удобно записать базисные сплайн-функции в виде единой матрицы-строки и

упорядочить коэффициенты сплайна в виде вектора СМС0;,0,С0;,\С0;,2,С;,0,С;,\.,С£2]Т. Тогда уравнение (1) переписывается как SL (и)=^ (и)С;. Аналогично обозначим базисные вейвлет-функции как М^‘=%^-2£+1), к=0,1,2, i=1,2,...,2L-1, и запишем их в виде матрицы-строки 1/=[М;,оМи,М;,2,...,М212]. Соответствующие коэффициенты вейвлет-разложения на уровне разрешения ; будем собирать в вектор

б;=[б;,0,б;д,б;,2,...,б£2]т.

Тогда для уровня разрешения ;-1 можно записать функции (р;-1 и ^;-1 в виде линейных комбинаций функций (^‘,^1-1=^‘Р; и \/‘-1=^‘О1‘, где блоки матрицы Р; составлены из коэффициентов соотношений (2), так как каждую широкую базисную функцию внутри отрезка аппроксимации можно построить из трех, а по краям интервала из двух, троек узких базисных функций, тогда как все элементы столбцов матрицы - нули, за исключением единственной единицы, так как каждый ленивый вейвлет - это одна узкая базисная функция.

Следовательно, справедлива цепочка равенств С;=^-1С ;-1+ ^1'-1В;-1=^'Р;С ;-1+(^О'Б1'-1.

Таким образом, процесс получения С из С;-1и Б;-1 может быть записан как С;=Р;С;-1+О;Б;-1или, используя обозначения для блочных матриц,

С1

(3)

Обратный процесс разбиения коэффициентов С на более грубую версию С;-1 и уточняющие коэффициенты Б;-1 состоит в решении системы линейных уравнений (3). Для данного случая блоки матрицы [А;|В;], обратной по отношению к [Р^], являются разреженными. Поэтому процесс создания версии с низшим разрешением, С1-1, характеризуемой меньшим количеством коэффициентов, можно выразить в явном виде с помощью матричного равенства С;-1=А;С;. При этом потерянные детали собираются в другой вектор Б;-1, определяемый выражением Б;-1=В;С;.

При выполнении анализа заданной функции в соответствии с полученным выше результатом грубое приближение получается из более точного путем исключения узлов, соответствующих нечетным числам. Следовательно, самое грубое прибли-

жение зависит только от нескольких начальных значений, и оно может оказаться очень плохим приближением исходной функции. Чтобы улучшить усредняющие свойства представленного метода анализа данных, мы предлагаем использовать в качестве вейвлетов для WL-1 функции NLk в У; с центрами в четных целых числах. Поскольку WL-1 должно являться дополнением У;-1 в У;, размерности этих пространств должны удовлетворять соотношению Dim(VL)=Dim(VL-1)+Dim(WL-1). Для выполнения этого условия мы предлагаем из исходных координат вычитать уравнение прямой, соединяющей первую и последнюю вершины при дополнительном условии, что вторая производная в последней точке обращается в нуль.

Будем обозначать базисные сплайн-функции и коэффициенты эрмитового сплайна 5-й степени с отсутствующими элементами по концам отрезка аппроксимации как (р0; и С/\ Аналогично обозначим базисные вейвлет-функции как М’г^у^^г), к=0,1,2, £=0,1,...,2;-1 и запишем их в виде матрицы-строки

м1 м1 м1 м1 м1

0,1 ’ 0,2 ’1,0 ’1,1 ’1,2 ’ ***’

мі , мі , мі , Мі

’ 2і-1,0’ 2і-1,1’ 21 -1,2’ 21,1

Соответствующие вейвлет-коэффициенты на уровне разрешения ; будем собирать в вектор

Д1 =

-'-'о

Дід Д1,2 Д1,0 Дід Дд2 0 ’ 0 ’ 1 ’ 1 ’ 1 ’

.,ДІ° , ДІ1 , ДІ2 , Д

’ 2і-Г 2і-1’ 2і-1’ 2

,1,1 21

В результате вейвлет-преобразование может быть записано как

С0 = [рМеТ ]

С1

А"

(4)

Нетрудно проверить, что матрицы, обратные по отношению к [Р0;|ОД теряют разреженную структуру. Поэтому систему линейных уравнений (4) приходится решать численно. При этом матрицу [Р0;|О0;] удобно сделать ленточной, изменив порядок неизвестных так, чтобы столбцы матриц Р0; и О; перемежались [3]. Тем не менее, хотя разрешимость полученной системы и гарантирована линейной независимостью базисных функций, вопрос ее хорошей обусловленности остается открытым. В [6] для частного случая кубических эрмитовых сплайн-вейвлетов с помощью метода неопределенных коэффициентов были впервые получены конечные неявные соотношения разложения. В матричном виде полученные результаты можно представить следующим равенством [P0L|Q0L]RL=GL, где матрица RL представляет собой простую ленточную матрицу, а матрица GL - трехдиагональную матрицу со строгим диагональным преобладанием. После этого решение системы уравнений типа (4) можно записать в матричном виде как [8]:

С1

= [ріа1 ]-1 Сі = я (С1 )-1 С0,

(5)

т. е. решение сводится к трехдиагональнои системе линейных уравнении. Для представленного выше типа мультивейвлетов, как нетрудно проверить непосредственным вычислением, справедливы аналогичные равенства, например,

0 0 0 0 0 0 1 0 0 1

0

0

-32

0

0

0 0 16

0 0 0

[р !Є2

0 0 0 0 0 0 0 0 00 0 0 1 0 01 0 0 0 0 00 0 0

0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -150 -

-1260 4

30

180

0

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

75

315

-4

-15

-45

0

0

0

0

0

1

0

0

0

-16

0

-1

-8

-48

0

8

0

1

4

12

0

32 0 0

0 64 0

0 0 0

-16 0 0

0 -16 0

0 0 1

0 0 0

0 0 0

0 0 0

111 16 0 0

0 99 2 0

0 0 1

0 0 0

0 0 0

0 0 0

_1_ 16 0 0

0 1 2 0

0 0 0

-168 0

-1960 0

2 0

14 0

76 0

0 0

84 0

490 0

-2 1

-7 0

0 0 0 0 0 0

-4 0 0 0 1

-19

0

1 0 0

0 1 0

0 0 -5 -5 1 0

0 0 7 -7 1 0

0 0 12 12

0 0 0

0 0 0

0 0 0

0 0 8

0 0 -15

0 0 0

0 0 0

1 0 0

0 1 0 0 0 8

0 0 15

0 0 0

0 0 0

0 0 0

0 0 0

0 0 16

0 0 0

0 0 0

0 0 0

0 0 0

0 0 0

0 0 -16

1 0 0

0 0 1 0

0 0 0 0

0 0 0

0 0 0

0 1 0

_1_ 16 -1 0

0 -4 0

0 0 0

0 0 0

0 0 0

0 1 0

111 16 1 0

0 -4 0

0 0 1

0 0

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

0 0

1 0

-8 0

48 64

-16 0

0 0

0 0

-1 0

4 0

-12 -16

8 0

(6)

Следующее утверждение дает последовательность вычисления коэффициентов вейвлет-анализа по известным коэффициентам сплайн-разложения на любом уровне разрешения Ь, Ь>2.

Теорема 1. Пусть значения сплайн-коэффициентов С^2, СЬ0 и СЬ1 в нечетных узлах пересчитаны из последовательного решения трех систем линейных уравнений, соответственно, вида

99 2 -4 1 2

1 2 -4 111 2

111

2

1 2

53 4 8 1 8

1 8 8 111 8

111

8

1 8

-4

-4

СІ'

Сзх'

сХ,(

" сХ,2" " С^2"

СХ ,2 С3 := СХ ,2 С3

СХ ,2 _ 2х-1_ СХ ,2 _ 2х-1_

с1і>0 - с^2

СХ ,0 С^ ,2

С3 - С1

с5 ,0

С"Х ,0

2х-3 тХ,0 "'і1 -1

с

С"Х ,2 2х-1 сХ ,2 2х-1 .

111 1

16 16

1 55

16 8

1 16

_1_

16 111 16 .

сХ

СХ

с

СХД + ИС^0 + СХ С3ХД - 15С1Х,0 - СХ СХД

С

С,

2х - 3 Х,1

-15С,і , + С,Г, 2х-1 2х-1

15С

С

-150 -16 -168

-1260 0 -1960

4 -1 2 16

30 -8 14 0

180 8 4 - 76 0

Здесь точки, расставленные по диагонали, означают, что предшествующий столбец повторяется соответствующее число раз, сдвигаясь при этом вниз на одну позицию.

Тогда вектор сплайн-коэффициентов на прореженной сетке АЬ-1 представляет собой результат умножения матрицы

1 0

- 8 0 48 64

-1 0 -8 0 - 48 0

на вектор [СЬ,0,СЬД,СЬ,2,...,СЙ]т сплайн-коэффициентов в нечетных узлах густой сетки АЬ, тогда как вектор вейвлет-коэффициентов равен такому же произведению с матрицей

75 8 84

315 0 490

-4 1 -2 -16 -1 0

-15 4 -7 0 4 0

-45 12 -19 0 -12 -1

4 0 2

-30 0 -14 '

180 0 76

20 1 24

0 -8 0

-360 48 -1048 '

4 -1 2' 16 1 0

30 -8 14 ' 0 -8 0

180 8 4 - 76 ' 0 48 64

0 -16 0

1

4

12

-4 0 -2

15 0 7 '

-45 0 -19 ''

-20 -1 -24 '

0 4 0 '

90 -12 262 '

-4 1 -2 ' -16 -1 0

-15 4 -7 ' 4 0 0

-45 12 -19 '' 0 -12 -16

0 8 0

при условии добавления к результату вектора [С0ь,1,Сь,2,СЬ,0,СЬ,1,СЬ,2,.,С2ьь1]т сплайн-коэффициентов в четных узлах густой сетки. Здесь точки, расставленные по диагонали, означают, что предшествующие три столбца повторяются соответствующее число раз, опускаясь при этом вниз на три ряда.

Доказательство. Согласно построению, на каждых т внутренних шагах сетки АЬ перекрываются по 3(т-і) широких базисных функций и вейвлетов и 3(2т-1) узких базисных функций. Поэтому после соответствующего изменения нумерации узлов на отрезке [0, т] можно записать следующие три конечных неявных соотношения разложения:

2т-2

£ а^к (2х - і) =

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

і=0

= £ Е (*ЛV (2X -1 - 2у) + 0к‘ V (х - Л), к = 0,1,2, (7)

1=0 У=0

где рк(2х-) - эрмитовы базисные сплайны на густой сетке; р(х-]) - эрмитовы базисные сплайны на прореженной сетке, р(2х-1-2]) - базисные вейвлеты на прореженной сетке.

Тогда для вычисления неопределенных коэффициентов соотношения (7) при к=0 с использованием табл. 1 и легко проверяемых равенств

1

—Гр,(2г) = 2', ? = -, к,/ = 0,1,2, соот-

ветственно, в точках

1

00 1 + с001 (- — 0 I 32

0 = с001!+с01 \-1.

0 8 + с° I 16

0 = с"! + с02 (-4) .

+ с0

+ сл’

64’

32,

х = 1: а0 = Ь000 + с™,

А £.01 О I .01

0 = Ьо ' 2 + со ,

0

02

0

02

0 = Ь - 4 + с0,

- + с„

32

- + с„

64

32

+ с/

64

+ с0015 + с01 (-— 1+ с02

1

16

П 01 I 3 | I 02 I 1

0 = с0 I — I + с I —

32

+ с

01 3 02 I 1

+ с, I —

'2 1 ( 4

х = т -1:

2т -1

= Ь

2 т—3 ‘-т-2

т_

0 = Ь01

+ ст

• 2 + с 2 2 + Ст-2 ,

2

0 = ьт2-2 • 4 + с“2,

а20 2 = с00 2 1 + с012 —

2т-2 т-2 2 т-2 32

64

0

0 = ст

0 = с0

15

— 1+с

- 3)+ст

— I+ст' 16) т

1

"32

Одно из решений полученной системы уравнений имеет вид: аот=Ь0=1, ./=0,...,т-2, при условии, что все остальные коэффициенты равны нулю. Это означает, что в нечетных узлах выполняются тождества для базисных сплайн-вейвлетов

%(2х-1-2;).

Попытаемся найти систему уравнений, связывающую коэффициенты разложения для четных узлов (случай, когда а01=а03=.=а02т-3=0). При т=2, 3 нетрудно убедиться, что система имеет только тривиальное решение. При т=4 решение имеет

вид: а01=а03=а05=Ь011=с011=0, Ь000=Ь002=-4, Ь00х=-20, Ь010=-15, Ь012=15, Ь020=Ь022=-45, Ь021=90, с000=с002= 4, с001=20, с010=30, с012=-30, с020=с022=180, с021=-360,

0 и 2 и 0 и 2

а00=а06=1/8, а02=а04=111/8.

Таким образом, из соотношения (7) выделяется неявное четырехчленное разложение вида р0(2х) +111ф0 (2 х - 2) + 111р0(2 х - 4) + р0(2 х - 6)

8

(>0(х) + 5^0(х -1) +^0(х - 2)-

1-р0(2х -1) - 5р0(2х - 3) - р0(2х - 5)) +15(2р1(х) - 2ф1( х - 2) - ф1(2 х -1) + ф1(2 х - 5)) + +180(р2(х) - 2р2(х -1) + р2(х - 2)) --45(р2(2 х -1) - 2р2(2 х - 3) + <р2(2 х - 5)).

Таблица. Значения базисных функций и их производных в точках отрезка [0,2]

X ро(х) р0(х) Ро(х) Р|(х) р1(х) р'1(х) рг(х) р2(х) р”(х)

0 0 0 0 0 0 0 0 0 0

1/2 1/2 15/8 0 -5/32 -7/16 3/2 1/64 1/32 -1/4

1 1 0 0 0 1 0 0 0 1

3/2 1/2 -15/8 0 5/32 -7/16 -3/2 1/64 -1/32 -1/4

2 0 0 0 0 0 0 0 0 0

Совершенно аналогично с помощью табл. 1 нетрудно убедиться, что для эрмитовых базисных сплайнов первых производных на густой сетке справедливо неявное трехчленное разложение вида (р1(2 х) +110р2(2 х - 2) + р2(2 х - 4)

16 =

Р0(2х - 1) - р0(2х - 3) - Р0(х) + Р0(х- 1) + +4(рД2х -1) + (р1(2 х - 3) - 2^1( х) - 2^1( х - 1)) + 12(^2(2х -1) -р2(2х - 3) - 4р2(х) + 4р2(х-1)),

а для базисных сплайнов вторых производных имеет место неявное четырехчленное разложение вида р2(2х) + 111р2(2 х - 2) + 111р2(2 х - 4) + р2(2 х - 6) = 8 = 2(Р0(х) + 12Р0(х - 1) + Р0(х - 2) - ^ +

1-р0 (2х -1) - 12р0(2х - 3) - р0(2х - 5)1 +7(2ф1(х) - 2ф1(х - 2) - (р1(2х -1) + (р1(2х - 5)) + +38р2 (х) - 524р2 (х -1) + 38р2(х - 2) -19р2 (2 х -1) + 262р2 (2 х - 3) - 19р2 (2 х - 5).

Соответствующие разложения по краям отрезка аппроксимации содержатся в матричном равенстве (6). Введем последовательности матриц Оь и Яь, блоки которых составлены из коэффициентов левых и правых частей полученных разложений, соответственно. В результате находим, что базисные функции пространства эрмитовых сплайнов 5-й степени на густой сетке, базисные функции на прореженной сетке и вейвлеты удовлетворяют равенству р0;Оь=[р!т1\у0г1] Нь, Ь>2.

Отсюда, используя свойство дополнения пространства вейвлетов, находим

х=

ivL-1]

cL

__ Ls~rL

= Co C0

= [Co

1 v0i-

1] rl ( gl )-1 c(

После этого решение системы уравнений (4) можно записать в виде (5), откуда после расщепления по четным и нечетным узлам вытекает утверждение Теоремы 1.

Пример. Рассмотрим в качестве тестовой функции функцию Хартена:

У$т(3жх), х < у,

/(х) =

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

|sin(4^x)|,

- "2sin(3^ x),

< x < f,

x >-

Начиная с верхнего уровня разрешения 1=5, то есть при числе разбиений и=2ь=32, на интервале 0<х<1 с длиной шага Ах=1/п находим последовательно, при Ь=5: Д)4=[-5,856, -88,1, 3,295, 5,904, -25,81, 14,12, 2,663, -173,8, 8,263, -1,359, -83,65, 11,54, -4,808, -144,9, 2,226, -19,38, 60,93, 66,25, 0,1275, -855,7, 6,406, -8,153, 0,7788, 23,93, -12,57, -352, 6,406, 8,153, 0,7788, 66,25, -0,1276, -855,7, 2,029, 10,15, 78,23, -8,348, -5,076, 95,12, -11,52, -1,357, 132,6, -10,86, 2,664, 124,9, -6,528, 5,788, 75,11, 6,961]т;

1=4: Д)3=[-407,3, -2522, -121,3, 15,31, 2087, 44,06, -29,11, -313,5, 17,2, -9,35, 2376, 127,6, -0,003728, -3,504, 17,1, 9,76, 2375, 74,52, -8,945, -885,9, 93,62, 10,86, -1204, 28,64]т; 1=3: Д)2=[-13980, -80860, -1099, 986, 23240, 459,8, -10,01, -5093, -1401, 235,3, 19260, -802,2] т; 1=2: £„‘=[-11760, 225000, 27040, 1022, -262700, 15090] т; Ь=1: на последнем шаге остается три коэффициента разложения производных на концах интервала С00=[-3,871-106, -4,62 5 107, 5,768 105]т и три вейвлет-коэффициента £00=[2,018-106,

1,2 1 4 107, -3,156 1 05] т.

При условии обнуления незначимых вейвлет-коэффициентов Д)2(6),ДД3), .,£04(48) итоговый коэффициент сжатия равняется 3*33/(96-32)=99/64=1,55.

Несложно предложить параллельную реализацию представленного в статье алгоритма вейвлет-преобразования эрмитовых сплайнов 5-й степени, в которой три прямых хода прогонки выполняются независимо, а три обратных хода выполняются с максимальным запаздыванием на два такта. Это позволит преодолеть некоторое отставание в вычислительной эффективности вейвлет-алгоритмов в сравнении с универсальными алгоритмами сжатия [9].

Работа выполнена при поддержке РФФИ 13-01-90900 молиннр, 13-07-90900 мол_ин_нр.

СПИСОК ЛИТЕРАТУРЫ

1. Добеши И. Десять лекций по вейвлетам / пер. с англ. -Ижевск: НИЦ «Регулярная и хаотическая динамика», 2001. -332 с.

2. Чуи Ч. Введение в вейвлеты / пер. с англ. - М.: Мир, 2001. -412 с.

3. Столниц Э., ДеРоуз Т., Салезин Д. Вейвлеты в компьютерной графике / пер. с англ. - Ижевск: НИЦ «Регулярная и хаотическая динамика», 2002. - 272 с.

4. Исаев Ю.Н. Конструирование биортогональных вэйвлет-бази-сов для оптимального представления сигналов // Известия Томского политехнического университета. - 2004. - Т. 307. -№ 1. - С. 37-42.

5. Strela V. Multiwavelets: Theory and Applications. - Cambridge: Massachusetts Institute of Technology, 1996. - 99 p.

6. Шумилов Б.М., Эшаров Э.А. Построение эрмитовых сплайн-вейвлетов // Вестник Томского государственного университе-

та. Сер. Математика. Кибернетика. Информатика. Приложение. - 2006. - № 19. - С. 260-266.

7. Стечкин С.Б., Субботин Ю.Н. Сплайны в вычислительной математике. - М.: Наука, 1976. - 248 с.

8. Шумилов Б.М. Алгоритм с расщеплением вейвлет-преобразования эрмитовых кубических сплайнов // Вестник Томского государственного университета. Сер. Математика. Механика. - 2010. - №4 (12). - С. 45-55.

9. Замятин А.В., То Динь Чыонг. Повышение эффективности трехэтапного алгоритма сжатия многозональных аэрокосмических изображений // Известия Томского политехнического университета. - 2008. - Т. 313. - № 5. - С. 24-28.

Поступила 09.01.2013 г.

i Надоели баннеры? Вы всегда можете отключить рекламу.