Файл: Обработка сигналов временное представление сигналов временное представление сигналов разложение сигналов по единичным импульсам.doc

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 10.01.2024

Просмотров: 31

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

ОБРАБОТКА СИГНАЛОВ 2. ВРЕМЕННОЕ ПРЕДСТАВЛЕНИЕ СИГНАЛОВ

2. ВРЕМЕННОЕ ПРЕДСТАВЛЕНИЕ СИГНАЛОВ

2.1. Разложение сигналов по единичным импульсам

Прямая координатная форма представления сигналов соответствует их естественной и привычной для нас форме математического описания в виде функций независимых переменных (аргументов или координат).

Единичные импульсы. Простейшим сигналом является единичный импульс. В качестве математической модели единичного импульса при анализе аналоговых сигналов используют дельта-функцию.

Дельта-функция. По определению, дельта-функция описывается следующими математическими выражениями (в совокупности):

(t-) = 0 при t  ,

(t-) dt = 1.

Функция (t-) равна нулю везде за исключением точки , в которой она бесконечно велика и не является дифференцируемой.

Дельта-функция является полезной математической абстракцией. На практике такие функции не могут быть реализованы с абсолютной точностью, так как невозможно реализовать значение, равное бесконечности, в точке t =  на временной шкале. Но во всех случаях, когда площадь импульса равна 1, а длительность импульса достаточно мала и за время его действия на входе сигнал на выходе системы практически не изменяется (реакция системы на импульс во много раз больше длительности самого импульса), входной сигнал можно считать дельта-функцией или, более точно, единичной импульсной функцией со свойствами дельта-функции.

При всей своей абстрактности дельта-функция имеет вполне определенный физический смысл. Представим себе импульсный сигнал прямоугольной формы П(t-) длительностью , амплитуда которого равна 1/, а площадь соответственно равна 1. При уменьшении значения  импульс, сокращаясь по длительности, сохраняет свою площадь, равную 1, и возрастает по амплитуде. Предел такой операции при   0 и носит название дельта-импульса. Этот сигнал (t-) сосредоточен в одной координатной точке t = , конкретное амплитудное значение его не определено, но площадь (интеграл) остается равной 1. Это не мгновенное значение функции в точке t = , а именно импульс – математическая модель короткого действия, значение которого равно 1.

Функция Кронекера. Для дискретных и цифровых систем разрешающая способность по аргументу сигнала определяется интервалом его дискретизации t. Это позволяет в качестве единичного импульса использовать дискретный интегральный аналог дельта-функции - функцию единичного отсчета (kt-nt), которая равна 1 в координатной точке k = n и нулю во всех остальных точках. При этом функция (kt-nt) может быть определена для любых значений t = const, но только для целых значений координат k и n.


Математические выражения (t-) и (kt-nt) называют также импульсами Дирака и Кронекера. Однако применяя такую терминологию, не будем забывать, что это не просто единичные импульсы в координатных точках  и nt, а импульсные функции, определяющие как значения импульсов в определенных координатных точках, так и нулевые значения по всем остальным координатам, в пределе от - до .

Разложение сигнала. Импульсы Дирака и Кронекера используются для разложения, соответственно, произвольных аналоговых сигналов s(t) и дискретных сигналов s(kt) в непрерывную последовательность неперекрывающихся импульсов:



Рис. 2.1.1.
s(t) = s()(t-) d. (2.1.1)

s(kt) = s(nt)(kt-nt). (2.1.1')

Для аналоговых сигналов разложение ( 2.1.1) в физическом представлении эквивалентно сканированию значений сигнала s(t) в моменты времени t =  бесконечно узкой щелью, бегущей вдоль оси t. Для цифровых сигналов эта щель равна одному отсчету. Пример разложения дискретного сигнала приведен на рис. 2.1.1.

С математических позиций единичные импульсные функции (t-), -<, и (kt-nt), -
Импульсный отклик линейной системы. Если на вход линейной системы в момент времени t = 0 подать единичный импульс, то на выходе мы получим реакцию системы на единичный входной сигнал. Эта реакция называется функцией импульсного отклика системы и однозначно определяется оператором преобразования h(..):

y(t) = T[(t-0)] = h(t). (2.1.2)

y(kt) = T[(kt-0)] = h(kt). (2.1.2')

Импульсный отклик аналоговой системы на входную дельта-функцию также в определенной степени представляет собой математическую абстракцию идеального преобразования. С практической точки зрения под импульсным откликом можно понимать отображение реакции системы на импульсный входной сигнал произвольной формы с единичной площадью, если длительность этого сигнала пренебрежимо мала по сравнению с временной (координатной) разрешающей способностью системы (например, с периодом ее собственных колебаний). Для цифровых систем импульсный отклик однозначно определяется реакцией системы на импульс Кронекера. Функцию импульсного отклика называют также

весовой функцией системы.

Очевидно, что в ЛИС-системах форма импульсного отклика не зависит от времени прихода входного сигнала и определяет только его положение на временной оси. Так, если входной импульс задержан (относительно 0) на время to, то соответствующий выходной сигнал будет определяться выражением:

y(t) = T[(t-to)] = h(t-to).

В любой системе, работающей в реальном масштабе времени, сигнала на выходе системы не может быть, если нет сигнала на ее входе. Отсюда следует односторонность импульсного отклика физических систем:

h(t-) = 0 при t<.

Для программных систем, работающих с массивами цифровых данных, импульсный отклик может быть и двусторонним, так как при обработке сигналов в любой текущей точке kt системе доступны как "прошлые" отсчеты kt-nt, так и "будущие" отсчеты kt+nt. Соответственно, сигнал на выходе системы может формироваться и по "будущим" отсчетам. Это расширяет возможности программной обработки сигналов по сравнению с физическими системами.

На рисунке 2.1.2 приведен пример импульсного отклика h(t) элементарной физической системы преобразования электрических сигналов – динамической интегрирующей RC-цепи.




Рис. 2.1.2.
При подаче на вход RC-цепи единичного и очень короткого (t << RC) импульса заряда q емкость С заряжается до напряжения Vо = q/C и начинает разряжаться через сопротивление R, при этом напряжение на емкости изменяется по закону v(t) = Voexp(-t/RC) = (q/C)exp(-t/RC). Отсюда, импульсный отклик RC-цепи на единичный входной сигнал с единичным значением заряда q = 1 равен: h(t) = (1/C)exp(-t/RC), где форма отклика определяется функцией экспоненты, а множитель (1/С) является масштабным преобразователем сигнала (заряда в напряжение). По существу, импульсным откликом системы определяется доля входного сигнала, которая действует на выходе системы по истечении времени t после поступления сигнала на вход (запаздывающая реакция системы).

Если функция импульсного отклика системы известна, то, с учетом принципа суперпозиции сигналов в линейной системе, можно выполнить расчет реакции системы в любой произвольный момент времени на любое количество входных сигналов с любыми моментами времени их прихода путем суммирования запаздывающих реакций системы на эти входные сигналы. На рис. 2.1.2 приведен пример входного сигнала s(t) для RC-цепи в виде последовательности импульсов и реакция системы y(t) на такой входной сигнал, образованная суммированием реакций системы на каждый импульс.


Допустим, что на вход RC-цепи в моменты времени t1=1 и t2=2 поступили очень короткие (по сравнению со значением RC) импульсы заряда величиной A и В. Математически это можно отобразить сигналом s(t) = q1(t)+q2(t), где q1(t) = At-t1) и q2 = B(t-t2). Выходной сигнал системы при известном импульсном отклике h(t) соответственно отобразится формулой:

y(t) = T[q1(t)+q2(t)] = T[A(t-t1)]+T[B(t-t2) = AT[(t-t1)+BT[(t-t2) = Ah(t-t1)+Bh(t-t2).

При расчете значений выходного сигнала в произвольный момент времени t после прихода на вход системы сигналов q1 и q2, например, для t = 5, для каждого из сигналов вычисляются значения их запаздывающих реакций: y1 = Ah(5-1) = Ah(4) и y2 = Bh(5-2) = Bh(3), после чего значения запаздывающих реакций суммируются у = у1+у2. Пример этой операции приведен на рис. 2.1.3, где для удобства графического представления приняты значения А = 1 и В = 1. Сущность операции не изменяется при любых значениях А и В, а в общем случае и для любого количества импульсов.




Рис. 2.1.3.
2.2. Свертка (конволюция)

Интеграл Дюамеля. Произвольный сигнал на входе системы с использованием выражений разложения сигнала может быть представлен в виде последовательной линейной комбинации взвешенных единичных импульсов:

y(t) = T[s(t)] = T[ s()(t-) d

На основании принципа суперпозиции линейный оператор Т может быть внесен под знак интеграла, т.к. последний представляет собой предельное значение суммы. При этом операция преобразования действует только по переменной t. Отсюда следует:

y(t) = s() Т[(t-)] d s() h(t-) d 2

Это выражение представляет собой интеграл Дюамеля или свертку (конволюцию) входного сигнала с импульсной характеристикой системы. Заменой переменных t-= можно убедиться в том, что свертка коммутативна:

s() h(t-) d  h() s(t-) d. (2.2.1')

Аналогично, для дискретных сигналов:


y(kt) = h(nt) s(kt-nt). (2.2.1'')

Выражения (2.2.1) имеют специальную форму упрощенной математической записи в символическом виде:

y(t) = s(t-) * h()  s(t) * h(t).

Техника свертки. Для вычисления свертки по выражению (2.2.1) функция импульсного отклика реверсируется по своей координате, т.е. строится в режиме обратного времени, и движется относительно функции входного сигнала в сторону возрастания значений t.




Рис. 2.2.1.
В каждый текущий момент времени значения обеих функций перемножаются (по всем значениям своих аргументов), а произведение интегрируется. Полученный результат относится к той координатной точке, против которой находится значение импульсного отклика h(0).

На рис. 2.2.1. приведен пример выполнения свертки прямоугольного импульса с импульсным откликом RC-цепи, площадь которого нормирована к 1. Если площадь импульсного отклика h(t) равна 1, то площадь выходного сигнала свертки всегда должна быть равна площади входного сигнала, что можно видеть на верхнем графике рисунка, при этом одномасштабное сравнение входного и выходного сигналов наглядно демонстрирует характер преобразования сигнала в данной системе. На последующих графиках рисунка демонстрируется вычисление результатов свертки в ряде последовательных точек ti = {3.5, 4, 5, 6, 7} временной оси. В силу отрицательного знака  в аргументах функции s(t-) интегрирование произведения h()s(t-) выполняется назад по времени и может ограничиваться только определенной длиной значимых значений импульсного отклика (которая в данном случае установлена равной r = 4), а результат относится к начальной точке h(0) импульсного отклика. Так как входной сигнал, рассмотренный на рисунке, представляет собой прямоугольный импульс с амплитудой 1, то интеграл свертки в каждой текущей точке расчета равен площади импульсного отклика в пределах границ входного прямоугольного импульса (заполнено точками).




Рис. 2.2.2.
Еще более наглядна техника выполнения цифровой свертки, приведенная на рис. 2.2.2. Для вычисления свертки массив одной из функций (s