Смекни!
smekni.com

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

Рис. 1.2.1. Примеры дискретной свертки.

Преобразование свертки однозначно определяет выходной сигнал для установленных значений входного сигнала при известном импульсном отклике системы. Обратная задача деконволюции - определение функции y(k) по функциям s(k) и h(n), имеет решение только при определенных условиях. Это объясняется тем, что свертка может существенно изменить частотный спектр сигнала s(k) и восстановление функции y(k) становится невозможным, если определенные частоты ее спектра в сигнале s(k) полностью утрачены.

Корреляция существует в двух формах: автокорреляции и взаимной корреляции.

Взаимно-корреляционная функция (ВКФ, cross-correlation function - CCF), и ее частный случай для центрированных сигналов функция взаимной ковариации (ФВК) – это показатель степени сходства формы и свойств двух сигналов. Для двух последовательностей x(k) и y(k) длиной К с нулевыми средними значениями оценка взаимной ковариации выполняется по формулам:

Kxy(n) = (1/(K-n+1))

x(k) y(k+n), n = 0, 1, 2, … (1.2.2)

Kxy(n) = (1/(K-n+1))

x(k-n) y(k), n = 0, -1, -2, … (1.2.2')

Рис. 1.2.2. Функция взаимной ковариации двух детерминированных сигналов.

Пример определения сдвига между двумя детерминированными сигналами, представленными радиоимпульсами, по максимуму ФВК приведен на рис. 1.2.2. По максимуму ФВК может определяться и сдвиг между сигналами, достаточно различными по форме.

На рис. 1.2.3 приведен аналогичный пример ФВК двух одинаковых по форме сигналов, на один из которых наложен шумовой сигнал. Мощность шума превышает мощность сигнала. Вычисление ФВК на рисунке выполнено в двух вариантах. Вариант 1 полностью соответствует формуле (1.2.2). Но в условиях присутствия в сигналах достаточно мощных шумов вычисление ФВК обычно выполняется по варианту 2 – с постоянным нормировочным множителем. Это определяется тем, что по мере увеличения сдвига n и уменьшения количества суммируемых членов в формуле (1.2.2) за счет шумовых сигналов существенно нарастает ошибка оценки ФВК, которая к тому же увеличивается за счет нелинейного увеличения значения нормировочного множителя, особенно при малом количестве отсчетов. Сохранение множителя постоянным в какой-то мере компенсирует этот эффект.

Рис. 1.2.3. ФВК двух сигналов, один из которых сильно зашумлен.

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

Рис. 1.2.4. ФВК двух зашумленных радиоимпульсов.

Относительный количественный показатель степени сходства двух сигналов x(k) и y(k) - функция взаимных корреляционных коэффициентов rxy(n). Она вычисляется через центрированные значения сигналов (для вычисления взаимной ковариации нецентрированных сигналов достаточно центрировать один из них), и нормируется на произведение значений стандартов (средних квадратических вариаций) функций x(k) и y(k):

rxy(n) = Kxy(n)/(sx sy). (1.2.3)

sx2 = Kxx(0) = (1/(K+1))

(x(k))2, sy2 = Kyy(0) = (1/(K+1))
(y(k))2. (1.2.4)

Интервал изменения значений корреляционных коэффициентов при сдвигах n может изменяться от –1 (полная обратная корреляция) до 1 (полное сходство или стопроцентная корреляция). При сдвигах n, на которых наблюдаются нулевые значения rxy(n), сигналы некоррелированны. Коэффициент взаимной корреляции позволяет устанавливать наличие определенной связи между сигналами вне зависимости от физических свойств сигналов и их величины.

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

Автокорреляционная функция (АКФ, correlation function, CF) является количественной интегральной характеристикой формы сигнала, дает информацию о структуре сигнала и его динамике во времени. Она, по существу, является частным случаем ВКФ для одного сигнала и представляет собой скалярное произведение сигнала и его копии в функциональной зависимости от переменной величины значения сдвига:

Bx(n) = (1/(K-n+1))

x(k) x(k+n), n = 0, 1, 2, … (1.2.5)

АКФ имеет максимальное значение при n=0 (умножение сигнала на самого себя), является четной функцией Bxy(-n)=Bxy(n), и значения АКФ для отрицательных координат обычно не вычисляются. АКФ центрированного сигнала Kx(n) представляет собой функцию автоковариации (ФАК). ФАК, нормированная на свое значение Kx(0)=sx2 в n=0:

rx(n) = Kx (n)/Kx(0) (1.2.6)

называется функцией автокорреляционных коэффициентов.

Рис. 1.2.5. Автокорреляционные функции.

В качестве примера на рис. 1.2.5 приведены два сигнала – прямоугольный импульс и радиоимпульс одинаковой длительности Т, и соответствующие данным сигналам формы их АКФ. Амплитуда колебаний радиоимпульса установлена равной

амплитуды прямоугольного импульса, при этом энергии сигналов будут одинаковыми, что подтверждается равными значениями максимумов АКФ. При конечной длительности импульсов длительности АКФ также конечны, и равны удвоенным значениям длительности импульсов (при сдвиге копии конечного импульса на интервал его длительности как влево, так и вправо, произведение импульса со своей копией становится равным нулю). Частота колебаний АКФ радиоимпульса равна частоте колебаний заполнения радиоимпульса (боковые минимумы и максимумы АКФ возникают каждый раз при последовательных сдвигах копии радиоимпульса на половину периода колебаний его заполнения).

Линейная цифровая фильтрация является одной из операций ЦОС, имеющих первостепенное значение, и определяется как

s(k) =

h(n) y(k-n), (1.2.7)

Рис. 1.2.6. Трансверсальный цифровой фильтр.

где: h(n), n=0, 1, 2, … , N – коэффициенты фильтра, y(k) и s(k) – вход и выход фильтра. Это по сути свертка сигнала с импульсной характеристикой фильтра.

На рис. 1.2.6 показана блок-схема фильтра, который в таком виде широко известен, как трансверсальный (z – задержка на один интервал дискретизации).

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

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

Самым распространенным преобразованием является дискретное преобразование Фурье. При K отсчетов функции:

S(n) =

s(k) exp(-j 2p kn/K). (1.2.8)

Напомним, что дискретизация функции по времени приводит к периодизации ее спектра, а дискретизация спектра по частоте - к периодизации функции. Для дискретных преобразований s(kDt) Û S(nDf), и функция, и ее спектр дискретны и периодичны, а числовые массивы их представления соответствуют заданию на главных периодах Т = KDt (от 0 до Т или от -Т/2 до Т/2), и 2fN = NDf (от -fN до fN), где K, N – количество отсчетов сигнала и его спектра соответственно. При этом:

Df = 1/T = 1/(KDt), Dt = 1/2fN = 1/(NDf), DtDf = 1/N, N = 2TfN = K. (1.2.9)

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

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

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