- •Учебно-методический комплекс дисциплины «введение в распознавание сигналов»
- •230201 «Информационные системы и технологии»
- •230102 «Автоматизированные системы обработки информации и управления»
- •Рабочая программа
- •230102 «Автоматизированные системы обработки информации и управления»
- •1. Цели и задачи дисциплины, ее место в учебном процессе
- •2. Содержание дисциплины
- •Наименование тем, их содержание, объем
- •1. Введение в теорию сигналов. (12 час)
- •2. Метрология сигналов. (14 час.)
- •3. Динамическое представление сигналов. (4 час.)
- •Самостоятельная работа студентов
- •Вопросы к зачету
- •Виды контроля
- •Требования к знаниям студентов, предъявляемые на зачете
- •3. Учебно-методические материалы по дисциплине Перечень обязательной (основной) литературы
- •Перечень дополнительной литературы
- •4. Учебно-методическая (технологическая) карта дисциплины
- •Практические работы
- •Тема 1: введение в теорию сигналов
- •1.1. Общие сведения и понятия [1,10, 15, 25]
- •1.2. Типы сигналов [1,10,15]
- •1.3. Системы преобразования сигналов [1,9,14,18]
- •Литература
- •Тема 2: метрология сигналов
- •2.1. Пространство сигналов [1,3,16,29].
- •2.2. Мощность и энергия сигналов [1,3,16].
- •2.3. Пространства функций [1,3,11,16,29].
- •2.4. Функции корреляции сигналов [1, 25, 29].
- •2.5. Математическое описание шумов и помех [1, 30].
- •Литература
- •Тема 3: динамическое представление сигналов
- •3.1. Разложение сигналов по единичным импульсам [1,11].
- •3.2. Свертка (конволюция) [1,11].
- •Литература
- •Самостоятельная работа
- •5.1. Задачи дискретизации функций [10, 21].
- •5.2. Равномерная дискретизация [16,21].
- •5.3. Дискретизация по критерию наибольшего отклонения [10].
- •5.4. Адаптивная дискретизация [10].
- •5.5. Квантование сигналов [5,21].
- •5.6. Децимация и интерполяция данных [4,5,17].
- •Литература
- •Карта обеспеченности дисциплины кадрами профессорско-преподавательского состава
3.2. Свертка (конволюция) [1,11].
Интеграл Дюамеля. Произвольный сигнал на входе системы с использованием выражений разложения сигнала может быть представлен в виде последовательной линейной комбинации взвешенных единичных импульсов:
y(t) =T[s(t)] =T[s()(t-)d
На основании принципа суперпозиции линейный оператор Т может быть внесен под знак интеграла, т.к. последний представляет собой предельное значение суммы. При этом операция преобразования действует только по переменной t. Отсюда следует:
y(t) = s() Т[(t-)] ds() h(t-) d
Это выражение представляет собой интеграл Дюамеля или свертку (конволюцию) входного сигнала с импульсной характеристикой системы. Заменой переменных t-=можно убедиться в том, что свертка коммутативна:
s() h(t-) d h() s(t-) d. (3.2.1')
Аналогично, для дискретных сигналов:
y(kt) =h(nt) s(kt-nt). (3.2.1'')
Выражения (3.2.1) имеют специальную форму упрощенной математической записи в символическом виде:
y(t) = s(t-) * h() s(t) * h(t).
Сравнением выражений (3.2.1' и 3.2.1'') с выражениями (3.1.3) нетрудно убедиться в их полной идентичности, за исключением нижнего предела интегрирования (суммирования). Это и понятно, так как выражения (3.1.3) были получены при рассмотрении реальной физической системы, работающей в реальном масштабе времени, импульсный отклик которых является односторонним (равен нулю при <0). Для таких систем интегрирование (и суммирование) от -до 0 не имеет смысла. Кроме того, в реальных физических системах импульсный отклик, как правило, отличен от нуля только на определенном интервале, при этом верхний интервал интегрирования (суммирования) в выражениях (3.2.1) также может устанавливаться равным значению этого интервала. В общем случае пределы интегрирования (суммирования) ограничиваются значениями, на которых функция h() существует или имеет значимые значения.
Рис. 3.2.1.
На рис. 3.2.1. приведен пример выполнения свертки прямоугольного импульса с импульсным откликом RC-цепи, площадь которого нормирована к 1. Если площадь импульсного отклика h(t) равна 1, то площадь выходного сигнала свертки всегда должна быть равна площади входного сигнала, что можно видеть на верхнем графике рисунка, при этом одномасштабное сравнение входного и выходного сигналов наглядно демонстрирует характер преобразования сигнала в данной системе. На последующих графиках рисунка демонстрируется вычисление результатов свертки в ряде последовательных точек ti= {3.5, 4, 5, 6, 7} временной оси. В силу отрицательного знакав аргументах функции s(t-) интегрирование произведения h()s(t-) выполняется назад по времени и может ограничиваться только определенной длиной значимых значений импульсного отклика (которая в данном случае установлена равной r = 4), а результат относится к начальной точке h(0) импульсного отклика. Так как входной сигнал, рассмотренный на рисунке, представляет собой прямоугольный импульс с амплитудой 1, то интеграл свертки в каждой текущей точке расчета равен площади импульсного отклика в пределах границ входного прямоугольного импульса (заполнено точками).
Рис. 3.2.2.
Свойства свертки. Для свертки характерны следующие свойства:
1. Дистрибутивность: h(t) * [a(t)+b(t)] = h(t) * a(t)+h(t) * b(t).
2. Коммутативность: h(t) * a(t) * b(t) = a(t) * b(t) * h(t).
3. Ассоциативность: [a(t) * b(t)] * h(t) =h(t) * a(t) * b(t).
Преобразование свертки однозначно определяет выходной сигнал y(t) для установленного значения входного сигнала s(t) при известном значении функции импульсного отклика системы h(t). Обратная задача деконволюции - определение функции s(t) по функциям y(t) и h(t), относится к разряду некорректных, и имеет решение только при вполне определенных условиях. Это объясняется тем, что свертка может существенно изменить частотный спектр сигнала y(t) относительно s(t) и восстановление функции s(t) становится невозможным, если определенные частоты ее спектра в сигнале y(t) полностью утрачены.
Любая практическая система должна быть устойчивой, т.е. для сигналов, конечных по энергии или средней мощности, выходные сигналы также должны быть конечными по этим параметрам. Устойчивость обеспечивается при выполнении условия абсолютной интегрируемости импульсного отклика системы:
|h(t)|dt<.
Для систем с m входами и n выходами аналогично определяются парциальные импульсные отклики hij(t), i = {1,2, ... ,n}, j = {1,2, ... ,m}, каждым из которых отображается сигнал на i-м выходе при поступлении сигналаt) на j-й вход. Полная совокупность импульсных откликов образует матрицу:
,
а выражение свертки приобретает вид:
(t) =()(t-) d
Здесь (и в дальнейшем тексте) жирным шрифтом с "крышкой" выделяются векторные величины.
Системы свертки. Свертка выполняется системой (физическим или программным устройством). Физические системы, работающие в реальном времени, вычисляют текущее значение выходного сигнала по всем прошлым значениям входного сигнала, и не могут иметь в своем распоряжении будущих значений входного сигнала. Операторы таких систем являются односторонними (каузальными). Вышеприведенная нормированная к 1 по площади функция RC-цепи h(t) = (1/RC)exp(-t/RC), принятая в качестве системного оператора на рис. 3.2.1, является именно таким односторонним каузальным оператором. При сравнении выходного сигнала такой системы с входным нетрудно заметить, что выходной сигнал сдвигается относительно входного сигнала. Для каузальных систем такой "сдвиг по фазе" существует всегда и не может быть исключен (сигнал на выходе системы не может быть раньше сигнала на ее входе).
Входным сигналом программных систем является сигнал в целом, записанный в память вычислительного устройства. При обработке таких данных в распоряжении системы при вычислении любой текущей точки выходного сигнала имеются как "прошлые" для данной точки, так и "будущие" значения входного сигнала. Это позволяет создавать системы без сдвига фазы выходного сигнала относительно входного. Для создания таких систем может использоваться два способа:
1. Первый способ иллюстрирует рис. 3.2.3. Задается система с односторонним каузальным оператором h(). Входной сигнал s(t) пропускается через систему в обычном порядке, и выполняется свертка g(t) = h()*s(t). Затем выходной сигнал g(t) реверсируется (g(t)=>g(-t), конец сигнала становится его началом в порядке возрастания t) и повторно пропускается через систему, т.е. выполняется свертка y(-t) = h()*g(-t) . Полученный сигнал снова реверсируется (y(-t) => y(t) = h(-)*g(t)) и является окончательным выходным сигналом y(t) системы.
Рис. 3.2.3.
Рис. 3.2.4.
Приведенное выше формирование двустороннего симметричного оператора свертки имеет чисто познавательный характер. На практике вполне естественным является расчет непосредственно симметричных двусторонних операторов под требуемые задачи обработки числовых данных (сигналов, зарегистрированных в дискретной числовой форме).
Рис. 3.2.5.