Оконное преобразование Фурье

Материал из Википедии — свободной энциклопедии
Перейти к: навигация, поиск

Оконное преобразование Фурье — это разновидность преобразования Фурье, определяемая следующим образом:

F(t,\omega) = \int\limits_{-\infty}^\infty f(\tau) W(\tau-t) e^{-i\omega \tau}\,d\tau,

где W(\tau-t) — некоторая оконная функция. В случае дискретного преобразования оконная функция используется аналогично:

 F(m,\omega) = \sum_{n=-\infty}^{\infty} f[n]w[n-m]e^{-j \omega n}

Существует множество математических формул визуально улучшающих частотный спектр на разрыве границ окна. Для этого применяются преобразования: прямоугольное (никакое), треугольное, сужающийся косинус, фрагмент синусоиды, синус в кубе, синус в 4-ой степени, преобразование Парзена, Велша, Гаусса, Хеннинга, Хэмминга, Чебышева, с пульсациями, Розенфилда, Блэкмана-Харриса, горизонтальное и с плоской вершиной. Также существует методика по взаимному перекрытию окон, при этом обычно можно выбрать сколько сэмплов из предыдущего окна будет усреднено с текущим окном.

Применение[править | править вики-текст]

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

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

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

Частотно-временное разрешение[править | править вики-текст]

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

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

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

Типы оконных функций[править | править вики-текст]

Прямоугольное окно[править | править вики-текст]

Прямоугольное окно; B=1.00
w(n)=\left\{ \begin{matrix}
   1 , & n\in[0,N-1] \\
   0, & n\notin[0,N-1]\\
\end{matrix} \right.

Получается автоматически при ограничении выборки N отсчетами. Максимальный уровень боковых лепестков частотной характеристики: -13 дБ.

Окно Ханна (Ханнинга)[править | править вики-текст]

Окно Ханна; B = 1.50
w(n) = 0.5\; \left(1 - \cos \left ( \frac{2 \pi n}{N-1} \right) \right)

где N — ширина окна. Уровень боковых лепестков: -31.5 дБ.

Окно Хэмминга[править | править вики-текст]

Окно Хэмминга
w(n)=0.53836 - 0.46164\; \cos \left ( \frac{2\pi n}{N-1} \right)

Уровень боковых лепестков: -42 дБ.

Окно Блэкмана[править | править вики-текст]

Окно Блэкмана; α = 0.16; B=1.73
w(n)=a_0 - a_1 \cos \left ( \frac{2 \pi n}{N-1} \right) + a_2 \cos \left ( \frac{4 \pi n}{N-1} \right)
a_0=\frac{1-\alpha}{2};\quad a_1=\frac{1}{2};\quad a_2=\frac{\alpha}{2}\,

Уровень боковых лепестков: -58 дБ (α=0.16).

Окно Кайзера[править | править вики-текст]

Окно Кайзера, α =2; B=1.5
Окно Кайзера, α =3; B=1.8
w(n) = \frac { | I _0  \left ( \beta \sqrt {1 - \left ( \frac{2n-N+1}{N-1} \right )^2} \right ) | } { | I _0 (\beta) | }

где I_0 — модифицированная функция Бесселя первого рода нулевого порядка; \beta — коэффициент определяющий долю энергии, сосредоточенной в главном лепестке спектра оконной функции. Чем больше \beta тем больше доля энергии, и шире главный лепесток, и меньше уровень боковых лепестков. На практике используются значения от 4 до 9.

См. также[править | править вики-текст]

Внешние ссылки[править | править вики-текст]