лекции / DSP_12
.pdf1
Лекции 12. Дискретное преобразование Фурье (ДПФ)
1.Вычисление спектральной плотности конечной последовательности.
2.Ряд Фурье непрерывной периодической функции.
3.ДПФ периодической последовательности.
4.ДПФ конечной последовательности.
5.Основные свойства ДПФ.
6.Вычисление свертки с помощью ДПФ.
12.1. Вычисление спектральной плотности конечной последовательности
Постановка задачи: вычислить СП конечной последовательности длины N
(см. (11.2)):
N 1 |
|
X(ej T ) x(nT)e j Tn . |
(12.1) |
n 0
СП — периодическая функция, поэтому она вычисляется на периоде д 2 .
T
При этом необходимое количество точек заранее неизвестно, т. е. вычисление невозможно выполнить методом прямой подстановки, а значит формула (12.1) не описывает алгоритм вычисления СП.
Таким образом, постановка задачи заключается в выводе формулы, описывающей алгоритм вычисления СП.
12.2. Ряд Фурье непрерывной периодической функции
Ряд Фурье непрерывной периодической функции xp(t) |
с периодом Ts : |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
xp(t) Xа (k)ejk t |
, |
|
(12.2) |
|||||||
|
k |
|
|
|
||||||
Разложение в ряд Фурье справедливо для функции |
xp(t) , удовлетворяющей |
|||||||||
условиям Дирихле (функция |
xp(t) |
имеет ограниченную вариацию: определена, |
||||||||
однозначна и имеет конечное число экстремумов и разрывов 1-го и 2-го рода). |
||||||||||
Xа (k) — коэффициенты Фурье: |
|
|
|
|
|
|
|
|||
|
|
1 Ts |
|
|
|
|||||
Xа (k) |
|
xp(t)e jk tdt , |
(12.3) |
|||||||
T |
||||||||||
|
|
s |
0 |
|
|
|
|
|
|
|
где индекс «а» означает «аналоговый»; |
|
|
|
|
|
|
|
|||
Δω — период дискретизации по частоте: |
|
|
|
|||||||
|
Δω 2π Ts ; |
|
|
(12.4) |
||||||
k — дискретная нормированная частота: |
|
|
|
|||||||
|
|
|
|
, |
|
|
|
|
||
|
|
k |
kΔω |
|
|
(12.5) |
||||
|
|
Δω |
|
|
||||||
|
|
|
|
|
|
|
|
|
||
подобно дискретному нормированному времени n |
nT |
|
(формально T 1). |
|||||||
T |
||||||||||
|
|
|
|
|
|
|
||||
Совокупность коэффициентов |
Фурье |
|
Xа (k), k , называют спектром |
|||||||
периодической функции xp(t) . |
|
|
|
|
|
|
|
|
|
Спектр аналогового периодического сигнала представляет собой бесконечную последовательность в частотной области.
Ряд Фурье (12.2) по определению должен быть абсолютно сходящимся:
|
|
|
2 |
||
|
|
|
|
||
|
Xа (k) |
. |
(12.6) |
12.3.ДПФ периодической последовательности
Рассмотрим возможность разложения в ряд Фурье периодической последовательности.
Выполнив замену t nT , получим:
Ts NT ; xp (t) xp (nT);
Δω |
2π |
, |
(12.7) |
|
NT |
||||
|
|
|
запишем формулу для коэффициентов Фурье — X(k) (индекс «а» снят; назначение «тильды» поясним позже):
|
|
|
|
|
1 |
|
N 1 |
j |
2π |
nTk |
|
|||||||
|
|
|
X(k) |
|
xp(nT)e NT |
. |
|
|||||||||||
|
|
|
|
|
||||||||||||||
|
|
|
|
|
NT n 0 |
|
|
|
|
|
|
|
|
|||||
В области дискретного нормированного времени (T = 1): |
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
1 |
|
N 1 |
j |
2π |
nk |
|
|||||
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
Xp(k) |
|
xp(n)e N |
. |
|
(12.8) |
|||||||||
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
N n 0 |
|
|
|
|
|
|
|
|
|||
Последовательность Xp (k) — периодическая, |
т. к. это функция |
периодического |
||||||||||||||||
j |
2π |
nk |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
аргумента e N с периодом N: |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
2π |
|
2π |
|
|
|
|
|
2π |
|
||||||
|
|
ej |
|
n(k mN) ej |
|
nke j2πmN ej |
|
nk . |
|
|||||||||
|
|
N |
N |
N |
|
Вывод: переход t nT привел к тому, что коэффициенты стали периодическими (поэтому добавлен индекс «p»).
Рассмотрим, к чему приведет подстановка коэффициентов Xp (k) в правую часть
ряда (12.2):
|
|
j |
2 |
nk |
|
|
|
Xp(k)e N . |
(12.9) |
k
Оказывается, что полученный ряд не удовлетворят условию абсолютной сходимости (абсолютно расходящийся):
|
|
Xp(k) |
, |
аследовательно, разложение периодической последовательности в ряд Фурье невозможно.
Определим, чему соответствует один период ряда (12.9):
N 1 |
j |
2 |
nk |
|
|
|
|||
... |
Xp(k) e N . |
(12.10) |
k0
Сэтой целью определим значения непрерывной периодической функции xp (t)
(12.2) в дискретных точках nT :
|
|
|
|
j |
2 |
nk |
xp(t) |
|
t nT xp(nT) xp(n) |
|
|||
|
Xа(k)e N . |
k
3
Заменим вычисление бесконечной суммы вычислением бесконечного числа конечных сумм из N слагаемых, сдвинутых друг относительно друга по оси частот k на N :
|
|
N 1 |
|
j |
2 |
n(k mN) |
|
N 1 |
j |
2 |
nk |
|||
xp(n) |
|
|
|
|||||||||||
Xа (k mN)e N |
|
Xа (k mN)e N . |
||||||||||||
|
m k 0 |
|
|
|
|
m k 0 |
|
|
|
|||||
Изменим порядок суммирования: |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
2 |
|
|
|
|
||
|
|
N 1 |
|
|
|
|
j |
nk |
|
|
|
|||
|
|
|
|
|
|
|
||||||||
|
|
xp(n) |
|
|
Xа (k mN) |
e |
N . |
(12.11) |
||||||
|
|
k 0 |
m |
|
|
|
|
|
|
|
|
|
||
При равенстве правых частей (12.10) и (12.11): |
|
|
|
|
|
|
|
|
||||||
|
|
Xp(k) |
|
|
|
|
|
|
|
|
|
|||
|
|
Xа (k mN) |
, |
|
(12.12) |
m
обеспечивается равенство левых частей, т.е. в (12.10) можем записать левую часть:
N 1 |
j |
2 |
nk |
|
|
|
|||
xp(n) Xp(k)e |
|
N . |
(12.13) |
|
k 0 |
|
|
|
|
Совокупность периодических коэффициентов |
|
Xp (k) |
называют спектром |
|
периодической последовательности xp(n). |
|
|
|
|
Выводы: |
|
|
|
|
1.Соотношение (12.12) описывает связь между спектрами периодических дискретного и аналогового сигналов подобно связи между СП непериодических сигналов, а именно: спектр периодического дискретного сигнала равен бесконечной сумме копий спектров аналогового сигнала, сдвинутых друг относительно друга по оси частот k на N .
2.Последовательность xp(n) и ее спектр Xp (k) имеют одинаковые периоды N
соответственно по n и по k .
Дискретным преобразованием Фурье периодической последовательности xp(n)
называют пару формул:
|
|
N 1 |
j |
2π |
nk |
|
|
|
|||||
|
|
|
|
|
|
|
|
||||||
ДПФ |
Xp(k) xp(n)e |
|
N |
, |
k 0,1,..., (N 1); |
(12.14) |
|||||||
|
|
n 0 |
|
|
|
|
|
|
|
|
|
||
|
|
1 |
N 1 |
|
j |
2π |
nk |
|
|
||||
|
|
|
|
|
|
|
|||||||
ОДПФ |
xp(n) |
Xp(k)e N |
, |
n 0,1,..., (N 1). |
(12.15) |
||||||||
|
|||||||||||||
|
|
N k 0 |
|
|
|
|
|
|
|
|
|
||
Здесь вместо Xp (k) принято записывать Xp (k) : |
|
||||||||||||
|
|
|
Xp(k) |
1 |
Xp (k) , |
(12.16) |
|||||||
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
N |
|
|
|
т. к. где постоянный множитель принято ставить в обратном преобразовании; xp(n) — N-точечная последовательность (один период);
Xp (k) — N-точечное ДПФ (коротко — ДПФ) (один период).
Вывод: ДПФ Xp (k) периодической последовательности xp(n) совпадает с ее |
|||
спектром Xp |
(k) с точностью до множителя |
1 |
. |
|
|||
|
|
N |
12.4. ДПФ конечной последовательности
Сделаем два предположения:
1. Конечная последовательность x(n) длины N — это__________________
4
_____________________________________________________________
x(n) xp (n), n 0,1,..., N 1.
2. ДПФ Xp (k) — это спектральная плотность |
X(ej T ) в N дискретных точках |
||||||
на периоде |
|
2 |
: |
|
|
|
|
|
|
|
|||||
д |
|
T |
|
|
|||
|
|
Xp(k) X(ej T ) |
|
k |
, |
k 0,1,..., N 1, |
|
|
|
|
где Δω 2π .
NT
Формулы ДПФ (12.14)—(12.15) примут общий вид для периодической последовательности с периодом N и конечной последовательности длины N :
|
|
|
|
|
N 1 |
|
|
|
|
|
|
|
|
|
|
ДПФ |
X(k) x(n)WNnk, |
k 0,1,..., N 1 |
; |
|
(12.17) |
||||||||
|
|
|
|
|
n 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
N 1 |
|
|
|
|
|
|
|
|
||
|
ОДПФ |
x(n) |
|
|
X(k)WNnk , |
n 0,1,..., N 1 |
, |
(12.18) |
||||||
|
|
|
||||||||||||
|
|
|
N k 0 |
|
|
|
|
|
|
|
|
|||
где WNnk — поворачивающий множитель: |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
j |
2 |
nk |
|
|
|
|
|
|
|
|
Wnk e |
|
N |
; |
|
|
|
(12.19) |
|||||
|
|
|
|
N |
|
|
|
|
|
|
|
|
||
x(n) — N-точечная последовательность (один период); |
|
|||||||||||||
X(k) — N-точечное ДПФ (один период). |
|
|
|
|
|
|
|
|
||||||
Таким образом, формулы ДПФ (12.14)—(12.15) имеют двойную трактовку: |
|
|||||||||||||
для периодической последовательности |
x(n) |
с периодом N ДПФ X(k) |
— это |
|||||||||||
ее спектр с точностью до множителя |
1 |
; |
|
|
|
|
|
|||||||
N |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
для конечной последовательности x(n) длины N ДПФ X(k) — это дискретные отсчеты ее СП в N дискретных точках на периоде д .
Поставленная задача — получить формулу, описывающую алгоритм вычисления СП конечной последовательности решена. Количество точек при вычислении СП совпадает с длиной конечной последовательности N .
Вопрос о достаточности точек N дает теорема Котельникова в частотной области, симметричная теореме Котельникова во временной области.
Теорема Котельникова во временной области: непрерывный сигнал с финитной спектральной плотностью «длительности» Fs 2 fв может быть точно восстановлен
по своим отсчетам с периодом дискретизации по времени T |
1 |
|
1 |
. |
|
|
|||
|
Fs |
2fв |
Теорема Котельникова в частотной области: непрерывная спектральная плотность
финитного сигнала длительности Ts NT может быть точно восстановлена по
своим отсчетам с периодом дискретизации по частоте f |
1 |
|
1 |
(или |
||
|
|
|||||
|
2π |
|
Ts |
NT |
||
Δω |
). |
|
|
|
|
|
|
|
|
|
|
||
|
NT |
|
|
|
|
12.5. Основные свойства ДПФ
Свойства ДПФ подобны свойствам СП:
1. ДПФ — периодическая последовательность с периодом N .
5
2.Линейность.
3.Сдвиг ДПФ на k0 отсчетов вправо:
x(n) X(k) x(n)WNk0n X(k k0)
4.ДПФ задержанной последовательности:
x(n) X(k)
|
|
x(n m) X(k)Wmk |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
N |
|
|
|
|
5. Равенство Парсеваля: |
|
|
|
|
|
|
|
|
|
|
|
||
N 1 |
|
|
|
2 |
|
1 N 1 |
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
x(n) |
|
|
|
|
|
|
X(k) |
|
. |
|
|
|
|
|
|
|
|
||||||||
n 0 |
|
|
|
|
|
N k 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
12.6. Вычисление свертки с помощью ДПФ и ОДПФ |
|
||||||||||||
Формула линейной свертки для воздействия длины N1 |
и ИХ длины N2 : |
|
|||||||||||
L 1 |
|
|
|
|
L 1 |
|
|
|
|
||||
y(n) h(m)x(n m) h(n m)x(m), |
(12.20) |
||||||||||||
m 0 |
|
|
|
|
m 0 |
|
|
|
|
где L N1 N2 1 — длина свертки.
С целью сокращения количества арифметических операций свертку вычисляют с использованием ДПФ и ОДПФ:
1.Воздействие и ИХ дополняют нулями до одинаковой длины L.
2.Переходят к круговой свертке для последовательностей одинаковых длин L:
L 1 |
|
L 1 |
|
(12.21) |
y(n) h(m)x(n m) h(n m)x(m). |
||||
m 0 |
|
m 0 |
|
|
Результат вычисления круговой свертки представляет собой периодическую последовательность y(n) с периодом L.
3. |
Вычисляют ДПФ y(n) (12.21): |
|
|
|
Y(k) H(k)X(k) , |
k 0,1,...,(L 1). |
(12.22) |
4. |
Вычисляют ОДПФ Y(k) — реакцию y(n) |
(12.20). |
|