2.4. Функции корреляции сигналов

Функции корреляции сигналов применяются для интегральных количественных оценок формы сигналов и степени их сходства друг с другом.

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

B_s() = s(t) s(t+) dt. (2.25)

Как следует из этого выражения, АКФ является скалярным произведением сигнала и его копии в функциональной зависимости от переменной величины значения сдвига . Соответственно, АКФ имеет физическую размерность энергии, а при  = 0 значение АКФ непосредственно равно энергии сигнала:

B_s(0) = s(t)² dt = E_s.

Функция АКФ является непрерывной и четной. В последнем нетрудно убедиться заменой переменной t = t- в выражении (2.25):

B_s() = s(t-) s(t) dt = s(t) s(t-) dt = B_s(-). (2.25')

С учетом четности, графическое представление АКФ производится только для положительных значений . На практике сигналы обычно задаются на интервале положительных значений аргументов от 0-Т. Знак + в выражении (2.25) означает, что при увеличении значений  копия сигнала s(t+) сдвигается влево по оси t и уходит за 0, что требует соответствующего продления сигнала в область отрицательных значений аргумента. А так как при вычислениях интервал задания как правило, много меньше интервала задания сигнала, то более практичным является сдвиг копии сигнала влево по оси аргументов, т.е. применение в выражении (2.25) функции s(t-) вместо s(t+).

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

Пример. На интервале (0,Т) задан прямоугольный импульс с амплитудным значением, равным А. Вычислить автокорреляционную функцию импульса.

П ри сдвиге копии импульса по оси t вправо, при 0≤≤T сигналы перекрываются на интервале от  до Т. Скалярное произведение:

B_s() = A² dt = A²(T-).

П ри сдвиге копии импульса влево, при -T≤<0 сигналы перекрываются на интервале от 0 до Т-. Скалярное произведение:

B_s() = A² dt = A²(T+).

При || > T сигнал и его копия не имеют точек пересечения и скалярное произведение сигналов равно нулю (сигнал и его сдвинутая копия становятся ортогональными).

О бобщая вычисления, можем записать:

B_s() = .

В случае периодических сигналов АКФ вычисляется по одному периоду Т, с усреднением скалярного произведения и его сдвинутой копии в пределах периода:

B_s() = (1/Т) s(t) s(t-) dt.

При =0 значение АКФ в этом случае равно не энергии, а средней мощности сигналов в пределах интервала Т. АКФ периодических сигналов также является периодической функцией с тем же периодом Т. Для однотонального гармонического сигнала это очевидно. Первое максимальное значение АКФ будет соответствовать =0. При сдвиге копии сигнала на четверть периода относительно оригинала подынтегральные функции становятся ортогональными друг другу (cos _o(t- = cos (_ot-/2)  sin _ot) и дают нулевое значение АКФ. При сдвиге на =T/2 копия сигнала по направлению становится противоположной самому сигналу и скалярное произведение достигает минимального значения. При дальнейшем увеличении сдвига начинается обратный процесс увеличения значений скалярного произведения с пересечением нуля при =3T/2 и повторением максимального значения при =T=2_o (cos _ot-2 копии  cos _ot сигнала). Аналогичный процесс имеет место и для периодических сигналов произвольной формы (рис. 2.11).

Рис. 2.11

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

Для сигналов, заданных на определенном интервале [a, b], вычисление АКФ производится с нормировкой на длину интервала [a, b]:

B_s() = s(t) s(t+) dt. (2.26)

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

r_s() = cos () = s(t), s(t+) /||s(t)||².

Взаимная корреляционная функция (ВКФ) сигналов (cross-correlation function, CCF) показывает как степень сходства формы двух сигналов, так и их взаимное расположение друг относительно друга по координате (независимой переменной), для чего используется та же формула (2.25), что и для АКФ, но под интегралом стоит произведение двух разных сигналов, один из которых сдвинут на время :

B₁₂() = s₁(t) s₂(t+) dt. (2.27)

При замене переменной t = t- в формуле (2.4.3), получаем:

B₁₂() = s₁(t-) s₂(t) dt = s₂(t) s₁(t-) dt = B₂₁(-)

Рис. 2.12. Сигналы и ВКФ

Отсюда следует, что для ВКФ не выполняется условие четности, а значения ВКФ не обязаны иметь максимум при  = 0. Это можно наглядно видеть на рис. 2.12, где заданы два одинаковых сигнала с центрами на точках 0.5 и 1.5. Вычисление по формуле (2.27) с постепенным увеличением значений  означает последовательные сдвиги сигнала s2(t) влево по оси времени (для каждого значения s1(t) для подынтегрального умножения берутся значения s2(t+)).

При =0 сигналы ортогональны и значение B₁₂()=0. Максимум В₁₂() будет наблюдаться при сдвиге сигнала s2(t) влево на значение =1, при котором происходит полное совмещение сигналов s1(t) и s2(t+). При вычислении значений B₂₁(-) аналогичный процесс выполняется последовательным сдвигом сигнала s1(t) вправо по временной оси с постепенным увеличением отрицательных значений , а соответственно значения B₂₁(-) являются зеркальным (относительно оси t=0) отображением значений B₁₂(), и наоборот. На рис. 2.13 это можно видеть наглядно.

Рис. 2.13. Сигналы и ВКФ

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

Для периодических сигналов понятие ВКФ обычно не применяется, за исключением сигналов с одинаковым периодом, например, сигналов входа и выхода систем при изучении характеристик систем.

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

r_sv() = cos () = s(t), v(t+) /||s(t)|| ||v(t)||. (2.28)

Значение коэффициентов взаимной корреляции может изменяться от -1 до 1.