- •Содержание лабораторной работы
- •Задание на лабораторную работу
- •Типовой script-файл для выполнения лабораторной работы
- •Задание на самостоятельную работу
- •Отчет и контрольные вопросы
- •Литература
- •Дискретное преобразование Фурье
- •Восстановление аналогового сигнала
- •Выделение дискретных гармоник полезного сигнала
- •Восстановление спектральной плотности
- •Код скрипта lr_04
Отчет и контрольные вопросы
Отчет
составляется в редакторе МS
Word
и содержит исходные данные и результаты
выполнения каждого пункта задания,
включая копируемые из окна Command
Window
результаты вычислений
(шрифт Courier
New),
созданные графики (копируются
по команде Edit
Сору
Figure
в
окне Figure)
и
ответы на поставленные вопросы
(шрифт Times
New
Roman).
Защита лабораторной работы проводится на основании представленного отчета и контрольных вопросов из следующего списка:
Запишите формулы ДПФ.
Что такое поворачивающий множитель?
Чему равно разрешение по частоте при вычислении ДПФ?
С чем связаны трудности прямого вычисления ДПФ по формуле (8)?
Что такое БПФ?
Каков порядок сложности алгоритмов ДПФ и БПФ Кули—Тьюки?
Назовите основные свойства ДПФ.
Дайте определение дискретной нормированной частоты.
Поясните смысл ДПФ для периодической последовательности.
Как с помощью ДПФ рассчитывается амплитудный и фазовый спектры периодической последовательности?
Поясните смысл ДПФ для конечной последовательности.
Как связаны значения абсолютных частот (в герцах [Гц] и радианах в секунду [рад/с]) и дискретных нормированных частот?
Поясните смысл приведенных критериев для выделения полезного сигнала из его аддитивной смеси с шумом.
Как задаются значения порогов в первом и втором критериях выделения полезного сигнала?
Как восстановить аналоговый периодический сигнал с финитным спектром по отсчетам ДПФ и на основе ряда Котельникова?
Как вычислить спектральную плотность в L точках на основе ДПФ при
?Как определить разрешение по частоте в ДПФ при добавлении нулей к исходной последовательности?
Литература
Солонина А. И., Арбузов С. М. Цифровая обработка сигналов. Моделирование в MATLAB. – СПб.: БХВ-Петербург, 2008. – Главы 11.
Сергиенко А. Б. Цифровая обработка сигналов.– 3-е изд.– СПб.: БХВ-Петербург, 2010. – Глава 5
Приложение А.
Дискретное преобразование Фурье
Дискретным преобразованием Фурье (ДПФ) называется пара взаимно однозначных образований:
прямое ДПФ (Discrete Fourier Transform — DFT):
8\* MERGEFORMAT ()
обратное ДПФ (ОДПФ) (Inverse Discrete Fourier Transform — IDFT):
9\* MERGEFORMAT ()
где
– дискретное нормированное время
;
– дискретная нормированная частота
;
– период дискретизации по частоте
(разрешение по частоте),
–
-точечная
последовательность, т. е. периодическая
последовательность во временной области
с периодом
;
–
-точечное
ДПФ, т. е. периодическая последовательность
в частотной области с периодом
,
– период последовательности и ДПФ;
– поворачивающий множитель;
–
-я
дискретная гармоника.
Значения абсолютных частот дискретных гармоник связаны со значениями дискретных нормированных частот соотношением:
10\* MERGEFORMAT ()
Дискретное преобразование Фурье (8) трактуется по-разному в зависимости от вида последовательности – периодическая с периодом или конечная длины N.
Для
периодической последовательности
с периодом
ДПФ
(8) представляет собой ее спектр
с точностью до множителя
.
Модуль
ДПФ
(с точностью до множителя
)
называют амплитудным спектром,
а аргумент
– фазовым спектром периодической
последовательности.
Амплитудный спектр вещественной периодической последовательности равен модулю ДПФ с точностью до множителя:
11\* MERGEFORMAT ()
При вычислении ДПФ (8) периодической последовательности она может задаваться на периоде или на целом числе периодов , что не меняет результата.
Для
конечной последовательности
длины
ДПФ
(8) представляет собой
дискретных равноотстоящих значений
ее спектральной плотности
на периоде
.
Для вещественных последовательностей, периодических и конечных, модуль ДПФ – четная, а аргумент – нечетная функция частоты .
Согласно определению, при вычислении ДПФ предполагается, что последовательность является периодической, и конечная последовательность представляет собой один период периодической последовательности.
При
этом точное выделение гармоник
последовательности
с частотами
гарантируется только в том случае, если
они кратны периоду дискретизации
частоте
:
12\* MERGEFORMAT ()
в
свою очередь, возможно только в том
случае, если на интервале
последовательности
укладывается целое число периодов
,
т. е. отношение
13\* MERGEFORMAT ()
является целым числом.
В случае, если условие (12) не выполняется, наблюдается эффект растекания спектра.
В MATLAB ДПФ (8) – (9) вычисляется с использованием алгоритмов БПФ2 и ОБПФ с помощью функций:
X = fft(x)
x = ifft (X)
где X и x – -точечные последовательность и ее ДПФ – векторы, нижняя граница индексов которых равна единице, в отличие от ДПФ (8) – (9), где она равна нулю.
