- •Продолжение
- •Теорема 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) получаем связь между спектром
Следствие из свойств 3, 4. Инвариантность амплитудного спектра относительно циклических сдвигов. При любом циклическом сдвиге амплитуда компонентов ДПФ не меняется.
Yn X n
Доказательство. Используем результаты циклических сдвигов, отмеченные в формулах (19), (23), (24). Объединим эти результаты в виде формулы.
Yn e |
i 2 |
n m |
X n , m Z |
(25) |
N |
|
|||
|
|
|
|
33
Теперь вычислим модуль от выражения (25).
|
Yn |
|
|
|
e |
i 2 |
n m |
X n |
|
|
|
e |
i |
2 n m |
|
X n |
|
|
|
X n |
|
, |
m Z |
(26) |
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
N |
|
|
|
|
N |
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
В формуле (26) учтено, модуль фазового множителя равен единице. Действительно по определению модуля комплексного числа имеем.
ei x 2 ei x e i x 1
Следствие доказано.
34
Определение. Под сверткой двух векторов a (a0 , a1, , aN 1 ) и
b (b , b , , b |
)с периодом N , будем понимать вектор |
|||||
c (c |
0 |
1 |
N 1 |
|
2N вдвое большим. Причем компоненты |
|
, c , с периодом, c ) |
||||||
0 |
1 |
N |
1 |
|
|
|
вектора-сверки определяются следующими формулами. |
||||||
|
|
N 1 |
|
|
N 1 |
|
|
ck |
anbk n |
ak nbn , |
k 1, , 2N 1 (27) |
||
|
|
n 0 |
|
|
n 0 |
|
5. ДПФ свертки векторов. Вектор являющийся сверткой |
|
двух других векторов имеет ДПФ равный произведению ДПФ |
|
исходных векторов. |
N 1 |
if a A, b B, |
ck anbk n , |
|
n 0 |
then |
|
c C, |
(28) |
where |
|
Ck Ak Bk |
35 |
Доказательство. Дополним векторы a и b нулями, чтобы они имели период 2N .
a (a |
, a , , a |
N 1 |
, 0, |
,0), |
0 |
1 |
|
|
|
|
|
|
|
N |
b (b , b , , b |
, 0, |
,0) |
|
0 1 |
N 1 |
|
|
N
Выпишем ДПФ для этих векторов.
2 N 1 |
|
2 |
|
|
i 2 N k n |
||
Ak an e |
|||
n 0 |
|
|
|
2 N 1 |
|
2 |
|
i 2 N k n |
|||
Bk bn e |
|||
n 0
, k 0, 1, , 2N 1,
(29)
, k 0, 1, , 2N 1
36
Для удобства введем обозначение
|
i |
2 |
(30) |
|
e |
2 N |
|||
|
|
Выпишем ДПФ вектора-свертки с использованием формул (28), (29), (30).
2 N 1 |
i |
2 |
|
k n |
2 N 1 |
N 1 |
|
|
|
|
||
|
|
|
|
|
||||||||
Ck cn e |
2 N |
|
|
|
|
|
|
k n |
|
|||
|
|
|
|
ambn m |
|
|||||||
|
n 0 |
|
|
|
|
|
n 0 |
m 0 |
|
|
|
|
N 1 |
2N 1 |
|
|
|
|
|
|
|
|
(31) |
||
|
|
bn m |
|
|
k n |
|
|
|
|
|
||
am |
|
|
|
|
|
|
|
|||||
m 0 |
|
n 0 |
|
|
|
|
|
|
|
|
|
|
37
В формуле (31) изменен порядок суммирования. В сумме в скобках (31) сделаем замену индексов
n m l,
n 0, 1, , 2N 1 l m, 1 m, , 2N 1 m
В результате получи выражение:
N 1 |
|
Ck am |
|
m 0 |
|
2N 1 m |
|
|
N 1 |
k m |
2N 1 m |
k l |
|
|
k (m l ) |
am |
(32) |
||||
bl |
|
|
bl |
|
|||
n m |
|
|
m 0 |
|
l m |
|
|
38
Теперь учтем, что компоненты векторов a и b отличны от нуля только для следующих значений индексов l.
al 0, bl 0, |
l 0, , N 1 |
Поэтому b 0 |
для l 0 или для l N - 1 . Таким образом, |
||||
l |
|
|
|
|
|
формула (32) принимает вид. |
|
|
|
|
|
|
N 1 |
|
N 1 |
|
|
Ck am |
k m |
bl |
k l |
(33) |
|
|
|
|
|||
|
m 0 |
|
l 0 |
|
|
39
Далее, так как al 0, bl 0 для l N – 1, то суммы в (33) не изменяться, если в них добавить нулевые члены. Поэтому в этих суммах верхние пределы можно увеличить до 2N – 1 . Выпишем эти суммы. При этом учтем, чему равен параметр (30), и явный вид ДПФ векторов a и b (29) .
N 1 |
2 N 1 |
2 |
|
|
|
|
|
||
am k m |
am e i 2 N k m |
Ak |
||
m 0 |
m 0 |
|
|
|
(34)
N 1 |
2 N 1 |
2 |
|
|
|
|
|
||
bl k l |
bl e i 2 N |
k l |
Bk |
|
l 0 |
l 0 |
|
|
|
40
С учетом формул (34) выражение (33) принимает искомый вид.
Ck Ak Bk
Свойство доказано.
искретное преобразование Фурье и спектр сигнало
Рассмотрим, какую роль играет дискретное преобразование Фурье в спектральном описании сигналов. Начнем с дискретного сигнала. Как мы знаем, спектр дискретного сигнала выражается формулой.
SD ( f ) |
1 |
|
i |
n |
f |
|
|||||
|
sn e |
|
F |
(35) |
|
|
2F n |
|
|||
|
|
|
|
||
41
Здесь F - частота Найквиста, а sn отсчеты дискретного сигнала. Предположим, что дискретный сигнал определен конечным набором отсчетов.
{sn }, n 0, , N 1
Другими словами для n 0 или для n N – 1 можно считать sn 0 . Поэтому ряд (35) заменяется конечной суммой.
SD ( f ) |
1 |
N 1 |
i |
n |
f |
|
|
(36) |
|||||
|
sn e |
|
F |
|||
|
|
|
|
|
||
2F n 0
42
