Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Лаб. работа - Корреляционный приемник

.pdf
Скачиваний:
64
Добавлен:
10.02.2015
Размер:
245.13 Кб
Скачать

Московский государственный университет имени Н.Э. Баумана

А.С. Косолапов, А.И. Сенин

ИССЛЕДОВАНИЕ КОРРЕЛЯЦИОННОГО ПРИЕМНИКА

Издательство МГТУ им. Н.Э. Баумана

МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ имени Н.Э. БАУМАНА

А.С. Косолапов, А.И. Сенин

ИССЛЕДОВАНИЕ КОРРЕЛЯЦИОННОГО ПРИЕМНИКА

Методические указания к лабораторной работе по курсу «Статистическая радиотехника»

Под редакцией А.И. Сенина

М о с к в а Издательство МГТУ им. Н.Э. Баумана

2 0 0 6

УДК 621.396 ББК 32.84 К71

Рецензент Е.А. Скороходов

Косолапов А.С., Сенин А.И.

К71 Исследование корреляционного приемника: Методические указания к лабораторной работе по курсу «Статистическая радиотехника» / Под ред. А.И. Сенина. – М.: Изд-во МГТУ им. Н.Э. Баумана. 2006. – 16 с.: ил.

ISBN 5-7038-2817-1

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

Для студентов 3-го курса специальности «Радиоэлектронные системы и устройства».

Ил. 5.

УДК 621.396 ББК 32.84

ISBN 5-7038-2817-1

© МГТУ им. Н.Э. Баумана, 2006

Цель работы – изучение корреляционного метода приема сигналов и исследование помехоустойчивости корреляционного приемника.

ЗАДАНИЕ ПО РАБОТЕ

1.Ознакомиться с теоретическим материалом по учебному пособию «Информационные технологии в радиотехнических системах» авторов В.А. Васина, И.Б. Власова, Ю.М. Егорова и других под редакцией И.Б. Федорова. – М.: Изд-во МГТУ им. Н.Э. Баумана, 2004. 766 с.: ил. (Сер. Информатика в техническом университете) и по данным методическим указаниям.

2.Изучить лабораторную установку по настоящему руковод-

ству.

3.Выполнить работу в указанном порядке.

КРАТКИЕ ТЕОРЕТИЧЕСКИЕ СВЕДЕНИЯ

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

Рассмотрим задачу обнаружения сигнала s(t), все параметры

которого известны, на фоне гауссовского аддитивного шума со спектральной плотностью в полосе пропускания приемника f ,

равной N0.

Пусть сигнал на входе приемника имеет вид u(t) = λs(t) + n(t),

где λ – случайная величина, принимающая значение 1 или 0 с вероятностями p и 1 p соответственно. Задача заключается в том,

чтобы по принятому сигналу u(t) решить, присутствует или нет полезный сигнал s(t) на входе приемника.

3

Можно показать, что оптимальный обнаружитель должен вы-

числять отношение функций правдоподобия l(u) =

W (u H1)

и

W (u H0 )

 

 

 

 

 

сравнивать его с некоторым порогом

l0 , значение которого за-

висит от выбранного критерия:

 

 

 

 

в случае использования критерия Байеса

 

 

l =

П01(1 p)

;

 

 

 

 

 

0

 

П10 p

 

 

 

 

 

 

 

 

в случае критерия идеального наблюдателя l0 = 1 p p ;

в случае критерия максимального правдоподобия

l0 =1;

в случае критерия Неймана–Пирсона

l0 = c0 ( p{l(u) c0 H0} = α0 ).

Здесь W (uH1) и W (uH0 ) – функции правдоподобия при справедливости гипотез H1 (полезный сигнал присутствует) и H0 (полезный сигнал отсутствует); П01 и П10 – функции потерь при принятии ошибочных решений в пользу гипотез H0 и H1 соответственно; запись p{l(u) c0 H0 } = α0 означает вероятность того,

что отношение функций правдоподобия при условии справедливости гипотезы H0 не меньше значения с0 ; α0 <1 – некоторая напе-

ред заданная величина.

Нетрудно показать, что для рассматриваемого случая логарифм отношения функций правдоподобия

 

E

 

2

T

 

ln l(u) =

+

u(t)s(t)dt,

(1)

N

 

N

 

 

0

 

0 0

 

 

 

 

 

 

4

где Е– энергия сигнала; Т – его длительность.

Из выражения (1) следует, что оптимальный обнаружитель должен вычислить интеграл

 

2

T

 

q =

u(t)s(t)dt

(2)

N

 

0 0

 

и сравнить его значение с порогом z0 = ln l0 + E N0 .

При превышении порога решение принимается в пользу гипотезы H1, в противном случае – в пользу H0.

Интеграл (2) называется корреляционным, так как является мерой взаимной корреляции между реализацией u(t) и ожидаемым сигналом s(t). Соответственно устройство, вычисляющее выраже-

ние (2), называется корреляционным приемником.

Таким образом, оптимальный обнаружитель состоит из коррелятора и порогового устройства – ПУ (рис. 1). Коррелятор, в свою очередь, состоит из перемножителя (X), интегратора и генератора опорного сигнала (ГОС), представляющего собой копию сигнала s(t).

Рис. 1. Оптимальный обнаружитель

Определим условные вероятности ложного обнаружения и пропуска сигнала:

 

 

 

E

 

 

 

 

α = p q ln l

+

H

 

 

=

N0

 

 

0

 

 

0

 

 

 

 

E

 

 

β = p q < ln l0

+

H1

 

N0

 

 

 

 

W (q H0 )dq,

z0

z0

=W (q H1 )dq,

−∞

(3)

(4)

5

где W (qH1) и W (qH0 ) – распределения напряжения на выходе корреляционного приемника при справедливости гипотез H1 и H0 соответственно.

Из (2) нетрудно понять, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W (q H0 ) =

1

 

 

 

 

 

 

 

 

 

q

2

 

 

 

 

 

 

 

 

 

 

 

 

exp

 

 

 

 

 

,

 

 

(5)

 

 

2E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2E

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

W (q H ) =

 

 

 

 

 

 

exp

 

 

 

 

 

N0

 

.

(6)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2E

 

 

 

 

 

 

 

 

 

2E

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

N0

 

 

 

 

 

N

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Используя (3)–(6), находим

 

 

 

α =

 

E

ln l

 

+

 

N0

0

 

 

ln l

+

 

E

 

N0

0

β =

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

E

 

 

 

 

 

 

 

 

 

 

 

 

ln l

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

1

 

 

 

 

q

 

− Φ

0

 

N0

 

 

 

 

 

exp

 

 

 

dq = 1

 

 

 

,

 

2E

 

 

 

 

 

 

 

2E

 

 

 

 

2

 

2E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N0

 

 

 

 

0

 

 

N0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2E

 

2

 

 

 

 

E

 

 

 

 

 

 

 

q

 

 

 

 

ln l0

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

N0

 

 

 

 

 

N0

 

 

 

 

 

 

 

 

 

exp

 

 

 

 

dq = Φ

 

 

 

 

 

 

 

.

 

2E

 

2E

 

 

 

2E

 

 

 

2

 

 

 

 

 

 

 

N0

 

 

N

 

 

 

 

N

 

 

 

 

0

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

x

 

 

z2

 

 

Здесь Φ(x) =

 

 

 

 

 

 

 

 

exp

 

dz

– интеграл вероятности.

 

−∞

 

 

2

 

 

На рис. 2 представлены графики распределения величин W (qH1), W (qH0 ) для случая наличия и отсутствия полезного

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

α и β.

6

Рассмотрим теперь задачу различения двух полностью известных сигналов s0 (t) и s1(t). Сигнал на входе приемника можно записать в виде

u(t) = λ s1(t) + (1−λ) s0 (t) + n(t),

где λ – случайная величина, принимающая значение 1 или 0 с вероятностями p и 1 p соответственно; n(t) – «белый» гауссов-

ский шум с односторонней спектральной плотностью H0. Необходимо оптимальным образом по принятой реализации u(t) решить, какой из сигналов передавался, или, что то же самое, произвести выбор между гипотезой H0 (присутствует сигнал s0 ) и альтернативной гипотезой H1 (присутствует сигнал s1 ).

Рис. 2. Графики распределения величин q при наличии и отсутствии полезного сигнала

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

 

W (u / H )

 

правдоподобия

l(u) = W (u / H01 )

и сравнить его значение с поро-

гом l

=

1p

.

 

0

p

 

 

7

Решение принимается в пользу гипотезы H1, если l(u) l0 , и в пользу гипотезы H0 , если l(u) < l0.

Для рассматриваемого случая нетрудно показать, что логарифм отношения правдоподобия определяется как

ln l(u) =

E1 E0

+

2

T u(t)[s (t) s (t)]dt,

 

 

 

 

 

 

 

 

N0

 

 

 

N0

 

1

0

 

 

 

 

 

 

0

 

 

 

 

где E1 и E0 – энергии сигналов s1(t)

и s0 (t)

соответственно (на

практике обычно E1 = E0 ,

 

p = 1/ 2).

 

 

 

 

При этом решение принимается в пользу сигнала s1(t),

если

q =

2

T u(t)[s (t)

s

(t)]dt 0.

(7)

 

 

 

N0

1

 

0

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Выражение (7) определяет структуру оптимального различителя. Он состоит из двух корреляционных приемников, вычитающего и порогового устройств (рис. 3). Один из корреляционных приемников настроен на прием сигнала s1(t), другой – на прием s0 (t).

Опорные сигналы, вырабатываемые на приемной стороне, должны быть синхронизированы с принятыми сигналами.

Рис. 3. Структурная схема оптимального различителя

Найдем среднюю вероятность ошибки для случая равновероятных сигналов

pош = 12 (α +β),

8

где α = W (qH0 )dq – условная вероятность принятия решения

ln l0

о наличии сигнала s1 в то время, когда в действительности пере-

ln l0

давался сигнал s0 ; β = W (qH1) dq – условная вероятность

−∞

принятия решения о наличии сигнала s0 в то время, когда в дейст-

вительности передавался сигнал s1.

 

 

 

 

Найдем вероятности α и β. Для

этого определим

сначала

плотности вероятности W (q H1) и W (q H0 ).

 

 

Рассмотрим случайную величину q

при наличии сигнала s1:

q = q =

2

T

[s (t) + n(t)][s (t) s

(t)]dt .

(8)

N

 

1

 

1

1

0

 

 

 

0 0

 

 

 

 

 

 

 

 

 

Так как n(t) – гауссовский шум, s1

и

s0

детерминированные

функции, а все операции, проводимые над n(t) в (8), – линейные,

значит,

случайная

величина q1 будет распределена по нор-

мальному закону. Ее среднее значение

 

 

 

 

 

 

 

 

M {q } =

2E(1Rs )

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

N0

 

 

 

 

 

 

 

 

 

 

 

а дисперсия

 

 

 

 

4E(1Rs )

 

 

 

 

 

 

 

 

 

σq21 =

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N0

 

 

 

E

 

 

 

 

 

 

 

где R

=

1

 

T s

(t)s (t)dt – коэффициент взаимной корреляции ме-

 

 

s

 

 

 

0

1

0

 

 

 

 

 

 

 

 

 

 

s1(t)

и s0 (t).

 

 

 

 

жду сигналами

 

 

 

 

Аналогично при наличии сигнала s0 величина q0 будет распределена по нормальному закону со средним значением

M {q0} = − 2E(1Rs )

N0

9