Дискретное преобразование Фурье
.docДискретное преобразование Фурье.
Преобразование Фурье служит для перехода из одной области, в которой мы рассматриваем некоторый объект (как правило это функция, описывающий этот объект), в другую, в которой те или иные свойства объекта могут быть наиболее полно и точно исследованы. Основой для преобразования Фурье является набор (базис) ортогональных функций, по которому производят разложение объекта (функции) в ряд Фурье.
Применительно к исследованию сигналов, преобразование Фурье используют для перехода из временной области в частотную, позволяющего получить частотный спектр сигнала. Суть и основные свойства метода изложены ниже.
Рассмотрим периодическую последовательность с периодом в N отсчетов. Ее можно записать следующим образом
xp(n)= Xp(k) , (1)
причем частоты спектральных составляющих, образующих xp(n) , могут принимать только значения k=2k/N, -<k<, поскольку периоды других частот не кратны N. В равенстве (1) коэффициенты Xp(k) представляют амплитуды синусоид с частотами k . Запись вида (1) избыточна вследствие периодичности функции , т.к. комплексные экспоненты с частотами k=2k/N и kmN=(2/N)(k mN), -<m<, не отличаются друг от друга, т.е.
= (2)
Следовательно, равенство (1) можно переписать в виде
xp(n)= Xp(k) (3)
подчеркивающем наличие всего N различных комплексных экспонент с периодом в N отсчетов. Для удобства перепишем равенство (3) в общепринятом виде
xp(n)= Xp(k) (4)
где деление на N не изменяет способа представления. Чтобы выразить коэффициенты Xp(k) через xp(n) умножим обе части равенства (4) на и просуммируем результаты по n
(5)
Меняя в правой части (5) порядок умножения и используя формулу
(6)
получим
(7)
или [после перестановки левой и правой частей равенства (7) и замены индекса m на k]
(8)
Соотношение (8) носит название дискретного преобразования Фурье (ДПФ), а (3) – обратного преобразования Фурье (ОДПФ).
Из определений (3) и (8) видно, что обе последовательности xp(n) и Xp(k) полностью определяются одним периодом xp(n) . Отсюда возникает интересный вопрос: как связаны z-преобразование конечной последовательности, образованной из одного периода периодической последовательности, и ДПФ всей периодической последовательности? Иначе говоря, рассмотрим последовательность конечной длины
(9)
причем последовательность xp(n) имеет период в N отсчетов, т.е. x(n) представляет собой один период последовательности xp(n). z-преобразование x(n) равно
(10)
Вычисляя сумму (10) при z=, т.е. в точке на единичной окружности с полярным углом 2k/N, находим
(11)
Сравнивая суммы (11) и (8) и учитывая, что xp(n) =x(n) на интервале 0 n N-1, получаем
(12)
Итак, коэффициенты ДПФ последовательности конечной длины равны значениям z-преобразования этой же последовательности в N точках, равномерно распределенных по единичной окружности. Еще более важный вывод состоит в том, что коэффициенты ДПФ последовательности конечной длины однозначно представляют саму последовательность, т.к. по ним можно точно восстановить исходную последовательность, используя ОДПФ. Итак, хотя ДПФ и ОДПФ вводятся для периодических последовательностей, важно, что через них можно представлять последовательности конечной длины.
Свойства ДПФ.
-
Линейность
Если xp(n) и yp(n) – периодические последовательности (с периодом в N отсчетов каждая), а Xp(k) и Yp(k) – их ДПФ, то дискретное преобразование Фурье последовательности xp(n)+ yp(n) равно Xp(k) + Yp(k). Это положение справедливо и для последовательностей конечной длины.
-
Сдвиг
Если последовательность xp(n) периодическая с периодом в N отсчетов, а ее ДПФ равно Xp(k), то ДПФ периодической последовательности xp(n-n0) будет равно
-
Свойства симметрии
Если периодическая последовательность xp(n) с периодом в N отсчетов является действительной, то ее ДПФ Xp(k) удовлетворяет следующим условиям симметрии
Re[Xp(k)]=Re[Xp(N-k)]
Im[Xp(k)]=-Im[Xp(N-k)]
| Xp(k)|=| Xp(N-k)|
arg Xp(k)=-arg Xp(N-k)
Аналогичные равенства справедливы и для конечной последовательности x(n), имеющей N-точечное ДПФ X(k).
-
Свертка последовательностей
Если xp(n) и hp(n) – две периодические последовательности с периодами по N отсчетов и ДПФ, равными
то N-точечное ДПФ последовательности yp(n) , являющейся круговой (или периодической) сверткой последовательностей xp(n) и hp(n), т.е.
,
равно
Yp(k)=Hp(k)Xp(k)