CC BY
158
ЭФФЕКТИВНЫЕ АЛГОРИТМЫ ВЫЧИСЛЕНИЯ ЛОКАЛЬНОГО ДИСКРЕТНОГО ВЕЙВЛЕТ-ПРЕОБРАЗОВАНИЯ
В. В. Мясников Институт систем обработки изображений РАН, Самарский государственный аэрокосмический университет
Аннотация
Основной задачей настоящей работы является разработка нового класса вычислительно эффективных алгоритмов локального вейвлет-преобразования. Предполагается замена иерархической (пирамидально-рекурсивной) вычислительной конструкции, присущей известным алгоритму «с дырами» («algorithme a trous», М.Холшнайдер и др, 1989) [12] и алгоритму быстрого ортогонального вейвлет-преобразования (алгоритм С.Малла, 1987) [3,6,11], на (горизонтально-) рекурсивную, в которой вычисление коэффициентов вейвлет-преобразования производится для всех позиций вейвлета заданного масштаба последовательно, то есть в режима «скользящего окна». Замена иерархической вычислительной конструкции на горизонтально-рекурсивную позволяет рассмотреть задачу построения базисных вейвлетов, которые удовлетворяют требованиям рекурсивности. Эта задача в общем случае включает в себя ряд подзадач, связанных с анализом основных классов вейвлетов: базисных (материнских), двухпараметрических, каркасов (фреймов), R-вейвлетов, полуортогональных и ортогональных. Настоящая работа посвящена вопросу построения эффективных алгоритмов для наиболее простого класса базисных (материнских) вейвлетов, а также определения условий/ограничений на вейвлеты, при которых такие алгоритмы существуют. Приводятся примеры новых и известных вейвлетов, для которых существуют эффективные (рекурсивные) алгоритмы вычисления локального дискретного вейвлет-преобразования.
Введение
Вейвлеты являются одним из эффективных инструментов пространственно-временного анализа данных. Во многих работах основной интерес представляют биортогональные и ортогональные вейвлеты [3,6,10,11,12]. Одной из причин этого является наличие для вейвлетов этих классов быстрого алгоритма вычисления вейвлет-разложения. Цель настоящего работы - показать, что существует подкласс вейвлетов, которые не обязательно являются ортогональными или биортогональными, но для которых существуют вычислительно эффективные алгоритмы расчета дискретного вейвлет преобразования (ДВП).
Работа организована следующим образом. В первом разделе приводятся известные понятия базисного (материнского) вейвлета, интегрального и дискретного вейвлет-разложения. Указываются условия существования обратного преобразования для вейвлет-преобразования сигнала.
Второй раздел содержит краткое изложение одного частного решения задачи построения эффективного алгоритма, рассмотренной ранее в работах автора [7,8]. В частности, в этом разделе приведен эффективный алгоритм вычисления дискретной линейной свертки, который порождается сплайн-представлением (дискретного) ядра свертки. Приведено выражение для сложности такого алгоритма как функции параметров дискретного обобщенного сплайна. Полное изложение результатов дано в работах [7,8].
В третьем разделе вводится понятие сплайн-вейвлетов как дискретных обобщенных сплайнов с конечным носителем, удовлетворяющих условию
обратимости интегрального вейвлет-преобразо-вания. Доказывается, что сплайн-вейвлеты являются базисными вейвлетами. Для введенного класса сплайн-вейвлетов оказываются применимы все результаты второго раздела работы. А именно, любой сплайн-вейвлет произвольного (не обязательно целого) масштаба порождает эффективный алгоритм вычисления локального ДВП с конкретной вычислительной сложностью. Очевидная связь между сплайн-вейвлетом и сложностью порождаемого им алгоритма позволяет рассмотреть задачу построения материнских вейвлетов, обладающих предельно низкой сложностью вычисления локального ДВП. Приводятся несколько примеров таких «эффективных» сплайн-вейвлетов. Также дано выражение для сложности эффективного алгоритма, порождаемого известными биортогональными полиномиальными сплайн-вейвлетами [3,6,11].
В заключение работы приводятся ее основные результаты, благодарности фондам, поддерживающим данную научную работу, а также список литературы.
1. Основные сведения из теории вейвлетов
В общем случае под вейвлетом понимают функцию у (г) е I) (К), для которой справедливо соотношение
) = 0. (1)
Ограничение (1) не является единственным. В зависимости от вида вейвлет-преобразования оно дополняется рядом других, обеспечивающих
обратимость преобразования. Различные типы вейвлет-преобразований, обратные к ним и условия обратимости приведены, например, в монографии [11].
Для определенности ниже принята следующая нотация для прямого и обратного непрерывного Фурье-преобразования функции / (г) е I) (К):
||/(г)|(1 + )Л <+» .
■ (8)
1.2. Дискретное вейвлет-преобразование
Дискретизация ИВП
При численном расчете ИВП его значения берутся на дискретном множестве точек:
/ Н=| / (гК"^ / (г) = ¿1 / (ю)Лю . (2) /-1 "Мм, = | / (Г)^ (Г)ЛГ (9)
—^ I I ' '
1.1. Интегральное вейвлет-преобразование Определение 1. [11] Если у (г)е I) (К) удовлетворяет условию допустимости в виде
|у (ю))
Лю < +<х> ,
(3)
то у(г) называется базисным (материнским) вейвлетом. Относительно каждого базисного вейвлета у (г) интегральное вейвлет-преобра-
зование (ИВП) функции / (г) е I2 (К) определяется формулой:
(Ш/)(а, Ь) = |ар1/2| / (г)у(—1 Лг , а (> 0), Ь е К . (4)
Условие восстановления исходной функции из (4) определяет следующая
Теорема 1. (Кальдерон, Гроссман, Морле) [6] Пусть у (г)е I} (К) - вещественная функция такая, что выполняется (3). Тогда любая /(г) е I) (К) удовлетворяет равенству:
/ (г )= £-Я(И/ )(а, Ь )(| а|-1'2 у
г — Ь 11 „ Ла
— II ЛЬ~Г-
а Л а
(5)
В [6] указано достаточное условие для выполнения условия допустимости (3):
у (о) = | у (г)Лг = 0 и
(6)
- у(ю) непрерывно дифференцируема, то есть производная существует и непрерывна.
Последнее условие выполняется, если у (г) достаточно быстро убывает во времени [6]:
I |у(г)(1 + |г|)Л <+ю , (7)
что является следствием следующего утверждения.
Предложение 1.[6] Функция / (ю) ограничена, "
раз непрерывно дифференцируема и имеет ограниченные производные, если
где
уь„.1." (г ) =
4т у( ^
ао ( ао
}," е 2 ,
(10)
а параметры а0,Ь0 е Я + выбирают из множества [6]:
а0
= 21/", V = 1,2,3...
Ь 1 1 1 I1 I1 13 Ь0 =...—,—,1,1—,1—,1—,..
0 4 2 4 2 4
(11)
Дискретный аналог выражения (9), в котором исходный сигнал является дискретным, называют дискретным вейвлет преобразованием (ДВП) [6].
Быстрый алгоритм ДВП на основе БПФ
Вычислительно эффективный способ расчета ДВП в общем случае заключается в использовании для каждой пары (], ") алгоритма быстрого преобразования Фурье (БПФ) [6,9]. В частности, если 8ирр у (г)=[—К /2, К /2], а обрабатываемый цифровой сигнал имеет N отсчетов, тогда:
• имеем V log2 (N / (2К)) масштабов,
• для каждого масштаба вычисление дискретного аналога «свертки» (9), использующей БПФ, требует О ((log2 N).
Поэтому общее число операций на вычисление дискретного аналога непрерывного вейвлет-преобразования равно:
О(тN)2).
(12)
Описанный алгоритм может быть использован для произвольных вейвлетов, в том числе и для всех рассматриваемых ниже. Однако, при некоторых дополнительных ограничениях на вид вейвлета и параметры а0, Ь0 е К+ существуют быстрые алгоритмы, сложность которых существенно ниже величины (12). Одним из наиболее известных алгоритмов является алгоритм быстрого ДВП С.Малла [6], который используется для вычисления ДВП для ортогональных и биортогональных вейвлетов.
Суть алгоритма быстрого ДВП состоит в использовании двух конечных наборов фильтров {к(и)},{я(и)} [6]. Обозначая {/0(и)} - входную последовательность значений цифрового сигнала, коэффициенты {Л]+1 (и)} ДВП для вейвлетов уровня
■
]+1 могут быть получены путем вычисления последовательных сверток сигнала {/ (и)} текущего уровня/разрешением у и отсчетов фильтров {к (и)},{g (и)}
с последующим прореживанием полученного сигнала. А именно [6]:
/+1 1»(/ * ¿ИМ ^ [и]=(/ * фи],
здесь обозначение «[2и]» означает прореживание
полученного сигнала, то есть выборку каждого второго отсчета. Таким образом, для входного сигнала длины N сложность алгоритма быстрого ДВП составляет О ^), то есть не зависит от ширины носителя вейвлета [6].
Наряду с очевидными преимуществами, быстрый алгоритм ДВП имеет несколько недостатков:
• алгоритм может быть использован только для ортогональных и биортогональных вейвлетов,
• список масштабов для вейвлетов ограничивается только степенями целых чисел (в каждом конкретном случае - это степени числа «2» или степени числа «3» и т.д.),
• расчет коэффициентов разложения производится «неравномерно»: для вейвлетов более высокого уровня (с большим носителем) число коэффициентов разложения меньше. Данный недостаток, несущественный в задачах сжатия или кодирования сигнала, оказывается важным в других задачах:
- при визуализации (скейлограммы) ДВП [6],
- при анализе скейлограмм,
- при локальном анализе/распознавании сигналов, когда вычисление ДВП производится локально, то есть в режиме «скользящего окна».
Заметим, что во многих приложениях (физика, медицина и др.) часто используются неортогональные (и не биортогональные) вейвлеты, которые не имеют быстрого алгоритма вычисления. Наиболее известным примером такого вейвлета является вейвлет «мексиканская шляпа», который представляет собой вторую производную функции Гаусса [3,6,10,11]. Очевидно, что для таких вейвлетов, а также для дробных масштабов процесс вычисления ДВП приходится производить с использования алгоритма БПФ. Такая реализация ДВП, с учетом выражения (12), является вычислительно очень трудоемкой вследствие больших размеров обрабатываемых сигналов и самих носителей вейвлетов.
Цель настоящей работы - показать, что существует подкласс базисных вейвлетов, которые не обязательно удовлетворяют условию ортогональности или биортогональности, и для которых существуют вычислительно эффективные алгоритмы расчета локального ДВП, свободные от указанных недостатков.
Ниже не просто показано, что такие вейвлеты существуют, но также приводятся несколько примеров
таких вейвлетов как новых, так и известных. Указывается сложность соответствующих алгоритмов вычисления локального ДВП.
Для построения таких специфических вейвлетов и алгоритмов быстрого вычисления локального ДВП используется подход, предложенный автором в работах [7,8]. Подход кратко изложен во втором разделе. А представление вейвлетов и эффективных алгоритмов ДВП, построенных с использованием этого подхода, производится в третьем разделе.
2. Эффективный алгоритм вычисления свертки,
порождаемый сплайн-представлением КИХ
Настоящий раздел содержит тезисное изложение одного частного решения задачи построения эффективного алгоритма, рассмотренной ранее в работах автора [7,8]. А именно, здесь приводится (эффективный) алгоритм вычисления линейной свертки цифрового сигнала с конечной импульсной характеристикой (КИХ), который порождается сплайн-представлением этой КИХ. Приведено выражение для сложности такого алгоритма как функции параметров сплайна. Более подробное изложение этих результатов с доказательствами дано в работах автора [7,8].
Под обобщенным (дискретным) сплайном порядка К, определенным на разбиении/сетке Д5 отрезка [а,Ь],
далее понимается функция БК (Дх;и) (К > 0), которая на частных отрезках [г*, ^+1) сетки Дх : а = Г0 < Г1 <... < = Ь , совпадает с функциями р*. (/), которые являются решениями единого линейного рекуррентного соотношения (ЛРС) порядка (К +1):
р* (и) = Ё акР* (и - к) 5 = 0 £-1.
Выбирая ак = (-1) СкК, к = 1, К +1, получаем
полиномиальный сплайн, где р* (г) - многочлены
степени не выше К.
Определим далее величины:
Д :
V* = 1, £ -1 ^ е Z л ^ - Г5-1 > К +1 л Г0 = 0 л = М -1 л +1 = М;
© =
т е 0, М + К +1:
К+1
(ДМ +1;т) * Ё а Бк (ДМ, х+1;т - к)
к =1 К+1
ф (т) = I Бк (ДМ +1;т )-Ё ак БК (ДМ +1;т - к)
к=1
т е©; г =
Величина гЕ далее называется суммарным дискретным дефектом сплайна.
к =1
В работе [7] показано, что если КИХ представи-ма в виде обобщенного дискретного сплайна, удовлетворяющего некоторым ограничениям, то алгоритм вычисления свертки с такой КИХ может быть представлен следующим образом.
Алгоритм ACR для обобщенного сплайна
Шаг 1. Предварительная обработка:
y (n) = ^ х(n - m)p(m), n = 0, N -1.
me©
Шаг 2. Постобработка:
K+1 _
y(n) = ^y(n-1)+ y(n), n = 0, M-1,
к=1
с инициирующими значениями y (t) = 0 .
□
Формально, этот алгоритм оказывается рекурсивным и реализует вычисления по единой формуле:
K+1
У (n) = E акУ (n -1) +
к =1
х (n - m )ф (m), n = 0, M-1.
me©
Алгоритм Acr для полиномиального сплайна Шаг 1. Предварительная обработка:
У (n) = ^ х(n - m)p(m), n = 0, N -1.
me©
Шаг 2. Постобработка:
У(-1)(n) = y (n), n = 0,N-1, yh (n )=yh (n -1)+yh+1 (n),
n = 0,N -1, к = Kh -2,0 (Kh > 1), y (n ) = yh (n), n = 0, N -1.
□
Для представленных алгоритмов справедливо
Предложение 2. [7] Вычислительная сложность алгоритма ACR , порождаемого сплайн-представлением sK (Mst) КИХ, имеет вид
add + 'mul = ^ 'add ' 'mul e R)
- для обобщенного сплайна :
u(ACR)N^T1 = K +1 + ^ , (13)
- для полиномиального сплайна:
u (ACR =r^+K'add. (14)
■
Существенным моментом является то, что вычислительная сложность алгоритма ACR , порождаемого сплайн-представлением sK (aMs+ii t) КИХ,
оказывается ограничена сверху. Это следует из того, что величина суммарного дискретного дефекта сплайна ограничена сверху величиной [7]:
гЕ<( 5 +1)( К +1) +1. (15)
Таким образом, основным результатом настоящего раздела и работ [7,8] является связь между дискретным обобщенным сплайном, порождаемым им алгоритмом Лся вычисления свертки и сложностью этого алгоритма.
3. Сплайн-вейвлеты с конечными носителями
и эффективный алгоритм локального ДВП
Для удобства дадим краткое изложение представленных ниже результатов. Сначала устанавливается, что ограниченная вещественная функция с конечным носителем и нулевым средним удовлетворяет условию допустимости (3), то есть является базисным вейвлетом. Далее доказывается, что масштабирование (по аргументу) вещественной функции, последовательность значений которой удовлетворяет однородному ЛРС порядка К, также приводит к функции, последовательность значений которой удовлетворяет однородному ЛРС порядка не выше К. Это позволяет рассматривать множество сплайн-функций, удовлетворяющих однородному или неоднородному ЛРС, с конечным носителем и нулевым средним, как замкнутое относительно операции масштабирования. С учетом того, что функции этого множества, содержащие различные масштабы материнских вейвлетов, являются сплайнами, при определенных условиях (малом числе узлов и малом порядке) для каждой такой функции существует эффективный алгоритм Лск, порождаемый сплайн-вейвлетом. Именно этот алгоритм может быть использован для построения вычислительно эффективного ДВП. В завершении приводятся примеры сплайн-вейвлетов, для которых существуют эффективные алгоритмы ДВП.
3.1. Сплайн-вейвлеты с конечным носителем
Предложение 3. Пусть вещественнозначная функция у(г) удовлетворяет условиям:
да
|у (г))<да, Бцрр (у (г))<да, I у (г)Лг = 0. (16)
—да
Тогда у (г) - базисный вейвлет.
Доказательство: Как было указано выше, достаточным условием для выполнения условия допустимости (3) является пара условий (6) и (7). Условие (6) выполняется в силу (16). Условие (7) для ограниченных функций с конечным носителем также выполняется. Следовательно, у (ю) непрерывно дифференцируема, у (г) удовлетворяет условию допустимости (3) и является базисным вейвлетом.
■
Известно, что если последовательность у(и)
удовлетворяет однородному ЛРС порядка К, то решением ЛРС является функция вида [1,2]:
1-1 я,-1
у (г ) = ЁК)' Ё ГГвг , (17)
1=0 г=0
где ю,, в1г е С . Набор корней {со, } 1, их количество I и кратности {Я, } 1 корней определяют «тип»
ЛРС. А значения {р,г} определяют начальные условия ЛРС.
Тогда справедливо
Предложение 4. Пусть функция у(г) допускает представление в виде (17) с некоторыми
I-1
коэффициентами ю,,и Куг)=ЁЯ . Тогда
I =0
у(аГ) (а е Я \ {0}) также допускают представление в виде (17) с другими значениями коэффициентов Ю,, Р,г и Ку(аг) < Ку(г) .
Доказательство: производится непосредственной проверкой. Понижение значения величины К (порядка соответствующего ЛРС) происходит при вы-
р* (г), определяющие поведение сплайна на частичных интервалах разбиения его области определения, имеют вид:
полнении
условия: 3l1,l2 : ( ) =(ю/2 ) •
Это предложение означает, что масштабирование (по аргументу) функции, последовательность значений которой удовлетворяет однородному ЛРС порядка K, также приводит к функции, последовательность значений которой удовлетворяет однородному ЛРС порядка не выше K. Этот вывод, очевидно, остается верным и в случае, когда первоначальная функция является обобщенным сплайном. Из предложений 3 и 4 имеем
Следствие 1. Пусть вещественная функция у ( t )
удовлетворяет условиям (16) и допускает представление в виде обобщенного сплайна порядка K с заданными на частичных интервалах функциями вида (17), различающимися коэффициентами
L-1
{ßfr}, и K = Ё R . Тогда y(at) (ае R \ {0}) -
l =0
базисный вейвлет. ■
Последнее утверждение позволяет выделить среди всего множества базисных вейвлетов подмножество, элементы которого попадают под
Определение 2. Вещественнозначная функция y(t) называется (базисным) сплайн-вейвлетом,
если она удовлетворяет условиям (16) и допускает представление в виде обобщенного сплайна порядка K с (5^-1)-м интервалом, а функции
)=£ )) "tr ßri
где
Vs = 0, S -1 Ls = L, roj = ral,
RS < R,, K = " R, •
(18)
(19)
Дополнительные ограничения (19) данного определения указывают на то, что последовательности
значений всех функций {ps (t)} ^ удовлетворяют
одному и тому же однородному ЛРС. Без этого ограничения получаемая функция, в силу следствия 1 также будет являться вейвлетом, однако она не будет (обобщенным) сплайном.
Способ получения сплайн-вейвлетов тривиален. Его фиксирует следующее
Предложение 5. (достаточное условие того, что обобщенный сплайн является базисным вейвлетом) Пусть y(t) - обобщенный сплайн с конечной
областью определения. Если для y(i) выполняется условие (1), тогда y(t) - базисный сплайн -вейвлет. ■
Ниже в таблице 1, представлены некоторые полиномиальные сплайн-вейвлеты. Очевидно, привести другие примеры как для полиномиальных, так и для обобщенных сплайн-вейвлеты не представляет сложности, поскольку, следуя условиям последнего предложения, от сплайна достаточно добиться выполнения единственного условия (1). Поскольку других условий на обобщенный сплайн не накладывается, то все множество степеней свободы сплайна с заданным порядком K и числом узлов S можно использовать для удовлетворения дополнительных (внешних) ограничений. Типичными для вейвлетов являются ограничения на число нулевых моментов и гладкость (регулярность) вейвлета [3,6,11]. Напомним, что для ортогональных вейвлетов число нулевых моментов и гладкость являются жестко взаимозависимыми [3,6,11]. Для введенных сплайн-вейвлетов, очевидно, эта зависимость существенно слабее, поскольку нет изначальных ограничений ортогональности.
Комментарий к понятию «сплайн-вейвлет» В теории вейвлетов давно используется понятие «сплайн-вейвлет» [3,6,11]. Наиболее известными являются сплайн-вейвлеты Баттла-Ламарье [3,6,11]. Эти вейвлеты являются ортогональными, но сплай-
l =0
□
■
нами они, в действительности, не являются. Название «сплайн-вейвлеты» они получили потому, что процесс их построения связан с цепочками кратномас-штабного анализа, состоящими из пространств функций-сплайнов: в качестве масштабирующей функции берутся 5-сплайны [3,4,6,11] с узлами в целых числах. Для масштабирующей функции в виде 5-сплайна нулевой степени получается сплайн-вейвлет в виде функции Хаара [3,4-6,11]. Для масштабирующих функций в виде B-сплайнов больших степеней в результате «ортонормализации» [3] вейвлеты перестают быть сплайнами, а их носители становятся бесконечными [3,6,11].
В отличие от ортогональных, биортогональные (полуортогональные) сплайн-вейвлеты с наименьшим компактным носителем действительно являются сплайнами и действительно имеют компактные носители [6,11]. Построение этих сплайн-вейвлетов также начинается c масштабирующей функции в виде 5-сплайна некоторой степени K, при этом результат - сплайн-вейвлет - оказывается сплайном той же степени K и имеет компактный носитель. Поэтому при вычислении скейлограмм для биорто-гональных сплайн-вейвлетов, наряду с известным быстрым алгоритмом (далее - БА) ДВП С.Малла, может быть использован и предложенный эффективный алгоритм.
В настоящей работе понятие «сплайн-вейвлет» задается определением 2. Заметим, что в качестве сплайна, в отличие от работ [6,11], не обязательно выступает полиномиальный сплайн.
3.2. Эффективный алгоритм вычисления
локального ДВП, порождаемый сплайн-вейвлетом с конечным носителем Естественным вопросом является вопрос о наличии для введенных сплайн-вейвлетов вычислительно эффективных алгоритмов ДВП. Учитывая, что сплайн-вейвлет является (обобщенным) сплайном и, следовательно, для него выполняются выводы второго раздела, имеем
Предложение 6. Пусть у (t) - сплайн-вейвлет
порядка K c (^+1)-м узлом, J - множество отличных от нуля масштабов, порождающих множество {y(at)} масштабированных функций исходного сплайн-вейвлета. Тогда для любой функции y(at) сплайн-вейвлета из множества
{у (at)} существует алгоритм AaR , порождаемый сплайн-представлением у (at), со сложностью u (ACR ), удовлетворяющей ограничениям ( M = supp (у (at)), N - длина обрабатываемого входного сигнала)
- для обобщенного сплайна:
и А )< ((5 + 2)(К +1) +1 — ^ ) (20)
- для полиномиального сплайна:
м А )< N-N7TT (5+!)(К+1)+1+К ^ ).(20')
Доказательство: следствие (13)-(15). ■
Естественно, алгоритм вычисления локального ДВП, порождаемый сплайн-вейвлетом, тем вычислительно эффективнее, чем ниже значения К и 5. Поэтому, в дополнении к традиционным для вейв-летов условиям, связанным с гладкостью и числом нулевых моментов, на конструируемые сплайн-вейвлеты можно наложить дополнительные требования. В частности, требование минимума сложности эффективного алгоритма вычисления локального ДВП для сплайна заданного порядка.
3.3. Примеры сплайн-вейвлетов с эффективными алгоритмами вычисления локального ДВП
Несколько примеров полиномиальных сплайн-вейвлетов с низкой вычислительной сложностью (до 10 операций на один отсчет рассчитываемой скейло-граммы при = ^^ =0,5) порождаемых ими алгоритмов вычисления локального ДВП приведены в табл. 1. Отметим, что в расчетах не учитывалось свойство симметрии или антисимметрии сплайн-вейвлета. Если это свойство учитывать (при аналитическом построении эффективного алгоритма), то приводимая вычислительная сложность алгоритма Лая , порождаемого сплайн-вейвлетом, будет ниже указанной в таблице примерно на %тиЛК /2 .
Из табл. 1 видно, что наименьшей вычислительной сложностью обладает алгоритм АО,, порождаемый сплайн-вейвлетом порядка К=0 с 5=2 интервалами. Этот сплайн-вейвлет является широко известным вейвлетом - функцией Хаара [3,5,6,11]. Вейвлеты на основе базиса Хаара являются орто-нормированными вейвлетами. Как следствие, для вычисления их ДВП можно использовать упомянутый в начале работы быстрый алгоритм ДВП С.Малла [6]. Заметим, что для базиса Хаара БА ДВП имеет наименьшую сложность среди всех БА ДВП и всех ортогональных вейвлетов. Сравнение БА ДВП для базиса Хаара с предложенным эффективным алгоритмом, порождаемым сплайн-представлением у (г), приведено в работе [5]. Как показало сравнение, эффективный алгоритм является вполне конкурентоспособным [5]: при определенных диапазонах требуемых (полезных) масштабов а ДВП порождаемый алгоритм Аа имеет меньшую вычислительную сложность, чем известный БА С.Малла.
Табл. 1 - Сплайн-вейвлеты, порождающие эффективные алгоритмы ДВП
K
<(Ar )
базисный вейвлет
K
'( Ar )
базисный вейвлет
4,5
6,5
О 0,5 1,0
-1,0 -О,
О 0,5 1,0
10
1,0 0,5
-0,5 -1,0
-1,0 -0,5 0 0,5 1,0
-1,0 -0,5
0,5 1,0
9,5
1,0 0,5
-0,5 -1,0
-1,0 -0,
0 0,5 1,0
-1,0 -0,5
0,5 1,0
1,0 0,5 0
-0,5 -1,0
-1,0 -0,5 0 0,5 1,0
-1,0 -0,
0 0,5 1,0
S
S
1
1
2
1
1
2
7
2
2
1
3
3
0
4
2
0
3
3
1
3.4. Основные свойства сплайн-вейвлетов и порождаемых ими эффективных алгоритмов ДВП Основные свойства введенных сплайн-вейвлетов и порождаемых ими эффективных алгоритмов вычисления локального ДВП суммированы в следующих положениях:
- для сплайн-вейвлетов с конечными носителями существуют (вычислительно) эффективные алгоритмы локального ДВП;
- сложность эффективного алгоритма вычисления локального ДВП для сплайн-вейвлета с конечным носителем не зависит от его масштаба;
- сплайн-вейвлеты обладают большим числом степеней свободы. Это число степеней свободы можно использовать для удовлетворения дополнительных (внешних) требований к сплайн-вейвлетам;
- при удовлетворении всех внешних требований можно построить сплайн-вейвлеты, которые порождают (для разных масштабов) алгоритмы c предельно низкой вычислительной сложностью (по отношению к сложности алгоритмов, порождаемых сплайн-вейвлетами того же порядка);
- эффективные алгоритмы вычисления локального ДВП, порождаемые сплайн-вейвлетами,
позволяют организовать расчет скеилограмм для масштабов, которые не обязательно соответствуют степеням целого числа. Это позволяет, при необходимости, дополнительно снизить сложность процесса вычисления скеилограмм.
3.5. Биортогональные полиномиальные сплайн-вейвлеты с конечным носителем
В завершении настоящего раздела приведем характеристики эффективного алгоритма вычисления локального ДВП, порождаемого широко известными биортогональными полиномиальными сплаин-вейвлетами с конечным носителем [3,6,11]. Как хорошо известно, для построения биортогонального полиномиального сплайн-вейвлета порядка к в качестве масштабирующей функции ф(/) следует использовать B-сплайн порядка к, определяемый выражением [4,11]:
B0 (t) =
[1, t е [0,1],
[0, t й [0,1];
1
Вк (Н Вк-1 ( -т)т, к > 1.
0
Размер носителя В-сплайна порядка к, очевидно, равен к+1. Изображения В-сплайнов порядков 0, 1 и 2 приведены на рисунках 1-3.
Далее, масштабирующие функции ф(), то есть В-сплайны, образуют базис пространства с текущим разрешением. Поэтому любой вейвлет /2) меньшего разрешения представим в виде линейной декомпозиции функций {ф( - п)} . Для вейвлета t)
текущего разрешения это представления имеет вид [3,6,10,11]:
(t ) = V* S(-1)"-1 h_ n_w(2t - n )
где
hn = | ф (t) ф (2t - n) dt.
(21)
(22)
Легко проверить, что при использовании в качестве масштабирующей функции ф() В-сплайна порядка к число ненулевых коэффициентов {кп} (22) в выражении (21) составит (3К+2). Следовательно, биортогональный сплайн-вейвлет ) является полиномиальный сплайном степени к с числом узлов (4К+3), определенным на носителе длиной (2К+1). Вид биортогональных полиномиальных сплайн-вейвлетов ) степеней 0,1 и 2 приведен на
рис. 1-3. Учитывая свойство гладкости для этих вейвлетов (исключая случай к=0), легко получить
выражение для максимального значения суммарного дискретного дефекта такого сплайна. А именно, учитывая число узлов и степень сплайна, имеем:
г Л3' к = 0,
Г = |к(4к+3), к > 0.
Подставляя этот результат в выражение (14) для сложности алгоритма АСК вычисления локального ДВП, порождаемого (полиномиальным) сплайн-представлением КИХ, имеем:
|3, к = 0,
i(Äcr)
N - M + 1
N [к(4к + 3) + к^, к > 0.
(23)
Для сплайн-вейвлетов до третьей степени сложность, очевидно, составит:
i(äcr)
N - M + 1 N
3, к = 0,
7 + ^add, к = 1,
22 + , к = 2
45 + , к = 3
(24)
операций на один отсчет вычисляемой скейлограм-мы. Существенно, что сложность эффективного алгоритма вычисления локального ДВП не зависит от масштаба вейвлета.
Выражения (23)-(24) можно улучшить «вручную», если учитывать равные по модулю множители, которые появляются в отсчетах неоднородности сплайн-вейвлета. Однако, такой детальный анализ эффективного алгоритма не является целью настоящей работы.
1,0 0,8 0,6 0,4 0,2 0
1.0 0,5 0
-0,5 -1,0
ф(0
-0 4 1 0 4 0 8 1 ,2
V(0
-1,2 -0,8
-0,4
0
0,4
Рис.1. Масштабирующая функция в виде В-сплайна нулевой степени и соответствующий ей (би)ортогональный сплайн-вейвлет
г
1,0 0,8 0,6 0,4 0,2 0
-0,4 0 0,4 0,8 1,2 1,6 2,0 2,4
1,0 0,6
0,2 -0,2 -0,6
-0,4 0 0,4 0,8 1,2 1,6 2,0 2,4
Рис.2. Масштабирующая функция в виде В-сплайна первой степени и соответствующий ей биортогональный сплайн-вейвлет
1,0 0,8
0,6
0,4
0,2 0
0,4 0,2 0 -0,2 -0,4
-3,0 -2,0 -1,0 0 1,0 2,0
Рис.3. Масштабирующая функция в виде В-сплайна второй степени и соответствующий ей биортогональный сплайн-вейвлет
Заключение
В работе предложен подкласс базисных вейвле-тов, для которых существует эффективный (рекурсивный) алгоритм вычисления локального ДВП. Показано, что в данный подкласс входят вейвлеты в виде обобщенных дискретных сплайнов. Приводится эффективный алгоритм вычисления локального ДВП, порождаемый соответствующим сплайн-вейвлетом, дано выражение для вычислительной сложности этого алгоритма.
Указаны преимущества предложенного эффективного алгоритма ДВП по сравнению с известным быстрым алгоритмом ДВП С.Малла. В частности, предложенный алгоритм для введенного подкласса вейвлетов не требует ортогональности или биортогональности вейвлетов и позволяет вычислять ДВП для произвольных (нецелых) масштабов. Приведены примеры некоторых сплайн-вейвлетов, порождающих эффективный алгоритм (локального) ДВП с малой сложностью (до 10 операций на отсчет). Приведено общее выражение для сложности эффективного алгоритма, порождаемого известными биорто-гональными полиномиальными сплайн-вейвлетами.
Благодарности Работа выполнена при поддержке:
• Российского фонда фундаментальных исследований (РФФИ), проект № 06-01-00616-а;
• Фонда содействия отечественной науке.
Литература
1 Андерсон, Дж.А. Дискретная математика и комбинаторика / Дж.А. Андерсон - М.: ИД «Вильямс», 2004.
2 Гельфонд, А. О. Исчисление конечных разностей /
A.О. Гельфонд - М.: Наука, 1967.
3 Добеши, И. Десять лекций по вейвлетам / И. Добеши
- Москва, Ижевск: НИЦ «Регулярная и хаотическая динамика», 2001.
4 Завьялов, Ю.С. Методы сплайн-функций / Ю.С. Завьялов, Б.И. Квасов, В.Л. Мирошниченко - М.: Наука, 1980.
5 Копенков, В.Н Быстрые алгоритмы локального дискретного вейвлет-преобразования с базисом Хаара /
B.Н. Копенков, В.В. Мясников // Труды научно-технической конференции с международным участием «Перспективные информационные технологии в научных исследованиях, проектировании и обучении "ПИТ-2006"», Самара, 2006. - С. 113-118.
6 Малла, С. Вейвлеты в обработке сигналов / С. Малла
- М.: Мир, 2005.
7 Мясников, В.В. Сплайны как средство построения эффективных алгоритмов локального линейного преобразования / В.В. Мясников // Компьютерная оптика, 2007. - Т. 31. - № 2. - С. 52-68.
8 Мясников, В.В. Эффективный алгоритм над множеством алгоритмов вычисления свертки / В. В. Мясников // Компьютерная оптика, 2006. - В. 29. - С. 78-117.
9 Нуссбаумер, Г. Быстрое преобразование Фурье и алгоритмы вычисления свёрток // Г. Нуссбаумер - М.: Радио и связь, 1985.
10 Петухов, А.П. Введение в теорию базисов всплесков. Учебное пособие // А.П. Петухов - СПб: Изд-во СПбГТУ, 1999.
11 Чуи, К. Введение в вейвлеты // К. Чуи - М.: Мир, 2001.
12 Holschneider, M. A real-time algorithm for signal analysis with help of the wavelet transform / M. Holschneider [and others] // Wavelets, Time-Frequency Methods and Phase Space, Chapter A. - Berlin: Springer-Verlag, 1989. P. 289-297.
-e /—s
\ л»
\ /
ф(0
0 0,5 1,0 1,5 2,0 2,5 3,0
