
книги из ГПНТБ / Вайнштейн Л.А. Выделение сигналов на фоне случайных помех
.pdfгде Р (т), Р (0), р[т(у)] и ртО имеют тот же смысл,
что и аналогичные выражения при сложном обнаружении.
При измерении нас |
интересует вопрос о том, |
присутст |
|
вует ли во входном сигнале f |
полезный сигнал т или нет, |
||
а если присутствует, |
то что |
можно сказать об |
его неиз |
вестном параметре т. В результате приема мы находим плот
ность апостериорной вероятности для всех возможных зна
чений |
измеряемого параметра |
|
|
|
|
Pf[fn (^)]=------ p(f} |
■ |
(29.27) |
|
С учетом формул (29.26) мы получаем следующее |
выра |
|||
жение для апостериорной плотности вероятности |
|
|||
|
pf[m(z)] = |
А(р(й’ |
|
(29.28) |
|
Л + Р(т) |
|
|
|
являющейся функцией т, причем |
|
|
|
|
|
А(х) = рто(т) |
P0(f) |
|
(29.29) |
|
А = J Л (т) Л. |
|
(29.30) |
|
4. |
Сложное измерение — сигнал с |
неизвестными |
пара |
|
метрами и и 6, из которых параметр т |
измеряется, а пара |
метр 0 не измеряется; возможные значения z и 0 образуют непрерывную совокупность. Вероятность появления функ
ции f (t) в |
этом случае равна |
|
|
|
Р (Л - П р |
(х,о)] рм (Л dxdh + р (0) р0 (/-) = j |
|||
= р и п рт |
рт^ |
м р (°)(/')’ |
I |
|
|
|
РИ + Р(0) = 1, |
I |
|
/7[^(т,0)] — Р (tri) рт (х,0), |
JJ>,„(^O)^rfO= 1, |
J |
||
причем обозначения понятны из предыдущего. |
интересует |
|||
Так как |
нас |
интересует |
параметр т и не |
параметр 6, мы должны образовать плотность апостериор ной вероятности для параметра т
Гр\т (т,6)] Pml, |
(f) d8 |
Р/ [т (т)] = j pf [т (т,0) ]</0 = |
, (29.32) |
12* |
179 |
произведя интегрирование |
по |
всем возможным |
значениям |
||||||||
параметра 0. Иначе можно написать |
|
|
|
||||||||
|
|
|
|
|
|
= —(29.33) |
|||||
|
|
|
|
|
|
|
A+?(m) |
|
|
|
|
где |
обозначено |
|
|
|
|
|
|
|
(29.34) |
||
|
|
|
Л ()* |
= ]■ Л (т,6)с?6, |
|
|
|||||
|
|
|
Л = JJ Л (т.б) fifx |
|
|
(29.35) |
|||||
|
|
|
л (’,(1) = .*/,„(!>) |
|
Р;У,1<П ■ |
|
|
(29.36) |
|||
три |
В данной задаче, таким образом, приходится |
различать |
|||||||||
величины: |
А, Л (-с) |
и |
Л (х,0), |
причем |
Л (т,6) |
есть коэф |
|||||
фициент |
правдоподобия |
при |
измерении |
параметров |
т и 6; |
||||||
А (т) — коэффициент |
правдоподобия при |
измерении |
одного |
||||||||
параметра т, а |
А — коэффициент правдоподобия |
при |
слож |
||||||||
ном обнаружении сигнала т (t, -г, |
0). |
|
|
|
|||||||
|
Если параметры т |
и |
0 |
независимы, то |
|
|
|
||||
|
|
|
/>„(’■ |
°) = />„(’) |
|
|
(29.37) |
||||
и формула (29.34) несколько |
упрощается: |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
<29-38) |
|
Вышеприведенные |
формулы |
дают возможность |
найти |
|||||||
|1вид оптимального приемника, |
если известны статистические |
||||||||||
свойства |
помехи n(t) |
и |
неизвестных параметров |
полезного |
|||||||
[.■сигнала |
т (/). |
Прежде чем приступить к |
соответствующим |
вычислениям, мы в следующем параграфе рассмотрим
вопрос об |
априорных вероятностях. |
|
иногда коэф |
||||
|
Заметим в заключение., |
что в |
литературе |
||||
фициенты |
правдоподобия |
Л (т), |
A (x,0) и |
т. |
д. |
вводятся |
|
по |
формулам |
|
|
|
|
|
|
|
|
|
|
|
|
|
<29'39> |
не |
включающим априорных вероятностей |
рт(р) |
и Pm(pfii). |
В этом случае выписанные выше формулы имеют несколько иной вид.
180
§ 30. |
АПРИОРНЫЕ ВЕРОЯТНОСТИ. РЕШАЮЩИЕ СХЕМЫ |
В |
статистической теории оптимальных приемников, |
основные понятия которой были рассмотрены в предыду щих параграфах, вопрос об априорных вероятностях по лезного сигнала связан с определенными трудностями. Дей ствительно, априорные вероятности нужны для вычисления апостериорных вероятностей, т. е. они необходимы для фактического осуществления оптимального приемника. Од нако априорные вероятности часто неизвестны. Так, Вуд ворд пишет: «Рассмотрим, например, априорную вероят ность обнаружения самолета некоторой радиолокационной установкой на расстоянии 16 км завтра в 9 ч утра. Если установка расположена на аэродроме с регулярным дви жением, статистический анализ прошлого может дать нам нужные вероятности в предположении, что движение са молетов представляет собой стационарный случайный про
цесс. Для большого класса задач, однако, мы не распо лагаем статистикой либо потому, что она не изучалась, либо вследствие более фундаментального обстоятельства: в прошлом не существовало совокупности сходных ситуа ций, из которой можно было бы вывести определенное суж дение».
Как |
мы показали в § |
29, плотности априорных вероят |
|
ностей |
р [т (т)], р \т (6)] и р [т (т, 0)] можно представить |
||
в виде |
двух множителей |
|
|
Р [т (0] = Р (т) рт (т), |
р [т. (0)] = Р (т) рт (6), |
(30.01) |
|
|
/?[/71(11,6)] = Р (т) рт(ъ, 0). |
||
|
|
||
Априорные вероятности |
Р(т) и Р(0) = 1—Р(т) являются |
соответственно вероятностями наличия и отсутствия по лезного сигнала на входе приемника. Эти вероятности наи более трудно оценить. Априорные вероятности рт (т), pmW)
и pn(t, 0) являются вероятностями распределения полез |
|
ных сигналов по неизвестным параметрам |
при условии, |
что полезный сигнал присутствует на входе |
приемника. |
Эти распределения в ряде случаев можно более или менее уверенно найти из теоретических соображений. Так, на пример, случайную высокочастотную фазу при некогерент ном приеме естественно предположить равномерно распре деленной по окружности, амплитуду флюктуирующего сиг
нала— по закону Релея. |
Дальность и азимут цели можно |
в некоторой небольшой |
области воздушного пространства |
181
предположить равномерно распределенными; при увеличе
нии размеров |
области это предположение может стать |
уже несправедливым. |
|
Учитывая выше приведенные рассуждения и предпола |
|
гая, что закон |
распределения априорных вероятностей по |
лезного сигнала по неизвестным параметрам т и 0 известен,
мы можем вычислить введенные выше для различных
случаев коэффициенты правдоподобия Л, Л (0), Л (х) и А (т;, 0).
Если далее образовать отношение апостериорных |
вероят |
|||||||||
ностей |
присутствия |
и |
отсутствия |
полезного сигнала, то |
||||||
получим |
при обнаружении |
|
|
|
|
|
||||
|
|
|
|
РЛт) Р (т) |
|
|
<30'02) |
|||
|
|
|
|
Д(и) = Р(0)А' |
|
|||||
а при |
измерении |
|
|
|
|
|
|
|
||
|
|
|
РДт(Ч] |
Р(т) |
|
|
(30М) |
|||
|
|
|
|
Л(О)-Р(О)Л(Х)- |
||||||
Эти формулы нетрудно |
вывести |
из |
выражений (29.09), |
|||||||
(29.22), |
(29.28), |
(29. |
33) |
и соотношений |
|
|
||||
|
|
|
|
|
|
|
|
Р(0) |
|
|
|
|
Р; (0) = 1 |
- Р |
|
= |
\ I |
, |
(30.04) |
||
|
|
|
|
|
|
|
|
_1_2_ |
|
|
|
|
|
|
|
|
|
|
|
Р (т) |
|
|
|
|
Pf (т) = СР{ [т (х)]сК |
(30.05) |
||||||
Формулы (30.02) |
и (30.03) показывают, что в |
отноше |
||||||||
ниях |
апостериорных |
вероятностей от априорных вероятно |
||||||||
стей |
Р(т) и |
Р (0) |
зависит |
лишь |
постоянный |
множи |
||||
тель |
Р(/ц)/7э(0), а |
принятая |
функция f (t) определяет коэф |
|||||||
фициенты правдоподобия Л и А (т). |
|
|
|
|||||||
Трудность, |
обусловленную |
незнанием отношения Ррп^, |
можно обойти, если изменить определение оптимального
приемника и назвать оптимальным приемник, образующий
коэффициенты правдоподобия (а не апостериорные вероят ности). В таком случае оптимальные приемники по опреде-
182
лению должны выдавать следующие математические ве личины:
1) |
при |
простом обнаружении |
А, |
| |
|
2) |
при |
сложном обнаружении |
А = Ja (G)d<), |
I |
|
3) |
при |
простом измерении |
А (т), |
| (30.0.’>) |
|
4) |
при |
сложном измерении |
А (т) = J А (т,0) М. |
I |
На основании входных данных и образованных с их
помощью величин (30.06) обычно приходится принимать
решения. Если решать должен человек, например ответить
«есть сигнал» или «нет сигнала», то оптимальный прием ник лишь помогает человеку, оставляя за ним операцию решения. Надо сказать, что в своих решениях человек всегда использует (часто ,не осознавая этого явно) апри орные знания о вероятности появления сигнала: в частно сти, если априорная вероятность появления сигнала доста точно мала, то для ответа «есть сигнал» потребуется бо
лее сильное превышение сигнала над шумами, т. е. боль шее значение Л.
Процесс решения нетрудно автоматизировать. Ограни
чиваясь |
задачей обнаружения |
(сложного или простого), |
|||||||
мы должны учесть, что вероятность наличия |
полезного |
||||||||
сигнала |
|
|
|
|
|
|
|
|
|
|
|
|
|
= |
|
|
|
|
(30.07) |
|
|
|
|
А + Р(т) |
|
|
|
||
есть монотонная функция |
коэффициента |
правдоподобия А. |
|||||||
Совершенно |
естественно |
считать, |
что |
сигнал tn |
присут |
||||
ствует, |
если |
вероятность |
Pf (т) достаточно велика (т. е. |
||||||
достаточно близка к единице), |
и |
что |
полезного |
сигнала |
|||||
нет, если вероятность РДт,') |
достаточно мала. |
Поэтому |
|||||||
простейшее правило решения имеет вид |
|
|
|
||||||
|
|
если |
Pf(m)^z Р*, |
то |
сигнал есть, |
) |
(30.08) |
||
|
|
если |
Pf>*(m)<ZP |
то |
сигнала нет, |
| |
где Р* —некоторое „пороговое" значение вероятности, ска
жем, Р* — 0,5; *Р |
= 0,9 или |
*Р = 0,99. |
|
|
Более сложное правило: |
|
|
||
если Pf (т) > Р*, |
то |
сигнал есть, |
1 |
|
если *Р > Pf (т) > ,Р* |
то имеем неопределенность, |
430.09) |
||
если *,РР |
то |
сигнала |
нет, |
| |
183
с двумя порогами Р* и Р# использует апостериорные
вероятности на выходе оптимального приемника более пол но, но при этом иногда дает неопределенный ответ. Если
сигнал |
f(t) принят, |
дальнейшая информация в |
приемник |
||
не поступает и |
на |
основании имеющихся сведений |
тре |
||
буется |
принять |
какое-то определенное решение, |
то |
един |
ственный выход заключается, очевидно, в применении пра вила (30.08) с одним порогом. Если же информация посту- ' пает в приемник постепенно, то на основании входных данных, накопившихся за фиксированный промежуток вре мени, можно принять и неопределенное решение, указы вающее на необходимость продолжать наблюдение. В этом случае можно применить «двухпороговое» правило (30.09);
в принципе можно было бы, вероятно, использовать и бо
лее сложные правила.
Рассмотрим более подробно правило (30.08). Коль скоро мы выберем одно из двух возможных решений, то мы всегда можем или принять правильное решение или ошибиться. Ошибки могут быть двух типов. Первый тип ошибки — принятие решения «да», когда на входе присут ствует только помеха. Эта ошибка называется ложной
тревогой, ее |
вероятность |
мы обозначим через F. Второй |
тип ошибки — принятие |
решения «нет», когда на входе |
|
присутствуют |
как помеха, |
так и полезный сигнал. Эта |
ошибка называется пропуском сигнала, вероятность этой
ошибки мы будем обозначать через Do. Вероятность лож ной тревоги F является вероятностью принять помеху за сумму сигнал + помеха; вероятность пропуска Do есть ве роятностью принять сумму сигнал + помеха за чистую по меху.
Правильные решения также могут быть двух типов: правильное обнаружение и правильное необнаружение. Вероятность правильного обнаружения, которую мы
обозначим через D, есть вероятность принять сумму сиг-
нал + помеха за сигнал + помеха, а вероятность правиль ного необнаружения, которую мы обозначим через Fo, есть вероятность принять помеху за помеху. Очевидно, что F, Dq, D и Fo — условные вероятности: Ео и F— вероятности принять правильное или неправильное решение при усло
вии, что полезного сигнала нет, D и Do — такие же вероят ности при условии, что полезный сигнал присутствует. По этому выполняются соотношения
Ов=1 —£); F0=l -F. |
(30.10) |
184
Полная вероятность |
принять |
правильное |
решение, оче |
||||||||||
видно, |
равна |
|
|
|
|
|
|
|
|
|
|
|
|
Г = р (т) D 4- Р (0) Fo = Р (т) D Д- Р (0) (1 |
— F), |
(30.11) |
|||||||||||
где Р(0) и Р (т) |
суть |
априорные вероятности отсутствия |
|||||||||||
и наличия сигнала т. |
|
|
|
|
|
|
|
|
|
||||
При использовании правила (30.08) необходимо задать, |
|||||||||||||
помимо порога Р , априорные |
вероятности |
|
Р (0) |
и Р (т). |
|||||||||
Если |
|
последние |
неизвестны, |
то |
можно |
воспользоваться, |
|||||||
как это было указано выше, |
коэффициентом |
правдоподо |
|||||||||||
бия, |
с |
помощью |
которого правило |
(30.08) |
перепишется |
||||||||
в виде |
|
если |
А > А *, |
то |
сигнал есть, |
|
|
|
|
||||
|
|
|
|
|
|
(30.12) |
|||||||
|
|
|
если |
А <С Л , |
то |
сигнала |
нет, |
|
|
|
|||
|
|
|
|
|
|
|
|||||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
_ Р(0) |
|
Р, |
|
|
|
|
(30.13) |
|
|
|
|
|
|
|
Р (т) ' I — Р |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|||
есть |
пороговое |
значение |
коэффициента |
правдоподобия. |
|||||||||
„Двухпороговое" правило (30.09) примет такой вид: |
|||||||||||||
|
|
если |
А > Л’, |
то |
сигнал есть, |
|
|
|
) |
|
|||
|
|
если |
. |
|
|
|
|
|
|
|
|
I |
(30.14) |
|
|
Л >• Л > Л# — неопределенность, |
j. |
||||||||||
|
|
если |
Л < Л#, то сигнала |
нет. |
|
|
|
j |
|
||||
Согласно этим правилам нетрудно построить схемы, |
|||||||||||||
автоматически |
принимающие |
решения. Таким |
образом, |
„решающий" оптимальный приемник должен образовывать коэффициент правдоподобия А и подавать его на вход ре
шающей схемы (30.12) или (30.14). Заметим, что вместо А можно использовать любую монотонно возрастающую функцию Л (например, In А), что часто упрощает схему оптимального приемника. Порог Л в формуле (30.12)
обычно находят из требования, чтобы вероятность ложных тревог равнялась заданному значению (часто весьма малому, например, /?=10-3, Е=10“5 или С=10~'°).
Остановимся в заключен е на терминологии, принятой в лите ратуре.
Наблюдателем Неймана-Пирсона (Neymann-Pearson) называют на блюдателя, который на основании принятых данных принимает реше ния о наличии сигнала m(t) по правилу, которое обеспечивает мак-
185
снмальную вероятность правильного обнаружения D при фиксиро ванной вероятности ложной тревоги F за данный промежуток вре
мени наблюдения Т. В математической статистике |
доказывается, что |
|||||
наблюдатель Неймана-Пирсона принимает решения |
как |
раз по „одно |
||||
пороговому “ правилу (30.12), |
причем |
величина |
порога |
А |
определяется |
|
фиксированным значением F. |
Любое |
другое |
правило |
решения приво |
дит к меньшим D (при заданных F и Г).
Идеальный наблюдатель Зигерта (Siegert) принимает решение, обеспечивающее максимальную вероятность W по формуле (30.11) при фиксированном времени наблюдения Т. Решение принимается
также по правилу (30.12), но величина порога |
выбирается |
равной |
|
. |
^(°) |
|
|
* * |
Р (m) ’ |
производит |
|
|
Последовательный наблюдатель Вальда (Wald) |
анализ |
данных, непрерывно поступающих на вход приемника. Последова тельный наблюдатель имеет возможность задержать решение до по ступления новых данных; правило решения для него имеет вид (30.14). Однако математическая теория последовательного наблюдения отли чается большей сложностью, и мы в дальнейшем будем исключи тельно применять схему решения (30.12) с одним порогом, интер претируя ее в духе наблюдателя Неймана-Пирсона.
Более глубокий подход к статистической теории приема дает
современная теория игр |
и статистических |
решений, |
использованная |
в теории оптимальных приемников Д. ван |
Метером и |
Д. Мидлтоном. |
|
Некоторые относящиеся |
сюда .вопросы рассмотрены в |
приложении I. |
§ 31. ПРОСТОЕ ОБНАРУЖЕНИЕ СИГНАЛА НА ФОНЕ КОРРЕЛИРОВАННЫХ НОРМАЛЬНЫХ ПОМЕХ
Согласно классификации |
§ 29 простое обнаружение — |
это обнаружение полностью |
известного сигнала на фоне |
.помех. В данной задаче полезный сигнал ,m(Z) может либо отсутствовать, либо быть вполне определенной функцией времени.
Что касается помехи, то в теории оптимальных прием ников мы будем считать, что она является стационарным
случайным процессом нормального (гауссова) |
типа со |
средним значением, равным нулю |
|
/Г(“0 = 0 |
(31.01) |
и произвольной функцией корреляции |
|
Rn ()* = n(t)n(t — z), |
(31.02) |
которая полностью определяет статистические свойства по мехи (см. гл. IX). Нормальный характер помехи позволяет сравнительно просто вычислять различные вероятности,
связанные с помехой, поскольку мы опираемся на распре-
деления Гаусса. Стационарность помехи не очень важна
186
в общих рассуждениях, однако лишь для стационарных помех удается в большинстве задач довести исследование
до выводов, которые могут быть эффективно использованы на практике.
Данные предположения о свойствах помехи позволяют охватить ряд помех, представляющих технический инте рес, в частности собственные шумы приемника и радиоло кационные помехи, обусловленные хаотическими отраже ниями (см. гл. IX и XI).
Начнем рассмотрение простого обнаружения со случая,
когда входной процесс f(t) нам известен только в дискрет ные моменты времени
=1)Аг\ /г=1,2, ...Н. (31.03)
Значения / (Д) мы будем обозначать через fh и называть
выборками входного процесса или элементами входной последовательности. Аналогично вводим выборки полезного сигнала т (t) и помехи п (t) по формулам
= n^n(th). (31.04)
Очевидно, что формула
f (t) = т (/) п (0 |
(31.05) |
примет вид |
|
^h = mh + nh- |
<31-06) |
Итак, если за время наблюдения мы сделали Н выбо рок, то корреляционная матрица помехи будет иметь эле менты
Rgh = ngnii = ( S ~ h ! A0; £- ^=1- 2, ... H.
(31.07)
Она обладает следующими особенностями: но главной диа гонали стоят элементы /?Д0), имеющие наибольшую аб
солютную величину; элемент, стоящий на h столбцов ле
вее или правее главной диагонали, равен Rn (hAt). Обратная матрица Qgh удовлетворяет соотношению
<3|08>
/-1
187
где Sgh — символ Кронекера, равный 1 при g = h и 0 при
g=£h. Элементы этой матрицы определяются соотноше нием
(31'09)
где Det || 7? |
д || есть |
детерминант, |
соответствующий корре |
|||||||||||
ляционной |
матрице, |
Ad (/? |
ft) — адьюнкта |
(алгебраическое |
||||||||||
дополнение) |
элемента |
R h. |
Обратная |
матрица |
Qgft, |
как |
и |
|||||||
прямая |
матрица, симметрична |
|
|
|
|
|
|
|
|
|||||
|
|
Q . =0 • R.=R. , |
|
|
|
(31.10) |
||||||||
|
|
^■gh |
|
^-kg’ |
gh |
|
xhg' |
|
|
|
' |
' |
||
однако |
ее элементы |
в |
общем |
случае зависят |
от |
g и |
h, |
|||||||
а не только от |g\—- h\. |
|
|
|
|
|
|
|
|
Считая |
|||||
Найдем коэффициент правдоподобия (29.10). |
||||||||||||||
помеху /z(Z) |
нормальной, |
мы можем сразу |
написать много |
|||||||||||
мерную |
плотность |
вероятности |
для |
выборок |
|
..., |
пн |
|||||||
в виде |
[см. формулу (59.13)] |
|
|
|
|
|
|
|
|
|
||||
|
|
/?(п) = р(п1, •••, |
пн) |
= |
|
|
|
|
|
|||||
— - |
1 |
|
= |
ехр / —4 V. Qahnjih 1 |
■ |
(31.11) |
||||||||
|
WDetliy |
Ч |
2 |
LiJsfl |
s Ч |
|
|
|
||||||
Это //-мерное распределение Гаусса. |
|
|
fH |
при |
отсут |
|||||||||
Плотность вероятности величин fr, |
|
|||||||||||||
ствии сигнала т, когда f =nh, |
равна |
|
|
|
|
|
||||||||
/’о(Л==/’(Л=-^=^== ехр/-4 |
У Qgd'gfA- |
|
||||||||||||
|
|
/(2^Det || ^ft|| |
|
I |
2 |
|
8 |
|
f |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
(31.12) |
|
Плотность вероятности |
тех |
же |
величин |
при наличии сиг |
нала, когда nh = fh — mk, равна
=— w) =
= —- --1 |
exp / —у V Q (j — m ) (Jh — mh) I. |
(31.13)
188