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

Salnikov_A_P_Teoria_elektricheskoy_svyazi_Kons

.pdf
Скачиваний:
535
Добавлен:
15.03.2015
Размер:
3.88 Mб
Скачать

m

 

ˆ

 

 

 

 

 

 

/ bi ) ,

 

 

 

Ri Li, j P(bj

 

 

 

j 1

 

 

 

 

 

 

а сам критерий требует

 

 

 

 

 

 

m

m

 

 

 

 

 

min R min P(bi

 

 

ˆ

/ bi )

 

 

) Li, j P(bj

,

(6.7)

i 1

j 1

 

 

Li ,i 0

 

 

 

 

 

 

 

 

 

где R – средний риск.

При использовании этого критерия оптимальной считается решающая схема, обеспечивающая наименьшее значение среднего риска

(6.7).

Из критерия минимального среднего риска, как наиболее общего, вытекают оба вышерассмотренных критерия:

1) критерий Котельникова при Li,i = 0 и Li,j = const (i ≠ j ),

2) критерий максимального правдоподобия при Li,i = 0 и

Li, j Li

 

1

.

 

P(b )

 

 

i

 

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

щений P(bi ) и сложностью объективного определения весовых коэффициентов Li,j.

6.2.4. Критерий Неймана-Пирсона

Критерий Неймана-Пирсона применяется в двоичных системах в ситуациях, когда невозможно определить априорные вероятности отдельных сообщений, а последствия ошибок разного рода неодинаковы. Такая ситуация типична для радиолокации, где осуществляется зондирование пространства узким радиолучом и прием отраженного от цели сигнала. При этом имеют место две ситуации: 1) наличие цели

– колебание на входе приемника содержит сигнал в аддитивной смеси с помехой (с неизвестной априорной вероятностью P(b1)), 2) отсутствие цели – на входе приемника действует одна помеха (с вероятностью P(b0) = 1 – P(b1)). Задача приема – обнаружение сигнала на фоне помех. При ее реализации возможны два вида ошибок:

1) пропуск цели (цель есть, но отраженный сигнал не обнаружен)

с условной вероятностью ˆ ;

P(b0 / b1)

2) ложная тревога (цель отсутствует, но принято решение о нали-

чии отраженного сигнала) с условной вероятностью ˆ .

P(b1 / b0 )

Очевидно, что последствия этих ошибок сильно различаются.

41

В таком случае целесообразно стремиться к уменьшению условной вероятности ошибки, вызывающей особо тяжелые последствия (пропуск цели), что можно сделать только за счет увеличения вероятности ошибки другого вида (ложной тревоги). Ясно, что это можно делать до определенной степени, т. к. слишком большая вероятность ложной тревоги приведет к ощутимым экономическим потерям и к подрыву доверия к системе в целом. Разумный выход – зафиксировать вероятность ложной тревоги на выбранном уровне ε

 

 

 

P(b1 / b0 ) w(z / b0 )dz ε ,

(6.8)

ˆ

B1

и затем минимизировать вероятность пропуска цели

min P(b0 / b1 ) min w(z / b1 )dz

ˆ

B0

min 1

w(z

ˆ

B1

/ b1 )dz . (6.9)

Минимизация (6.9) при заданной величине (6.8) достигается, если решение о наличии цели принимается при выполнении неравенства

 

 

 

1

w(z / b )

λ(ε) ,

1

 

w(z / 0)

 

где λ(ε) – пороговый уровень, определяемый заданной вероятностью ложной тревоги.

Контрольные вопросы

1.Сформулируйте задачу оптимального приема дискретных сообщений.

2.Дайте геометрическую трактовку задаче оптимального приема дискретных сообщений.

3.Что называют правилом решения (решающей схемой) демодулятора?

4.Что такое идеальный (оптимальный) приемник дискретных сообщений?

5.Что понимают под потенциальной помехоустойчивостью приема дискретных сообщений?

6.В чем суть теории потенциальной помехоустойчивости? Когда и кем были заложены ее основы?

7.Какой смысл вкладывают в понятие критерия качества приема дискретных сообщений? Перечислите известные Вам критерии.

8.В чем суть критерия идеального наблюдателя (критерия Котельникова)?

9.Укажите особенности критерия Котельникова.

42

10.Что представляет собой критерий максимального правдоподобия? Как он соотносится с критерием Котельникова?

11.Расскажите о критерии минимального среднего риска. В чем его общность?

12.При каких условиях критерий минимального среднего риска совпадает с критерием Котельникова?

13.При каких условиях критерий минимального среднего риска совпадает с критерием максимального правдоподобия?

14.В чем сущность критерия Неймана-Пирсона? В каких случаях целесообразно его использование?

6.3.Синтез оптимального демодулятора при известном ансамбле сигналов (когерентный прием)

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

 

Постановка задачи:

Известны:

 

1.

Ансамбль сигналов на выходе модулятора

 

{si(t)}; i = 1, 2,…, m; t (0, T).

2.

Непрерывный канал

 

 

Z t si (t) N (t);

i 1, 2,..., m ,

где N(t) – квазибелый нормальный шум, т. е.

 

N

 

const f (0, F )

 

GN

( f )

0

О

f (0, F )

.

 

 

 

 

3.В качестве критерия качества приема задан критерий максимального правдоподобия (6.6)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w(z / b )

ˆ

arg

 

max

arg

 

 

 

i

 

i

 

 

max

 

 

 

 

 

i

 

 

 

 

 

i

 

 

 

i

w(z / 0)

 

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

Решение

В основу решения положим выражение заданного критерия качества приема, для чего рассмотрим входящие в него функции правдоподобия гипотез:

1) о наличии во входном колебании z(t) i-го сигнала [z(t) = si(t) + n(t)]

43

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w(z / bi ) w(z1, z2

,..., zn

;t1, t2 ,..., tn / bi ) ,

 

2) об отсутствии в нем какого-либо сигнала [z(t) =

n(t)]

 

 

 

 

 

 

w(z1, z2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w(z / 0)

,..., zn

;t1, t2 ,..., tn / 0) ,

 

где n

T

,

t

1

,

n 2FT .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

2F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Начнем с последней. Учитывая, что сечения квазибелого шума,

разделенные

интервалами

t

 

1

 

,

 

не коррелированны, а в силу

 

 

 

2F

нормального распределения шума и независимы, получим

 

 

 

 

 

n

 

 

 

1

 

 

 

 

1

 

n

 

 

w(z / 0) w(zk

/ 0)

 

 

 

 

 

exp

 

 

 

zk2 .

 

 

 

 

 

 

n

2

 

 

 

 

 

 

 

 

 

 

k 1

 

 

 

2πζ

 

 

 

k 1

 

Поскольку СП Z(t) = si(t)+ N(t) отличается от шума N(t) только известным, а потому неслучайным сигналом si(t), играющим роль математического ожидания Z(t), то

 

 

1

 

 

1

 

n

2

 

w(z / bi )

 

 

 

exp

 

 

zk

si,k

,

 

 

n

2

 

 

2πζ

 

 

k 1

 

 

где использовано обозначение si,k = si(tk).

В итоге отношение правдоподобия гипотез о наличии и отсутствии сигнала принимает вид

i

или с учетом

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

w(z / b )

 

 

1

 

 

1

 

 

 

 

 

zk2

 

 

 

 

 

i

exp

 

 

 

 

 

 

 

 

2

2

 

 

w(z / 0)

 

 

 

 

k 1

 

 

2

2NO F

NO

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

n

 

 

i

 

exp

 

1

zk2 t

1

zk

 

 

 

 

 

 

 

NO k 1

 

 

 

 

NO k 1

 

n

 

 

zk

si,k

2

k 1

 

 

si,k 2 t .

Перейдем к белому шуму, сняв ограничение на ширину его спектра (F ). Иначе говоря, от евклидова пространства перейдем к гильбертовому. При этом

t

1

 

dt,

si,k si (tk ) si (t),

zk z(tk ) z(t),

 

2F

 

 

 

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

T

 

 

 

 

 

 

 

1

1

 

 

 

 

i

exp

z2 (t)dt

z(t) si

(t) 2 dt .

(6.10)

 

 

 

 

 

 

 

NO 0

NO 0

 

 

 

44

Синтезируемый демодулятор должен принимать решение в поль-

зу ˆ , обеспечивающего максимум выражения (6.10), или, что эквива- bi

лентно, максимум показателя экспоненты в нем

max

 

max

T

z

2

(t)dt

T

2

dt

 

ˆ

 

i

 

 

 

z(t) s (t)

 

b

(6.11)

 

 

 

 

 

i

 

i .

i

 

i

0

 

 

 

0

 

 

 

 

 

Нетрудно видеть, что максимум (6.11) достигается при минимуме вычитаемого

max

i

i

min

T

z(t) s

(t) 2 dt

 

 

i

 

 

i

0

 

 

 

mind

2

 

ˆ

 

 

z, s

b

(6.12)

 

 

i

i .

i

Демодулятор оптимальный по критерию максимального правдо-

подобия принимает решение в пользу того символа ˆ , сигнал si(t) ко- bi

торого отстоит от принятого колебания z(t) на меньшее расстояние. Рассматривая выражение (6.12) как алгоритм обработки принято-

го колебания z(t) приходим к схеме демодулятора, представленной на рис. 6.2.

Σ

-s1(t)

Г1

Σ

z(t)

-s2(t)

Г2

Σ

-sm(t)

Гm

КВ

 

Т

 

 

 

 

 

0

 

 

 

 

КВ

 

 

Т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

РУ

ˆ

 

 

 

 

0

 

min

bi

 

 

 

 

 

 

 

 

 

 

 

 

 

(по i)

 

j-ая ветвь

 

 

 

 

 

 

 

 

 

 

 

 

КВ

 

Т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kT

 

Рис. 6.2. Схема демодулятора, реализующего алгоритм (6.12)

45

Другую форму алгоритма можно получить из выражения (6.11)

max T

z2 (t)dt T

z(t) s

(t) 2 dt

 

 

 

 

i

 

 

 

i

0

0

 

 

 

 

max T

z2 (t) z2 (t) 2z(t)s (t) s2

(t) dt

 

 

 

i

i

 

 

 

i 0

 

 

 

 

 

 

T

 

1 T

 

 

 

 

 

1

 

 

 

 

 

 

 

2

 

ˆ

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

max (z, s )

 

 

 

 

 

 

 

 

 

 

 

max

 

z(t)s (t)dt

 

s

(t)dt

 

 

 

 

 

 

s

 

 

 

 

b

 

i

 

i

 

 

 

i

 

 

 

 

 

i

 

 

 

 

 

i ,

i

0

 

2

0

 

 

 

i

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ei

 

 

или

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

1

 

 

 

 

i

arg

 

max

(z, s )

 

E

 

,

(6.13)

 

 

 

i

2

i

 

 

 

i

 

 

 

 

 

где Ei – энергия i-го сигнала.

Схема оптимального демодулятора, реализующего алгоритм (6.13), приведена на рис. 6.3. Поскольку в каждой ветви такого демо-

дулятора присутствует вычислитель скалярного произведения (z, si )

коррелятор, то его называют демодулятором на корреляторах (активных фильтрах).

Х

 

 

Т

 

 

 

 

 

 

 

 

 

 

0

 

s1(t)

Г1

z(t) Х

s2(t)

Г2

Х

sm(t)

Гm

Т

0

j-ая ветвь

Т

0

Σ

12 E1

Σ

 

ˆ

 

РУ bi

 

1

max

 

E2

2

 

 

 

 

(по i)

Σ

 

1

kT

2

Em

 

 

Рис. 6.3. Оптимальный демодулятор на корреляторах (6.13)

46

Если использовать сигналы равных энергий, то алгоритм (6.13) и схема демодулятора (рис. 6.3) существенно упрощаются (рис. 6.4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

ˆ

 

 

arg

max

z(t)s (t)dt

 

 

 

 

i arg

 

max(z, s )

 

 

.

(6.14)

 

 

 

 

 

 

 

 

i

 

 

 

i

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

i 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Х

 

 

 

 

 

Т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s1(t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Г1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Х

 

 

 

 

 

Т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z(t)

 

 

 

 

 

 

 

 

 

 

 

 

 

РУ

 

bi

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

max

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s2(t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(по i)

 

 

 

 

 

 

 

Г2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j-ая

ветвь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Х

 

 

 

 

 

Т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sm

 

(t)

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Гm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 6.4. Оптимальный демодулятор на корреляторах (6.14)

 

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

6.3.2. Синтез оптимального когерентного демодулятора на согласованных фильтрах

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

Как известно, реакция линейного фильтра на воздействие z(t) вычисляется с помощью интеграла Дюамеля

t

y(t) z(η)g(t η)dη

0

47

Потребуем, чтобы в заранее выбранный момент времени t0 значение этой реакции y(t0) с точностью до коэффициента совпало со скалярным произведением (6.14)

t0

T

yi (t0 ) z(η)gi (t0 η)dη a z(η)si (η)dη .

0

0

Как видно, это достигается при gi (t0 η) аsi (η) и t0 T . После

замены переменных t0 η t получаем

 

gi (t) аsi (t0 t) .

(6.15)

Фильтры, обладающие такими импульсными характеристиками, называют согласованными (СФ) с соответствующими сигналами.

На рис. 6.5 изображены сигнал длительностью Т и импульсные характеристики согласованных с ним фильтров для t0 = Т и t0 > Т, из которых видно, что импульсная характеристика согласованного

фильтра является «зеркальным

отражением»

сигнала относительно

 

s(t)

 

момента времени 0,5t0.

 

 

Таким образом, фильтры с

 

 

 

 

 

 

импульсными характеристиками

 

 

 

(6.15) вполне могут заменить кор-

 

 

 

реляторы в

ветвях

оптимального

0

T

t

демодулятора (рис. 6.3 и 6.4), если

 

gСФ(t)

 

решения принимать

по отсчетам

 

 

 

 

 

 

их реакции yi(kT) (рис. 6.6).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

СФ1

 

 

 

0

t0 =T

t

 

 

 

 

 

 

 

 

 

 

 

 

 

gСФ(t)

 

z(t)

 

СФ2

 

 

ˆ

 

 

 

 

 

 

bi

 

 

 

 

 

 

 

РУ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

t0 >T

t

 

 

СФm

 

 

 

 

 

 

 

 

 

 

 

 

 

kT

Рис. 6.5. Сигнал и импульсные харак-

 

 

 

 

теристики согласованных фильтров

Рис. 6.6. Демодулятор на СФ

Свойства согласованных фильтров

1. Импульсная характеристика СФ является «зеркальным от-

ражением» сигнала, с которым он согласован, относительно момента времени 0,5t0 (с точностью до постоянного коэффициента)

gСФi (t) аsi (t0 t) .

48

Это свойство было положено в основу определения СФ (6.15).

2. Передаточная функция СФ

HСФ ( jω)

После замены t0 t = ,

 

 

gСФ (t)e jωt dt a s(t0 t)e jωt dt .

 

 

t = t0 – ,

dt = d , при t -

 

 

HСФ ( jω) a s(η)e jω t0 η dη ae jωt0 s(η)e jωηdη ae jωt0 S* jω .

 

 

Таким образом, передаточная функция СФ с точностью до мно-

жителя ae jωt0 совпадает с сопряженной спектральной функцией сигнала, с которым он согласован

HСФ ( jω) ae jωt0 S* jω .

Амплитудно-частотная характеристика СФ

HСФ (ω) HСФ ( jω) a S jω aS ω

с точностью до коэффициента а повторяет амплитудный спектр сигнала, с которым он согласован

Фазо-частотная характеристика СФ

СФ (ω) arg HСФ ( jω) arg S* jω ωt0 arg S jω ωt0

отличается знаком от фазового спектра сигнала, с которым он согласован (без учета слагаемого –ωt0).

3. Форма отклика СФ на «свой» сигнал (сигнал с которым он согласован)

yСФ (t) s(η)gСФ (t η)dη .

Учитывая, что из (6.15) вытекает gСФ (t) as(t0 t) , получим

yСФ (t) a s(η)st t0 )dη aBs (t t0 ) .

Таким образом, отклик СФ на «свой» сигнал с точностью до коэффициента совпадает с его корреляционной функцией, смещенной по оси времени на интервал t0 (рис. 6.7)

yСФ (t) aBs (t t0 ) .

49

s(t)

 

 

 

 

0

T

t

Bs(t)

0

T

t

yСФ(t)

0 t0 = T t

Рис. 6.7. Сигнал s(t), его корреляционная функция Bs(t) и отклик СФ yСФ(t)

Из полученного результата вытекают следующие выводы:

Отклик СФ на «свой» сигнал с точностью до постоянного коэффициента совпадает с его корреляционной функцией.

Длительность отклика на «свой» сигнал всегда равна 2Т.

СФ не восстанавливает форму сигнала, искаженного шумом. Его

задача создать один отсчет y(t0), по которому можно наилучшим образом судить о присутствии на входе «своего» сигнала.

4.СФ обеспечивает наибольшее отношение сигнал/шум (с/ш)

на своем выходе при действии на входе аддитивной смеси «своего»

сигнала

и

центрированного

нормального

белого

шума

Z(t) s(t) N (t) со спектральной плотностью мощности N=NО/2. Докажем это, уточнив предварительно, что под отношением с/ш

на выходе СФ понимают отношение математического ожидания отсчета случайной реакции СФ Y(t) в момент времени t0 = T к корню из ее дисперсии

 

с

 

 

 

Y (T )

 

 

 

 

 

 

 

 

 

 

 

 

.

(6.16)

 

 

 

 

 

 

 

 

 

 

ш

ВЫХ

 

ζY

 

Рассмотрим произвольный линейный фильтр с передаточной

функцией H ( jf ) . Поскольку Y (T ) представляет собой отсчет реакции ys(T) на математическое ожидание воздействия, каковым является сигнал s(t), то

50