§ Обратное вейвлет–преобразование. Основные свойства вейвлета
Теоретически
доказано, что если сигнал представлен своим вейвлет–преобразованием
, где
, то этот сигнал может
быть восстановлен по формуле:
, (1)
где – нормализующий коэффициент,
подобный
в
преобразовании Фурье;
- спектр Фурье вейвлета
.
Требуемое
свойство конечности коэффициента ограничивает возможный класс функций
, которые можно использовать в качестве вейвлетов.
Следует отметить, что прямое вейвлет преобразование
(2)
и обратное вейвлет
преобразование (1) правомерны, если вейвлеты образуют ортогональную систему базисных
функций
,
для выбранных разных значений
.
Определим основные свойства, которыми должен обладать вейвлет:
1)
Локализация во времени и по частоте, то есть вейвлет определен на
конечном временном интервале (например, ) и имеет ограниченный по частоте спектр
Фурье
.
2) Нулевое среднее значение
.
Замечание. Считается, что вейвлет обладает лучшими свойствами по
выделению небольших кратковременных изменений сигнала, если он имеет несколько
нулевых моментов m-го порядка ()
.
3) Ограниченность по энергии
.
4) Автомодельность базиса функций , составленного из
вейвлета
. Самоподобие
(похожесть) функций
одна
на другую, так как имеют одинаковое число осцилляций, что и
- базисный вейвлет. Это
определяется тем, что функции
получены из
путем масштабирования и
сдвига.
Помимо
рассмотренного вейвлета Хаара часто используются и иные более гладкие
вейвлеты:
a)
Wave-вейвлет , который
получается как первая производная от функции Гаусса:
.
Спектр Фурье wave-вейвлета
.
b)
Mhat-вейвлет или
«мексиканская шляпа» , получаемый как вторая производная функции
Гаусса
.
Спектр Фурье Mhat-вейвлета
.
В настоящее время используется несколько десятков различных вейвлетов.
Выбор типа (формы) анализирующего вейвлета, определяется тем, какую информацию необходимо более полно выявить из сигнала.
§ Дискретное вейвлет–преобразование
По аналогии с преобразованием Фурье, для непрерывного прямого (2) и обратного (1) вейвлет преобразования определяют их дискретные варианты.
Дискретное прямое вейвлет преобразование определяется формулой:
. (3)
В
случае перехода к дискретным отсчетам функции
, то в формуле (3) интеграл
заменяется суммой
.
Обратное дискретное вейвлет преобразование определяется формулой
, (4)
где - дискретный вейвлет
.
Для дискретных вейвлетов выполняется свойство ортогональности
и биортогональности
,
где - функция Дирака.
Замечание.
Прямое и обратное вейвлет
преобразование реализованы в модуле Wavelet-анализа пакета Matlab.
При этом пользователю представляется широкий набор используемых вейвлетов .
Замечание. Для прямого и обратного вейвлет преобразования разработаны быстрые алгоритмы вычислений, подобные БПФ, что позволяет значительно ускорить расчет вейвлет-спектров.
Замечание.
Свойство ортогональности вейвлетов может не выполняться
для выбранной в качестве вейвлета
функции. В этом случае прямое
преобразование будет осуществляться с использованием
, а обратное с помощью
специальных реконструкционных вейвлетов
, которые отличаются от
и могут быть найдены
по ним.
Такие вейвлеты иногда называют парными (dyadic wavelet) и он имеет следующую реконструкционную формулу (обратного преобразования)
.