Материал: конспект-лекций

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

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.

 

Рис. 6.2.2. Взаимноковариационные функции сиг-
налов.

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)

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

Источник: https://studfile.net/preview/16436030/