Связь ДПФ и ДВПФ. Ядро Дирихле
![]() DSPL-2.0 — свободная библиотека алгоритмов цифровой обработки сигналов Распространяется под лицензией LGPL v3
Страница проекта на SourceForge
|
В предыдущих параграфах мы рассмотрели
ДВПФ
дискретного сигнала
, в общем случае бесконечного во
времени,
а также ДПФ
,
, ограниченной во времени выборки сигнала
,
. Мы специально разделим обозначения сигналов
и
, а также их спектров
и
добавлением соответствующих аббревиатур.
Необходимо напомнить, что при вычислении ДВПФ ограниченного во времени сигнала при
и
, в то время как ДПФ предполагает периодическое повторение
во времени. Наглядно сигналы
,
, а также их спектры показаны на рисунке 1.

a — временное представление


б — спектры


Заметим, что при
, что обведено пунктиром на рисунке 1а.
ДВПФ непериодического сигнала возвращает
—непрерывную,
-периодическую функцию частоты
, в то время как ДПФ возвращает вектор отсчетов спектра
(см. рисунок 1б).
В данном параграфе, мы покажем, ДПФ и ДВПФ имеют однозначную взаимосвязь между собой.
Рассмотрим преобразование ДВПФ и ДПФ:


1805–1859




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

сигнала


а —временно́е представление; б —спектры
Тогда можно записать для любого целого :































Ядро Дирихле представляет собой интерполянт, позволяющий восстановить непрерывную периодическую спектральную плотность дискретного сигнала по
спектральным отсчетам ДПФ:


При программном расчете ДВПФ мы не можем получить значения для континуального множества непрерывной частоты
. Мы можем лишь рассчитать значения спектра для конечного набора фиксированных частот. ДПФ позволяет это делать для сетки
,
.
Однако если требуется произвести расчет для произвольной сетки частот
, где
количество частотных точек ДВПФ,
, то мы можем воспользоваться выражением ДВПФ (1):









На практике длину часто выбирают кратной
, т.е.
. В этом случае производится интерполяция спектра ДПФ, как это показано на рисунке 4.

а —добавление нулей во времени; б — спектры ДПФ
Из рисунка 4 можно видеть, что при расчете ДПФ по отсчетам, мы получаем только
спектральных отсчетов на одном периоде повторения
, показанной пунктирной линией. При дополнении нулями сигнала
до длины
(в приведенном примере
), мы получим
спектральных отсчетов ДПФ на одном периоде повторения
.
Таким образом, мы можем производить численный расчет одного периода на сколь угодно частой сетке частот, при дополнении исходного сигнала нулями.
Дополнение нулями можно также пояснить с позиций ДПФ следующим образом. При выводе выражения ДПФ мы раскладывали в ряд Фурье сигнал, периодически повторенный во времени бесконечное число раз. Период повторения сигнала равен отсчетов, и в результате мы получили дискретные гармоники периодического сигнала на сетке частот
,
, т.е. отсчеты ДПФ
.
Если мы дополним сигнал нулями до длины отсчетов,
, то это равносильно увеличению периода повторения сигнала до
отсчетов. Что в свою очередь, делает сетку частот при разложении в ряд Фурье более частой, и в результате мы получаем более детальную спектральную картину сигнала, потому что берем дискретные значения ДВПФ на более частотой сетке частот
,
.
Если мы увеличим период повоторения сигнала в раз, т.е. возьмем
, то сетка частот будет также в
раз чаще, что мы и наблюдаем на рисунке 4.
Необходимо сделать важное замечание. Дополнение нулями не добавляет информации в исходный сигнал , а значит не улучшает разрешающей способности спектрального анализа. И хотя, мы получили дополнительные спектральные отсчеты ДПФ после добавления нулей, но значения этих спектральных отсчетов не могут быть произвольными, а однозначно связаны с
отсчетами ДПФ
интерполяционной формулой (12).
Рассмотренная взаимосвязь дискретного спектра ДПФ и непрерывного спектра ДВПФ через интерполяционную формулу (12) позволяет сделать вывод, что отсчеты ДПФ не только задают амплитуду и фазу гармоники дискретного сигнала на заданной частоте, а в них заложена информация о поведении непрерывного спектра ДВПФ во всем диапазоне нормированной частоты от до
. Это можно понять из рисунка 3, поскольку ядро Дирихле — непрерывная функция, которая носит периодический затухающе-колебательный характер, а значит, при интерполяции (12), каждый спектральный отсчет ДПФ будет влиять на поведение спектра ДВПФ во всем частотном диапазоне от
до
.
Для подчеркивания этого спектральные отсчеты ДПФ
называют бинами ДПФ[1] [2, стр. 83].
Таким образом, каждый бин ДПФ несет информацию о спектре дискретного сигнала в окрестности соответствующей частоты. При этом частота гармонического синусоидального сигнала в общем случае не попадает точно на частоту бина ДПФ, в результате чего гармонический сигнал проявляется сразу во всех бинах ДПФ с разной степенью интенсивности. Данный эффект будет рассмотрен в следующем параграфе и носит название эффекта растекания спектра.
Итак мы рассмотрели однозначную связь спетра ДПФ и ДВПФ. Было показано, что ДПФ получается как результат дискретизации ДВПФ
на одном периоде повторения спектра.
Аналогично была получена интерполяционная формула, позволяющая восстановить ДВПФ спектр
как функцию непрерывной частоты
используя
спектральных отсчетов ДПФ
. В качестве интерполянта выступает периодическая функция (ядро Дирихле).
Также был проанализирован способ численного расчета ДВПФ, и применением алгоритма ДПФ на произвольной сетке частот, который получается в результате дополнения нулями вектора отсчетов исходного сигнала.
[1] от англ. bin (ларец или закрома)