Научная статья на тему 'Вейвлеты с конечной замкнутой формой представления'

Вейвлеты с конечной замкнутой формой представления Текст научной статьи по специальности «Математика»

CC BY
69
13
i Надоели баннеры? Вы всегда можете отключить рекламу.
Журнал
Научное приборостроение
ВАК
RSCI
Область наук

Аннотация научной статьи по математике, автор научной работы — Малыхина Г. Ф., Меркушева А. В.

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

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

WAVELETS WITH A FINITE CLOSED FORM OF REPRESENTATION

Transformation of non-stationary signals in data systems for dynamic spectrum analysis is based on the use of wavelets (a set of scaling and wavelet functions forming the basis for different scale subspaces, group of filters and algorithms for organizing this set of components into a system). Creating the ordinary wavelet components involves infinite products and sums, hence they have no closed representation and need fundamental calculations for their implementation, including approximations. We consider the method, procedures and applied constructions for wavelets having a finite closed form of representation.

Текст научной работы на тему «Вейвлеты с конечной замкнутой формой представления»

ISSNÜ868-5886

НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2ÜÜ4, том 14, № 3, с. 86-92

ОРИГИНАЛЬНЫЕ СТАТЬИ

УДК 621.391;519.21 © Г. Ф. Малыхина, А. В. Меркушева

ВЕЙВЛЕТЫ С КОНЕЧНОЙ ЗАМКНУТОЙ ФОРМОЙ ПРЕДСТАВЛЕНИЯ

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

ВВЕДЕНИЕ

Почти все известные ортонормальные вейвлеты, кроме простейших (и сейчас практически не применяемых) вейвлетов Хаара и Шеннона, не могут быть выражены в конечной замкнутой форме. Под таким видом выражения имеется в виду представление вейвлетов в виде явной зависимости от элементарных (тригонометрических, экспоненциальных, дробно-рациональных) функций или от более специальных функций математической физики (гамма, Бесселя). Суть проблемы в том, что основная часть вейвлет-преобразований, которые позволяют анализировать структуру нестационарных сигналов в информационно-измерительных системах (ИИС), базируется на вейвлет-базисах (вейвлетах), получение которых связано с построением бесконечных произведений, их преобразованием Фурье, конечными аппроксимациями, табулированием передаточных функций фильтров и разработкой программных средств, которые включены в некоторые пакеты (в виде так называемых Шо1Ьох’ов) [1]. Однако существуют процедуры для получения вейвлетов в замкнутой форме (в упомянутом выше смысле), и эта форма может быть использована для их непосредственного вычисления. Процедуры позволяют получить как ортогональные вейвлеты, так и неортогональные (но интерполирующие) вейвлеты.

Вейвлет-преобразования (ВП) находят применение при решении широкого круга задач: от обработки нестационарных сигналов в ИИС (для получения динамического спектра) до ряда областей математики (уравнения в частных производных и операторные), физики (область элементарных частиц высоких энергий) и биологии (диагностика на основе анализа структуры процессов) [2-4]. Это связано с тем, что вейвлеты образуют базисы сигналов или пространственных и временных функ-

ций обширного класса, обладают уникальным свойством экономного представления и локализации любых особенностей сигнала, сжатия информации, воспроизведения исходного сигнала без искажений. Тем не менее почти все ортонормиро-ванные вейвлеты не могут быть выражены в замкнутой форме [5, 6]. Так, построение наиболее часто используемых вейвлетов Добеши (Daubechies), имеющих важное свойство ортонормальности и компактной области определения, получается на основе поиска тригонометрического полинома m(t) с определенными свойствами. Полином m(t) определяет частотную форму масштабирующей функции ВП, а степень m(t) — меру регулярности (гладкость, определяемую числом непрерывных производных) масштабирующей функции и связанной с ней вейвлет-функции. А именно эти две функции ВП определяют структуру так называемого многоразрешающего анализа (multi-resolution analysis) сигналов и функций, т. е. представление на разных уровнях масштаба в виде аппроксимаций и их остатка (в форме деталей). Масштабирующая функция (МФ) ^(t) получается как обратное преобразование Фурье (ПФ) от ее частотного представления ф :

ф(т) =П m(2_k т). (1)

•\/2п k=0

После того как масштабирующая функция найдена, вейвлет-функция (ВФ), порождающая вейвлет-базис, определяется по одному из выражений: (2) — для частотной формы ВФ или (3) — для временной формы:

у(<ю) = e ~т /2 уп(ю /2 + п )<р(ю /2), (2)

V(t) =2 (-1)n h-n+i (2t - n), (3)

П=-<ж

где \кп }—ТО в (3) являются коэффициентами в пред- °б°значений и функций, которые вхдцят в замкну-

п “ , ч _ , ч тые выражения аналитических вейвлетов.

ставлении т(ш) рядом Фурье: т(ю) =

: (1/2) £ hne

jnm

Поскольку явный вид обратного ПФ от бесконечного произведения (1), приведенного выше, получен быть не может, то замкнутой формы выражений для вейвлетов Добеши неизвестно.

Близкий (с вейвлетами Добеши) по значимости — это вейвлет Мейера (Meyer), образующей которого служит единственная функция (4), т. к. логика конструкции Мейера требует от нее бесконечной гладкости (дифференцируемости) и компактности области определения:

ф(ю) = jexp(-1/(1 - ю2)), |ю|<1;

I 0, иначе.

(4)

Обратное ПФ от ф(ю) в виде (4) тоже не представляется в замкнутой форме, поэтому и вейвлет Мейера не имеет конечной замкнутой формы.

Существование полноценного ортонормиро-ванного базиса ВП с аналитическим представлением МФ и ВФ, которые могут вычисляться явным образом через элементарные или хорошо известные специальные функции, было установлено Вальтером (Walter) [7, 8]. Им найдено два класса таких вейвлетов, названных вейвлетами приподнятого косинуса.

(Можно отметить, что один из видов этих вейвлетов использовался ранее в цифровой связи в виде импульсов, имеющих спектр приподнятого косинуса, для удаления межсимвольной интерференции)

Прототип масштабирующей функции первого класса определяется выражением

Фі(t) =

sin(1 - в)t + 4в t cos п(1 + в)t nt[1 - (4 в t )2]

(5)

где 0 < в < 1/3. Связанная с этой МФ основная вейвлет-функция имеет вид

^(t + 1/2) = sin п(1 + вУ - 4в t cos п(1 - в)t

(6)

ВСПОМОГАТЕЛЬНЫЕ ОБОЗНАЧЕНИЯ И ЭЛЕМЕНТЫ ПРИКЛАДНОЙ ТЕОРИИ ВЕЙВЛЕТОВ

Будем считать, что отсутствие пределов у интегралов соответствует их значениям — ~, + то;

^к означает суммирование по к в пределах

— ^, + то; символ ф — это ПФ от функции ф (т. е. ее представление в частотной форме); все сигналы и функции имеют ограниченную Ь2-норму (конечную энергию, интеграл квадрата модуля); принадлежность сигнала или функции к Ь1 означает интегрируемость ее модуля; Ъ — множество всех целых чисел (<0, 0 и >0).

ПФ 5(ю) сигнала 5(0 и его обратное преобразование Фурье связаны соотношениями

5(ш) =[¿'(¿)в]М^ и 5(0 = — [5(т)в 10*d« . (7)

* 2п 1

Непрерывная функция Б считается функцией отсчетов (или интерполирующей), если она удовлетворяет двум условиям:

1) S (n) = 5

n,0

n = 0,± 1,± 2, ... и

2) Существуют постоянные А и В >0 такие, что для всех значений ш

0 < A <£ | S(rn + 2nk )|2 < B <<

пt[(4вt)2 -1] sin2n(1 - в)t + 8в t cos 2п(1 + в)t пt[(4вt)2 - 1] .

Другой класс вейвлетов приподнятого косинуса рассмотрен ниже.

Для анализа использованных Вальтером методов и процедур целесообразно сначала остановиться на элементах ВП, ряде используемых ниже

В частности, эти условия выполняются для широко известной в области обработки сигналов 8тс-функции (функции отсчетов), которая имеет вид: 8т(я?)/(я?) и представляет из себя масштабирующую функцию вейвлета Шеннона. (Заметим, что условие 2) делает {(¿-п)}п=-ш +« так называемым базисом Рисса).

В основе вейвлет-преобразования лежат ВФ [6, 9], которые образуют универсальные базисы {ут,п(0}т,п , хорошо приспособленные к обработке нестационарных сигналов и анализу динамических спектров:

Щт,п(() = 2т/2у (2mt - п); т, п — целые, (8)

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

зиса), либо принадлежать к классу базисов Рисса

[5].

Обычно ВФ формируются посредством первоначального построения МФ. Эти МФ должны удовлетворять определенным свойствам, которые позволяют выполнять многоразрешающий анализ, т.е. представление сигнала на нескольких уровнях масштаба, являющихся мерой разрешающей способности наблюдения частотной структуры нестационарного сигнала. Строго говоря — это структура сигнала на нескольких уровнях масштаба его представления. Однако, поскольку параметр масштаба а = 1/ю, то это является и частотным представлением с гиперболической (1/ю) шкалой частот.

Свойства, которым должны удовлетворять МФ вейвлет-преобразования, могут быть представлены как во временной области (9):

| ф(0 ф(г - п)& = 30'„ ,

п — целое (п е Ъ, Ъ — решетка целых на вещественной оси);

ф(г) = 2 СкФ(2г - к) {Ск }kеZ е 12;

к

для любого сигнала / (г)

/т (г) = ЕпаМпФ(2 тг - п), /т (г) /(г) (9)

— так и в виде эквивалентного набора (10) свойств в частотной области:

^ |ф(ю + 2пк)|2 = 1,

ф(ю) =

,jkw/2

ф(ю/2) =

(10)

= т(ю/2) ф(ю/2), где т(ю/2) = 2 ^ к

с ке

jkm / 2

ф непрерывна при

2к ю = 0 и р(0) = 1.

По полученной МФ строится основная ВФ. ВФ тоже может быть определена как во временной области, так и в частотной. Во временной области для этого используется соотношение

V(г) = 2к (-1)к+1 С1-к (2г - к). (11)

В связи с формой зависимости ВФ от МФ это выражение называют уравнением растяжения. Частотная форма основной вейвлет-функции имеет вид

у(ю) = е]т12т(ю/2 + п)ф(ю) .

(12)

Частотная форма МФ вейвлета Шеннона является характеристической функцией интервала [-п, п]: ф(т) = Х[-ж п](т). Согласно выражению (8),

т(а/2) является 4п периодичным расширением от ф(ю), и, следовательно, выполняется соотношение

т

ю ю

1 + ф

2 2

0,

п < \ю | < 2п; иначе.

УЧУ

Это приводит к ВФ Шеннона в форме выражения sin 2п(t -1/2) - sin n(t -1/2)

V (t) = -

п(t -1/2)

Кроме приведенных элементов связи МФ и ВФ в теории вейвлетов и иллюстрирующего примера ВФ Шеннона, до выполнения процедуры получения вейвлетов с замкнутой формой мы предварительно сформулируем положение, которое определяет необходимое и достаточное условие того, чтобы некоторая функция Б была функцией отсчетов и как из частотной формы такой функции может быть получена МФ другого вейвлета. Более строго: для Б и Б еЬ1 Б будет являться функцией отсчетов (т.е. Б(п) = 5п0), если (и только если)

выполняется условие 2 Б(ю + 2кп) = 1. Причем

к =-^

при Б(ю) > 0 функция Б порождает ортонормаль-ное семейство функций { фп (г) = (г - п)}^, (рп (г) е Ь2, где (рп (г) является временной формой

масштабирующей функции ф(о) = ТБй.

ПРОЦЕДУРА, ПОРОЖДАЮЩАЯ ЗАМКНУТУЮ АНАЛИТИЧЕСКУЮ ФОРМУ ВЕЙВЛЕТОВ

В основе процедуры лежит функция к со свойствами: к > 0, четная (к(-х) = к(х)), нормированная

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

(| к( x)dx = 1) и имеет в качестве области определения интервал [-п/3, я/3]. (В связи с особой ролью функции к используют понятие к-процедуры). Функция к, проинтегрированная на интервале 2п, эквивалентна ПФ (или частотной форме) некоторой масштабирующей функции р(ю), т.е. имеет место соотношение

ШТ/1

ф (ю) = J h( x)dx,

(13)

где р(ю) >0 — четная, непрерывная функция, имеющая в качестве области определения интервал [-2я/3, 2п/3] и удовлетворяющая условию

Ш-П

2 Ф (ю + 2kn ) = 1.

(14)

k =-~

Обратное ПФ (временная форма ф(^) является масштабирующей функцией и одновременно функцией отсчетов. Соответствующая ей ВФ представляется с помощью выражения

у(t + 1/2) = 2 ф(2t) - ф(t).

(15)

Эта ВФ в силу (15) и приведенного выше условия является функцией отсчетов на полуцелых значениях аргумента

у(п + 1/2) = Ô о .

(16)

Описанная процедура для образования ф(а) может быть использована для получения частотной формы некоторой новой ортогональной МФ

фо(а):

Фо (ю) =VФ(ю) ,

(17)

Ф2 (t) =

sin п(1 - ß)t + sin п(1 + ß)t 2п t [1 + 2 ß t ]

(18)

Этот принцип порождения новых видов вейвлетов из уже полученных МФ можно реализовать путем использования в качестве первичных такие вейвлеты, как sinc-функцию и два вида приподнятого косинуса (5) и (18).

Первичный вейвлет — 8те-функция

Частотная форма МФ для sinc-функции представляет характеристическую функцию интервала [-п, п], а в(ю) может быть принята в виде кусочнопостоянной функции. Такая в(ю) не удовлетворяет условию непрерывности в нуле, но она (как легко проверить) может быть использована аналогично функции к, включая формулу (13). Соответствующая функция ш^ю) тоже будет кусочнопостоянной, 2п-периодичной и (при условии

| ш(т) |2 + | ш(т + п) |2 = 1) дополнительно удовлетворяет соотношению

| тв (ю)|2 + | тв (ю + п) |2 = 1.

(19)

где ф0(ю) удовлетворяет условиям, которые необходимы для МФ вейвлет-преобразования и выражены в частотной форме (10).

Метод такого подхода введен Лиу (Liu) (см. [10] и [11]) и позволил ему получить масштабирующие функции, которые являются одновременно и ортогональными, и функциями отсчетов, т.е. удовлетворяют условию второго соотношения в (10) ("уравнению сжатия") и условию (14). Одна из полученных таким образом масштабирующих функций в замкнутой форме [7] приведена выше (5), вторая может быть выражена соотношением

Причем имеет место встречная импликация (тогда и только тогда):

| шв (т) |2 + | шв (т + п) |2 = 1 о

о | ш(т) |2 + | ш(т + п) |2 = 1.

Это позволяет вывести ряд вейвлетов, имеющих замкнутую форму выражения.

Общий набор процедур

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

• Для целого п > 0 вводится набор полуоткрытых (справа) интервалов:

ОРТОГОНАЛЬНЫЕ ВЕЙВЛЕТЫ В ЗАМКНУТОЙ ФОРМЕ

Ортогональные вейвлеты более общей формы могут быть получены на основе метода Вальтера [10], который ввел способ модификации ортогональной МФ ф (t), позволяющий получить новый ортогональный вейвлет-базис масштабирующей функции. Им показано, что, используя МФ ф (t), ее ПФ ф(ю) и нечетную, непрерывную в нуле функцию в(а>), для которой евсо является 2п-периодичной, можно образовать новую МФ фв (а) = е]всоф(а). Эта частотная форма фв (а) новой МФ будет удовлетворять условиям (8), а обратное ПФ от нее ф e(t) является также ортогональной МФ. Т.е. фo(t) образует ортогональный базисный набор { фв (t - n)}„eZ .

I о =

п ^

0,— , Ik =

s 2

п

2

n-k+1

2

n-k

(2о)

к = 1,2,..., п.

• Определяется д1(ю) на интервале [0, п]:

в(ю) = £к для юе 1к, к = 1,2,...,п.

• В качестве в (ю) берется нечетное продолжение в^ю) на интервал [-п, п] и его периодическое продолжение на всю вещественную ось (-да, да).

Определенная таким образом функция в(ю) удовлетворяет условию (19) и порождает фв (т) и фв(0. Причем {фв ^ - п)}пе2 является ортонор-мальным базисным набором на Ь2 (с элементами базиса, интегрируемыми с квадратом), а фв (0 удовлетворяет уравнению растяжения (второе со-

ж

отношение в условиях (9)).

Если в качестве ек используются значения, кратные п (например, 0 и п), то фв (т) — вещественная. И как можно показать [12], в этом случае функция швю/2) — четная, и представляет собой расширение на 4п-периодичную модификацию, получаемую из функции тв (т / 2) вида

mд (ю /2) =

j(st -sk-i)

Цg (ю) = e

-0}^/2

m.

ю ю

+ п Фд —

2 2

V J

= еЛт/2-6(ю+2п)+6(т/2+п)] m

ю л ю

+ п Фв —

2 2

J

= e j[^/2-0(m+2n)+e(m/2+n)] e^l^ уп

f .

Ю —+ П

v 2

v j

2

- j(®/2)

Ц(ю) =

eStak, если A;

/7 ^‘/''ТТТТ

0,

если иначе;

A:

B:

2--

2--

V

где ao

<| ю |< п

<| ю |< 2п ,

= ej (st -st -i) k =

k = 1,2,..., n;

1, ak = e (k k-1), k = 1, 2, ... , n.

+

2 (ak - ak-1) sin ak (t - 1/2) -

k =1

n+1

(t -1/2)},

где ak = п

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

f 2 n-k+2 ___1 ^

^ n-k+1

k = 1,2,..., n + 1.

на Io;

на Ik, k = 1,2,..., n; (21) на [п, 2п).

Частотное представление вейвлет-функции (ее ПФ) может быть получено с помощью (12) и цепочки соотношений:

(22)

С учетом указанных выражений преобразование Фурье ВФ определяется формулой

(23)

Процедура получения основной вейвлет-функции заканчивается выполнением обратного преобразования Фурье от у(ю). После проведения ряда необходимых для этого преобразований получается ВФ в замкнутой форме в виде выражения

, ч ± 1

V ({) =-----------х

п(Г -1/2)

х 2п^ - 1/2) +

ПРИМЕНЕНИЕ ПРОЦЕДУР ВАЛЬТЕРА ДЛЯ ПОЛУЧЕНИЯ ВЕЙВЛЕТ-ФУНКЦИЙ В ЗАМКНУТОЙ ФОРМЕ

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

сложное построение. Его общность определяется возможностью широкого выбора для величины n и значений параметров ek. В связи с недавним появлением самого метода совокупность ВФ этого типа пока не исследована; объем всего класса обширен, но строго не установлен. Тем не менее представляется полезным проследить схему построения вейвлетов в замкнутой форме с помощью рассмотренных процедур для небольших значений n и величин параметров ek, которые наиболее удобны при расчете.

Для n = 1, 10 = [0, я/2), 11 = [я/2, я) положим £0 = 0, тогда, если выбрать значение £1 = 0, то можно получить sinc-функцию. Значение е1 = п при прежней величине остальных параметров процедуры (в соответствии с выражением ak = e1 (Sk -Sk-1) и соотношением из (24)) дает a0 = 1; a1 = -1; а1 = = 3п/2 и а2 =п Поэтому вейвлет-функция для таких условий будет иметь вид

Ц (t) = -^ х п(t -1/2)

х {sin 2п(t -1/2) + sin п(t -1/2) - 2 sin[(3n/2)(t -1/2)] }

Для n = 2, 10 = [0, я/4), 11 = [п/4, я/2) и 12 = [п/2,

я) положим £0 = 0, тогда имеется три возможных случая.

— При значениях е1 = 0 и е2 = п получается рассмотренная выше ВФ.

— При е1 = е2 = п получаются оценки a0 = 1, a1 = -1, a2 = 1 и тогда, согласно (23), получается выражение

e

-1 (®/2)

-1 при 7п/4 < | ю | < 2п;

1 при 3п/2 < | ю |< 7п /4;

-1 при п < | ю |< 3п /2;

0 в других случаях.

(24)

Поскольку для этого варианта параметров получаются значения а = 7п/4; а = 3п/2 и а3 = п, то непосредственным преобразованием может быть определена ВФ в виде выражения

e

0

1

1

ж

n-k

n - k +1

2

2

п

Ц (t) =----------х

п(t -1/2)

х {sin п(t -1/2) - 2 si^^ / 2)(t - 1/2)] +

+ 2 ял[(7п / 4)(t -1/2)] - sin 2п(t -1/2)}.

— При e1 = п; e2 = 0 получаются оценки a0 =1, a1 = a2 = -1.Тогда на основе (23) определяется выражение для частотной формы ВФ

e

- j (œ/2)

V (ш) =

і

при 7п /4 <| ш |< 2п;

-1 при

П<| ш |< 7п/4;

0

в других случаях.

И аналогично предыдущему варианту путем обратного ПФ получаем ВФ

V (t ) =

1

п (t -1/2)

х {sin 2п(Y -1/2) + sin п(t -1/2) -

- 2sin[(7n / 4)(t - 1/2)]} .

ЗАКЛЮЧЕНИЕ

По принятой в области обработки сигналов (Signal processing) терминологии вейвлетами называют совокупность компонент вейвлет-преобразования: масштабирующую функцию и основную (материнскую) вейвлет-функцию (порождающую вейвлет-базис с помощью аффинного преобразования временного или пространственного аргумента ВФ) и концепцию мультимасштабно-го анализа — анализа частотной структуры исследуемого сигнала на различных уровнях масштабов. По логике построения ВФ (т.е. последовательности МФ — га-функция — ВФ) структура вложенных подпространств различного масштаба формируется путем дискретного бинарного растяжения оси времени у МФ и создания локального базиса в каждом из таких подпространств набором МФ с определенным шагом временного смещения. Технология построения включает математические преобразования с бесконечными произведениями и суммами, использование аппроксимаций, вычислительной фазы, организации алгоритмов с программной их реализацией и использование в Toolbox прикладных систем Matlab, Matead, Mathematica. В связи со столь сложной методологией построения у всех вейвлетов (кроме двух примитивных видов) отсутствуют конечные замкнутые выражения.

В отличие от вейвлетов, доступных для практического использования только в форме пакетов указанных выше систем, сравнительно недавно (94-98 г.) разработаны вейвлеты с замкнутой фор-

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

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

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

1. Дьяконов В.П., Абраменкова И.О. Matlab 5.0/5.3. Система символьной математики. СПб.: Питер, 1999. 633 с.

2. Кратиров Д.В., Меркушева А.В. Алгоритм, основанный на вейвлет-преобразовании и нейронной сети, для бесконтактного измерения параметров газожидкостного потока // Сборник докладов Международной конференции "Датчики и системы". СПб.: Изд. СПбГТУ, 2002. Т. 3. С.51-55.

3. Малыхина Г.Ф., Меркушева А.В. Вейвлет-фильтрация нестационарного сигнала с адаптацией на основе нейронной сети // Труды Международной конференции по мягким вычислениям и измерениям SCM-2001. СПб.: СПбГЭТУ, 2001. Т. 1. С. 239-242.

4. Дремин И.М., Иванов О.В., Нечитайло В.А. Вейвлеты и их использование // Успехи физических наук. 2001. Т. 171, № 5. С. 465-501.

5. Daubechies I. The wavelet transformation, time-frequency localization, and signal analysis // IEEE, Transactions on Information Theory. 1990. V. 36. P.961-1005.

6. Чуи К. Введение в вейвлет. М.: Мир, 2001. 412 с.

7. Walter G., Zhang J. Orthonormal wavelets with simple closed-form expressions // IEEE, Transactions on Signal Processing. 1998. V. 46. P. 22482251.

8. Walter G. Translation and diliation invariance in orthogonal wavelets // Applied computations and harmonic analysis. 1994. V. 10. P. 344-349.

9. Меркушева А.В. Классы преобразований нестационарного сигнала в информационноизмерительных системах. Время-масштабные (вейвлет-) преобразования для спектральновременного анализа // Научное приборостроение. 2002. Т. 12, № 3. С. 68-82.

10. Walter G. Wavelets and other orthogonal systems

X

with applications. CRC Press, Boca Raton, FL,

1996. 230 p. Санкт-Петербург

11. Hernandez E., Weiss G. A Course on Wavelets.

CRC рге^ Boca Б'Ц 1996 189 p. Материал поступил в редакцию 11.03.2004.

WAVELETS WITH A FINITE CLOSED FORM OF REPRESENTATION

G. F. Malykhina, A. V. Merkusheva

Saint-Petersburg

Transformation of non-stationary signals in data systems for dynamic spectrum analysis is based on the use of wavelets (a set of scaling and wavelet functions forming the basis for different scale subspaces, group of filters and algorithms for organizing this set of components into a system). Creating the ordinary wavelet components involves infinite products and sums, hence they have no closed representation and need fundamental calculations for their implementation, including approximations. We consider the method, procedures and applied constructions for wavelets having a finite closed form of representation.

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