- •Лекция 7
- •Vector x
- •Получим дискретное преобразование Фурье (ДПФ) для этого вектора, что символически выражается следующим образом.
- •Последовательность X n , является, вообще говоря, последовательностью комплексных чисел. Поэтому интересно построить
- •Если ДПФ (2) рассматривать для любых значений индекса k, то обнаруживаются интересные свойства
- •На этом рисунке можно увидеть еще одно свойство. График АЧХ симметричен относительно 17-ого
- •Теперь обратимся к свойствам ДПФ, рассмотренным на предыдущей лекции.
- •2. Линейность. Линейной комбинации векторов соответствуетx1,x2 линейная комбинация ДПФ.
- •Vector x1
- •Vector x2
- •3. Циклический сдвиг влево на m позиций. Циклическому сдвигу влево на m позиций
- •Рассмотрим для примера, последовательность, рассмотренную выше (1) у которой компоненты убывают по экспоненциальному
- •Vector x
- •Определение. Под сверткой двух векторов a (a0 , a1, , aNи 1 )
- •Обратим внимание на следующие свойства свертки.
- •Vector b
- •Vector c
- •Третье, если периоды векторов a и b разные и равны N1 и N2,
- •Vector a
- •Vector b
- •Vector c
- •Свойство ДПФ свертки векторов. Вектор являющийся сверткой двух других векторов имеет ДПФ равный
- •Поэтому ДПФ (2) для всех трех векторов a, b, c строится не по
- •Пусть компоненты второго вектора имеют вид единичной ступеньки.
- •Vector a
- •Vector b
- •Vector c
- •Дискретное преобразование Фурье и пакет MATLAB
- •Обратное дискретное преобразование Фурье ОДПФ, выполняет функция ifft(X) . Вызов этой функции осуществляется
- •Эти формулы в пакете MATLAB выглядят следующим образом.
- •Дискретное преобразование Фурье и спектр сигналов
- •Здесь F - частота Найквиста, а sn отсчеты дискретного сигнала. Предположим, что дискретный
- •Спектр дискретного сигнала (36) является непрерывной периодической функцией с периодом 2F . На
- •Мы видим, что сумма, стоящая в формуле (38) является ДПФ для последовательности sn
- •Таким образом, дискретный спектр SD ( fk )
- •Поэтому можно сказать, что дискретное преобразование Фурье является по сути дела, дискретным спектром
- •Если преобразование Фурье сводится к аналитическим выражениям, то задача нахождения спектра сигнала решается
- •Finite signal
- •Тогда спектр S( f ) такого сигнала будет равен интегралу с конечными пределами.
- •где N для удобства четное число. Отсчеты сигнала берем в дискретные моменты времени
- •Выразим шаг дискретизации через частоту Найквиста
- •Сравним формулу (47) с выражением (36) для спектра дискретного сигнала. Мы видим, что
- •Это связано с тем, что в интеграле (45) мы рассматриваем непрерывный сигнал для
- •Подставляем дискретные частоты (50) в формулу (47) и получаем
- •Пусть несущая частота в этом импульсе равна Гц, а
- •Spectrum | S |
- •Spectrum arg( S ) 3
- •Изменим нумерацию компонент дискретного сигнала , чтобы это было похоже на последовательность, для
- •Вместо спектра заданного формулой (51) рассмотрим спектр, полученный с помощью ДПФ (54)
- •Spectrum | S |
- •Spectrum arg( S ) 3
- •Простой анализ показывает, что связь между компонентами дискретного спектра (51) непрерывного сигнала и
- •Быстрое преобразование Фурье (БПФ)
- •В этом алгоритме важным моментом является число членов N суммы (57). Рассмотрим ДПФ
- •Для вычисления каждого коэффициента y(k) , как легко видеть, требуется около 2n комплексных
- •Введем вектор:
- •Далее разобьем сумму (59) на два слагаемых с четными и нечетными членами. При
- •Сравнивая соотношения (64) с формулой (59), мы видим, что
- •Объединим формулы (65), (67)
- •через коэффициенты ДПФ размерности 2n 1 .
- •Итого, реализация ДПФ размерности
- •В общем случае для ДПФ размерности 2n операция сведения к двум ДПФ меньшей
3. Циклический сдвиг влево на m позиций. Циклическому сдвигу влево на m позиций компонент вектора-сигнала, соответствует умножение компонент вектора-сигнала на фазовый множитель.
if
x (x0 , x1, , xN 1 ) y (xm , xm 1, , xN 1, x0 , x1, , xm 1 ),
x X, |
y Y, |
|||
then |
|
|
|
(11) |
|
|
|
|
|
Yn e |
i |
2 |
n m |
X n |
|
N |
|
Следствие. Инвариантность амплитудного спектра относительно циклических сдвигов. При любом циклическом сдвиге амплитуда компонентов ДПФ не меняется.
Yn X n
23
Рассмотрим для примера, последовательность, рассмотренную выше (1) у которой компоненты убывают по экспоненциальному закону. Возьмем период равный N = 32 и осуществим сдвиг влево на три позиции m = 3 .
Результат показан на следующих трех рисунках.
24
Vector x
1 |
|
|
|
|
0.8 |
|
|
|
|
0.6 |
|
|
|
|
0.4 |
|
|
|
|
0.2 |
|
|
|
|
00 |
10 |
20 |
30 |
40 |
25
12 |
|
| X | |
|
|
|
|
|
|
|
10 |
|
|
|
|
8 |
|
|
|
|
6 |
|
|
|
|
4 |
|
|
|
|
2 |
|
|
|
|
00 |
10 |
20 |
30 |
40 |
26
4 |
|
arg(X) |
|
|
|
|
|
|
|
3 |
|
|
|
|
2 |
|
|
|
|
1 |
|
|
|
|
0 |
|
|
|
|
1 |
|
|
|
|
2 |
|
|
|
|
30 |
10 |
20 |
30 |
40 |
27
Определение. Под сверткой двух векторов a (a0 , a1, , aNи 1 )
b (b0 , b1, , bN 1 )с периодом N , будем понимать вектор
c (c0 , c1, , cN 1 ) с периодом 2N вдвое большим. Компоненты вектора-сверки определяются следующими формулами.
N 1 |
N 1 |
|
ck anbk n ak nbn , |
k 0, , 2N 1 (12) |
|
n 0 |
n 0 |
|
Здесь предполагается, что компоненты векторов a и b равны нулю для следующих значений индексов.
am bm 0, if |
m 0 or |
m N 1 |
(13)
28
Обратим внимание на следующие свойства свертки.
Первое, если периоды векторов a и b одинаковые и равны N , то период свертки будет в два раза больше 2N . На рисунках рассматриваются вектора с периодом N = 4 .
Второе, как легко показать из формул (12), (13) последняя компонента свертки всегда равна нулю.
c2N 1 |
0 |
(14) |
! Доказать самостоятельно свойство (14)
29
6 |
|
|
Vector a |
|
|
|
|
|
|
|
|
|
|
5 |
|
|
|
|
|
|
4 |
|
|
|
|
|
|
3 |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
00 |
2 |
4 |
6 |
8 |
10 |
12 |
30
Vector b
6
5
4
3
2 |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
00 |
2 |
4 |
6 |
8 |
10 |
12 |
31
Vector c
6 |
|
|
|
|
|
|
5 |
|
|
|
|
|
|
4 |
|
|
|
|
|
|
3 |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
00 |
2 |
4 |
6 |
8 |
10 |
12 |
32