Скачиваний:
5
Добавлен:
28.04.2023
Размер:
203.63 Кб
Скачать

1

Лекции 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);

Δω

,

(12.7)

NT

 

 

 

запишем формулу для коэффициентов Фурье — X(k) (индекс «а» снят; назначение «тильды» поясним позже):

 

 

 

 

 

1

 

N 1

j

nTk

 

 

 

 

X(k)

 

xp(nT)e NT

.

 

 

 

 

 

 

 

 

 

 

 

NT n 0

 

 

 

 

 

 

 

 

В области дискретного нормированного времени (T = 1):

 

 

 

 

 

 

 

 

 

 

1

 

N 1

j

nk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Xp(k)

 

xp(n)e N

.

 

(12.8)

 

 

 

 

 

 

 

 

 

 

 

 

N n 0

 

 

 

 

 

 

 

 

Последовательность Xp (k) — периодическая,

т. к. это функция

периодического

j

nk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

аргумента e N с периодом N:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ej

 

n(k mN) ej

 

nke jmN 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

nk

 

 

 

 

 

 

 

 

 

 

 

ДПФ

Xp(k) xp(n)e

 

N

,

k 0,1,..., (N 1);

(12.14)

 

 

n 0

 

 

 

 

 

 

 

 

 

 

 

1

N 1

 

j

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,

 

 

 

где Δω .

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

(или

 

 

 

 

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).

 

Соседние файлы в папке лекции