86102 (Частотно-временной анализ сигналов)

Череповецкий военный инженерный институт радиоэлектроники.

Кафедра №8

Курсовая работа по математике

Тема: «Частотно-временной анализ сигналов»

Выполнили:

Плотников Е.А.

г. Череповец-2008.

Содержание

1. Плоскость частота-время

2. Базисные функции частотно-временного анализа

3. Прямое и обратное преобразование Фурье

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

4.1 Дискретизация масштаба

4.2 Дискретизация масштаба и сдвига. Фреймы

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

Литература

1. Плоскость частота-время

Для анализа и сравнения частотно-временных локализационных свойств различных базисов используют плоскость частота-время. Любая функция может характеризоваться интервалом I_t на временной оси и интервалом в Фурье области, в которых содержится 90% ее энергии, сосредоточенной около центра тяжести функции Тогда в этой плоскости функцию можно изобразить в виде прямоугольника, как показано на рис. 3.1.

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

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

Энергия такой функции E:

Следовательно, ширина функции равна . В соответствии со свойством масштабирования Фурье-преобразования ( ) Влияние масштабирования на положение функции в плоскости время-частота показано на рис. 3.3.

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

2. Базисные функции частотно-временного анализа

Итак, частотно-временной анализ предназначен для выявления локальных частотно-временных возмущений сигнала. Вследствие кратковременности таких возмущений, сам сигнал может рассматриваться как заданный в L₂ т.е. для одномерных сигналов – на всей действительной оси с нормой . Следовательно, базисные функции, которые получили название вейвлетов, также должны принадлежать L₂ и быстро убывать при Тогда, чтобы перекрыть такими базисными функциями все возможные временные положения сигнала, необходимо, чтобы базисные функции представляли собой набор смещенных во времени функций. Удобнее всего, если этот набор образуется из одной и той же "материнской" функции (прототипа), сдвинутой по оси t т.е. Чтобы обеспечить частотный анализ, базисная функция должна иметь еще один аргумент – масштабный коэффициент, который является аналогом частоты в Фурье-анализе. Тогда базисные функции для частотно-временного анализа будут иметь вид

Где масштабный коэффициент введен как делитель t, причем масштабированию подвергается также и сдвиг b. Это позволяет сохранить относительную "плотность" расположения базисных функций по оси t при расширении или сжатии самой функции и при (рис 3.6)

Таким образом, базисные функции для частотно-временного анализа должны обладать следующими свойствами. Ограниченность, т.е. принадлежность L₂ . Локализация. Базисные функции вейвлет - анализа, в отличие от преобразования Фурье, должны быть локализованы, т.е. определены на конечном интервале как во временной, так и в частотной областях. Для этого достаточно, чтобы выполнялись условия:

Нулевое среднее. Равенство нулю нулевого момента

или, что иногда необходимо – равенство нулю момента m-го порядка

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

3. Прямое и обратное преобразование Фурье

При

- прямое преобразование Фурье

- обратное преобразование Фурье.

Комплексная функция имеет смысл спектральной плотности, ее иногда называют непрерывным спектром Фурье-функции f(t).

Также как и в случае периодической функции, предполагается, что f(t) удовлетворяет условиям Дирихле или, что эквивалентно, абсолютно интегрируема и удовлетворяет условию Дини.

Отметим также, что:

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

Представление функции f(t) через ее непрерывное вейвлет – преобразование является избыточным. В задачах обработки информации, встречающихся на практике, сигнал, во-первых, имеет ограниченную полосу и, во-вторых, допускаются те или иные погрешности в получаемых результатах. Поэтому используют дискретное представление непрерывных сигналов, при которых параметры преобразования, в данном случае a и b, приобретают дискретные значения. Вейвлет-преобразование, при котором значения a и b дискретны, называют дискретным вейвлет-преобразованием (DWT - Discrete Wavelet Transform).

4.1 Дискретизация масштаба

Рассмотрим сначала случай дискретного масштаба a и положим . Это равноценно разбиению частотной оси на поддиапазоны (частотные полосы). Предположим, что (это можно сделать всегда, умножив функцию ψ на некоторый модуляционный множитель (см. ). Тогда частотное окно будет равно:

а центральная частота m-го вейвлета:

.

Базисом для DWT является функция, полученная из

( )

при :

.

Если справедливо и если

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

(3.5.2.)

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

, (3.5.3)

где константы А и В такие, что . Условие (3.5.3.) в терминах радиотехники имеет довольно прозрачное толкование. Действительно, так как при каждом значении масштаба вейвлет представляет собой полосовой фильтр, то набор (сумма) этих фильтров (блок фильтров) является некоторым устройством с неравномерной частотной характеристикой, определяемой константами A и B (рис. 3.12). Сигнал, например звуковой, на выходе такого устройства при сильной неравномерности частотной характеристики претерпевает существенные искажения. Поэтому для его восстановления принимают специальные меры, в частности, устанавливают фильтр, компенсирующий искажения частотной характеристики. В вейвлет-преобразовании таким фильтром является дуальный (или двойственный) вейвлет , Фурье-образ которого имеет вид:

. (3.5.4.).

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

( )

и формулу получим (3.5.4.):

Из (3.5.4.) и (3.5.3.) можно показать, что

4.2 Дискретизация масштаба и сдвига. Фреймы

В этом случае полагают дискретными величины a и b, т.е. Частотное окно для анализа сохраняется прежним. Ширина временного окна

равна , а среднее значение изменяется дискретно пропорционально m -ой степени a₀ - масштабу вейвлета. Чем уже функция ψ, т.е. меньше величина , тем меньше (на ту же величину) шаг сдвига этой функции. Базисными функциями для дискретного вейвлет-преобразования будут функции, получаемые из ,при и

Коэффициенты разложения любой функции из L₂ могут быть получены как

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

для всех если константы A и B такие, что В этом случае формула для восстановления функции f(t) по коэффициентам будет иметь вид

(3.5.8)

где ошибку восстановления 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₀ = 2. Такое разбиение называют логарифмическим, так как отношение верхней и нижней границ диапазонов постоянно и равно 2. Такое разбиение является еще и идеальным, так как оно реализуется идеальными полосовыми фильтрами. Подобная идеализация нужна для исследования свойств частотного разложения с помощью идеализированных вейвлетов, что позволит в дальнейшем перейти к более сложным разложениям. Любой сигнал со спектром может занимать полосу частот, охватывающую несколько таких поддиапазонов.

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