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

Главные признаки вейвлет-преобразования одномерного сигнала - его представления в виде обобщенного ряда или интеграла Фурье по системе базисных функций. Построение вейвлет-преобразования в системе Matlab, график функции и ее дискретное преобразование.

Рубрика Программирование, компьютеры и кибернетика
Вид курсовая работа
Язык русский
Дата добавления 18.02.2019
Размер файла 2,4 M

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

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

Размещено на http://www.allbest.ru//

Размещено на http://www.allbest.ru//

Федеральное агентство железнодорожного транспорта

Омский государственный университет путей сообщения (ОмГУПС (ОМИИТ)

Кафедра «Автоматика и системы управления»

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

Пояснительная записка к курсовому проекту

по дисциплине «Методы и алгоритмы цифровой обработки информации»

Омск 2016

Задание

Необходимо построить сигнал , его спектр (ДПФ) и Вейвлет-преобразование.

Реферат

Курсовой проект содержит 15 страниц, 6 рисунков, 6 формул, 4 источника.

Сигнал, спектр, ДПФ, Вейвлет-преобразование, функция.

Объектом исследования является Вейвлет-преобразование и его моделирование в среде Matlab.

Пояснительная записка выполнена в программе MicrosoftWord 2010, расчеты выполнены в среде Matlab R2013b.

Введение

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

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

Вейвлеты. Общие замечания

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

(1)

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

Рисунок 1 - Вейвлет мексиканская шляпа (а) и модуль его спектральной плотности (б)

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

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

Рисунок 2 - Спектры вейвлетов

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

Главные признаки вейвлета

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

Ограниченность. Квадрат нормы функции должен быть конечным:

(2)

Локализация. ВП в отличие от преобразования Фурье использует локализованную исходную функцию и во времени, и по частоте. Для этого достаточно, чтобы выполнялись условия:

(3)

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

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

(4)

Из этого условия становится понятным выбор названия «вейвлет» - маленькая волна.

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

Часто для приложений бывает необходимо, чтобы не только нулевой, но и все первые моментов были равны нулю

(5)

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

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

Сопоставление с преобразованием Фурье

Классическое преобразование Фурье (ПФ) является традиционным математическим аппаратом для анализа стационарных процессов. При этом сигналы разлагаются в базисе косинусов и синусов или комплексных экспонент. Эти базисные функции простираются вдоль всей оси времени. С практической точки зрения и с позиций точного представления произвольных сигналов ПФ имеет ряд ограничений и недостатков. Обладая хорошей локализацией по частоте, оно не обладает временным разрешением. ПФ даже для одной заданной частоты требует знания сигнала не только в прошлом, но и в будущем, а это - теоретическая абстракция. Обусловлено это тем, что базисной функцией при разложении в ряд Фурье является гармоническое колебание, которое математически определено на временном интервале от - ? до + ?. ПФ не учитывает, что частота колебания может изменяться во времени. Локальные особенности сигнала (разрывы, ступеньки, пики и т.п.) дают едва заметные составляющие спектра, по которым обнаружить эти особенности, и тем более их место и характер, практически невозможно. В этом случае невозможно и точное восстановление сигнала из-за появления эффекта Гиббса. Для получения о сигнале высокочастотной информации с хорошей точностью следует извлекать ее из относительно малых временных интервалов, а не из всего сигнала, а для низкочастотной спектральной информации - наоборот. Кроме того, на практике не все сигналы стационарны, а для нестационарных сигналов трудности ПФ возрастают многократно. Часть указанных трудностей преодолевается при использовании оконного ПФ:

(6)

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

Рисунок 3 - Локальная во времени функция, перемещаемая вдоль оси времени

Если окно, показанное на рисунке 3, перемещать скачками (через ф) вдоль всего времени существования сигнала , то за некоторое число таких перемещений возможен «просмотр» всего сигнала. Так что вместо обычной спектрограммы получится набор спектрограмм, схематично представленный в виде прямоугольников на рисунке 4, а. Такой спектральный анализ равносилен анализу с помощью набора фильтров с постоянной шириной полосы пропускания, равной .

Рисунок 4 - Набор спектрограмм в виде прямоугольников

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

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

Построение вейвлет-преобразования в системе Matlab

Необходимо построить сигнал, его спектр(ДПФ) и Вейвлет-преобразование.

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

Для построения спектра будем использовать стандартную функцию Matlab fft, учитывая все вышесказанное, получим:

t = 0:0.01:10;

w=t;

y = sin(w.*t);

subplot(2,1,1),plot(y), grid ;

Y=fft(y);

subplot(2,1,2),plot(t,abs(fft(y))), grid ;

На рисунке 5 изображена функция и ее дискретное преобразование.

Рисунок 5 - Функция (а) и ее дискретное преобразование (б)

Для построения Вейвлет преобразования будем использовать функцию cwt. вейвлет преобразование сигнал

Функция выполняет континуальное преобразование сигнала и возвращает массив спектральных коэффициентов преобразования.

Входными параметрами функции являются следующие:

«wavelet» - идентификатор имени материнской вейвлет-функции (ядра преобразования);

scales - вектор значений масштабирующей переменной базисного вейвлета.

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

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

Функция, записанная в третьей форме, обладает более широкими возможностями визуализации, поскольку содержит расширение mode - опцию, позволяющую указать вид «подкрашивания» спектра. В частности, значение «lvl» позволяет улучшить контрастность изображения спектральных коэффициентов сигнала, значение «glb» - затенить спектральный шум коэффициентов, сглаживая, таким образом, цветность общей спектральной картины, значения «abslvl» или «lvlabs», а также «absglb» или «glbabs» - отконтрастировать и сгладить цветовую гамму спектра соответственно. При этом использование расширения «absglb» является эквивалентным включению в функцию cwt опции «plot».

Спектральную картину вейвлет-преобразования можно представить также аксонометрически, указав в опции mode префикс «3D», например, «3Dplot», «3Dlvl» и т.д.

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

scales = [1: 10, 25: 35];

cwt(y,1:40, 'mexh','abslvl',[100 400]);

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

Рисунок 6 - Сигнал и его вейвлет-преобразование

Заключение

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

Библиографический список

Matlab.Exponenta[Электронный ресурс] / Список функций Wavelet Toolbox. - Режим доступа: http://matlab.exponenta.ru/wavelet/book1/4/cwt.php

Смоленцев Н.К. Основы теории вейвлетов. Вейвлеты в MATLAB. - М.: ДМК Пресс, 2005. -304 с. 

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

Новиков Л.В. Основы вейвлет-анализа сигналов. Учеб.пособие. - СПб.: Изд-во 000 «МОДУС», 1999. - 152 с.

Размещено на Allbest.ru

...

Подобные документы

  • Применение вейвлет-преобразования для сжатия и обработки медицинских сигналов и изображений. Разработка алгоритма автоматизированного выделения PQRST-признаков в сигнале электрокардиограмм с помощью вейвлет-инструментария математического пакета Matlab.

    дипломная работа [4,6 M], добавлен 16.07.2013

  • Характеристика сигнала и его представление в виде математического ряда. Условия ортогональности двух базисных функций. Ряд Фурье, его интегральное преобразование и практическое использование в цифровой технике для обработки дискретной информации.

    реферат [69,9 K], добавлен 14.07.2009

  • Разработка функции вычисления дискретного преобразования Фурье от входного вектора. Исследование свойств симметрии ДПФ при мнимых, четных и нечетных входных сигналах. Применение обратного преобразования Фурье для генерации периодической функции косинуса.

    лабораторная работа [228,8 K], добавлен 13.11.2010

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

    курсовая работа [1021,4 K], добавлен 29.10.2017

  • Исследование простейших радиотехнических сигналов, разложение их в ряд Фурье. Построение амплитудных спектров синуса, суммы синусов и синка. Создание в среде программирования Matlab программ с параметрами: длина сигнала, амплитуда, частота дискретизации.

    лабораторная работа [990,4 K], добавлен 23.11.2014

  • Анализ таблиц сопряженности и коэффициента сопряженности Крамера. Выявление структуры нечисловых данных. Определение эмпирического среднего с помощью медианы Кемени. Очистка тестового сигнала от шума с использованием дискретного вейвлет-преобразования.

    контрольная работа [408,8 K], добавлен 23.12.2016

  • Общие сведения о языке программирования Matlab. Функции работы с векторами и матрицами. Операторы условных переходов. Построение двумерных графиков. Построение гистограммы изображения. Функции его преобразования и зашумления, метрики определения качества.

    лабораторная работа [853,5 K], добавлен 25.10.2015

  • Графики вещественнозначных функций вещественного переменного. Построение двумерных графиков. Пример построения графика синусоиды. Пример использования функции subplot. Оформление двумерных графиков в системе MatLab. Основные функции оформления графиков.

    курсовая работа [826,3 K], добавлен 30.04.2019

  • Сигнал как некоторое средство для передачи информации. Знакомство с параллельными алгоритмами двумерного быстрого преобразования Фурье, анализ способов вычисления. Общая характеристика процессора Power5 64-bit RISC. Рассмотрение функций библиотеки MPI.

    дипломная работа [1,6 M], добавлен 09.10.2013

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

    дипломная работа [3,0 M], добавлен 01.07.2008

  • Анализ проблем, возникающих при совмещении изображений в корреляционно-экстремальных навигационных системах. Использование двумерного дискретного преобразования Фурье. Нахождение корреляционной функции радиолокационного и моделируемого изображений.

    дипломная работа [3,6 M], добавлен 07.07.2012

  • Генерирование и сохранение мелодии в виде звукового файла формата wav. Проведение частотного анализа полученного сигнала. Зависимость объема wav-файлов от разрядности кодирования сигнала. Спектр нот записанного wav-файла с заданной разрядностью.

    лабораторная работа [191,0 K], добавлен 30.03.2015

  • Особенности работы в режиме командной строки в системе Matlab. Переменные и присваивание им значений. Комплексные числа и вычисления в системе Matlab. Вычисления с использованием функции sqrt. Неправильное использование функций с комплексными аргументами.

    дипломная работа [1,9 M], добавлен 30.07.2015

  • Основные понятия стеганографии. Атаки на стегосистемы, стегосистемы водяных знаков. Применение дискретных вейвлет преобразований в кодировании цифровых зображений. Алгоритмы стеганографического встраивания информации в изображения формата JPEG2000.

    дипломная работа [3,5 M], добавлен 09.06.2013

  • Описание архитектуры процессора TMS320C25. Моделирование фильтра в модуле FDATool программной среды Matlab. Алгоритм нерекурсивной фильтрации сигнала. Расчет массива отсчетов входного сигнала. Моделирование фильтра при различных частотах входного сигнала.

    курсовая работа [119,2 K], добавлен 14.06.2015

  • Проектирование информационной системы (ИС) преобразования данных с помощью математических и алгоритмических подходов. Автоматизированная ИС преобразования измеренных значений сил и моментов в расчетные случаи для виртуальной модели автомобиля для ОММиР.

    курсовая работа [2,6 M], добавлен 25.12.2011

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

    курсовая работа [191,3 K], добавлен 07.02.2010

  • Обработка детерминированного сигнала. Классификация измерительных сигналов. Формула исходного сигнала. Построение спектра амплитуд и спектра фаз. Точность спектрального анализа. Нормальный закон распределения. Спектральный анализ случайного сигнала.

    курсовая работа [616,8 K], добавлен 07.07.2013

  • Використання методів обробки сигналів, які базуються на використанні малохвильової теорії. Вимоги до алгоритмів компресії та критерії порівняння алгоритмів. Застосування вейвлет-перетворень. Критерії оцінювання оптимальності вибору малохвильових функцій.

    реферат [1,1 M], добавлен 26.05.2019

  • Математическое описание, алгоритм и программа вычисления определенного интеграла методом трапеций. Расчет n-значений исследуемой функции и вывод их в виде таблицы. Технические и программные средства. Входные и выходные данные, функциональное назначение.

    курсовая работа [21,0 K], добавлен 03.01.2010

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.