- •Продолжение
- •Теорема 2. (Теорема отсчетов). Если сигнал s(t) имеет спектр ограниченный полосы, т.е. S(f)
- •Доказательство. Выразим сигнал через его спектр
- •Спектр S(f) , принимающий ненулевые значения только на отрезке f [-F, F] ,
- •Тогда получаем
- •Берем формулу для сигнала (59) и заменяем в ней спектр S(f) на спектр
- •Отсюда, учитывая, что интеграл в (65) равен
- •Итак, шаг дискретизации t позволяет точно восстановить аналоговый сигнал, если его спектр ограничен
- •Это условие можно записать через частоту Найквиста F
- •Лекция 6
- •Соотношение (2) носит название дискретного преобразования Фурье (ДПФ). Иногда говорят, что вектор
- •Теорема 1. Вектор x (x0 , x1, , xN 1 ) можно восстановить
- •В выражении (4) поменяем порядок суммирования.
- •Используем формулу суммы геометрической прогрессии (7) для вычисления суммы (6).
- •В результате формула (8) принимает вид.
- •Таким образом, если n m , сумма в выражении (6) равна нулю.
- •Подставляем (12) в формулу (5), и используем свойство символа Кронекера, сворачивать сумму. В
- •Можно сказать, что вектор x (x0 , x1, , xN 1 ) и
- •Фрагмент кода программы показывает, как можно получить ДПФ для заданного вектора x ,
- •Vector x
- •Обратное дискретное преобразование Фурье ОДПФ, выполняет функция ifft(X) . Вызов этой функции осуществляется
- •Свойства дискретное преобразование Фурье
- •1. Свойства симметрии. Если вещественный вектор x с
- •Из соотношений (16) вытекают условия симметрии для действительных и мнимых частей компонент вектора
- •2. Линейность. Линейной комбинации векторов соответствуетx1, x2 линейная комбинация ДПФ.
- •3. Циклический сдвиг влево. Циклическому сдвигу влево компонент вектора-сигнала, соответствует умножение компонент вектора-сигнала
- •Заменим в формуле (20) компоненты вектора y на компоненты вектора x с помощью
- •Далее выполняем в выражении (21) следующие преобразования.
- •Следствие. Циклический сдвиг влево на m позиций. Циклическому сдвигу влево на m позиций
- •4. Циклический сдвиг вправо на m позиций. Циклическому сдвигу вправо на m позиций
- •Следствие из свойств 3, 4. Инвариантность амплитудного спектра относительно циклических сдвигов. При любом
- •Теперь вычислим модуль от выражения (25).
- •Определение. Под сверткой двух векторов a (a0 , a1, , aN 1 )
- •Доказательство. Дополним векторы a и b нулями, чтобы они имели период 2N .
- •Для удобства введем обозначение
- •В формуле (31) изменен порядок суммирования. В сумме в скобках (31) сделаем замену
- •Теперь учтем, что компоненты векторов a и b отличны от нуля только для
- •Далее, так как al 0, bl 0 для l N – 1, то
- •С учетом формул (34) выражение (33) принимает искомый вид.
- •Здесь F - частота Найквиста, а sn отсчеты дискретного сигнала. Предположим, что дискретный
- •Спектр дискретного сигнала (36) является непрерывной периодической функцией с периодом 2F . На
- •Мы видим, что сумма, стоящая в формуле (38) является ДПФ для последовательности sn
- •Таким образом, дискретный спектр
- •Теперь рассмотрим, как связано дискретное преобразование Фурье со спектром непрерывных сигналов. Спектр непрерывного
- •Тогда спектр S( f ) такого сигнала будет равен интегралу с конечными пределами.
- •В методе прямоугольников интеграл (3) заменяем следующей суммой
- •Сравним формулу (47) с выражением (36) для спектра дискретного сигнала. Мы видим, что
- •Обычно нас интересует спектр, как для положительных значений частоты f 0 ,так и
- •Чтобы сумме в выражении (51) придать вид ДПФ, сделаем
- •Заменим последовательность sn другой
- •Разобьем этот интервал на две части. Рассмотрим сначала правую часть этого интервала
- •Таким образом, мы получили формулу для вычисления половины спектра.
- •Сумма (55) в этом случае примет вид
- •Подставляя (64) в формулы (55), (54) для интервала (61) получаем связь между спектром
Заменим последовательность sn другой |
|
последовательностью xn |
по правилу. |
s |
|
N |
xm , |
m 0, 1, , N 1 |
(53) |
|
m |
|
1 |
|
|
|
2 |
|
|
||
|
|
|
|
|
В этом случае спектр сигнала (52) принимает вид
|
1 |
|
|
2 |
|
N |
1 |
|
N 1 |
|
2 |
|
|
|
|
i |
|
k |
2 |
|
i |
k m |
(54) |
||||
|
|
N |
|||||||||||
|
|
N |
|||||||||||
S( fk ) |
e |
|
|
|
|
xm e |
|||||||
2F |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
m 0 |
|
|
|
|
Сумма в (54) внешне похожа на ДПФ для последовательности xn .
Выпишем эту сумму.
53
~ |
N 1 |
i |
2 |
k m |
(55) |
|
N |
||||
X k xm e |
|
|
|
|
m 0
Отличие суммы (55) от настоящего ДПФ, приведенного в формуле (2), состоит в следующем. В настоящем ДПФ индекс k принимает следующие значения.
k 0, 1, |
, N 1 |
(56) |
|
В сумме (55) индекс k принимает другие значения.
|
N |
|
N |
(57) |
||
|
|
|
|
|
||
k |
2 1, |
, 1, 0, 1, , |
2 |
|||
|
54
Разобьем этот интервал на две части. Рассмотрим сначала правую часть этого интервала
k 0, 1, , |
N |
(58) |
|
2 |
|||
|
|
Интервал (58) содержится в интервале (56) , поэтому для значений (58) сумма (55) совпадает с ДПФ.
~ |
(59) |
X k X k |
Итак, для интервала (58) получаем связь между спектром и ДПФ
S( fk ) 1 e |
i |
2 |
|
N |
|
|
k 0, 1, , N |
|||
N |
2 |
|
Xk , |
|||||||
|
|
|
|
k |
|
1 |
|
|
|
|
|
2F |
|
|
|
|
|
|
2 |
(60)
55
Таким образом, мы получили формулу для вычисления половины спектра.
Теперь рассмотрим левую часть интервала (57)
k |
N |
1, |
, |
1 |
(61) |
2 |
|
Сделаем замену индексов в уравнении (55).
|
|
k N l, |
|
|
|
|
|
k |
N |
1, , 1 |
l |
N |
1, |
, N 1 |
(62) |
2 |
2 |
|
56
Сумма (55) в этом случае примет вид
~ |
N 1 |
|
|
i |
2 |
(l N )m |
|
|
|
|
|
|
|
|
|
N |
|
|
|
|
|||
Xl N xm e |
|
|
|
|
|
|
|||||
|
m 0 |
|
|
|
|
|
|
|
|
(63) |
|
|
|
2 |
|
|
|
|
|
|
|
2 |
|
N 1 |
i |
l m |
|
|
|
N 1 |
|
i |
l m |
||
xm e |
N |
e |
i 2 m |
xm e |
N |
||||||
|
|
|
|
|
|
|
|
||||
m 0 |
|
|
|
|
|
|
m 0 |
|
|
|
|
Последняя сумма в (63) является ДПФ, так как интервал изменения индекса l (62) содержится в интервале (56). Поэтому формула (63) принимает вид
~ |
|
N |
(64) |
|
|
|
|||
X l N X l , |
l |
|
1, , N 1 |
|
2 |
||||
|
|
|
57
Подставляя (64) в формулы (55), (54) для интервала (61) получаем связь между спектром и ДПФ
|
1 |
|
|
2 |
|
N |
1 |
|
|
|
N |
|
|
|
|
|
i |
|
l |
|
|
|
|
|
, N 1 (65) |
||||
S( fl N ) |
e |
N |
2 |
X l , |
l |
1, |
||||||||
|
|
|
|
|||||||||||
2F |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
2 |
|
|
Таким образом, мы получили формулу для вычисления второй половины спектра.
! Доказать самим соотношение (65).
Итак формулы (60) и (65) устанавливают связь между дискретным спектром финитного непрерывного сигнала и ДПФ.
58