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

Salnikov_A_P_Teoria_elektricheskoy_svyazi_Kons

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

 

 

 

 

A

 

 

 

 

 

 

,

 

 

Bc (t) U A(t) cos (t),

 

 

Bc (t) U 2 Bs2 (t)

 

 

Bs (t) A(t) sin (t),

 

 

 

 

(t) arctg

 

Bs (t)

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Bc (t) U

 

 

Отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w( A, ) w(Bc , Bs )

 

J

 

,

 

 

(5.10)

 

 

 

 

 

 

 

 

 

Bc (t) U

Bc (t) U

 

 

cos

Asin

 

 

 

 

 

 

 

 

 

 

 

 

 

J

 

A

 

 

 

 

 

 

 

 

A . (5.11)

 

Bs (t)

 

Bs (t)

 

 

sin

 

Acos

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для нахождения w(Bc , Bs )

обратимся к аналитическому СП

 

 

 

j ωct θ(t )

B(t)e

jθ(t )

e

jωct

 

B(t) cosθ(t) jB(t) sin θ(t) e

jωct

 

X (t) B(t)e

 

 

 

 

 

 

 

Bc (t) jBs (t) e jωct .

 

 

 

 

 

 

 

Из его выражения видно,

что Bc (t) и Bs (t) являются линейными пре-

образованиями центрированного нормального СП X(t):

B(t) Bc (t) jBs (t) X (t)e jωct ,

Bc (t) Re X (t)e jωct ,

Bs (t) Im X (t)e jωct

и, следовательно, имеют нормальное распределение с дисперсиями

ζ2Bc ζ2Bs ζ2X .

Докажем их некоррелированность (а следовательно и независимость) в совпадающие моменты времени

 

 

 

 

 

 

 

 

 

B

 

(t ,t ) B (t )B (t ) B2

(t ) cosθ(t )sin θ(t )

B ,B

1 1

c 1

s 1

1

 

1

1

 

 

c

s

 

 

 

 

 

 

 

 

 

 

 

 

1

B2

(t )sin 2θ(t )w(B,θ)dBdθ

1

B2

(t )w(B)dB

 

 

2

 

1

1

 

 

2

 

1

 

0

0

 

 

w( B)w(θ)

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

12 B2 (t1)sin 2θ(t1)

1

sin 2θ(t1)dθ 0 .

0

w(θ)

0

Здесь учтено, что B(t) и θ(t) – огибающая и фаза нормального СП являются, как выше установлено, независимыми.

Таким образом,

 

 

 

 

B2

B2

 

 

 

 

Acos U 2

A2 sin2

 

 

1

 

 

c

s

 

1

 

 

 

 

 

w(Bc , Bs ) w(Bc )w(Bs )

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

e

 

X

 

e

 

X

 

2πζ2X

 

 

 

2πζ2X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и с учетом (5.10) и (5.11) получаем

31

 

A

 

 

A2 U 2 2 AU cos

 

 

w( A, )

 

 

 

2

 

 

 

 

e

 

X

.

(5.12)

2πζ

2

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку выражение (5.12) невозможно представить в виде про-

изведения одномерных функций

w( A, ) w( A)w( ) , то можно сде-

лать вывод о зависимости процессов A(t) и (t) .

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

A

 

A2

U 2

AU cos

 

w( A) w( A, )d

 

 

 

2

e

 

2

d .

 

 

e X

 

ζ X

2πζ

2

0

X

 

0

 

 

 

 

 

 

 

 

Интеграл вида

1 ea cos d I0 a

0

известен в математике как модифицированная функция Бесселя нулевого порядка. С его учетом окончательно имеем

 

A

A2

U 2

 

 

AU

 

 

 

 

 

2

 

 

 

w( A)

 

 

 

 

 

 

 

 

( A 0) .

 

 

e

 

X

I0

 

 

 

 

 

 

 

(5.13)

2

 

 

 

2

 

ζ X

 

 

 

 

 

ζ X

 

 

 

Выражение (5.13) называют обобщенным распределением Рэлея или

распределением Райса. Графики этого выражения приведены на рис. 5.10 для следующих частных случаев:

 

 

 

 

A

 

A2

 

 

 

 

 

 

1)

U = 0

w( A)

e 2X – обычное распределение Рэлея,

2

 

 

 

 

ζ X

2)

ζ X

0

w( A) δ A U – случай отсутствия в Y(t) СП X(t),

3)

ζ X

U X

– обобщенное распределение Рэлея (Райса).

w(A)

(A-U)

1

2

3

0

 

2

3

4

5

6

7

8

9 U 10

A

 

Рис. 5.10. Обобщенное распределение Рэлея для частных случаев

 

32

Из

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

U / ζ X

тем правее смещен максимум плотности вероятности и тем

симметричнее (ближе к нормальному распределению) кривая w( A) .

Выводы

1. Если мгновенные значения центрированного СП X(t) имеют нормальное распределение, то его огибающая A(t) распределена по закону Релея

 

A

 

A2

 

 

 

w( A)

e 2X ,

2

 

ζ X

а фаза (t) равномерно

w( ) 21π .

2. Распределение огибающей аддитивной смеси центрированного нормального СП и гармонического сигнала подчиняется обобщенному распределению Рэлея (оно же распределение Райса)

 

A

A2

U 2

 

 

AU cos

 

 

 

 

2

 

 

w( A)

 

 

 

 

 

 

 

 

( A 0) .

 

e

 

X

I0

 

 

 

 

 

 

2

 

 

 

2

 

ζ X

 

 

 

 

 

ζ X

 

 

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

1.Сформулируйте задачу анализа прохождения СП через заданный функциональный узел.

2.Как вычисляют плотность вероятности w(y) реакции безынерционной цепи по известной плотности вероятности w(x) воздействия?

3.Как вычисляют математическое ожидание реакции безынерционной цепи на случайное воздействие X(t)?

4.Как вычисляют дисперсию реакции безынерционной цепи на случайное воздействие X(t)?

5.Как вычисляют функцию корреляции реакции безынерционной цепи на случайное воздействие X(t)?

6.Как вычисляют совместную плотность вероятности w(у1, у2; t) двух

СП Y1(t) и Y2(t), связанных известными функциональными зависимостями y1 f1 x1, x2 и y2 f2 x1, x2 с двумя другими СП X1(t) и

X2(t)?

33

7.Как меняется распределение нормального СП при его прохождении через линейную цепь?

8.Как меняется произвольное распределение СП при его прохождении через узкополосный фильтр?

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

10.Опишите процедуру корреляционного анализа прохождения СП через линейную цепь.

11.Дайте определение огибающей и фазы СП.

12.Дайте определения аналитическому СП, его математическому ожиданию, дисперсии и функции корреляции.

13.Каким условиям удовлетворяет стационарный аналитический СП?

14.Каково распределение огибающей центрированного нормального СП?

15.Каково распределение фазы центрированного нормального СП?

16.Каково распределение огибающей суммы центрированного нормального СП и гармонического сигнала?

17.Напишите аналитическое выражение закона Рэлея. Распределение какого СП он характеризует?

18.Напишите аналитическое выражение обобщенного закона Рэлея (закона Райса). Распределение какого СП он характеризует?

Рекомендации по проведению экспериментальных исследований прохождения случайных процессов через различные ФУ

Для закрепления знаний, полученных при изучении данного раздела рекомендуется выполнить в рамках виртуальной лаборатории работу № 20 «Прохождение случайных процессов через различные функциональные узлы» в полном объеме (рис. 5.11).

Обратите внимание на характер распределения СП на выходах одностороннего и двустороннего ограничителей - реальное проявление -функций в виде выбросов на гистограммах плотности вероятности распределения, соответствующих порогам ограничения. Убедитесь в нормализации СП с произвольными распределениями после их прохождения через ФНЧ и ПФ и в отсутствии нормализации после прохождения СП через ФВЧ (объясните почему?).

34

Рис. 5.11. Исследование прохождения СП через различные ФУ

6.Оптимальный прием дискретных сообщений

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

Дано:

1. Источник дискретных сообщений. Это значит, что известен ансамбль передаваемых сообщений

B bi m i 1, 2,..., m , где m – объем алфавита источника

иих статистика (распределение вероятностей) P(bi ) .

2.Модулятор. Это значит, что известны правила преобразования каждого сообщения в непрерывный сигнал и длительность сигнала T

bi si(t); i = 1, 2,…, m;

t (0, T).

3.Непрерывный канал. Канал задается своей математической моде-

лью, описывающей связь его реакции Z(t) с воздействием si(t) и канальными помехами N(t), например

35

Z(t) si(t) N (t);

i 1, 2,..., m

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

сигналами si(t) в приемнике определяются точно, иначе говоря, в нем осуществляется дискретизация времени функцией (t-kT), при которой границы тактов совпадают с границами сигналов.

Требуется:

Определить правило решения (решающую схему) вида

ˆ

t T, j 1, 2,...,m ,

z(t) bj

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

ˆ

о

z(t) СП Z(t) на каждом интервале Т следует принимать решение bj

переданном символе bi (при j = i имеет место правильный прием, иначе (при j ≠ i) – ошибочный).

 

Дадим геометрическую трактовку этой постановке задачи (рис.

 

 

 

6.1). Совокупность всех возможных

 

 

 

реализаций z(t) образует пространство

ˆ

 

 

принимаемых колебаний (обычно бес-

Bi

n2

 

конечномерное пространство Гильберта

 

ˆ

 

 

 

L2(T)) в котором присутствуют m раз-

 

Bj

ˆ

 

 

si

z2

Bk

личных векторов si передаваемых сиг-

 

 

 

 

 

sk

n1

налов si(t) (i = 1, 2,…, m). Выбор прави-

 

 

 

ла решения таким образом сводится к

 

 

z1

 

 

 

разбиению этого пространства на m не-

 

ˆ

Рис. 6.1. Геометрическая

пересекающихся областей Bi , каждая из

которых соответствует принятию реше-

трактовка постановки задачи

 

ния о передаче конкретного сообщения

bi (сигналом si(t)). На рис. 6.1. показаны две ситуации: 1) конец векто-

ра колебания

z (t) s (t) n (t)

попадает в область

ˆ

отведенную под

B

 

1

k

1

 

k

 

решение о передаче сообщения bk сигналом sk(t), что соответствует правильному приему; 2) конец вектора колебания z2(t) si(t) n2 (t)

попадает в область ˆ , отведенную под решение о передаче сообще-

Bj

ния bj сигналом sj(t), что соответствует ошибочному приему.

Разные правила решения (разные приемные устройства) различаются способом разбиения пространства принимаемых колебаний на

области ˆ . В этой связи возникает задача наилучшего разбиения, ко-

Bi

торое, очевидно, всегда существует в определенном смысле. Напри-

36

мер, если сообщение bi передается чаще сообщения bj и важно , чтобы как можно меньше передаваемых символов принимались ошибочно,

ˆ

расширить за счет области

ˆ

то следует область Bi

Bj . Наилучшее раз-

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

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

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

6.2. Критерии качества приема дискретных сообщений

6.2.1. Критерий идеального наблюдателя (критерий Котельникова)

Этот критерий требует обеспечения минимума средней вероятности ошибочного приема.

Для двоичной системы

min P min P(0)P(1/ 0) P(1)P(0 /1) ,

для m-ичной системы

m

m

 

 

 

 

min P min P(bi

 

ˆ

/ bi

)

 

) P(bj

,

i 1

j 1

 

j i

 

 

 

 

 

 

где

ˆ – условная вероятность j-ой ошибки при передаче

P(bj / bi )

i-го сообщения,

37

m

ˆ

P(bj / bi ) – условная вероятность любой ошибки при передаче j 1

j i

i-го сообщения,

Р – безусловная вероятность любой ошибки. Вычислим условную вероятность конкретной ошибки

ˆ

/ bi )

P(bj

 

ˆ

 

B j

 

 

w(z / bi

)dz ,

 

 

 

 

 

где w(z / bi ) w(z1

, z2

,..., zn

;t1, t2 ,...tn

/ bi ) – n-мерная условная плотность

вероятности (при разложении z в n-мерном евклидовом пространстве

по любому базису), а интеграл, вычисляемый по векторной перемен-

 

 

 

 

 

 

 

 

ной z , очевидно, n-кратный. Таким образом, критерий Котельникова

приобретает вид

 

 

 

 

 

 

 

 

m

 

m

 

 

 

 

min P min

 

P(b )

 

 

 

 

w(z / b )dz

,

(6.1)

 

i

i

 

ˆ

i 1

 

j 1 ˆ

 

 

 

 

B j

 

 

B j

 

 

 

 

j i

где находится варьированием областей ˆ . min P Bj

Минимуму средней вероятности ошибок соответствует максимум средней вероятности правильного приема (иная эквивалентная форма записи критерия Котельникова)

max(1 P) max

m

 

 

 

 

 

 

 

 

P(b )

 

 

 

 

w(z / b )dz

.

(6.2)

 

i

 

i

 

ˆ

i 1

 

ˆ

 

 

 

 

 

Bi

 

 

Bi

 

 

 

 

 

Учитывая, что демодулятор должен реализовать критерий (6.1)

ˆ

на основе анализа единственной реа-

или (6.2), принимая решение bj

 

рассмотрим апостериорную вероят-

лизации z на интервале 0 – Т,

 

 

 

 

 

 

 

 

 

ность вида P(bi / z) , т.е. вероятность того,

что при приеме сигнала

передавалось сообщение bi . Очевидно, что максимум сред- z(t) z

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

реализацию принятого колебания z(t) относить к той области ˆ , для

Bi

которой апостериорная вероятность P(bi / z) максимальна, т.е. реше-

ˆ

принимается при совместном выполнении совокупно-

ние в пользу bi

сти неравенств

 

 

 

 

j i .

 

P(bi / z) P(bj

/ z), j 1, 2,..., m

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

38

ˆ

arg

 

 

 

 

 

 

i

 

 

max P(b / z)

 

(6.3)

 

 

 

 

i

.

 

 

 

 

 

i

 

 

 

 

Для выполнения анализа (6.3) воспользуемся известной формулой

Байеса

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P(bi / z)w(z) P(bi )w(z / bi

) .

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w(z / b )

 

max P(bi / z) max P(bi )

 

i

,

 

 

 

 

 

 

 

w(z)

 

а выражение (6.3) принимает вид

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

i arg

 

 

max P(b )w(z / b )

(6.4)

 

 

 

i

i

 

 

 

 

 

i

 

 

 

 

(безусловная плотность вероятности

здесь исключена, т. к. она

w(z)

не зависит от i и, следовательно, не влияет на решение).

В развернутом виде критерий (6.4) можно записать в виде системы из m-1 неравенств

P(bi )w(z

/bi ) P(bj )w(z / bj )

j 1,2,...,m; j i

ˆ bi ,

или

 

 

 

 

 

 

 

P(bj )

 

 

 

w(z / b )

 

 

ˆ

 

i

 

 

 

 

 

 

 

( j 1,2,...,m;

j i) bi .

 

P(bi )

 

w(z / bj )

 

 

 

 

 

 

 

 

Условную плотность вероятности w(z / bi ) , рассматриваемую при

 

 

 

 

 

 

известном после приема векторе z как функцию аргумента bi, назы-

вают функцией правдоподобия гипотезы о передаче сообщения bi, а

 

 

 

 

 

 

 

 

 

 

w(z / b )

 

i, j

 

i

)

- отношением правдоподобия двух гипотез о переда-

 

 

 

w(z / b

 

 

 

 

j

 

 

че сообщений bi и bj. С учетом этого критерий Котельникова можно записать в виде:

 

 

P(bj

)

 

ˆ

 

если i, j

 

 

 

( j 1,2,..., m;

j i) , то решение b .

(6.5)

 

 

P(bi )

 

i

 

 

 

 

 

 

Рассмотренный критерий Котельникова обладает следующими особенностями:

1)требует знания априорных безусловных вероятностей отдельных сообщений P(bi ) ;

39

2) безразличен к виду ошибок ˆ (все виды ошибок одинаково

P(bj / bi )

нежелательны), что приводит к росту ошибок при приеме менее вероятных сообщений, а они являются более информативными.

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

Полагая, что все передаваемые сообщения равновероятны

P(b1 ) P(b2 ) ... P(bm ) m1 ,

из (6.5) получим

если i, j

1

( j 1,2,..., m;

ˆ

j i) , то решение bi .

Удобно помимо гипотез о передаче сообщений bi (i = 1, 2,…, m) ввести еще одну «нулевую» гипотезу о том, что никакое сообщение

(сигнал) не передавалось, т. е. принятое колебание является реализа-

цией только помехи

z(t) n(t)

 

 

 

 

 

(z n) . Обозначим отношение прав-

доподобия

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w(z / b )

i ,

 

 

 

 

i

 

 

 

 

 

 

 

w(z / 0)

 

 

 

тогда правило решения можно записать в виде

 

если i

j , при всех

 

ˆ

 

j i , то решение bi

 

или

 

 

 

 

 

 

 

 

ˆ

 

max

 

 

 

i arg

 

 

(6.6)

 

 

 

i

.

 

 

 

 

i

 

 

 

Критерий (6.6) называют критерием максимального правдоподобия. Он совпадает с критерием Котельникова при равных вероятностях передаваемых сообщений.

6.2.3. Критерий минимального среднего риска (байесовский критерий)

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

сумму условных вероятностей ˆ в его выражении (6.1) с зара-

P(bj / bi )

нее назначенными весами (ценой, платой) Li,j. Средневзвешенная

сумма условных вероятностей ˆ при передаче сообщения bi

P(bj / bi )

(обычно называемая условным риском) имеет вид

40