- •Содержание лабораторной работы
- •Задание на лабораторную работу
- •Типовой script-файл для выполнения лабораторной работы
- •Задание на самостоятельную работу
- •Отчет и контрольные вопросы
- •Литература
- •Дискретное преобразование Фурье
- •Восстановление аналогового сигнала
- •Выделение дискретных гармоник полезного сигнала
- •Восстановление спектральной плотности
- •Код скрипта lr_04
Типовой script-файл для выполнения лабораторной работы
Перед
выполнением работы должна быть
представлена табл.1 исходных данных
для
своего номера бригады
.
Для запуска лабораторной работы необходимо обратиться к script-файлу 1r_04 по его имени:
» 1r_04
Для принудительного снятия script-файла с выполнения следует нажать комбинацию клавиш <Ctrl>+<Break>.
При выполнении script-файла текущие окна с графиками не закрывать.
Используемые внешние функции
В script-файле 1r_04 используется две внешние функции.
Внешняя функция fft_e1, предназначенная для определения значений модуля ДПФ (MODm) и дискретных нормированных частот (m) гармоник, которые согласно первому критерию (14) при заданном пороге (17) (e1) относят к полезному сигналу:
function [MODm,m] = fft_e1(MODX,e1)
% Определение значений модуля ДПФ и частот полезного сигнала
%
% MODX — вектор значений модуля ДПФ смеси сигнала с шумом
% e1 — заданный порог
%
% MODm — вектор значений модуля ДПФ полезного сигнала
% m — вектор значений частот полезного сигнала
i = 1;
МАХ = max(MODX);
for k = 1: length(MODX)
if (MODX(k)/MAX)>el
MODm(i) = MODX(k);
m(i) = k-1;
i = i+1;
end
end
Внешняя функция fft_e2, предназначенная для определения модуля значений ДПФ (MODm) и дискретных нормированных частот (m) гармоник, которые согласно второму критерию (15) при заданном пороге (18) (е2) относят к полезному сигналу:
function [MODm,m] = fft_e2 (MODX,e2)
% Определение значений модуля ДПФ и частот полезного сигнала
%
% MODX — вектор значений модуля ДПФ смеси сигнала с шумом
% e2 — заданный порог
% MODm — вектор значений модуля ДПФ полезного сигнала
% m — вектор значений частот полезного сигнала
i = 1;
% Р — средняя мощность смеси сигнала с шумом
P = sum(MODX.^2)/length(MODX);
for k = 1:length(MODX)
if ((MODX(k).^2)/Р)>е2
MODm(i) = MODX(k);
m(i) = к-1;
i = i+1;
end
end
Задание на самостоятельную работу
Задание на самостоятельную работу заключается в создании function-файлов для вычисления ДПФ последовательностей с использованием исходных данных из табл. 1 для своего номера бригады .
Последовательности выбираются из представленного далее списка:
1C.
Периодическая последовательность с
периодом
:
5\* MERGEFORMAT ()
Вывести графики амплитудного и фазового спектра периодической последовательности.
Определить амплитуды и частоты дискретных гармоник, используя function-файл fft_e1.
2С. Конечная последовательность (5) длины .
Вывести графики модуля и аргумента ДПФ конечной последовательности.
3С. Конечная последовательность длины :
6\* MERGEFORMAT ()
где
Вывести
графики конечных последовательностей
(2) и
(6) и модулей их ДПФ.
4С.
Цифровой единичный импульс (7) на
интервале
.
7\* MERGEFORMAT ()
Вывести графики цифрового единичного импульса и модуля его ДПФ.
5С. Последовательность с однотональной амплитудной модуляцией:
Задать
значения
,
,
,
,
,
и период последовательности
.
Вывести графики последовательности и ее амплитудного спектра.
6С.
Последовательность
:
Задать
частоту дискретизации
Гц.
Вывести
графики последовательности на интервале
с шагом
и модуля её ДПФ.
7C. Гауссов радиоимпульс:
.
Задать
и
.
Вывести
графики последовательности на интервале
и модуля ее ДПФ.
