
Литература / Гришин Ю.П., Ипатов В.П., Казаринов Ю.М. Радиотехнические системы (1990)
.pdfPm = A0FB— средняя |
мощность |
помехи |
«(/) в полосе |
частот [ — FB, FB], |
Аналогично, |
в силу |
независимости |
слагаемых (3.38) |
дисперсия |
D {^} = 2FBTD {i2(/)}, где |
последний сомножитель есть дисперсия квадрата гаус совской случайной величины с нулевым средним. Из теории вероятностей известно, что для нормальных
величин |
Z){5y(z)}=5'y(z)-[5y(z)]2 = 2D2{5},(z)} |
и, |
а |
сле |
|||||
довательно, |
при гипотезе |
Но |
D {Q =2(2FBF)P2, |
при |
|||||
гипотезе |
Я, |
D {Q=2(2FBT)(Pm + Pc)2. |
Интегрируя |
и |
нор |
||||
мальные |
ПВ |
и |
^(^Я^ |
от |
до |
оо |
от |
||
— оо до |
£п |
соответственно, нетрудно |
прийти |
к |
выра |
||||
жениям |
|
|
|
|
|
|
|
|
|
Рлт=1-Ф(Л); рпс=1-Ф |
|
|
|
h |
|
|
|
||
* |
|
|
PJPw+\ |
|
|
||||
|
|
__ |
С/ 111* |
“Г * |
|
|
|||
в которых h = |
|
|
|
|
|
|
|
||
Как |
видно, при рлт = const |
вероятность пропуска |
рпс |
в энергетическом приемнике убывает с- ростом числа
независимых |
отсчетов |
2FBT на интервале наблюдения |
|
[О, |
Т ] и с увеличением |
отношения мощности сигнала Рс |
|
к |
мощности |
помехи Рш |
в полосе [ —FB, FB], |
§3.6. СТРУКТУРЫ И ПОКАЗАТЕЛИ РАЗЛИЧИТЕЛЕН ДЕТЕРМИНИРОВАННЫХ СИГНАЛОВ
Предшествующий материал относился только к обна ружению, представляющему собой частный случай раз личения двух сигналов, один из которых равен нулю. В настоящем и следующем параграфах будет рассмот рено различение ненулевых сигналов одинаковой энер гии. При этом за основу будет принято правило МП (2.10), оптимальное в том случае, когда критерием качества служит сумма условных вероятностей ошибок Рошусл (2-3), либо полная вероятность ошибки Рош при равных априорных вероятностях всех сигналов р^Х/М,
либо средний риск |
(2.1) |
при равновероятности |
сигналов |
и равной опасности |
всех |
ошибок П;1 = П, i^k. |
Подобные |
допущения упростят рассмотрение и придадут ему мак симальную физическую наглядность, обобщения же об суждаемых правил и структур на случаи неравных энергий и байесовских критериев с произвольными Пй, очевидны и у читателя, постигшего суть вопроса, затруднений не вызовут. Начнем с различения двух детерминированных сигналов.
60
Различение двух детерминированных сигналов. Согласно (2.10), (2.17), действующий по правилу МП различитель двух детерминированных сигналов s0(i) и $,(/) равной энергии Ео = £\ = Е должен принимать решение о присут ствии в колебании y(t) сигнала, имеющего с у(г) большую корреляцию:
Я, |
г |
|
1 = 0, |
|
|
Z1^ZO, 2;= |у(ф;(/)(1/, |
1. |
|
|||
Яо |
о |
|
|
|
|
Перепишем это |
правило в |
следующем |
виде: |
||
|
т |
й, |
|
|
|
z = zi-^o = f>'(0's'(0d/% о, |
|
(3.39) |
|||
|
0 |
я» |
|
|
|
где |
s(t) = sl(t) — s0(t). |
Таким |
образом, |
чтобы различить |
два детерминированных сигнала, достаточно вычислить единственную корреляцию z с разностным сигналом s(t), причем какое из решений будет вынесено — зависит только от знака z. Структура, реализующая правило (3.39) и ис пользующая коррелятор, показана на рис. 3.14. Как всегда, коррелятор можно заменить фильтром, который должен быть согласованным с разностным сигналом s(t).
Рассмотрим геометрическое толкование процедуры ра зличения, проясняющее попутно вопросы о вероятностях ошибок и о том, как лучше выбрать сами различаемые сигналы. Будем интерпретировать s0(t) и sl(t) как векторы
в некотором |
евклидовом пространстве |
(см. § 2.4), учтя, |
||
что длины |
последних x/pf(/)d/ = v/£), |
(=0, 1, |
в |
силу |
равенства энергий одинаковы: х/£)=ч/£, |
i = 0, 1. |
Так |
как |
в любом пространстве через два вектора с общим началом всегда проходит плоскость, векторы s0(P) и sl(t) можно считать расположенными на плоскости Р (рис. 3.15). До пустим, что к присутствующему на входе различителя сигналу s0(t) добавилась помеха и(/), которую можно также интерпретировать как вектор, но уже не обязательно
лежащий в |
плоскости |
Р. Тогда результирующий |
вектор |
||
у(1) = .?0 (/) + «(/) и решение будет |
выноситься |
в |
пользу |
||
сигнального |
вектора, |
ближайшего |
к у(/), т. е. |
имеющего |
Рис. .3.14
61
Рис. 3.15
от y(t) меньшее евклидово расстояние dyi - [y(z) — s((z)]2 dt
(cm. § 2.4). Понятно, что вектор y(z) ближе к тому из векторов 50(z), ^(О, к которому ближе егопроекция y'(z) на плоскость Р. Таким образом, если через середину
вектора j,(z)=j1(z) —s0(f) |
провести |
перпендикуляр (аа' |
на |
рис. 3.15), разделяющий |
Р на |
две полуплоскости, |
то |
получится следующая модель принятия решения: считается, что в y(z) присутствует тот сигнал, который лежит в той же полуплоскости, что и проекция y(t) на Р.
Отсюда уже нетрудно прийти к выражениям для вероятностей перепутывания сигналов. Действительно, как видно из рис. 3.15, сигнал s0(z) будет ошибочно принят за сигнал st(z) тогда и только тогда, когда проекция шумового вектора «(z) на направление s(z) окажется
больше половины |
длины |
d0l разностного |
вектора |
s(t) = |
|
= s1(t) — s0(f). Так |
как |
проекция |
n(z) на |
направление |
|
|
|
г |
|
|
|
j(z) равна 2ш/</01, |
где |
= f |
—скалярное |
про- |
|
|
|
о |
|
|
|
изведение (корреляция) n(z) и s(z), то для вероятности перепутывания s0(f) с Sj(z) имеем
со
p0l = P(Hl\H0)=P^ZJd0l>^ = J W(zm)dzm.
^о./2
Воспользовавшись’нормальностью гш как линейного пре
образования h(z) и |
тем, что среднее гш = 0, |
а дисперсия |
|
D {зш} =2V0£s/2 |
(см. |
§3.1), где Es — энергия |
разностного |
сигнала, причем |
d0г |
получим |
|
62
1 |
7 |
I |
I j |
, b / ^0 1 |
|
Poi= r=—=3= |
J |
exp---- <Ьш = 1-ф “r= |
|
||
2TiN°d°l diJ2 |
' 2 N°^ |
° V |
|
V^/lNo, |
|
Вероятность |
перепутывания |
уменьшается |
(3.40) |
||
с ростом |
длины разностного вектора s(t), т. е. с увеличением евкли
дова |
расстояния |
между з0(г) и зД/) |
|
do i = |
|
= '/2E-2$s0(f)sl(t)dt = ч/2£'(1 — р), |
|
в |
котором |
|
|
|
1 |
з0(г)зД/)(1/ |
(3.41) |
Р |
Е |
— коэффициент корреляции сигналов s0(t) и s^t). Теперь выражение (3.40) можно представить в более традиционном виде:
р01=1-Ф
где q = -J2EINo — уже известный параметр обнаружения, равный отношению сигнал/шум на выходе фильтра, согла
сованного |
с 3,(z), ПРИ гипотезе |
//,. |
|
В силу симметрии векторов з0(?) и гДО относительно |
|||
аа’ на рис. 3.15 |
вероятности |
перепутывания р01 и р10 |
|
одинаковы |
и полная вероятность ошибки (2.2) при равно |
||
вероятных |
сигналах (р0 =pY = 1 /2) |
||
|
|
|
(3.42) |
Выражение |
(3.42), как и |
рис. 3.15, отвечают на |
вопрос о наилучшей в смысле минимума Рош паре сигналов фиксированной энергии Е. Действительно, ко эффициент корреляции (3.41) как частное от деления скалярного произведения s0(t) и зДД на произведение длин есть косинус угла между векторами s0(t) и зД/) и, следовательно, минимален (р= —1), когда эти векторы противоположны. Таким образом, вероятность ошибки
(3.42) |
минимальна |
и равна |
|
Л>ш=1-Ф(<7) |
|
(3-43) |
|
для противоположных сигналов зД/)= — з0(/), |
почему такая |
||
пара |
и считается |
оптимальной в любых |
приложениях, |
63
где требуется различение двух детерминированных сигналов равной энергии. Тем не менее на практике по разным причинам нередко используют и неоптимальные, например
ортогональные, сигналы, для которых’ р = 0 и |
[см. (3.42)] |
Рош = 1-Ф(д1^). |
(3.44) |
Сравнивая выражение (3.44) с (3.43), нетрудно видеть, что применение ортогональных сигналов вместо противо положных требует для сохранения значения Раш в ч/2 раз большего значения q, т. е. двукратного увеличения энергии сигналов Е. Поэтому в классе детерминированных сигналов ортогональная пара имеет энергетические потери по отно шению к противоположной £дЕ = 3 дБ.
Легко видеть, что противоположную пару образуют
два радиосигнала, |
отличающиеся сдвигом |
на угол л фазы |
|
несущей частоты. |
Характерных примеров |
ортогональных |
|
пар значительно больше, и среди них |
такие, как |
два |
|
отрезка длительностью Т гармонического |
колебания |
ча |
стоты f = kfT (к — натуральное число), сдвинутые по фазе на угол +п/2; любые два сигнала, не перекрывающиеся по времени или по спектру; разнообразные фазоманипу-
лированные |
сигналы и |
пр. |
|
|
|
|
|
|
|
||
Различение М сигналов. При произвольном М раз- |
|||||||||||
личитель, |
придерживающийся правила |
МП, |
|
по-преж |
|||||||
нему |
считает |
присутствующим |
в |
у(1) |
|
сигнал, |
наи |
||||
менее удаленный от у(г) в смысле евклидова |
расстояния |
||||||||||
^ = Л[Я')-^')]2<1/ |
[СМ. (2.10), |
(2.16)] |
или, |
что |
при |
||||||
одинаковых |
энергиях |
сигналов |
равносильно, |
имеющий |
|||||||
с у(/) |
максимальную корреляцию zt [см. (2.17)]. Используя |
||||||||||
соглашение о символике решающих правил в |
(2.6) — (2.10), |
||||||||||
алгоритм МП можно записать в |
виде |
|
|
|
|
|
|||||
й, |
i=Q, 1, ..., M-l. |
|
|
|
|
|
|
|
(3.45) |
||
zk^Zi, |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
Структура, |
в |
осно |
|||
|
|
|
|
|
ву |
которой |
положено |
||||
|
|
|
|
|
правило |
(3.45), |
изобра |
||||
|
|
|
|
|
жена на рис. 3.16. Она |
||||||
|
|
|
|
|
содержит |
М каналов |
|||||
|
|
|
|
|
(по числу сигналов), вы |
||||||
|
|
|
|
|
числяющих М корреля |
||||||
|
|
|
|
|
ций |
наблюдаемой |
реа |
||||
|
|
Рис. |
3.16 |
|
лизации |
у(/) |
со всеми |
||||
|
|
|
сигналами. |
|
Решающий |
64
блок РБ отыскивает наибольшую из корреляций. Сигнал, поданный в качестве опорного на коррелятор, на выходе которого зафиксирован максимум z;, объявляется присут ствующим в y(t). Вместо корреляторов можно использовать Л/ фильтров, каждый из которых согласован со своим сигналом Число каналов в схеме можно уменьшить на единицу, если в каждом из них вычислять корреляцию не с 5f(z), а с Si(t) — s0(t), i=l, 2, ..., М-1. Тогда с выходов
всех М— 1 каналов |
будут |
сниматься |
разности zt — z0, |
|||
i= 1, 2, ..., М— 1, |
и логика |
работы РБ |
должна |
быть |
из |
|
менена: решение |
Но |
будет |
выноситься |
только |
если |
все |
z; —2о^0 (это значит, что при всех i>0 z0>z;). Если же среди разностей zi — z0 хоть одна положительна, то будет решено, что в y(t) содержится тот сигнал, для которого zk — z^Zj —z0, i= 1, 2, ..., M— 1.
Геометрические |
образы, |
использованные |
ранее, |
|||
можно |
распространить и на случай Л/>2, только те |
|||||
перь |
пучок |
векторов |
s0(t), sk(/),..., |
располо |
||
жится |
не на |
плоскости, |
а в |
ЛГ-мерном пространстве. |
Для того чтобы по возможности уменьшить вероятность перепутывания /-го сигнала с к-м, следует максима льно «раздвинуть» /-Й и к-й векторы. Однако если М>2, то, раздвигая два вектора, необходимо следить, чтобы каждый из них не приближался к какому-то третьему, поскольку все сигналы равноправны и каждый из них может с равной вероятностью присутствовать в y(z). Таким образом, оптимальный выбор М детерминирован ных сигналов сводится к поиску такой конфигурации пучка М векторов, в которой минимальное евклидово расстояние между парой векторов было бы максимальным:
min <7,k = max. Так как при равенстве энергий, т. е. длин
//к
векторов, dik = x/f[^(/)-^(0]2d/ = V2£(l-pjt), где pit — ко
эффициент корреляции /-го и к-го сигналов, то требова ние максимума минимального расстония тождественно условию минимума максимального коэффициента кор реляции в множестве (ансамбле) сигналов S= ...
..., sM_, (/)}: max p,t = min. Предельно достижимый минимум
*ki
максимального |
коэффициента корреляции устанавливается |
|||||
довольно легко. |
Просуммировав р,к |
по всем i и к, получим |
||||
|
|
1 С |
1 (Тм-1 |
I2 |
(3.46) |
|
Zp.fc=T S’.('U(')d'=T |
Е |
dz>0’ |
||||
I. |
к |
Ji. к |
L J L i = 0 |
|
|
|
3 |
Заказ |
31/3 |
|
|
|
оэ |
где неравенство следует из неотрицательности квадрата под интегралом. Кроме того, в сумме слева М слагаемых при i=k равны единице, а остальные-Л/(Л/—1) не больше
Рмакс» где Рмакс = тах pit. Поэтому М + М(М-1)рмакс^0
i/k
ИРмакс^-1/^-!)-
Конфигурацию из М векторов, в которой косинус угла между любой парой векторов равен — 1/(Л7—1), называют правильным симплексом. Если эти векторы взять в качестве М сигналов, то полученный детерминированный ансамбль (симплексный код) при равновероятности всех
.s,(z) обеспечит минимум полной вероятности ошибки Рош (2.2), что и решает вопрос об оптимальном выборе М сигналов. Примерами правильных симплексов служат при М=2 два противоположных вектора; при Л/=3— три
вектора, лежащих в |
одной плоскости и разбивающих ее |
на три равных угла |
по 120°; при М=4 — четыре вектора, |
концы которых являются вершинами правильного тетра эдра, и т. д. Попятно, что при Л/» 1 выполняется соотношение — 1/(Л7—1)~0, и поэтому при большом числе различаемых сигналов ортогональный ансамбль, в котором
Pik-O, i^k, |
практически нс проигрывает симплексному |
в значении |
Рош. |
Последовательность вывода точного выражения для вероятности ошибки различения М сигналов с произволь ными p,t такова [9]. ПВ системы случайных величин 20, zt, ..., zM~i есть Л/-мерный нормальный закон, для
задания которого достаточно знать средние всех |
и их |
|
корреляционную матрицу [см. |
(1.4)]. Для средних |
при |
истинности гипотезы Ht имеем |
2, = £р1(. Корреляционный |
же момент ьй и k-й корреляций равен N0Epik/2. После того как M-мерная ПВ найдена, ее Л/-кратнып интеграл по области i—Q, 1, ..., М— 1, позволяет получить вероятность правильного решения Р(Я(|Я() при условии истинности Я(. Сумма таких вероятностей, деленная на Л/ (с учетом равновероятности сигналов), будет полной вероятностью правильного решения Рпр, связанной с Рош (2.2) очевидным равенством Рош=1 — Рпр.
Получаемый таким образом Л/-кратный интеграл в ря де важных случаев удается свести к однократному. Так, для любых равнокоррелированных (равноудаленных) сигна лов (pik = p, *ik)
СО
*V \ |
|
|
( |
(-V + <7 л/1-р) -*d |
(3.47) |
- у) фМ ’1 |
— со
66
В практических расчетах выражение (3.47) используют редко из-за необходимости численного интегрирования. Чрезвычайно полезна оценка (3.47) сверху, для вывода которой будем считать, что истинна гипотеза Н,. При этом ошибка происходит всегда, когда истинно хотя бы одно из событий Zi>zh i^l. Вероятность ее Рош1, равная
вероятности объединения (J (г, >2,) событий 2, >2,, ///,
по теореме сложения вероятностей,
=(2,.>2()|яХ £ Р& >2,| Н,)~
1*
- Е >-’/) Л (-4>-i)|Ar()+
i<k l*i,k
и в силу известного неравенства Буля не больше первой суммы справа. Так как каждое слагаемое этой суммы есть вероятность перепутывания двух сигналов [$,(/) с •$■,(/)], то, согласно (3.42), для равноудаленных сигналов
При равновероятных сигналах (р,= 1/Л/) приходим к так называемой аддитивной границе полной вероятности ошибки
Г |
/ /1 _ р\ |
(3-48) |
Р°Ш=М ,?0 |
) |
При М—2. как видно из сравнения с (3.42), последнее соотношение является строгим равенством. Использование (3.48) при М>2 оправдывается, с одной стороны, асимп тотическим сближением его правой части и Рош по мере роста требований к качеству различения (Рот->0), а с дру гой—тем, что, выбирая необходимую энергию сигналов (минимальное значение q) исходя из правой части (3.48), разработчик всегда действует с известной перестраховкой, гарантируя удержание фактической вероятности ошибки ниже цифры, принятой им при расчете.
§3.7. РАЗЛИЧЕНИЕ СИГНАЛОВ СО СЛУЧАЙНЫМИ НАЧАЛЬНЫМИ ФАЗАМИ
Если различаемые сигналы известны, за исключением
начальной фазы (р, т. е. |
описываются моделью (3.9) |
3* |
67 |
s,(r, (p) = Re {S,(z)exp [/(2n/oz + (p)]}, то применению правила
МП (2.10) |
должно |
предшествовать |
вычисление |
ФП |
|
f/Mz) | |
согласно |
(2.19), |
т. е. ' |
усреднение |
ФП |
^(д/) | Ht, (р), |
построенной для |
детерминированных |
сиг |
налов с фиксированной фазой (р по всем ее возможным значениям с учетом априорной ПВ lT0((p). При равномер ной ПВ фазы И'о (ср) = 1/(2тт), |(р|^тг, такое усреденение уже проводилось в § 3.2 — там усреднялось ОП, отличающееся от ФП постоянным относительно г* и <р сомножителем.
Поэтому ФП |
повторяет результат (3.13) и с уче |
||
том равенства энергий всех различаемых сигналов |
|||
|
|
(17\ |
|
fV(y(t)\Hi) = cI0^ |
|
||
где с —коэффициент, |
содержащий сомножители, не зави |
||
сящие |
от |
/, а |
|
1 |
Г» |
• |
(3.49) |
2,= 2 |
Y(f)Si(t)dt |
О
— модуль корреляции комплексных огибающих принятого колебания y(t) и z'-го сигнала. Монотонность функции /0( ) на положительной полуоси позволяет перейти к достаточ
ной статистике |
Z, и записать правило МП в |
виде |
Zk^Zi, i=Q, 1, |
..., М-1. |
(3.50) |
Таким образом, оптимальный различитель М сигналов равной энергии со случайными начальными фазами должен вычислить все М величин (3.49) и, если максимальной из них является Zk, принять решите о присутствии в y(t) к-го сигнала. Это означает, что содержащимся в наблюда емом колебании y(t) считается тот сигнал, комплексная огибающая которого имеет наибольшую по модулю кор реляцию с комплексной огибающей y(z).
Стандартные способы нахождения Z, рассматривались в § 3.2. Один из них состоит в вычислении Z; как длины двумерного вектора, компонентами которого служат кор
реляции |
и z2i колебания |
y(t) с квадратурными составля |
ющими |
Sj(t) и 5,1(0 г’го |
сигнала 5,(/; (р). Техническим |
воплощением этого варианта было бы М-канальное устрой ство, каждый канал которого состоял бы из двух корреля торов и схемы объединения zlt, z2i по правилу (3.14).
68
Другая возможность ис пользована в схеме раз-
личителя на |
рис. 3.17, |
|
в которой Z; находится |
|
|
как огибающая на вы |
|
|
ходе фильтра, |
согласо |
|
ванного с z'-м сигналом |
|
|
[с 5,(/; (р) при некотором |
|
|
фиксированном |
значе |
рис |
нии <р, например с квад- |
||
ратурной составляющей |
|
Вместо линейных детекторов ЛД в этой схеме можно применить детекторы огибающей с любыми мо нотонными характеристиками, так как логика работы решающего блока РБ [см. (3.50)] сводится к отбору максимальной из величин Z, и вынесению решения о при сутствии на входе того сигнала, с которым согласован фильтр, заметнее других прореагировавший на у(/). Отме тим, что уменьшить число каналов различителя до М— 1 в отличие от случая детерминированных сигналов нельзя, так как сравнению в РБ подлежат уже не сами реакции СФ, а результаты их нелинейных преобразований (оги бающие).
Точные формулы для вероятностей ошибок различения М произвольных сигналов достаточно громоздки даже при М = 2, однако в приложениях чаще других встречаются ансамбли сигналов, ортогональных в усиленном смысле.
Последнее означает, |
что любые два несовпадающих сигна |
||||
ла s,(z; (р,), |
sk(t; <pfc) |
ортогональны |
при |
любых |
значениях |
начальных |
фаз: |
|
|
|
|
JsJz; <p,)xt(G <pk)dz = O |
при любых |
<р;, <pt |
и z/A', |
(3.51) |
или, что эквивалентно, ортогональны детерминированные
комплексные огибающие этих |
сигналов: |
fsf(z)s;(z)d/ = o, *ki . |
(3.52) |
Условия (3.51), (3.52) жестче обычного требования ортогональности, фигурировавшего ранее в применении к детерминированным сигналам. Так, два отрезка косинусо иды, сдвинутые на угол + л/2, являясь ортогональными в обычном смысле, не ортогональны при изменении сдвига фаз, т. е. в усиленном смысле. В то же время сигналы, не перекрывающиеся по времени или по спектру, ор тогональны и в усиленном смысле.
69