16.6. ОПРЕДЕЛЕНИЕ ЗАДЕРЖКИ СИГНАЛА

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

(16.23)

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

Пусть исходному сигналу соответствует -преобразование .

Тогда -преобразование последовательности будет а суммарного сигнала

(16.24)

Из (16.24) следует, что определяемый выражением (16.23) сигнал можно трактовать как свертку где — сигнал, z-преобразование которого равно . Таким сигналом является сумма двух дельтафункций (рис. 16.12):

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

Рис. 16.12. Сигнал

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

Обратимся к (16.24) и, учитывая, что контур интегрирования в (16.20) совпадает с окружностью единичного радиуса на z-плоскости, подставим в множитель вместо z переменную .

Тогда

(16.25)

Таким образом, выражение (16.24) позволяет составить следующее соотношение:

(16.26)

Из (16.26) видно, что наложение задержанной копии на исходный сигнал создает эффект модуляции спектра энергии по закону . Глубина модуляции определяется коэффициентом а период модуляции равен .

С аналогичным явлением мы встретились в примере предыдущего тараграфа, где спектр энергии исходного сигнала также являлся периодической функцией со.

Прологарифмировав выражение (16.26), получим

(16.27)

Примерный вид слагаемых правой части показан на рис. 16.13.

Вычислим кепстр мощности по выражению (16.20), которое на основании (16.25) -(16.27) можно записать в форме

(16.28)

Как видно, информация о задержке содержится в кепстре , поэтому вычисление начнем именно с , не уточняя пока структуры игнала и кепстра .

Рис. 16.13. К выражению (16.27);

Основываясь на выражении (16.25), получаем

Так как , можно воспользоваться разложением

Тогда

Подставив этот результат во второе слагаемое в правой части (16.28), получим

Очевидно, что отлично от нуля только в точках и т. д., причем

Кепстр представлен на рис. 16.14. Истинная задержка определяется по положению первого пика.

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

Рис. 16.14. Кепстр при

Рис. 16.15. Сигнал нормированный к величине