
Литература / бакулев радиолокация распозн
.pdfнеблагоприятно распределение 1110 =w , а байесово решение ь• - мини
0
максное. Это позволяет облегчить отыскание минимаксных значений и наименее благоприятных априорных распределений, которые часто ока
зываются равномерными.
Критерий последовательной проверки гипотез Вальда. В рас смотренных критериях ограничивалось (фиксировалось) время приня
тия решений Т.шбл или объем выборки у1, у2,.. ·J'k· Однако можно заранее объем выборки не фиксировать. При критерии Вальда область Г делится на три подобласти Г1, ГO и Г2 нижним Т.1 и верхним Тв порогами, как по
казано на рис. 3.3 ,6:
б(у) = d0 , если Лкs Т.1 - в этом случае справедлива гипотеза Но и у
принадлежит области Г0 ; |
· |
|
|
|
б(у) = d1,если Лк~Тв - |
в этом случае справедлива гипотеза Н1 и у |
|||
принадлежит области Г1 ; |
|
|
|
|
б(у) = d2 , если |
Т.1sЛкs Тв - в этом случае |
принимается решение |
||
продолжить наблюдение. |
|
|
|
|
_ |
1v(y1...yk /0 ~ l) |
|
|
|
здесь Лк-Л(у1,У2,, .. ,у") -~-"-----. |
|
|
||
|
w(y1... yk /0 ~ 0) |
|
|
|
Таким образом критерий Вальда двухпороговый: |
|
|||
lтн~Лк~тв.1 |
|
|
|
(3.8) |
Пороги определяются |
вероятностями D и |
F: Тн~ l - D , |
Тв~ D . |
|
Длительность наблюдений - |
|
1-F |
F |
|
величина случайная. Критерий Вальда яв |
ляется оптимальным в смысле минимизации среднего времени наблю
дения (обнаружения) по большему ансамблю экспериментов.
Сведение сложной гипотезы к простой. Если кроме параметров 0=1 и 0 = О, имеются другие:µ в пространстве .Q 1 для случая 0 = l и v в пространстве .Q2 для случая 0 = О с распределениями w(y/µ,l) и w(y/v,O) при известных w0(µ) и w0(v), то можно сформировать отношение прав доподобия, не зависящее от параметров ~L, v:
J1v(y / µ, l)w0 (p)dp
|
|
|
|
== 1-v(y / е=l) |
n, |
(3 .9) |
Л( |
|
/ 0 l |
v)Jl,1' |
|||
|
у |
,1- |
' |
w(y / 0 == О) |
Jw(y / 11, O)w (v)dv . |
|
|
|
|
|
|
0 |
|
|
|
|
|
|
П2 |
|
Структура обнаружителя. В соответствии с полученными алго ритмами обнаружения можно представить их структуру при различных
критериях оптимальности:
51

- однопороrовые критерии с фиксированным временем Т.,абл(п)
(рис. 3.4,а);
- двухпороrовый критерий с переменным временем Тнабл(п)
(рис. 3.4,6).
Устройство Л(у) преобразует
распределение
w(y/0) в распре
деление w(Л). По роговые устрой ства (ПУ), назы
|
|
|
|
|
|
|
ваемые |
реле |
или |
||||||
|
|
|
|
|
|
|
компараторам и, |
||||||||
|
|
|
|
|
|
|
осуществляют |
||||||||
|
|
|
|
|
|
|
сравнение |
|
|
Л с |
|||||
|
6) |
|
|
|
|
порогом |
|
Т. |
Гра |
||||||
Р11с. 3.4. Струк,урные схемы однопороговых (а) и двухпороговых фическая |
интер |
||||||||||||||
(6) обнаруж11телей |
|
|
|
|
|
|
претация |
различ- |
|||||||
ных критериев приведена в таблице. |
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
||||
Критерий |
|
|
Порог |
Гщн\>Ическая интерпретация |
|||||||||||
Однопороrовый |
|
q Со1 -Соо |
|
wrл,tш, |
|
||||||||||
Байеса |
|
|
|
|
|
|
|
|
|
|
|
||||
|
Р С10-С11 |
|
d(I |
|
|
f |
|
,'\ |
|
||||||
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
Т di |
|
• |
|
|
|
||
Однопороrовый |
|
|
|
!.!... |
|
1'1Л!~ |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
максимума |
|
|
|
р |
|
[J |
|
1' |
|
л |
|
||||
апостериорной |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
da |
!!. |
|
J, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
BCf)OЯTIIOCTИ |
|
|
|
|
|
|
|
р |
|
|
|
|
|
|
|
Од11011ороrовый |
|
|
|
1 |
|
W(ЛI~ |
|
|
|
|
|||||
максимума |
|
|
|
|
|
|
|
|
../,.. |
- л |
|
||||
правдоподобия |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d«1 |
1 d1 |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
||||
НейманаПирсона |
f,v(y / O)dy =F |
W(Л!t'>{;h~ |
|
||||||||||||
Одно1юроrовый |
U) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
т |
|
|
|
|
dQ |
|
|
r d, |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Однопоро1·овый |
(_<L) |
с., -с•• |
W(Л)ro,\ |
|
|||||||||||
|
|||||||||||||||
МИIIИМi\КСНЫЙ |
|
|
|
|
|
|
|
|
|
|
|||||
Р |
мм С10-С11 |
|
|
|
|
|
|
|
|
|
|
||||
|
|
<lo |
T.w.,. |
d1 |
|
|
|
|
|||||||
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||
Двухпоро1·овый |
|
|
D |
1-D |
W(Л)оо |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||
Вальда |
Т=- Т=- |
|
|
|
|
|
|
|
|
|
|
||||
|
в |
F |
" 1-F |
|
D |
|
f |
|
|
Л |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
do |
Tw |
dJ ·тw |
tl1 |
' |
|
|
52
Задача обнаружения решается в каждом элементе разрешения. Способ просмотра элементов разрешения определяется выбранным ме тодом обзора пространства. Количество элементов разрешения зависит от величины области обзора или пространства обнаружения. Вероят ность ложной тревоги F и вероятность пропуска цели 1-D обыLJно за даются на все пространство обнаружения. Вероятность правильного не
обнаружения во всем пространстве 1- F = F равна произведению веро
ятностей правильного необнаружения во всех т элементах:
тt/1
F =ПF, =П() -F, ).
1=1 |
i=I |
|
111 |
Если F,=~= ... =f'i=const, то F'=П(I-F,)=(1-Fi)m. |
|
|
i=I |
При F 1<<1 |
можно считать F=(1-Fi)m~l-mFi, или 1-F=F= |
=mF1• Таким образом, если задано F, то
IFi=Flm.j
Пример. Задано F = 10-3, число элементов разрешения по дально
сти mu = 150км/150м = 1ООО, число элементов разрешения по азимуту
та=360°/0,36°= 1ООО, число элементов разрешения по углу места
тр=90°/О,9°=100. Общее число элементов разрешения т = mumam/J = 108
и Fi =1о-з 1108 = 10- 11 •
Что касается D, то вероятность правильного обнаружения равна
1- D, а эта последняя равна произведению вероятности пропуска цели в
одном элементе на вероятность правильного необнаружения во всех ос тальных т -1 элемен:rах:
D=] -[D Fi(m-l)] =] -(1- D )[1-(m- l)F ] ~ 1-(] - D ) =D , |
||||
1 |
1 |
1 |
1 |
1 |
т.е.
ID=D1-I
Таким образом, D = D 1, а F= mF1.
3.3.Модели радиолокационных сигналов
Детер.11инирован11ый сиг11ш1, или сигна1 с полностью известны.wи
параметрами имеет вид
y(t) = 0и(/) + n(t), u(t) = Иш(l) cos[ro 0/ + \/f(/) - <р],
где U111(t) - амплитуда сигнала; ro 0 - несущая частота сигнала; \/f(/) - функция угловой модуляции; <р - начальная фаза колебаний.
53

Все это точно известные неслучайные величины. Считается, что мы знаем время запаздывания сигнала, эффективную площадь рассея ния цели, форму сигнала и все параметры его модуляции. Эта модель наиболее идеализирована.
Квазидетерминированные сигналы, или сигналы со случайными
параметрами имеют вид
y(t) = eu/1/(t,µ) + n(t)'
где µ - вектор случайных параметров сигнала.
Возможны два случая:
а) сигнет со случайной начальной фазой <р:
u(t,<p) = u,,,(t) cos[юоt + ч1(t) - <р],
где <р - неизвестная начальная фаза, распределенная равномерно от О до
2n, т.е. w0(<p) = 1/(2л);
б) сигнал со случайной фазой <р и флуктуирующей аиmитудой аиш(t), где <р - неизвестная начальная фаза с распределением }Vо(<р)=1/(2л);
а - коэффициент флуктуации амплитуды с распределением
}V(a}=~exp{- а\}.
(Уа 2а-а
Таким образом, u(t,<p,a) = aU111(t) cos[юof + ЧJ(f}- <р].
Модели сигналов охватывают случаи одиночных импу,пьсов и па чек импульсов. Пачки импульсов разделяют на пачки когерентных и не когерентных импульсов. Вид пачек показан на рис. 3.5.
,,,,,- -......
а) |
6) |
Рис. 3.5. Форма пачек импульсов при плавном (а) и ступенчатом (б) обзоре пространства
В пачке когерентных радиоимпульсов начальные фазы импульсов
коррелированы: (<pi<pk) -:t:- О . Пачка некогерентных радиоимпульсов со-
стоит из импульсов с независимыми начальными фазами (<p,<pk) = О .
54

Формирование
когерентных радиоим
пульсов и особенности
их спектров показаны на
рис. 3.6. Стабильный за
дающий генератор (син
тезатор частот) форми рует колебание на часто
те /j. После умножителя
частоты (УмЧ) колеба
ния несущей частоты
fo=n.fj попадают на усиР11с. 3.6. Формирование коrере•m{ЫХ радиоимпульсов
литель мощности (УМ), где усиливаются и модулируются импульсами. Формирование некогерентных радиоимпульсов показано на рис. 3.7.
На генератор колебаний |
ра- |
м |
f. |
|
|
||||
|
V |
Or~ |
||
диочастоты, работающии в |
~ |
|
||
режиме самовозбуждения, по- |
м6 |
П ГJ • |
||
даются модулирующие |
им- |
_ |
||
пульсы и формируются мощ- |
f 0 |
1 |
||
|
|
ные радиоимпульсы частотой
/о со случайнои начальной
фазой, что объясняется слуР11с. 3.7. Форм11рование неко,·ерентных р,щионм чайным характером началь пульсов
ных условий самовозбуждения.
Что касается помех: то в дальнейшем рассматриваются модели не
коррелированной помехи - "белый" шум, коррелированной помехи -
"пассивная" помеха и "неrауссовой" активной помехи.
3.4.Синтез оптимальных обнаружителей
Обнаружение одиночного радиоимпульса с полностью извест ными параметрами на фоне «белого» шума. «Белый» шум складыва
ется с сигналом аддитивно: y(t) = 0u(t) + n(t), 0 = О, 1, при этом
M{n(t)} =О,
R(т)= M{n(t)n(t- т)}= No б(т),
2
G(ш)= } R(т)ехр{- jшт}dт= : 0 ,
o-l = R(O) = cFJfG(ш)dш =оо
-мощность шума на входе приемника.
55
В пределах полосы пропускания приемника ci=N0Лf= =КшкТ0ЛJ;
где Кш - коэффициент шума приемника (Кш=ТJТ0+ +Тпрм!Тu-1, То=290К, kT0=4· 1о-21 Вт/Гц). Для прикидочных расчетов можно ориентироваться
на следующие значения коэффициентов шума.
Тип приемника: |
Кш • дБ |
|
с параметрическим усилителем |
..................................................... 2 - |
4 |
сЛБВ ................................................................................................ |
4-8 |
|
с туннельным диодом ...................................................................... |
5 - |
6 |
с балансным смесителем ................................................................. |
6 - |
9 |
Пусть наблюдение входной реализации ведется дискретно через
интервалы времени Лt = Тна6/п в моменты времени tk:{t1,/2, •.. ,t1.,...,tn}. По скольку п = 2.fвТ, отсчеты в моменты tk независимы: y(tk) = Yk = 0щ + nk, k = 1,2,3,4, ... ,п. В этом случае «белый» шум имеет нормальный закон распределения вероятностей (0 = О) и плотность распределения вероят ностей (ПРВ) имеет вид
w(yk /0 = 0) =iv(yk /0) = d-:- |
ехр{- Yi,,}, |
|
|
|
|
|
|
-v21UY |
20-- |
|
|
|
|
|
|
8 силу независимости отсчетов совместная ПРВ |
|
|
|
|
|
|
|
И\У1,У2,Уз,···,У11!0=O)=w(у/O)=п,v(yk/0)= |
г;:- |
ехр |
{ |
|
11 |
. |
|
- ~ LYk |
|||||||
|
,, |
1 |
|
l |
2 } |
|
|
|
k=I |
"ШО- |
|
|
20- |
k=I |
|
При 0 =1 отсчеты входной выборки Yk = щ+пk, поэтому щ =Yгlik и тогда
w(vtf0 = 1) = w(v/1) = г::l-- ехр{ |
(Yk - |
Щ) |
2 |
} • |
|
||||||||
|
|
|
|
|
"2,ro- |
|
2о-2 |
|
|
|
|
||
с учетом независимости отсчетов |
|
|
|
|
|
|
|
||||||
1 |
|
|
|
п |
|
1 |
ехр |
{ |
1 |
|
п |
|
2 } |
w(i !0 =t)= Пw(yk/1)= |
r--- |
-- L<Yk-щ) |
• |
||||||||||
|
|
|
|
|
|
'\J 21ro- |
|
2 |
|
|
|
|
|
|
|
|
|
k=I |
|
20- |
k=I |
|
|
||||
Подставляя w(v! l) и w(v/0) в выражение для Л, получаем |
|||||||||||||
Л=ехр |
{ |
2 |
1 |
Lukyk --- 2LИi . |
|
|
|
|
|
|
|||
, |
|
п |
1 |
п |
'} |
|
|
|
|
|
|
||
|
|
О- |
k=I |
20- |
k=I |
|
|
|
|
|
|
|
Можно сравнивать Л с порогом решения Т, можно также сравни вать монотонную функцию от Л с такой же функцией порога. Напри
мер, часто используют сравнение lnЛ~ lnT. Учитывая, что
1 |
11 |
1 |
п |
111 л =-2 LЩYk --- 2 LИi |
|||
О- |
k=I |
20- |
k=I |
56

и обозначая а-2 |
1 п |
получаем алгоритм оптимального |
|
ln Т+ |
Lиf =ипор , |
||
|
2k=I |
|
tl
обнаружения LllkYk ~ llnop•
k=I
Если перейти к непрерывному времени, то нужно устремить ин тервал дискретизации к нулю Лt О, при этом T,iaбл=const, n oo,
Nт•=~., I тнаб.,
-0 а-2 • Тогда lnЛ -2 Ju(t)y(t)dt-- f u2ct)dt.
2Лt |
|
|
|
|
|
N0 |
о |
|
N0 |
о |
||
|
|
|
|
|
1;.;.б.1 |
|
|
|
|
|||
|
|
|
Обозначим |
f u(t)y(t)dt = =(t) и заметим, что это корреляцион |
||||||||
|
|
|
|
|
|
о |
|
|
|
|
||
|
|
|
|
|
1;136_, |
|
|
|
|
|
|
|
ный интеграл, а |
Ju2 (t)dt = Е - |
энергия сигнала. Объединяя : 0 ln Т+ |
||||||||||
|
|
|
|
|
|
о |
|
|
|
|
||
1 |
т |
|
|
|
N |
|
|
Е |
|
|
|
|
11:16., |
|
|
|
|
|
|
|
|||||
+- |
|
f и2(t)dt |
=- 0 ln Т+- =ипор, получаем алгоритм обнаружения: |
|||||||||
2 |
|
о |
|
2 |
2 |
|
|
|
|
|||
|
|
|
.,. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
z = набf., |
u(t)y(t)dt ;u110P • |
|
|
|
(3.1 О) |
||||
|
|
|
о |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
Структурная |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|||
схема корреляционно |
|
|
|
|
|
|||||||
го обнаружителя (КО) |
|
|
|
|
|
|||||||
показана на рис. 3.8,а. |
|
|
|
|
|
|||||||
|
|
|
Этой |
структуре |
|
|
|
|
|
|||
эквивалентна |
структу |
|
|
|
|
|
||||||
ра |
фильтрового обна |
|
|
|
|
а) |
||||||
ружителя |
|
(ФО) |
|
|
|
|
z |
|||||
(рис.3.8,6). |
|
|
|
|
|
У(~, ОФ |
|
|||||
|
|
|
|
|
|
--1 |
||||||
|
|
|
Импульсная |
ха |
|
|
|
|
||||
рактеристика |
фильт |
|
|
|
|
|
||||||
ра, |
максимизирующе |
|
|
|
|
|
||||||
го отношение мощно- |
|
|
|
|
б |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
сти сигнала к мощно- |
|
|
|
|
||||||||
сти |
|
шума q, |
|
Р11с. 3.8. СтруКl)'ры карреляцнаннаrа (КО) (а) и фильтра |
||||||||
|
является вага (ФО) |
(6) абнаружrт:лей |
||||||||||
зеркальным |
|
отобра- |
|
|
|
|
||||||
жением входного сигнала ГJ(t) |
= и(Тс - t), поэтому выходной сигнал |
имеет форму
57
т т
up1,,x (Те)= J17(~ - t)y(t)dt = fu(t - ~)y(t)dt =z .
оо
Соотношение для коэффициента передачи согласованного фильтра получаем по формуле
ш
k(jш) = f!]{t)exp{ - jcot}dt =ехр{-jшt}S*(jш),
при этом формируется максимальное отношение сигнала к шуму:
q111ax = 2EIN0 (при реШ!ьных шумах в диапазоне частот О оо q = EIN0 ).
Поскольку выходной сигнал обнаружителей, описываемый корреляци
онным интегралом, зависит от времени запаздывания tu и расстройки по
частоте сод, корреляционный обнаружитель оказывается многоканаль ным по дальности и скорости. Фильтровой обнаружитель многоканален
только по скорости.
В КО на выходе инерционного фильтра (ИнФ) будет нарастающее
напряжение в момент fu+t'и, равное z(т11), в то время как в ФО на выходе
возникает радиоимпульс, по форме совпадающий с корреляционной функ
цией входного сигнала. Для устранения колебаний внутри огибающей ра диоимпульса на выходе ОФ ставят детектор огибающей (ДО) (рис. 3.9, а-г).
Для определения величины вероятностей правильного обнаруже
ния D и ложной тревоги F надо знать плотности распределения вероят ностей величины z на входе порогового устройства при 0 = О и 0 =). Ес
ли 0 = О, то на входе - толы<о шум n(t), поэтому y(t)=n(t). Операция ин
тегрирования является линейной, и(t) - детерминированная величина.
Следовательно z(t) будет иметь то же распределение вероятностей, что и n(t), т.е. нормальное распределение с параметрами
M{z/0 = О}=О,
o/=M{z2}-(M{z} )2=M{z2 }= 7~абf., 1;r.абf., u(t)u(t')M {n(t)n(t')}dtdt' =
|
|
|
|
|
|
о |
|
о |
|
|
|
Т |
, |
7' ~ |
N |
|
|
|
|
N |
Т |
|
N Е |
нао., |
нао., |
|
|
|
|
наб., |
|||||
= |
f |
f u(t)u(t')-fб(t-t')dtdt'=-f |
|
f |
и2ct')dt=T· |
||||||
|
о |
о |
|
|
|
|
|
|
|
о |
|
Таким образом, |
|
|
|
|
|
|
|
|
|||
w(x/0 =0) = |
] |
ехр |
{ |
.,.2 |
} |
. |
|
|
|
||
г;::- |
|
- ~ |
|
|
|
||||||
|
|
|
v2tr6z |
|
|
26= |
|
|
|
|
|
58

При 0 =1 меняется лишь среднее значение из-за того, что u(t) - де
терминированная величина:
М{z 10 =1} =М{rJu(t)[u(t) +n(t))dt}} =7'и'(l)dt=Е. |
|||
|
|
|
1 |
Следовательно, |
|
|
|
|
J'(t)/ tl)'(t)U(tч ИвФ 1 |
||
W(z/0 = 1) = |
|
|
|
=-1-ехр{ (z--E} |
2 }· |
||
'Гша-= |
~ |
|
U(t) |
|
|
|
|
Отсюда можно |
вы |
а) |
|
числить искомые |
вероят |
|
|
ности Fи D: |
|
|
|
в
Рис. 3.9. Форма С11ГШ1ЛОВ в коррСЛЯЦИОIIНОМ (а).
фильтровом (6) и фильтровом с детектором оп1баю щей (в) обнаружителях
где h = u110/z; Ф(h)- интеграл вероятности.
Аналогично вычисляется вероятность правильного обнаружения: |
||||||||||||
D = |
rл |
|
|
|
{ ( |
=- |
Е2)2 |
} d== |
а, |
{ |
,, |
r-y- |
|
f Б |
ехр |
|
|
|
~ fехр |
|
_У~ |
|
|||
|
|
21ra-? |
|
|
2а-.,. |
|
"27r |
|
2 |
|
||
|
11пор |
- |
|
|
|
- |
|
_,, |
|
|
|
|
-;,_,, |
h-,Jц;;;; |
ехр |
{ |
21,1 |
|
|
|
|
|
|
||
1 |
|
_J; |
JY=1-Ф(h-.jq.,~), |
|
|
|
где 'lmax=E/No.
Обнаружение квазидетерминированного радиоимпульса. «Бе
лый» шум аддитивно складывается с квазидетерминированным сигна
лом, у которого начальная фаза неизвестна: u(t) = U,u(t)cos[ro 0!+\/f(l)-<p].
Здесь <р - неизвестная начальная фаза, распределенная равномерно в пре делах от О до 21r, с плотностью распределения вероятностей w(<p) = 1/21r.
При этом отношение правдоподобия имеет вид
59
