Курсовая работа: Частотно-временной анализ сигналов
![](/file/7d5aacc146456be8f3642.jpg)
⚡ 👉🏻👉🏻👉🏻 ИНФОРМАЦИЯ ДОСТУПНА ЗДЕСЬ ЖМИТЕ 👈🏻👈🏻👈🏻
Череповецкий военный инженерный институт радиоэлектроники.
Тема: «Частотно-временной анализ сигналов»
2. Базисные функции частотно-временного анализа
3. Прямое и обратное преобразование Фурье
4. Дискретное вейвлет-преобразование
4.2 Дискретизация масштаба и сдвига. Фреймы
4.3 Примеры дискретного вейвлет-преобразования
Для анализа и сравнения частотно-временных локализационных свойств различных базисов используют плоскость частота-время. Любая функция может характеризоваться интервалом I t
на временной оси и интервалом в Фурье области, в которых содержится 90% ее энергии, сосредоточенной около центра тяжести функции Тогда в этой плоскости функцию можно изобразить в виде прямоугольника, как показано на рис. 3.1.
Очевидно, что смещение функции на от исходного состояния вызовет перемещение прямоугольника параллельно оси t. Модуляция этой функции комплексной экспонентой сдвигает прямоугольник параллельно оси (рис.3.2.).
Масштабирование функции (ее сжатие или растяжение) приводит к развороту прямоугольника. Действительно, получим новую функцию масштабированием функции на коэффициент a:
Следовательно, ширина функции равна . В соответствии со свойством масштабирования Фурье-преобразования ( ) Влияние масштабирования на положение функции в плоскости время-частота показано на рис. 3.3.
В качестве примеров функций, иллюстрирующих эффективность их представления в плоскости время-частота, рассмотрим функцию Дирака и Фурье - базис. Известно, что - функция является идеальным базисом для временного анализа сигналов. Результатом такого анализа являются отсчеты, которые можно рассматривать как временной спектр сигнала. На плоскости время - частота функция выглядит как показано на рис. 3.4а, т.е. эта функция обладает свойством хорошей временной локализации, но плохой локализацией в спектральной области (она имеет равномерный спектр на всех частотах). Базисные функции Фурье-анализа, наоборот, обладают хорошей частотной локализацией в то время, как во временной области они имеют бесконечную протяженность (см. рис. 3.4б).
2. Базисные функции частотно-временного анализа
Итак, частотно-временной анализ предназначен для выявления локальных частотно-временных возмущений сигнала. Вследствие кратковременности таких возмущений, сам сигнал может рассматриваться как заданный в L 2
т.е. для одномерных сигналов – на всей действительной оси с нормой . Следовательно, базисные функции, которые получили название вейвлетов, также должны принадлежать L 2
и быстро убывать при Тогда, чтобы перекрыть такими базисными функциями все возможные временные положения сигнала, необходимо, чтобы базисные функции представляли собой набор смещенных во времени функций. Удобнее всего, если этот набор образуется из одной и той же "материнской" функции (прототипа), сдвинутой по оси t т.е. Чтобы обеспечить частотный анализ, базисная функция должна иметь еще один аргумент – масштабный коэффициент, который является аналогом частоты в Фурье-анализе. Тогда базисные функции для частотно-временного анализа будут иметь вид
Где масштабный коэффициент введен как делитель t, причем масштабированию подвергается также и сдвиг b. Это позволяет сохранить относительную "плотность" расположения базисных функций по оси t при расширении или сжатии самой функции и при (рис 3.6)
Таким образом, базисные функции для частотно-временного анализа должны обладать следующими свойствами. Ограниченность
, т.е. принадлежность L 2
. Локализация
. Базисные функции вейвлет - анализа, в отличие от преобразования Фурье, должны быть локализованы, т.е. определены на конечном интервале как во временной, так и в частотной областях. Для этого достаточно, чтобы выполнялись условия:
Нулевое среднее.
Равенство нулю нулевого момента
или, что иногда необходимо – равенство нулю момента m-го порядка
Это – вейвлеты m-го порядка, позволяющие анализировать более тонкую структуру сигнала, подавляя медленно изменяющиеся его составляющие.
3. Прямое и обратное преобразование Фурье
Комплексная функция имеет смысл спектральной плотности, ее иногда называют непрерывным спектром Фурье-функции f(t).
Также как и в случае периодической функции, предполагается, что f(t) удовлетворяет условиям Дирихле или, что эквивалентно, абсолютно интегрируема и удовлетворяет условию Дини.
4. Дискретное вейвлет-преобразование
Представление функции f(t) через ее непрерывное вейвлет – преобразование является избыточным. В задачах обработки информации, встречающихся на практике, сигнал, во-первых, имеет ограниченную полосу и, во-вторых, допускаются те или иные погрешности в получаемых результатах. Поэтому используют дискретное представление непрерывных сигналов, при которых параметры преобразования, в данном случае a и b, приобретают дискретные значения. Вейвлет-преобразование, при котором значения a и b дискретны, называют дискретным вейвлет-преобразованием (DWT - Discrete Wavelet Transform).
Рассмотрим сначала случай дискретного масштаба a и положим . Это равноценно разбиению частотной оси на поддиапазоны (частотные полосы). Предположим, что (это можно сделать всегда, умножив функцию ψ на некоторый модуляционный множитель (см. ). Тогда частотное окно будет равно:
а центральная частота m-го вейвлета:
Базисом для DWT является функция, полученная из
Если справедливо и если достаточно быстро затухает, то любая функция из L 2
может быть представлена в виде дискретной по последовательности
Для восстановления f(t) по дискретным значениям (3.5.2.) на базис (t) налагаются дополнительные ограничения, а именно, образ Фурье вейвлета (t) должен удовлетворять соотношению
где константы А и В такие, что . Условие (3.5.3.) в терминах радиотехники имеет довольно прозрачное толкование. Действительно, так как при каждом значении масштаба вейвлет представляет собой полосовой фильтр, то набор (сумма) этих фильтров (блок фильтров) является некоторым устройством с неравномерной частотной характеристикой, определяемой константами A и B (рис. 3.12). Сигнал, например звуковой, на выходе такого устройства при сильной неравномерности частотной характеристики претерпевает существенные искажения. Поэтому для его восстановления принимают специальные меры, в частности, устанавливают фильтр, компенсирующий искажения частотной характеристики. В вейвлет-преобразовании таким фильтром является дуальный (или двойственный) вейвлет , Фурье-образ которого имеет вид:
Покажем, что с помощью такого вейвлета по коэффициентам DWT полностью восстанавливается сигнал. Действительно, используя соотношение Парсеваля
Из (3.5.4.) и (3.5.3.) можно показать, что
4.2 Дискретизация масштаба и сдвига. Фреймы
В этом случае полагают дискретными величины a и b, т.е. Частотное окно для анализа сохраняется прежним. Ширина временного окна
равна , а среднее значение изменяется дискретно пропорционально m -ой степени a 0
- масштабу вейвлета. Чем уже функция ψ, т.е. меньше величина , тем меньше (на ту же величину) шаг сдвига этой функции. Базисными функциями для дискретного вейвлет-преобразования будут функции, получаемые из ,при и
Коэффициенты разложения любой функции из L 2
могут быть получены как
Выражение (3.5.6) является дискретным вейвлет-преобразованием функции . Чтобы обратное преобразование во временную область было справедливым, должно выполняться следующее условие:
для всех если константы A и B такие, что В этом случае формула для восстановления функции f(t) по коэффициентам будет иметь вид
где ошибку восстановления R можно оценить как Разделив все члены неравенства (3.5.7) на , можно видеть, что константы A и B являются границами нормированной на энергии – скалярного произведения . Они (эти константы) как бы "обрамляют" нормированную энергию коэффициентов Отсюда произошел термин фрейм (frame), которым называют множество функций при которых условие (3.5.7) выполняется. Если A= B , то и множество называют плотным фреймом. При этом выражение вытекающее из (3.5.7), является обобщением теоремы Парсеваля на плотные фреймы. Для плотных фреймов из (3.5.8) получаем
Если A=B=1, то плотный фрейм становится ортогональным базисом. Заметим, что для вейвлетов, образованных материнским вейвлетом (3.3.6), хорошие результаты при восстановлении сигналов получаются при так как . Для больших величин, например будет т.е. восстановление приводит к большим искажениям.
4.3 Примеры вейвлетов для дискретного преобразования
Как было отмечено выше, функции вейвлет обладают свойством частотно-временной локализации, т.е. они ограничены как в частотной, так и во временной областях. Ниже рассмотрим два примера: первый – спектр вейвлетов в частотной области представляет собой идеальный полосовой фильтр, второй – сами функции вейвлет представляют собой прямоугольники. Все вейвлеты, с точки зрения частотно-временных свойств, занимают промежуточное положение между этими крайними случаями.
Sinc-базис. Разобьем ось частот на интервалы (поддиапазоны), как показано на рис. 3.13 при a 0
= 2. Такое разбиение называют логарифмическим, так как отношение верхней и нижней границ диапазонов постоянно и равно 2. Такое разбиение является еще и идеальным, так как оно реализуется идеальными полосовыми фильтрами. Подобная идеализация нужна для исследования свойств частотного разложения с помощью идеализированных вейвлетов, что позволит в дальнейшем перейти к более сложным разложениям. Любой сигнал со спектром может занимать полосу частот, охватывающую несколько таких поддиапазонов.
Тогда и т.е. сигнал представляет собой сумму некоторого числа элементарных сигналов. В рассматриваемом идеальном случае частотные каналы не перекрываются, поэтому имеет место ортогональность этих элементарных сигналов, т.е.
Выберем из всего множества сигналов такие, которые ограничены полосой частот 2 I
, т.е. имеющие спектр . Рассмотрим периодическую функцию такую, что: , т.е. полученную периодизацией F 1
(ω) (рис. 3.14)
Тогда спектр функции: F i
(ω) при произвольном I можно представить в виде:
Посмотрим, как при этих условиях можно представить функцию f (t) во временной области. Для этого разложим периодическую функцию с периодом , в ряд Фурье (см. ):
Где, подставляя (3.5.10а) в (3.5.9) и выполняя обратное преобразование Фурье, получим:
Вычислим первый интеграл. Переставляя операции суммирования и интегрирования и ограничивая пределы интегрирования с учетом функции окна, получим:
Выражение (3.5.13) является представлением функции f (t) в базисе вейвлет. В рассматриваемом частном случае идеальной полосовой фильтрации вейвлетом является функция (3.5.14), образованная из материнской функции по (3.5.15) с учетом (3.5.12). Такой вейвлет называется sinc –вейвлетом по имени функции (3.5.12), которая его образует, а функция (3.5.12) получила название масштабной функции.
Множитель при необходим для сохранения нормы вне зависимости от величины масштаба, так как:
Покажем, что в рассматриваемом частном случае т.е. определяется отсчетами функции при . Рассмотрим интеграл Фурье ( ) при дискретных значениях функции , заданной на интервале Имеем, с учетом (3.5.10б):
Последнее равенство справедливо при и вещественных
Выполнив преобразование Фурье выражения (3.5.14), можно видеть, что спектр Фурье sinc -вейвлета представляет собой идеальный полосовой фильтр, в общем случае занимающий полосу частот от до
Вейвлет Хаара. Разобьем теперь временную ось на интервалы, как показано на рис. 3.17 и определим на единичном интервале функцию
Эта функция является материнским вейвлетом, так как она удовлетворяет условию ( ). Система сдвигов таких функций образует ортонормальный базис, так как их взаимная энергия равна нулю при и равна единице при
Преобразование Фурье ( ) вейвлета Хаара имеет вид и показано на рис. 3.17б.
Функции Хаара, также как sinc -вейвлет, могут быть получены с помощью масштабной функции
Из приведенных примеров следует ряд интересных выводов:
1. Представление вейвлет-функции в виде прямоугольников в любой из областей (частотной или временной) ведет к бесконечному расширению в противоположной области. Следовательно, для того, чтобы функции вейвлет были локализованы одновременно во временной и частотной областях, они должны убывать с ростом аргумента, по крайней мере, по закону обратной пропорциональности (см.( и )).
2. Вейвлеты ψ(t), спектры Фурье которых представляют собой полосовые фильтры, могут быть выражены через масштабные функции (t), спектры Фурье которых представляют собой фильтры нижних частот (см. формулы (3.5.15) и (3.5.19)).
3. Базисные функции для DWT могут быть получены из одной материнской функции путем ее масштабирования и сдвига (см. формулы (3.5.14) и (3.5.15)).
4. Любой сигнал f(t) из L 2
может быть представлен своим вейвлет- разложением (3.5.13), если число компонентов f i
(t) таково, что они занимают полосу частот большую, чем полоса сигнала.
1. Новиков И.Я., Стечкин СБ. Основы теории всплесков // Успехи математических наук. 1998. V. 53. № 6. С.9-13.
2.Петухов А.П. Введение в теорию базисов всплесков. СПб.: Изд. СПбГТУ, 1999. 131 с.
3.Воробьев В.И., Грибунин В.Г. Теория и практика вейвлет-преобразования. СПб.: ВУС, 1999. 203 с.
4. Астафьева Н.М. Вейвлет-анализ: основы теории и примеры применения// УФН . 1996. Т. 166, № 11. С. 1145-1170.
5. Martin Vatterli, Jelena Kovačevic. Wavelets and Subband Coding. Prentice Hall, New Jersey, 1995.
Название: Частотно-временной анализ сигналов
Раздел: Рефераты по математике
Тип: курсовая работа
Добавлен 15:05:24 16 ноября 2010 Похожие работы
Просмотров: 350
Комментариев: 15
Оценило: 3 человек
Средний балл: 5
Оценка: неизвестно Скачать
Срочная помощь учащимся в написании различных работ. Бесплатные корректировки! Круглосуточная поддержка! Узнай стоимость твоей работы на сайте 64362.ru
Привет студентам) если возникают трудности с любой работой (от реферата и контрольных до диплома), можете обратиться на FAST-REFERAT.RU , я там обычно заказываю, все качественно и в срок) в любом случае попробуйте, за спрос денег не берут)
Да, но только в случае крайней необходимости.
Курсовая работа: Частотно-временной анализ сигналов
Курсовая Применение Производной К Исследованию Корней Уравнений
Контрольная работа по теме Классификация объектов контроллинга. Роль и место анализа в системе стратегического контроллинга
Реферат по теме Искусство древних цивилизаций Востока. Общие черты и особенности
Курсовая работа: Политика государственных расходов в переходной экономике. Скачать бесплатно и без регистрации
Реферат: Форсирование Днепра
Эссе Социальная Психология
Реферат: Value Of Literary Studies Essay Research Paper
Контрольная работа по теме Основы управления инновационными процессами в школе
Дипломная работа по теме Техническое обслуживание ГАЗ-33075
Реферат На Тему Вселение Камчатского Краба В Баренцево Море: Итоги И Перспективы
Моих Целей Мне Мешает Достичь Эссе
Хромосомная Теория Наследственности Реферат
Курсовая Работа На Тему Методика Вивчення Теми "Короткі Історичні Відомості Про Обчислювальні Прилади Та З Історії Створення Еом"
Реферат по теме Лечение пролежней
Дубровский Разбойник Или Жертва Сочинение План
Курсовая работа: Проектирование микропроцессорной системы на основе микроконтроллера К1816ВЕ31
Реферат по теме Философские воззрения Дж. Д. Сэлинджера
Лермонтов Бородино План Для Сочинения 5 Класс
Мини Сочинение О Красоте Растений
Темы Дипломных Работ По Плаванию
Реферат: Поведенческий маркетинг
Реферат: Программное обеспечение персональных компьютеров
Реферат: Шпаргалки