книги из ГПНТБ / Гинзбург, В. В. Теория синхронизации демодуляторов
.pdf6
ВЛИЯНИЕ СИНХРОНИЗАЦИИ НА ПОМЕХОУСТОЙЧИВОСТЬ
*8.1. Условная вероятность ошибки в оптимальных приемниках
Исследование вероятности ошибки Рош при стационарной ра боте приемника сводится, как видно ив {1 .6 ), ж нахождению ус ловной вероятности ошибки р(ср) при фиксированной фс и к ус реднению р(ф) по множеству значений фс.
Вероятность р(ф) определяется ф-лой (1.4). Для удобства ра счетов эту формулу запишем в виде
р(ф )= Г р (ф \ jc) ay (х) dX, |
(6.1) |
где х — вектор, описывающий вариант сигнала, |
а р (ф|х) — вероят |
ность ошибки ири действии нормальной помехи при условии, что х зафиксирован.
Условная вероятность ошибки при когерентном приеме. Бу дем считать, что в реальном когерентном приемнике (демодуля торе) начало интегрирования не связано с опорными колебания ми, определяющими варианты сигнала. Тогда РУ на l-и посылке принимает решение в пользу г-го варианта сигнала, если при всех i^bj выполняется неравенство
<|+Г |
ti~\-T |
|
Г x{t)st (t, y ) d t > |
[ x(t)Sj(t, ty)dt. |
(6.2) |
*i |
'U |
|
Здесь ti — момент формирования синхроимпульса, который счи тается совпадающим с началом интервала интегрирования и вы ражается через фс по ф-ле ( 1 Л); sj{t, ф) — i-ый вариант сигнала, зависящий от текущей фазы ф опорного колебания; x(t) — при
нятый сигнал, равный |
(6.3) |
x{t) = z(t) +n(t), |
■где n(t) — слагаемое, обусловленное помехой; z(t) — слагаемое, обусловленное полезным сигналом.
С учетом (6.3) |
неравенство (6.2) принимает вид |
|
|
% i ( ф. |
Ф ) — 2 / ( ф , ф ) > 7 У / ( ф , ф ) — |
ф ), |
(6.4) |
150
|
ti+T |
|
где |
2<(<р, ij))= j z(t) s( (t, t|i)dt\ |
(6.5) |
U
</+r
N t (<p,$)= J n(f)st (t, ф)Л. (6 .6) U
Если помеха n(t) — нормальный стационарный процесс с ну левым средним, то распределение величины JV,- нормально, имеет нулевое среднее и не зависит от ф. Дисперсию ее обозначим че рез а2,(ф). Тогда дисперсия разности АЛ^(ф) =М,-(ф, t|>)—УДф, ф) равна
о], (ф) = О» (Ф) + |
О* (Ф) - 2К „ (ф), |
(6.7> |
где /Сгj ('ф) — корреляция величин |
и N+ |
причем |
Процесс z(t) определяется передаваемой информацией, |
при ограниченной полосе канала связи на текущее значение z(ty влияет не только данная, но и соседние посылки сигнала. В реаль ных условиях можно ограничиться учетом, помимо данной, пред шествующей и последующей посылок. Тогда величины Z* при фик сированных ф и ф являются детерминированными функциями ин формационных символов r/_1, Г[ и гi+i, переданных соответственно*
на |
(/—1)-й, l-й и (7+1)щ |
посылках, т. е. Zf(<p, ф )= 2;(ф , |
ф, гг_,, |
г и |
r l + i ) . |
|
|
|
Обозначим разность в правой части (6.4) через |
|
|
|
А гг/(ф, ф , г , , |
ri+1) = Z ,( ф, ф) — г / (ф, ф). |
(6 .8) |
Решение в пользу г-ro варианта сигнала будет правильным, ес ли /= г/. Вероятность такого события равна
Я ц { Ф. Ф. r i - v П+i) = Р 1А 2 </(Ф’ |
Ф- r i - v *. H+1) > A ^ / W ] |
= |
|
= F[ht,( ф, ф, |
rt_ v |
г|+1)], |
(6.9> |
где F(x) — функция Лапласа, |
|
ф , rt_ v i, r[+1) |
|
д г ,/( ф , |
(6. 10) |
||
М Ф- Ф- ri - 1- ri+0 = |
|
оц (Ф) |
|
|
|
Усредняя величину 1—qi, по всем значениям символов г, г;_ц гi+i и предполагая, что символы на соседних посылках независимы, найдем условную вероятность ошибки при заданных ф и ф 1):
м м |
м м |
Р(ф- Ф) = X 2 |
X ^PiPmPnV — Яи (ф, Ф, П_,, Гж )], (6 .1 1 ) |
(=1 /= 1 |
m=l п= 1 |
1*1
•) В дальнейшем обе вероятности р(<р,ф) и р(<р) называются одинаково* «условной вероятностью ошибки», так как они легко отличаются по своим аргу ментам.
151
тде М — число вариантов сигнала на посылке, a qij определяется
(ф-лами (6.9), (6.10) при г/_1 = т , Г(+1 —п.
Заметим, что (6.11) и последующие соотношения определяют вероятность ошибки в М-позиционной системе. Вероятность ошиб ки, приведенная к одному двоичному подканалу [137], вычисляет ся иначе.
Усреднив вероятность р(ф, ф) по всем значениям фазы опор ного колебания, получим выражение для условной вероятности ошибки при условии, что фс равна <р:
Я |
(дс)dx, |
(6.12) |
Р(ф)= JР(ф. х) |
—Л
где w^(х) — плотность вероятности фазы ф.
Для приближенных расчетов ф-лу (6.11) можно упростить, за менив функцию Лапласа в (6.9) ее разложением в ряд Тейлора в окрестности условного математического ожидания
мм
Мф , ф) = 2 X РпРтКЛф- Ф. Ш, п). n=l m==1
Обозначив |
|
м |
м |
н а (ф. Ф) = |
X РпРт [h ‘> (ф’ ф’ т ’ ^ ~ h il (ф> 'М 2 |
п—\ т = 1
и ограничиваясь квадратичным членом указанного разложения, получим после преобразований
м |
м |
|
|
|
|
р (ф, ф) = 2 |
£ |
Pi\ 1 — F [hii ^ ^ |
+ |
— 7 = г #</(ф. Ф) ha(<p, |
Ф )Х |
i=i /=1 |
|
|
2 / 2 я |
|
|
|
|
|
|
||
|
|
X ехр ["-----Ф) ]} • |
|
(6.13) |
|
Расчеты |
по |
ф-ле (6.13) могут |
оказаться проще, чем по |
ф-ле |
(6 .1 1 ), особенно при больших М.
Условная вероятность ошибки при оптимальном некогерентном
приеме. В данном случае РУ на |
I-й посылке принимает решение |
в пользу t-го варианта сигнала, |
если при всех i ^ j выполняется |
неравенство |
|
х (t) s, (t) dt |
х (t) Sj (t) dt + |
|
(6.14) |
где S j ( t ) — преобразованный по |
Гильберту j-й вариант сигнала. |
152
Обозначив
U+т
|
Х , - ( ф . |
r t_ v r h rz+I) = |
j |
z(t) s M d t - , |
(6.15a) |
||||
|
|
|
|
|
|
4 |
|
|
|
|
у I(ф. г1- 1, ri> rl+i) = |
*‘+т |
Л |
(6.156) |
|||||
|
j |
2 ( 0 |
Sc(t) di- |
||||||
|
TV. = |
U+T |
|
Mi |
= |
U+T |
A |
(6.16) |
|
|
J |
n(0st (t)dty |
J |
n(0s,(0^ ; |
|||||
|
|
и |
|
|
|
и |
|
|
|
v?(cp, |
|
r /f /-/+1) = 1/2 = |
[Хг(ф, |
r,_„ |
гь гж )+ЛГ(]* + |
||||
|
|
+ |
|
(ф , r f _ l t |
r , , |
r l + l ) + M i f , |
( 6 . 1 7 ) |
||
представим неравенство |
( 6 . 1 4 ) в виде |
|
|
|
|||||
|
|
|
|
I/* > 1 / 2 . |
|
|
( 6 . 1 8 ) |
||
Решение в пользу t-го варианта сигнала будет правильным, |
|||||||||
если i— ri. Вероятность такого события |
|
|
|
||||||
|
|
|
) - |
О*, |
|
|
|
|
|
Ян (ф . |
|
|
оИо |
|
|
|
j, rl+x)dx2dxv (6.19) |
||
i—v r/+i> = |
v r i ’ |
ф . |
/ - |
1’ |
|
||||
где Wv (x1, |
jc2 | qp, |
r*_1, i, |
j, n+i) — условная совместная |
плотность |
|||||
вероятности величин V i |
и Vj. |
|
|
|
|
|
|||
Вычислим условную вероятность ошибки, определяемую вы |
|||||||||
ражением, аналогичным (6 .1 1 ), |
для системы связи с однократной |
||||||||
модуляцией |
с ортогональными |
в усиленном смысле |
сигналами |
||||||
1132], для которых при i=£j |
Л |
|
“ |
|
|
||||
U+T |
|
|
|
и+т |
|
*l+TA |
|
||
U |
si (0 s, (t)dt = |
h |
|
|
I * |
|
|||
|
Г st (t) Sj (t) dt |
|
(t)Sj(t)di = |
||||||
|
|
|
|
Ч+г |
Л |
|
|
|
|
=| Si(t)s{(t)dt = 0 ,
авеличины У, и Vj независимы и, так как TV* и М{ нормальны,, распределены по закону Райса [132].
При фиксированных ф, /7, i, j вероятность ошибки 1—фр (ф,
П- i, n+i) равна
о - |
= |
|
) |
Х |
|
О х, |
|
||
|
|
|
|
|
X ехр |
А + А |
*2 + Zl2 |
dxYdx2, |
(6.20) |
|
2ста |
2аг |
|
|
153
где о2 — дисперсия величин Ni и Л1*, h (x) — функция Бесселя мнимого аргумента нулевого порядка;
Z? |
■Z?(«P. г/ - р |
ri+1) = |
X?(cpi, |
rt_v rh |
rl+l) + |
|
|
+ YU Ф- 0 -i> О, гWJ |
|
(6.21) |
|||
Проинтегрировав (6.20), |
получаем |
|
|
|
||
•1 Qa (ф. |
ri—v ri+1) = |
ехР |
Z? |
z l ( ф . |
г-p |
'■г-н) |
4 а 2 |
|
|
|
|||
|
|
|
|
|
|
|
|
2 ,-(Ф. гг-1 >'■ |
r/+i) |
|
|
(6.22) |
|
|
|
|
|
|
|
Здесь г1г(х, у) — специальная функция, вычисленная в приложе нии 3.
Заметим, что при Zj —0 '(например, при точной синхронизации) гиз (6 .2 2 ) получается известное выражение для помехоустойчиво сти [132]
1— <7г/(ф- 'W |
ri+1 |
1 |
£-(ф, r,_lt », rl+l) |
|
= — ехр |
4 а 2 |
(6.23) |
||
|
|
2 |
|
|
Проделав выкладки, |
аналогичные |
проделанным |
при выводе |
|
б(6 .1 1 ), и приняв во |
внимание, что величины Vi при |
некогерент- |
яом приеме не зависят от начальной фазы опорного колебания,
найдем с учетом (6 .2 2 ) |
условную вероятность ошибки |
|
||
2 |
2 |
iPmPnехр |
(ф, m, /, n) |
|
1 |
|
|
||
p{(p)==iEi=\£т—1 n=l |
4a2 |
X |
||
|
|
|
|
|
Zy(cp, |
m , i, n) Z? (ф , |
m, i, n) |
(6.24) |
|
X V |
|
|
|
О вычислении функций Z;. Условные вероятности ошибки в -ф-лах (6.12) и (6.24) определяются через величины Zit несущие информацию о вариантах принимаемого сигнала с учетом влия ния рассинхронизации. При когерентном приеме и неограниченной полосе пропускания канала связи эти величины отличаются от вариантов переданного на трех соседних посылках сигнала только -сдвигом, определяемым фазой синхросигнала <р. Обозначим эти варианты через 5,(ф , ф, гм, п, П+0- Если соответствующий тако му значению сигнал z(t) в ф-ле (6.5) обозначить через y(t), то для канала с ограниченной полосой пропускания в установившем уся режиме имеем
г (0 = J У (t — td g Vi) dtu |
(6.25) |
о |
|
.где g ( t ) — импульсная реакция канала. Подставив (6.25) |
в (6.5) |
154
и поменяв в полученном выражении порядок интегрирования,, имеем
|
|
|
|
со |
|
|
|
|
2;(ф> |
Ф> |
rt_ v г,, |
г,+1) = |
j‘ St |
ф, |
r/—I, rh |
г<+1) g (tj div (6.26> |
|
|
|
|
|
0 |
|
|
|
|
Аналогичные |
соотношения |
нетрудно |
установить для |
величин |
||||
Xi и |
Уi, |
определяющих |
в соответствии с (6 .2 1 ) |
вариант |
сигнала |
для оптимального цегогерентного приема.
Заметим, что при неограниченной полосе канала величины Z: определяются только двумя посылками сигнала, попадающими в интервал обработки, так как взаимного влияния посылок при этом нет.
6.2.Вероятность ошибки при когерентном приеме сигналов
соднократной ФМ и ЧМ
Вероятность ошибки при фазовой модуляции. Условные веро ятности ошибки р(<р, ф) три ФРМ и ФМ связанны известным со отношением [56, 137]
Р ф р м (Ф . Ф) = 2 Р ф м (Ф . Ф) [ 1 — Р ф м (Ф. Ф)]- |
(6 ' 2 7 > |
которое при малых вероятностях ошибки можно заменить более простым
РфРм(ф> Ф) ~ 2 РФМ(ф| ф)- |
(6-28> |
Для нахождения условной вероятности ошибки при ФМ заме тим, что варианты опорного колебания, отличающиеся только зна
ком, равны (при 0<t<T) |
|
h (t, ф0) = — s2(t, фо) = a cos (и0t + фц), |
(6.29> |
и воспользуемся (6.5) —(6.11)
При неограниченной полосе частот интегралы Zj, которые в. этом случае обозначаются через Si, зависят только от сигнала на (7—1 )-й и /-й посылках, попавших в интервал интегрирования. При этом «условные» отношения сигнал/помеха йц не зависят от
Гi+i, т. е.
|
л</(ф. Ф- ri-v n+i) |
= М ф- |
ri-i)> 1 ф ] - |
|
||
При однократной модуляции h\i{ф, ф, ги ) = /г21 (ф, ф, r;-i), по |
||||||
этому |
|
|
|
|
|
|
|
р(ф, ф) = |
1 — 0,5^[Л1 2 (ф, ф, 1) — 0,5/^ [/г1 2 (ф, ijj, 2)]. |
(6.30)* |
|||
Отношения h:j выражаются через S* и N{. Подставив |
(6.29) в |
|||||
(6.5), |
найдем St: |
|
|
|
|
|
|
ЗДф, |
ф, |
1 , 1 ) = — 52 (ф, |
г|), 1 1) = |
0 ,5a2Tcos(i|5— 1(50), (6.31а> |
|
5Дф, |
i|j, 1 , |
2) = |
— S2(ф, 1(5, 1 , |
2) — 0,5а?Т |^ 1.-----5 -jco s^ — ф0) + |
155
|
|
— sin -5- cos (Ф + 2 лр0 — ф)j . |
|
(6.316) |
|||
|
|
л k |
2 |
|
I |
|
|
Здесь &=(йо/сйг= |
ио7’/2л; |
предполагается для простоты целым чис |
|||||
лом. |
|
|
|
|
|
|
|
Дисперсии величин АЛ(ф, лр) и -Vsj(<р, |
ф) одинаковы и равны |
||||||
|
|
‘i+t |
|
|
|
|
|
v- (Ф) = < |
j |
п (t) a cos (w01+ |
тр)dt |
> |
|
|
|
= | |
j В (tx— t) a2 cos (ш01 -f- ф) cos (co0 tl + |
ф) dtdtv |
(6.32) |
||||
где B(ti—t) — корреляционная функция помехи n(t). |
|
||||||
При белом шуме |
|
|
|
|
|
||
|
|
|
о? |
- 0 , |
|
|
(6.33) |
|
|
B (tv /) = - f 6 & |
|
|
|||
где a2/ — спектральная плотность шума, поэтому |
|
|
|||||
|
|
|
a2f а‘ |
|
|
|
(6.34) |
|
|
|
^(ф ) = - ^ - Г . |
|
|
|
|
Подставив |
(6.7), |
(6 .8 ), |
(6.31) и (6.34) |
в (6.10) |
и учитывая, что |
||
А 12(Ф) = Д 21(Ф) = |
—оМф), получаем |
|
|
|
|
||
|
|
*12 (ф. Ф. О = V 2 йсоэ(ф — ф0); |
|
(6.35а) |
|||
Ли (ф, |
Ф, 2 )= 1^2”* [(^ 1 -----cos (гр _ |
aPo) + |
-j- |
sin-2 -X |
|
||
|
|
х COS (ф + 2ф0 —ф)J , |
|
|
(6.356) |
||
где h2= a 2T/2of2, откуда на основании (6.30) |
|
|
|
РФМ(ф, Ф) = 1 — О.б/ЧуТЛсовф — тро)] — 0,5F{l/'2 h X
X [(1 — ф/зх) cos(гр — ф0) + (1 /Я&) sin (p/2)cos(p + 2ф0 — ф)]}. (6.36)
Как видно, из рис. 6.1, построенного по ф-ле (6.36), при ма лых k ( k = \, 2 ) вероятность ошибки зависит не только от пог решностей тактовой синхронизации ф и опорного колебания ф—яро, но и от формы сигнала, на которую влияет ф0. Из сопоставления рис. 6Ла и рис. 6.46, на котором показан сигнал S\(t, фо) при двух значениях ф0, следует, что чувствительность приемника к величи не фс зависит от распределения энергии сигнала по длине посыл ки. Если плотность энергии сигнала ib окрестности границ мала (ф0=90°), то приемник нечувствителен к ошибкам тактовой син хронизации. При больших k(k>3) энергия сигнала почти равно мерно распределена по длине посылки при любой начальной фа-
156
зе tyo, поэтому условная вероятность ошибки практически зависит только от ф и г|з фо и не зависит от ф0. Зависимость вероятности ошибки от отношения сигнал/помеха при разных <р и ф—ф0 пред ставлена на |рис. 6.1в. Как видно из ipwc. 6.1в, изменение этик вели-
Рис. 6.1. Условная (Вероятность ошибки три когерентном приеме:
J — Ф ="ф —гр0= 0 ; 2 — 1 |> -0 ,1 я , Ш— гр0“ 0; 3 — ф-0,1я, —tpo“ 0,1я; 4 — 0 ,2 л , ф—фо—0;
о — ф - 0 ,2 п , г|)—г|)о -0,1 л ;
--------- неограниченная п о л о са ;---------ограниченная полоса (Ло»3)
чин на 0 ,1 л (5% периода) эквивалентно энергетическим потерям
(0,5-1,0) дБ.
При ограниченной полосе канала условная вероятность ошиб ки зависит от формы частотной характеристики канала. При мем ее прямоугольной с (полосой ДF, пропускающей колебания с частотами <о0—й < ш < ю о+ й, где Q=nAF. Будем считать, что на текущую посылку сигнала влияют в основном только предшест-
157
вующая и последующая посылки *) и рассмотрим все варианты сигнала, соответствующие «тройке» посылок, а вне этих посылок продолжим варианты периодически. Последнее позволяет при на хождении функций Zi вместо (6.26) определить для каждого ва рианта сигнала Si(iр, ф, п-ь i, rl+l) =Si(t), раосматриваемого как
функция начала интервала обработки t— <pTI2n, его ряд Фурье и отбросить члены ряда с частотами, лежащими вне полосы канала
связи.
Первый из четырех вариантов функции Si(t), соответствующих различным сочетаниям информационных символов 4, г* и г;+ 1 (остальные варианты отличаются только знаком) является посто янной (рис. 6.3а), а при другие (рис. 6.36, в, г) получаются ин вертированием и сдвигом на величину, кратную Т, функции So(t) (рис. 6.36). Поэтому первый вариант функции Zi(t) то же является постоянной, а три другие получаются сдвигом и инвертированием функции Za(t), ранной сумме первых n0 — E(SAFT/>2)2) членов раз ложения So(t) в ряд Фурье (53]:
(6.37)
ft=1
*) Влияние последующей посылки объясняется тем, что вследствие задержки в канале начало отклика канала на эту посылку может существенно опережать интервал времени, на котором сосредоточена основная часть энергии посылки.
2) Е(х) означает целую часть х.
158
Таким образом, варианты функций Z* равны: |
|
||||||||
|
|
Zi(q>, |
ф, 1, 1, |
1) = — а*Т cos (ф — ф„); |
(6.38а) |
||||
z i(ф. Ф- 2- 1. 2) = y а2Тc o s — ’М — Т + 7 ^ х |
|||||||||
|
|
|
|
|
по |
, |
пл |
N |
|
|
|
|
|
|
sm2 — |
(6.386) |
|||
|
|
|
|
|
П=1 |
|
|
|
|
|
|
|
|
|
|
|
7 |
|
|
ZJqj, |
ф, |
2 , |
1 , |
1 ) = |
-----a2/ |
1cos (ф — ф0) ( ------- 1-+-Щ -Х |
|||
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
«о |
п |
П |
|
|
|
|
|
|
|
sin2 — |
п (ф + 2л) |
|
|||
|
|
|
X |
|
|
3 |
(6.38в) |
||
|
|
|
л=1 |
|
|
•cos — |
1— - |
||
|
|
|
|
|
|
|
|
|
|
z i(y, |
ф, |
1 , |
1 , |
2 ) = |
-----2-а2Т cos(ф— ф0) [ ---- L _ ) _ i i x |
||||
|
|
|
|
|
|
^ |
|
\ О |
Л* |
а п п sin2 —
хпS= 1 "" |
я« 3 со5 т ( ф ~ 2 я ) 1 • |
(6.38г) |
|
Для нахождения дисперсии величин jV,-(ф, ф) при ограничен ной полосе канала учтем, что корреляционная функция процесса Ni (как и раньше, имеется в виду изменение Nt в зависимости от t — q>T/2n) при неограниченной полосе имеет форму симметрично го треугольника [85, 136]. Поэтому корреляционную функцию про
цесса Ni(ф, |
ф) |
при т <Т можно представить с помощью функции |
|||||||||
а д . - |
|
|
|
х) _ 2 ±•2 а?Т [1 |
|
|
|
|
|
|
|
|
|
|
V |
+ Z0 (т)], |
|
|
|
|
|
||
откуда с учетом |
(6.37) |
|
|
л п |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
а2, а2Т |
|
„ п9 |
sin2 — |
, |
t ' = l , |
2. |
(6.39) |
||
ст2 (ф) = _1 |
_ |
J - + - S - V |
------ 2_ |
||||||||
|
|
|
|
3 |
я* Ц |
п* |
|
|
|
|
|
|
|
|
|
|
Л=1 |
|
|
|
|
|
|
Подставляя |
(6.38), |
(6.39) в |
(6 .1 0 ) и учитывая, |
как |
и выше, |
что |
|||||
А^)г(Ф) =/Сг1 (ф) = —o2i (ф), .получим: |
|
п* |
|
|
\ —1/2 |
|
|||||
|
|
|
|
|
/ |
|
|
|
|
||
М ф . Ф. |
1, 1) = ЛУТсоз(ф — Ф о )(^ -+ ^ S |
^ - s i n ’ T |
) |
’ |
|||||||
|
|
|
|
|
|
П=1 |
|
|
(6.40а) |
159