126
АКФ зашумленных сигналов. Зашумленный сигнал записывается в виде суммы v(k) = s(k)+q(k). В общем случае, шум не обязательно должен иметь нулевое среднее значение, и норми-
рованная по мощности автокорреляционная функция цифрового сигнала, содержащая N – отсче-
тов, записывается в следующем виде:
Bv(n) = (1/N) s(k)+q(k), s(k-n)+q(k-n) =
=(1/N) [ s(k), s(k-n) + s(k), q(k-n) + q(k), s(k-n) + q(k), q(k-n) ] =
=Bs(n) + M{sk qk-n} + M{qk sk-n} + M{qk qk-n}.
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Bv(n) = Bs(n) + skqk n + qksk n + qkqk n . |
(6.1.13) |
|||||||||||||
При статистической независимости полезного сигнала s(k) и шума q(k) с учетом разложе- |
|||||||||||||||
ния математического ожидания |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
M{sk qk-n} = M{sk} M{qk-n} = sk qk |
|
|
|
|||||||||||
может использоваться следующая формула: |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
2 . |
|
|
|
||||
|
Bv(n) = Bs(n) + 2 sk qk |
+ qk |
|
|
(6.1.13') |
||||||||||
|
|
|
Пример зашумленного сигнала и его АКФ в со- |
||||||||||||
|
|
|
|||||||||||||
|
|
поставлении с незашумленным сигналом приведен на |
|||||||||||||
|
|
рис. 6.1.5. |
|
|
|
|
|
|
|
||||||
|
|
|
Из формул (6.1.13) следует, что АКФ зашум- |
||||||||||||
|
|
ленного сигнала состоит из АКФ сигнальной компо- |
|||||||||||||
|
|
ненты полезного сигнала с наложенной затухающей |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
до |
значения |
2 sk qk + qk |
шумовой функцией. При |
||||||||||
|
|
больших значениях |
K, когда |
qk → 0, имеет место |
|||||||||||
|
|
Bv(n) Bs(n). Это дает возможность не только выде- |
|||||||||||||
Рис. 6.1.5. |
|
лять |
по АКФ |
периодические |
сигналы, практически |
||||||||||
|
|
полностью скрытые в шуме (мощность шумов много |
|||||||||||||
|
|
||||||||||||||
больше мощности сигнала), но и с высокой точностью определять их период и форму в пределах периода, а для одночастотных гармонических сигналов – и их амплитуду с использованием выра-
жения (6.1.6).
|
|
127 |
|
|
|
|
Таблица 6.1. |
Кодовые сигналы являются разновид- |
|
|
|
ностью дискретных сигналов. На определен- |
||
|
|
|
||
M |
Сигнал Баркера |
АКФ сигнала |
|
|
|
|
|
ном интервале кодового слова М t они мо- |
|
2 |
1, -1 |
2, -1 |
||
гут иметь только два амплитудных значения: |
||||
|
|
|
||
3 |
1, 1, -1 |
3, 0, -1 |
||
0 и 1 или 1 и –1. При выделении кодов на |
||||
|
|
|
||
4 |
1, 1, 1, -1 |
4, 1, 0, -1 |
существенном уровне шумов форма АКФ ко- |
|
|
|
|
||
|
1, 1, -1, 1 |
4, -1, 0, 1 |
дового слова имеет особое значение. С этой |
|
|
|
|
||
5 |
1, 1, 1, -1, 1 |
5, 0, 1, 0, 1 |
позиции наилучшими считаются такие коды, |
|
|
|
|
||
7 |
1, 1, 1, -1, -1, 1, -1 |
7, 0, -1, 0, -1, 0, -1 |
значения боковых лепестков АКФ которых |
|
|
|
|
||
11 |
1,1,1,-1,-1,-1,1,-1,-1,1,- |
11,0,-1,0,-1,0,-1,0,- |
минимальны по всей длине интервала кодо- |
|
|
|
|
||
|
1 |
1,0,-1 |
вого слова при максимальном значении цен- |
|
|
|
|
||
|
|
|
|
|
13 |
1,1,1,1,1,-1,-1,1,1-1,1,- |
13,0,1,0,1,0,1,0,1,0,1,0 |
трального пика. К числу таких кодов отно- |
|
|
|
|
||
|
1,1 |
,1 |
сится код Баркера, приведенный в таблице |
|
|
|
|
||
|
|
|
|
6.1. Как видно из таблицы, амплитуда центрального пика кода численно равна значению М, при этом амплитуда боковых осцилляций при n 0 не превышает 1.
6.2. ВЗАИМНЫЕ КОРРЕЛЯЦИОННЫЕ ФУНКЦИИ СИГНАЛОВ [1,19].
Взаимная корреляционная функция (ВКФ) разных сигналов (cross-correlation function,
CCF) описывает как степень сходства формы двух сигналов, так и их взаимное расположение друг относительно друга по координате (независимой переменной). Обобщая формулу (6.1.1) автокор-
реляционной функции на два различных сигнала s(t) и u(t), получаем следующее скалярное про-
изведение сигналов:
Bsu( ) =
s(t) u(t+ ) dt.
(6.2.1)
Взаимная корреляция сигналов характеризует определенную корреляцию явлений и физи-
ческих процессов, отображаемых данными сигналами, и может служить мерой “устойчивости” данной взаимосвязи при раздельной обработке сигналов в различных устройствах. Для конечных
по энергии сигналов ВКФ также конечна, при этом:
|Bsu( )| ||s(t)|| ||u(t)||,
что следует из неравенства Коши-Буняковского и независимости норм сигналов от сдвига по ко-
ординатам.
При замене переменной t = t- в формуле (6.2.1), получаем:
|
|
|
Bsu( ) = s(t- ) u(t) dt = |
u(t) s(t- ) dt = Bus(- ). |
|
Отсюда следует, что для ВКФ не выполняется условие четности, Bsu( ) Bsu(- ), и значе- |
||
ния ВКФ не обязаны иметь максимум при |
= 0. |
|
128
Это можно наглядно видеть на рис. 6.2.1, где заданы два одинаковых сигнала с центрами на точках 0.5 и 1.5. Вычисление по формуле (6.2.1) с постепенным увеличением значений озна-
чает последовательные сдвиги сигнала s2(t) влево по оси времени
(для каждого значения s1(t) для подынтегрального умножения бе-
рутся значения s2(t+ )). При =0 сигналы ортогональны и значение B12( )=0. Максимум В12( )
будет наблюдаться при сдвиге сигнала s2(t) влево на значение =1, при котором происходит пол-
ное совмещение сигналов s1(t) и s2(t+ ).
Одни и те же значения ВКФ по формулам (6.2.1) и (6.2.1') наблюдаются при одном и том же взаимном положении сигналов: при сдвиге на интервал сигнала u(t) относительно s(t) вправо по оси ординат и сигнала s(t) относительно сигнала u(t) влево, т.е. Bsu( ) = Bus(-
На рис. 6.2.2 приведены примеры ВКФ для прямоугольного сигнала s(t) и двух одинаковых треугольных сигналов u(t) и v(t).
Все сигналы имеют одинаковую длитель-
ность Т, при этом сигнал v(t) сдвинут вперед на интервал Т/2.
Сигналы s(t) и u(t) одинаковы по вре-
менному расположению и площадь "пере-
крытия" сигналов максимальна при =0, что и фиксируется функцией Bsu. Вместе с тем функция
Bsu резко асимметрична, так как при асимметричной форме сигнала u(t) для симметричной формы s(t) (относительно центра сигналов) площадь "перекрытия" сигналов изменяется по разному в за-
висимости от направления сдвига (знака при увеличения значения от нуля). При смещении исходного положения сигнала u(t) влево по оси ординат (на опережение сигнала s(t) - сигнал v(t))
форма ВКФ остается без изменения и сдвигается вправо на такое же значение величины сдвига – функция Bsv на рис. 6.2.2. Если поменять местами выражения функций в (6.2.1), то новая функция
Bvs будет зеркально повернутой относительно =0 функцией Bsv.
С учетом этих особенностей полное ВКФ вычисляется, как правило, отдельно для положи-
тельных и отрицательных запаздываний:
Bsu( ) =
s(t) u(t+
) dt. Bus(
) =
u(t) s(t+ ) dt.
(6.2.1')
Взаимная корреляция зашумленных сигналов. Для двух зашумленных сигналов u(t) = s1(t)+q1(t) и v(t) = s2(t)+q2(t), применяя методику вывода формул (6.1.13) с заменой копии сигнала
s(t) на сигнал s2(t), нетрудно вывести формулу взаимной корреляции в следующем виде:
Buv( ) = Bs1s2( ) + Bs1q2( ) + Bq1s2( ) + Bq1q2( ). |
(6.2.2) |
Последние три члена в правой части (6.2.2) затухают до нуля при увеличении |
. При боль- |
ших интервалах задания сигналов выражение может быть записано в следующей форме:
129
Buv( ) = Bs1s2( ) + s1( ) q2( ) + q1( ) s2(
При нулевых средних значениях шумов и статистической
) + |
q1( ) q2( ) . |
(6.2.3) |
независимости от сигналов имеет
место:
Buv( ) → Bs1s2( ).
ВКФ дискретных сигналов. Все свойства ВКФ аналоговых сигналов действительны и для ВКФ дискретных сигналов, при этом для них действительны и особенности дискретных сигналов,
изложенные выше для дискретных АКФ (формулы 6.1.9-6.1.12). В частности, при t = const =1
для сигналов x(k) и y(k) с числом отсчетов К:
Bxy(n) = |
K |
|
K n |
||
|
При нормировании в единицах мощности:
|
1 |
K |
|
Bxy(n) = |
xk |
||
K |
|||
|
k 0 |
||
|
|
K-n |
|
xk yk-n. |
|
k 0 |
|
yk-n |
xk yk n |
.
(6.2.4)
(6.2.5)
Оценка периодических сигналов в шуме. Зашумленный сигнал можно оценить по взаим-
ной корреляции с "эталонным" сигналом методом проб и ошибок с настройкой функции взаимной
корреляции до максимального значения.
Для сигнала u(k)=s(k)+q(k) при статистической независимости шума и qk |
→ 0 функция вза- |
имной корреляции (6.2.2) с шаблоном сигнала p(k) при q2(k)=0 принимает вид: |
|
Bup(k) = Bsp(k) + Bqp(k) = Bsp(k) + qk pk . |
|
А поскольку qk → 0 при увеличении N, то Bup(k) → Bsp(k). Очевидно, что функция Bup(k)
будет иметь максимум, когда p(k) = s(k). Меняя форму шаблона p(k) и добиваясь максимизации
функции Bup(k), можно получить оценку s(k) в виде оптимальной формы p(k).
Функция взаимных корреляционных коэффициентов (ВКФ) является количественным показателем степени сходства сигналов s(t) и u(t). Аналогично функции автокорреляционных ко-
эффициентов, она вычисляется через центрированные значения функций (для вычисления взаим-
ной ковариации достаточно центрировать только одну из функций), и нормируется на произведе-
ние значений стандартов функций s(t) и v(t): |
|
su( ) = Csu( )/ s v. |
(6.2.6) |
Интервал изменения значений корреляционных коэффициентов при сдвигах |
может изме- |
няться от –1 (полная обратная корреляция) до 1 (полное сходство или стопроцентная корреляция).
При сдвигах , на которых наблюдаются нулевые значения su( ), сигналы независимы друг от друга (некоррелированны). Коэффициент взаимной корреляции позволяет устанавливать наличие связи между сигналами вне зависимости от физических свойств сигналов и их величины.
При вычислении ВКФ зашумленных дискретных сигналов ограниченной длины с исполь-
зованием формулы (6.2.4) имеется вероятность появления значений
su(n)| > 1.
Для периодических сигналов понятие ВКФ обычно не применяется, за исключением сигна-
130
лов с одинаковым периодом, например, сигналов входа и выхода при изучении характеристик си-
стем.
6.3. СПЕКТРАЛЬНЫЕ ПЛОТНОСТИ КОРРЕЛЯЦИОННЫХ ФУНКЦИЙ [1,25].
Спектральная плотность АКФ может быть определена из следующих простых сообра-
жений.
В соответствии с выражением (6.1.1) АКФ представляет собой функцию скалярного произ-
ведения сигнала и его копии, сдвинутой на интервал |
при - < |
< : |
Bs( ) = s(t), s(t- |
) . |
|
Скалярное произведение может быть определено через спектральные плотности сигнала и
его копии, произведение которых представляет собой спектральную плотность взаимной мощно-
сти:
s(t), s(t- |
) = (1/2 |
)
S(
) S *( ) d
|
Смещение сигнала по оси абсцисс на интервал |
отображается в спектральном представле- |
||
нии |
умножением спектра сигнала на exp(-j |
), а |
для |
сопряженного спектра на множитель |
exp(j |
): |
|
|
|
|
S *( ) = S*( |
) exp(j |
). |
|
|
С учетом этого получаем: |
|
|
|
s( ) = (1/2
= (1/2
)
)
S(
|S(
) S*( ) exp(j
) d
)|2 exp(j
) d
(6.3.1)
Но последнее выражение представляет собой обратное преобразование Фурье энергетиче-
ского спектра сигнала (спектральной плотности энергии). Следовательно, энергетический спектр сигнала и его автокорреляционная функция связаны преобразованием Фурье:
Bs( ) |S( )|2 = Ws( ). (6.3.2)
Таким образом, спектральная плотность АКФ есть не что иное, как спектральная плотность мощности сигнала, которая, в свою очередь, может определяться прямым преобразованием Фурье через АКФ:
|S( )|2 |
|
|
|
|
= |
Bs( ) exp(-j |
) d . |
(6.3.3) |
Последние выражение накладывает определенные ограничения на форму АКФ и методику их ограничения по длительности.