- •Лабораторная работа № 3 спектральный анализ сигналов
- •1.1 Изучить разложение в ряд Фурье последовательности прямоугольных импульсов. Проанализировать явление Гиббса.
- •1.2 Реализовать разложение в ряд Фурье пилообразного сигнала и последовательности треугольных импульсов.
- •1.3 Сделать общие выводы.
- •Часть 2
- •1. Смоделировать различные виды сигналов и получить их спектры в spTool (Signal Processing Tool).
- •2.Создать многочастотный сигнал на фоне сильного шума, создаваемого генератором случайных чисел согласно своему варианту.
- •3. Построить график спектральной плотности полученного сигнала с помощью прямого преобразования Фурье.
- •Часть 3
- •1. Смоделировать звуковой сигнал.
- •2 Выполнить дискретное преобразование Фурье для этого сигнала получить спектр звукового сигнала.
- •3. Восстановить исходный сигнал при помощи обратного преобразования Фурье.
Лабораторная работа № 3 спектральный анализ сигналов
Цель работы: ознакомление со спектральным описанием периодических функций с помощью рядов Фурье.
Необходимые теоретические сведения. Разложение в ряд Фурье
Первым рассматриваемым сигналом будет последовательность прямоугольных импульсов с амплитудой А, длительностью и периодом повторенияТ. Начало отсчета времени примем расположенным в середине импульса (рис.1).
Рис 1. - Периодическая последовательность прямоугольных импульсов
Данный сигнал является четной функцией, поэтому для его представления удобнее использовать синусно-косинусную форму ряда Фурье— в ней будут присутствовать только косинусные слагаемые , равные
Введем скважность в полученную формулу для коэффициентов ряда Фурье, а затем приведем формулу к виду .
Представление последовательности прямоугольных импульсов в виде ряда Фурье имеет вид:
Амплитуды гармонических слагаемых ряда зависят от номера гармоники по закону (см. рис. 2).График функции имеет лепестковый характер. Итак, ширина лепестков, измеренная в количестве гармоник, равна скважности последовательности (при имеем, если). Отсюда следует важное свойство спектра последовательности прямоугольных импульсов — в нем отсутствуют (имеют нулевые амплитуды) гармоники с номерами, кратными скважности.
Рис. 2 - Коэффициенты ряда Фурье для последовательности прямоугольных импульсов.
Расстояние по частоте между соседними гармониками равно частоте следования импульсов — . Ширина лепестков спектра, измеренная в единицах частоты, равна , то есть обратно пропорциональна длительности импульсов, т.е. чем короче сигнал, тем шире его спектр.
Важным частным случаем предыдущего сигнала является меандр (рис. 3) — последовательность прямоугольных импульсов со скважностью, равной , когда длительности импульсов и промежутков между ними становятся равными.
Рис. 3 - Меандр
Тогда,
,
где m – произвольное целое число.
Таким образом, в спектре меандра присутствуют только нечетные гармоники. Представление меандра в виде ряда Фурье с учетом этого может быть записано следующим образом:
.
Гармонические составляющие, из которых складывается меандр, имеют амплитуды, обратно пропорциональные номерам гармоник, и чередующиеся знаки. На примыкающих к разрыву участках сумма ряда Фурье дает заметные пульсации. Это явление, присущее рядам Фурье для любых сигналов с разрывами первого рода (скачками), называется эффектом Гиббса. Можно показать, что амплитуда первого (самого большого) выброса составляет примерно 9 % от величины скачка.
Рисунок 4. Эффект Гиббса.
Пилообразный сигнал (рис. 5). в пределах периода описывается линейной функцией:
, .
Данный сигнал является нечетной функцией, поэтому его ряд Фурье в синусно-косинусной форме будет содержать только синусные слагаемые:
.
Сам ряд Фурье для пилообразного сигнала выглядит следующим образом:
Рис. 5 - Пилообразный сигнал.
Периодическая последовательность треугольных импульсов имеет симметричную форму (рис. 6):
, .
Рис. 6 - Последовательность треугольных импульсов.
Ряд Фурье имеет следующий вид:
.
Рассмотрим программу, реализующую разложение в ряд Фурье прямоугольной последовательности импульсов.
» N = 8; % число ненулевых гармоник » t = -1:0.01:1; % вектор моментов времени » А = 1; % амплитуда » Т = 1; % период » nh = (1:N)*2-1; % номера ненулевых гармоник » % строки - гармоники » harmonics = cos(2*pi*nh'*t/T); » Am = 2/pi./nh; % амплитуды гармоник » Am(2:2:end) = -Am(2:2:end); % чередование знаков » s1 =harmonics .* repmat(Am', 1, length(t)); » % строки - частичные суммы гармоник » s2 = cumsum(s1); » for k=1:N, subplot(ceil(N/2),2, k), plot(t,s2(k,:)),end
|
ЗАДАНИЕ1.