Скачать

Вейвлет-анализ сигналов и его применение


Курсовая работа

на тему:

«Вейвлет-анализ сигналов и его применение»



1. Идея и возможности вейвлет-преобразования

вейвлет преобразование редактирование дискретный

Вейвлет-технологии начали серьёзно развиваться в 80–90 годы прошлого века, хотя первый тип вейвлета был описан ещё в 1909 году учёным Хааром. Многие типы и семейства вейвлетов были названы именами учёных, которые внесли большой вклад в разработку теоретических основ вейвлетов: Мейер, Добеши, Маллат.

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

Ряд Фурье использует в качестве базиса синусоиды, которые предельно локализованы в частотной области (вырождаются на спектрограмме в вертикальную линию), и вообще не локализованы во временной области.

Противоположный пример – импульсная базисная дельта-функция d(t).Она чётко локализована во временной области и потому идеально подходит для представления разрывов сигнала. Но она не несёт информации о частоте сигнала и потому плохо приспособлена для представления сигналов на заданном отрезке времени.

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

Вейвлеты характеризуются своим временным и частотным образом.

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


Сравнение представления сигналов в различных областях

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

График коэффициентов Фурье (например, полученный с помощью команды fft) этого сигнала не показывает ничего особенно интересного: плоский спектр с двумя пиками, представляющими одну частоту. Однако график вейвлет коэффициентов ясно показывает точное расположение во времени рассмотренного выше разрыва.


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

2. Свойства вейвлетов

Вейвлет («короткая волна», «всплеск») – это волновая форма сигнала эффективно ограниченной длительности, которая имеет среднее значение ноль.

Сравним вейвлет с синусоидальной волной, которая является основой анализа Фурье. Синусоиды не имеют ограниченной длительности – они продолжаются от минус до плюс бесконечности. И где синусоиды гладкие и предсказуемые, вейвлеты стремятся быть неровными и асимметричными.


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

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

Математически процесс анализа Фурье представлен преобразованием Фурье:

,

которое является суммой по всему времени сигнала f(t) умноженного на комплексную экспоненту.

Результатами этого преобразования являются коэффициенты ФурьеF(w), умножение которых на синусоиду соответствующей частоты даст синусную компоненту исходного сигнала. Графически этот процесс выглядит так:


(Сигнал) (Преобразование Фурье) (Синусные компоненты исходного сигнала)

Аналогично, непрерывное прямое Wavelet-преобразование определяется как сумма по всему времени сигнала, умноженного на масштабируемые, сдвинутые версии вейвлет функции:

,

где y(t) Wavelet-функция, f(t) – сигнал.

Результатом НВП будет вейвлет коэффициенты С(t, a), которые являются функцией позиции tи масштаба a.

Умножением каждого коэффициента С на соответственно масштабируемый и сдвинутый вейвлет получают непосредственные вейвлеты исходного сигнала:


Сигнал Вейвлет Вейвлеты

преобразование исходного сигнала

Амплитудно-временное представление нестаціонарного сигнала и его результат непрерывного вейвлет преобразования

Масштабирование вейвлета просто означает его растяжение (или сжатие).

Вводится понятие – масштабный коэффициент, который обозначают буквой а. Если речь идет, например, о синусоидах, то эффект от масштабного коэффициента очень легко увидеть:


Чем больше частота, тем более сжатая синусоида.

Масштабный коэффициент действует и на вейвлеты. Чем меньше масштаб, тем более «сжатым» будет вейвлет.

Из диаграмм видно, что для синусоиды sin(wt) масштаб аобратно пропорционален частоте w. Аналогично, с вейвлет-анализом, масштаб обратно пропорционален частоте сигнала.

Сдвиг вейвлета просто означает задержку или ускорение его фронта. Математически задержка функции на k представляется в виде:


Вейвлет функция  Сдвинутая вейвлет функция

2. Семейства вейвлет-функций

Можно привести несколько ярких представителей семейств вейвлет-функций

Haar

Daubechies


Morlet-преобразование

MexicanHat-преобразование

Непрерывное прямое вейвлет-преобразование

Для создания НВП необходимо выполнить следующих пять шагов:

1. Взять вейвлет и установить его на начальный интервал исходного сигнала.

2. Вычислить коэффициент С, который показывает как тесно коррелированны вейвлет и сигнал на этом интервале. Высокое значение С означает большую схожесть. Заметьте, что результаты будут зависеть от формы вейвлета, выбранного Вами.

3. Сдвинуть вейвлет вправо и повторять шаг 1 и 2 до тех пор, пока Вы не исследуете весь сигнал.

4. Масштабировать (растянуть) вейвлет и повторить шаги 1 – 3.

5. Повторить шаги 1 – 4 для всех масштабов.

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

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

Трёхмерное представление результатов расчёта – графики коэффициентов напоминают вид сверху на неровную (ухабистую) поверхность.


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

Из графиков видно, что чем выше масштаб, тем «протяженнее» вейвлет. Чем протяженнее вейвлет, тем длиннее часть сигнала, с которой он сравнивается, и более крупные черты сигнала будут измерены вейвлет коэффициентами.

Таким образом, есть связь между масштабом вейвлет и частотой, как показано вейвлет анализом:

Малый масштаб а Þ Сжатый вейвлет Þ быстро изменяющиеся составляющие Þ высокая частота w.

Большой масштаб а Þ Растянутый вейвлет Þ медленно изменяющиеся, крупные черты Þ низкая частота w.

Непрерывное обратное вейвлет-преобразование

Обратное непрерывное вейвлет-преобразование осуществляется по формуле реконструкции во временной области. Одна из форм может быть представлена


где f(t) – восстановленный сигнал, y(t)– вейвлет-функция, С(t, a) – вейвлет коэффициенты, которые являются функцией позиции tи масштаба a, Kyкоэффициент, зависящий от выбора вейвлет-функции,R область ограничения сигнала.

3. Дискретное вейвлет преобразование

Дискретное вейвлет преобразование одномерного сигнала

Цифровая обработка сигнала требует его дискретизации. Как и в случае преобразования Фурье существует дискретная форма вейвлет преобразования. Выше было отмечена определенная степень свободы в выборе базиса вейвлет преобразования. В данном разделе нами будет использоваться один из самых простых вейлвет базисов – базис Хаара.

Рассмотрим дискретизированный и квантованный сигнал (сигнал 2.1) – рисунок 2.1 (а). Будем постепенно усреднять данный сигнал, усредняя попарно его отсчеты. Таким образом, каждый шаг усреднения будет сокращать разрешение сигнала в 2 раза (т.е. для его представления будет требоваться в два раза меньшее число отсчетов). Однако при таком усреднении мы теряем часть информации о сигнале, для того чтобы восстановить сигнал после усреднения нам потребуется дополнительная информация. Будем сохранять разности между усредненным отсчетом и отсчетами, из которых усредненный отсчет состоит при более высоком разрешении. Данные разности показывают детали сигнала – его флуктуации вокруг среднего при данном уровне разрешения. На рисунке 2.1 детализирующее коэффициенты показаны в правой части рисунков 2.1 (б, в, г, д). Теперь воспользовавшись детализирующими коэффициентами мы сможем восстановить прежнею форму сигнала.

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

Заметим, что сигнал 2.1, изображенный на рисунке 2.1 (а), может быть представлен следующим образом – http://gmdidro.googlepages.com/formula8.gif, где -http://gmdidro.googlepages.com/formula8_1.gif некоторые базисные функции, а http://gmdidro.googlepages.com/formula8_2.gif – координаты сигнала 2.1 в этом базисе. Очевидно, что если мы выберем в качестве http://gmdidro.googlepages.com/formula8_3.gif единичную ступеньку, изображенную на рисунке 2.2 (а), то, сдвигая http://gmdidro.googlepages.com/formula8_3.gif необходимое число раз, мы сможем представить сигнал 2.1 с помощью суммы таких единичных ступенек. Таким образом, мы ввели базис, в котором мы можем представить сигнал 2.1. Отметим, что поскольку функцииhttp://gmdidro.googlepages.com/formula8_1.gif, изображенные на рисунке 2.2, не пересекаются между собой, то построенный нами базис является ортогональным. Функции http://gmdidro.googlepages.com/formula8_1.gif называются масштабирующими функциями.

Рисунок 2.1 – Усреднение дискретизированного сигнала 2.1


Рисунок 2.2

Теперь необходимо ввести некоторый базис для представления детализирующих коэффициентов. Такой базис был введен Хааром и его базисные функции, названные вейвлетами, изображены на рисунке 2.2 (б).

Рассмотрим теперь процедуру усреднения сигнала, проиллюстрированную на рисунке 2.1, с точки зрения только что введенных базисов. Рассмотрим конкретный сигнал, заданный следующим вектором значений – (9 7 3 5). С помощью масштабирующей функции Хаара мы можем представить сигнал так, как это изображено на рисунке 2.3.

Рисунок 2.3 – Представление исходного сигнала в базисе Хаара

Проведем процедуру декомпозиции сигнала на две части – усредненный сигнал с двое уменьшенным разрешением и детализирующие коэффициенты. Получим следующий вектор – http://gmdidro.googlepages.com/formula9.gif = (8 4 | 1 –1), представление которого в базисе Хаара с помощью масштабирующей функции и вейвлетов изображено на рисунке 2.4.

Рисунок 2.4 – Представление усредненного сигнала в базисе Хаара

Выделим в векторе (8 4 | 1 –1) часть, представляющую усредненный сигнал (первая половина вектора), и проведем относительно неё повторное усреднение и нахождение детализирующих коэффициентов. Получим следующий вектор – (6 | 2 1 –1) представление которого в базисе Хаара с помощью масштабирующей функции и вейвлетов изображено на рисунке 2.5.

Рисунок 2.5 – Представление дважды усредненного сигнала в базисе Хаара

Таким образом, мы представили исходный сигнал с помощью его усредненной части (среднего по сигналу) и детализирующих коэффициентов. Отметим, что размерность исходного и преобразованного векторов совпадают, это говорит о том, что при преобразовании не было потерь информации и, следовательно, возможно полное восстановление исходного вектора. Шаги описанной процедуры ещё раз проиллюстрированы на рисунке 2.6.

Рисунок 2.6 – Представление сигнала в базисе Хаара

Отметим, что если мы будем восстанавливать сигнал после его разложения в базисе Хаара, то мы можем остановить процесс восстановления «на полпути» и получить представление сигнала с заданным разрешением. Другими словами нами получен математический инструмент изменения разрешения сигнала.

Дискретное вейвлет преобразование изображения

Рассмотрим вейвлет преобразование изображения. Общая идея вейвлет преобразования многомерных сигналов заключается в декомпозиции многомерного сигнала до одномерных сигналов и, последующего их вейвлет преобразования с композицией результатов. Существуют два метода такого преобразования – стандартное и нестандартное вейвлет преобразование. Этими методы различаются порядком применения вейвлет преобразования я к декомпозированным одномерным сигналам. Ниже будет рассмотрено нестандартное вейвлет преобразование.

При нестандартном вейвлет преобразовании изображения вейвлет преобразование попеременно применяется то к строкам, то к столбцам изображения. Иллюстрация этого метода представлена на рисунке 2.7.

Рисунок 2.7 – Нестандартное вейвлет преобразование изображения

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

На рисунке 2.7 показан также псевдокод рекурсивного применения DWT к изображению. При этом на каждом шаге преобразования удобно представлять изображение никак матрицу, а как вектор.


4. Применения дискретного вейвлет преобразования

Сжатие изображений. JPEG 2000

Основная идея, используемая при сжатии сигналов с помощью вейвлет преобразования, заключается в том, чтобы отбрасывать детализирующие коэффициенты, значения которых близки к нулю. В 1999 году был разработан новый стандарт сжатия изображений, названный JPEG 2000 и призванный заменить стандартный алгоритм сжатия JPEG. Одним из основных отличий JPEG 2000 от JPEG является изменение основной процедуры преобразования изображения. В то время как в JPEG использовалось преобразование Фурье и в JPEG 2000 используется вейвлет преобразование, что позволило не только улучшить визуальное качество изображения, но и добавить некоторую интересную функциональность, принципиально не достижимую в JPEG.

Одной из таких дополнительных функции является ROI (Region of Interest). ROI позволяет динамически в пространстве и во времени повышать разрешение изображения. Под динамическим повышением разрешения изображения в пространстве понимается то, что мы можем повысить разрешение только выделенной области изображения. Под динамическим повышением разрешения изображения во времени понимается то, что мы можем повышать разрешение выделенной области изображения постепенно, шаг за шагом.

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

Поиск изображений по образцу

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

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

Основная идея метода заключается в описании каждого изображения с помощью 20 наибольших детализирующих коэффициентов его вейвлет разложения. Эти двадцать коэффициентов называются ярлыком изображения («ключевыми словами» изображения в базисе вейвлетов) и именно по ним ведется поиск в базе данных.

Рисунок 3.1. Иллюстрация работы алгоритма поиска изображения по образцу


Многомасштабное редактирование

Как уже неоднократно подчеркивалось, основу различных применений вейлетов составляет возможность простого и быстрого изменения разрешения сигнала, преобразованного с помощью DWT. Но эта черта вейвлетов негде так не очевидна как при многомасштабном редактировании изображений и трехмерных моделей. Дело в том, что при многомасштабном редактировании изменение разрешения редактируемого объекта происходит интерактивно, что особенно хорошо выявляет описанные преимущества вейвлетов. Рисунок 3.6 иллюстрирует идею многомасштабного редактирования.

вейвлет преобразование редактирование дискретный

Рисунок 3.2 – Многомасштабное редактирование трехмерной модели

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


Заключение

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

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


Перечень использованных источников

1. И.М. Дремин, О.В. Иванов, В.А. Нечитайло. Вейвлеты и их использование. – Успехи физических наук, 2001

2. Wavelet Digest – www.wavelet.org

3. Br. Vidakovic, P. Mueller. Wavelets for kids – Duke University.

4. А. Переберин. Многомасшабные методы синтеза и анализа изображений – Москва, 2001.

5. А. Петров. Вейвлеты и их приложения – Рыбинск, РГАТА 2007