Индексация и перестановка спектральных отсчетов дискретного преобразования Фурье

Content

DSPL-2.0 is free digital signal processing algorithms library

Distributed under LGPL v3 license

GitHub project page.

Found a mistake? Select it with the mouse and press ctrl+enter
Введение

В предыдущем разделе мы получили выражения для прямого и обратного ДПФ:

(1)

Расчет ДПФ и ОДПФ ведется на основе индексов временных и частотных отсчетов без учета частоты дискретизации. Таким образом, можно использовать выражения для ДПФ и ОДПФ при любой частоте дискретизации, не меняя вычислительную программу. В данном разделе мы рассмотрим, как привязать индексы ДПФ к значениям частоты , выраженной в Герц, или к значениям циклической частоты рад/c.

Индексация спектральных отсчетов

При рассмотрении ДПФ мы говорили, что спектр дискретного сигнала , , это периодическая функция с периодом рад/c, где частота дискретизации исходного сигнала (Гц). Соответственно период повторения спектра дискретного сигнала , , для частоты , выраженной в Герц, равен частоте дискретизации Гц.

ДПФ получается путем дискретизации периодической функции на одном периоде повторения с шагом

(2)

Таким образом, ый спектральный отсчет соответствует частоте

(3)
или
(4)

Пример 1. При частоте дискретизации , при , первый спектральный отсчет соответствует частоте .

Пример 2. При частоте дискретизации , при , спектральный отсчет с номером , соответствует частоте .

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

Пусть входной сигнал , , является комплексным и состоит из одной комплексной экспоненты с частотой Гц:

(5)

Зададим частоту дискретизации равной Гц и возьмем отсчетов исходного сигнала.

Амплитудный спектр сигнала
Figure 1. Амплитудный спектр сигнала

Рассчитаем ДПФ , , данного сигнала и получим амплитудный спектр, как это показано на рисунке 1.

Как можно заметить на рисунке 1, в спектре сигнала присутствует только одна компонента с индексом что согласно (4) соответствует частоте

(6)
что может показаться странным, потому что мы задавали частоту исходного сигнала Гц.

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

Периодический спектр  сигнала
Figure 2. Периодический спектр сигнала

Из рисунка 2 следует, что при дискретизации спектра на одном периоде повторения от 0 до Гц (на рисунке 1 обозначено как период ДПФ) в выборку попадает периодическая гармоника на частоте 100 Гц.

При этом спектр в интервале частот от Гц до Гц периодически повторяет спектр в интервале частот от Гц до Гц.

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

Сделаем важное замечание. Частотная компонента, соответствующая частоте Гц, в силу периодичности спектра дискретного сигнала, также соответствует частоте Гц.

При перестановке мы будем относить эту компоненту к частоте Гц.

Перестановка спектральных отсчетов ДПФ для четного

В случае четного спектральный отсчет , согласно (4), соответствует частоте Гц. Как мы отметили выше, этот же отсчет соответствует частоте Гц. Тогда можно записать спектр после перестановки для четного в виде:

(7)
ый спектральный отсчет после перестановки соответствует частоте
(8)
или
(9)

Таким образом, отсчет соответствует частоте Гц, отсчет соответствует частоте Гц и отсчет соответствует частоте Гц.

Перестановка спектральных отсчетов для четного показана на рисунке 3.

Перестановка спектральных отсчетов для четного
Figure 3. Перестановка спектральных отсчетов для четного

Перестановка спектральных отсчетов ДПФ для нечетного

В случае нечетного спектральный отсчет , согласно (4), соответствует частоте Гц, а спектральный отсчет соответствует частоте Гц.

Спектр после перестановки для нечетного :

(10)

ый спектральный отсчет после перестановки соответствует частоте

(11)
или
(12)

После перестановки спектральных отсчетов, в случае нечетного , спектральный отсчет , согласно (12), соответствует частоте Гц, спектральный отсчет соответствует частоте Гц и последний отсчет соответствует частоте Гц.

Перестановка спектральных отсчетов для нечетного показана на рисунке 4.

Перестановка спектральных отсчетов для нечетного
Figure 4. Перестановка спектральных отсчетов для нечетного

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

Произведем перестановку спектральных отсчетов для корректного отображения отрицательных частот для приведенного выше примера (рисунок 1).

Количество отсчетов в приведенном примере , значит, мы можем воспользоваться выражением (7) для получения . Тогда после перестановки спектральные отсчеты будут соответствовать частотам (9).

На рисунке 5 приведен амплитудный спектр после перестановки, а также значения частоты (Гц), которой соответствуют спектральные отсчеты .

Амплитудный  спектр  после   перестановки  и  соответствующие значения  частоты 
Figure 5. Амплитудный спектр после перестановки и соответствующие значения частоты 

Как можно видеть на рисунке 5, после перестановки спектральных отсчетов компонента соответствует частоте Гц, как мы и задавали для исходного сигнала.

Выводы

В данном разделе мы рассмотрели вопрос индексации и перестановки спектральных отсчетов на выходе ДПФ.

Были приведены выражения для перестановки спектральных отсчетов для четного и нечетного для корректного отображения отрицательных частот после ДПФ.

Более детально свойства ДПФ будут рассмотрены в следующих разделах.

Page update: 23.06.2020 (23:45:52)